Integrated analysis of patient networks and plasmid genomes reveals a regional, multi-species outbreak of carbapenemase-producing Enterobacterales carrying both blaIMP and mcr-9 genes.

BACKGROUND: Carbapenemase-producing Enterobacterales (CPE) are challenging in healthcare, with resistance to multiple classes of antibiotics. This study describes the emergence of IMP-encoding CPE amongst diverse Enterobacterales species between 2016 and 2019 across a London regional network. METHODS: We performed a network analysis of patient pathways, using electronic health records, to identify contacts between IMP-encoding CPE positive patients. Genomes of IMP-encoding CPE isolates were overlayed with patient contacts to imply potential transmission events. RESULTS: Genomic analysis of 84 Enterobacterales isolates revealed diverse species (predominantly Klebsiella spp, Enterobacter spp, E. coli); 86% (72/84) harboured an IncHI2 plasmid carrying blaIMP and colistin resistance gene mcr-9 (68/72). Phylogenetic analysis of IncHI2 plasmids identified three lineages showing significant association with patient contacts and movements between four hospital sites and across medical specialities, which was missed on initial investigations. CONCLUSIONS: Combined, our patient network and plasmid analyses demonstrate an interspecies, plasmid-mediated outbreak of blaIMPCPE, which remained unidentified during standard investigations. With DNA sequencing and multi-modal data incorporation, the outbreak investigation approach proposed here provides a framework for real-time identification of key factors causing pathogen spread. Plasmid-level outbreak analysis reveals that resistance spread may be wider than suspected, allowing more interventions to stop transmission within hospital networks.

RamanSPy: An Open-Source Python Package for Integrative Raman Spectroscopy Data Analysis.

Raman spectroscopy is a nondestructive and label-free chemical analysis technique, which plays a key role in the analysis and discovery cycle of various branches of science. Nonetheless, progress in Raman spectroscopic analysis is still impeded by the lack of software, methodological and data standardization, and the ensuing fragmentation and lack of reproducibility of analysis workflows thereof. To address these issues, we introduce RamanSPy, an open-source Python package for Raman spectroscopic research and analysis. RamanSPy provides a comprehensive library of tools for spectroscopic analysis that supports day-to-day tasks, integrative analyses, the development of methods and protocols, and the integration of advanced data analytics. RamanSPy is modular and open source, not tied to a particular technology or data format, and can be readily interfaced with the burgeoning ecosystem for data science, statistical analysis, and machine learning in Python. RamanSPy is hosted at https://github.com/barahona-research-group/RamanSPy, supplemented with extended online documentation, available at https://ramanspy.readthedocs.io, that includes tutorials, example applications, and details about the real-world research applications presented in this paper.

Systems level profiling of chemotherapy-induced stress resolution in cancer cells reveals druggable trade-offs.

Cancer cells can survive chemotherapy-induced stress, but how they recover from it is not known. Using a temporal multiomics approach, we delineate the global mechanisms of proteotoxic stress resolution in multiple myeloma cells recovering from proteasome inhibition. Our observations define layered and protracted programs for stress resolution that encompass extensive changes across the transcriptome, proteome, and metabolome. Cellular recovery from proteasome inhibition involved protracted and dynamic changes of glucose and lipid metabolism and suppression of mitochondrial function. We demonstrate that recovering cells are more vulnerable to specific insults than acutely stressed cells and identify the general control nonderepressable 2 (GCN2)-driven cellular response to amino acid scarcity as a key recovery-associated vulnerability. Using a transcriptome analysis pipeline, we further show that GCN2 is also a stress-independent bona fide target in transcriptional signature-defined subsets of solid cancers that share molecular characteristics. Thus, identifying cellular trade-offs tied to the resolution of chemotherapy-induced stress in tumor cells may reveal new therapeutic targets and routes for cancer therapy optimization.

ProteinLens: a web-based application for the analysis of allosteric signalling on atomistic graphs of biomolecules.

The investigation of allosteric effects in biomolecular structures is of great current interest in diverse areas, from fundamental biological enquiry to drug discovery. Here we present ProteinLens, a user-friendly and interactive web application for the investigation of allosteric signalling based on atomistic graph-theoretical methods. Starting from the PDB file of a biomolecule (or a biomolecular complex) ProteinLens obtains an atomistic, energy-weighted graph description of the structure of the biomolecule, and subsequently provides a systematic analysis of allosteric signalling and communication across the structure using two computationally efficient methods: Markov Transients and bond-to-bond propensities. ProteinLens scores and ranks every bond and residue according to the speed and magnitude of the propagation of fluctuations emanating from any site of choice (e.g. the active site). The results are presented through statistical quantile scores visualised with interactive plots and adjustable 3D structure viewers, which can also be downloaded. ProteinLens thus allows the investigation of signalling in biomolecular structures of interest to aid the detection of allosteric sites and pathways. ProteinLens is implemented in Python/SQL and freely available to use at: www.proteinlens.io.

Kernel-Based Independence Tests for Causal Structure Learning on Functional Data.

Measurements of systems taken along a continuous functional dimension, such as time or space, are ubiquitous in many fields, from the physical and biological sciences to economics and engineering. Such measurements can be viewed as realisations of an underlying smooth process sampled over the continuum. However, traditional methods for independence testing and causal learning are not directly applicable to such data, as they do not take into account the dependence along the functional dimension. By using specifically designed kernels, we introduce statistical tests for bivariate, joint, and conditional independence for functional variables. Our method not only extends the applicability to functional data of the Hilbert-Schmidt independence criterion (hsic) and its d-variate version (d-hsic), but also allows us to introduce a test for conditional independence by defining a novel statistic for the conditional permutation test (cpt) based on the Hilbert-Schmidt conditional independence criterion (hscic), with optimised regularisation strength estimated through an evaluation rejection rate. Our empirical results of the size and power of these tests on synthetic functional data show good performance, and we then exemplify their application to several constraint- and regression-based causal structure learning problems, including both synthetic examples and real socioeconomic data.

Learning compositional sequences with multiple time scales through a hierarchical network of spiking neurons.

Sequential behaviour is often compositional and organised across multiple time scales: a set of individual elements developing on short time scales (motifs) are combined to form longer functional sequences (syntax). Such organisation leads to a natural hierarchy that can be used advantageously for learning, since the motifs and the syntax can be acquired independently. Despite mounting experimental evidence for hierarchical structures in neuroscience, models for temporal learning based on neuronal networks have mostly focused on serial methods. Here, we introduce a network model of spiking neurons with a hierarchical organisation aimed at sequence learning on multiple time scales. Using biophysically motivated neuron dynamics and local plasticity rules, the model can learn motifs and syntax independently. Furthermore, the model can relearn sequences efficiently and store multiple sequences. Compared to serial learning, the hierarchical model displays faster learning, more flexible relearning, increased capacity, and higher robustness to perturbations. The hierarchical model redistributes the variability: it achieves high motif fidelity at the cost of higher variability in the between-motif timings.

Development and Delivery of a Real-time Hospital-onset COVID-19 Surveillance System Using Network Analysis.

BACKGROUND: Understanding nosocomial acquisition, outbreaks, and transmission chains in real time will be fundamental to ensuring infection-prevention measures are effective in controlling coronavirus disease 2019 (COVID-19) in healthcare. We report the design and implementation of a hospital-onset COVID-19 infection (HOCI) surveillance system for an acute healthcare setting to target prevention interventions. METHODS: The study took place in a large teaching hospital group in London, United Kingdom. All patients tested for SARS-CoV-2 between 4 March and 14 April 2020 were included. Utilizing data routinely collected through electronic healthcare systems we developed a novel surveillance system for determining and reporting HOCI incidence and providing real-time network analysis. We provided daily reports on incidence and trends over time to support HOCI investigation and generated geotemporal reports using network analysis to interrogate admission pathways for common epidemiological links to infer transmission chains. By working with stakeholders the reports were co-designed for end users. RESULTS: Real-time surveillance reports revealed changing rates of HOCI throughout the course of the COVID-19 epidemic, key wards fueling probable transmission events, HOCIs overrepresented in particular specialties managing high-risk patients, the importance of integrating analysis of individual prior pathways, and the value of co-design in producing data visualization. Our surveillance system can effectively support national surveillance. CONCLUSIONS: Through early analysis of the novel surveillance system we have provided a description of HOCI rates and trends over time using real-time shifting denominator data. We demonstrate the importance of including the analysis of patient pathways and networks in characterizing risk of transmission and targeting infection-control interventions.

Learning spatiotemporal signals using a recurrent spiking network that discretizes time.

Learning to produce spatiotemporal sequences is a common task that the brain has to solve. The same neurons may be used to produce different sequential behaviours. The way the brain learns and encodes such tasks remains unknown as current computational models do not typically use realistic biologically-plausible learning. Here, we propose a model where a spiking recurrent network of excitatory and inhibitory spiking neurons drives a read-out layer: the dynamics of the driver recurrent network is trained to encode time which is then mapped through the read-out neurons to encode another dimension, such as space or a phase. Different spatiotemporal patterns can be learned and encoded through the synaptic weights to the read-out neurons that follow common Hebbian learning rules. We demonstrate that the model is able to learn spatiotemporal dynamics on time scales that are behaviourally relevant and we show that the learned sequences are robustly replayed during a regime of spontaneous activity.

Informing antimicrobial management in the context of COVID-19: understanding the longitudinal dynamics of C-reactive protein and procalcitonin.

BACKGROUND: To characterise the longitudinal dynamics of C-reactive protein (CRP) and Procalcitonin (PCT) in a cohort of hospitalised patients with COVID-19 and support antimicrobial decision-making. METHODS: Longitudinal CRP and PCT concentrations and trajectories of 237 hospitalised patients with COVID-19 were modelled. The dataset comprised of 2,021 data points for CRP and 284 points for PCT. Pairwise comparisons were performed between: (i) those with or without significant bacterial growth from cultures, and (ii) those who survived or died in hospital. RESULTS: CRP concentrations were higher over time in COVID-19 patients with positive microbiology (day 9: 236 vs 123 mg/L, p < 0.0001) and in those who died (day 8: 226 vs 152 mg/L, p < 0.0001) but only after day 7 of COVID-related symptom onset. Failure for CRP to reduce in the first week of hospital admission was associated with significantly higher odds of death. PCT concentrations were higher in patients with COVID-19 and positive microbiology or in those who died, although these differences were not statistically significant. CONCLUSIONS: Both the absolute CRP concentration and the trajectory during the first week of hospital admission are important factors predicting microbiology culture positivity and outcome in patients hospitalised with COVID-19. Further work is needed to describe the role of PCT for co-infection. Understanding relationships of these biomarkers can support development of risk models and inform optimal antimicrobial strategies.

SC3: consensus clustering of single-cell RNA-seq data.

Single-cell RNA-seq enables the quantitative characterization of cell types based on global transcriptome profiles. We present single-cell consensus clustering (SC3), a user-friendly tool for unsupervised clustering, which achieves high accuracy and robustness by combining multiple clustering solutions through a consensus approach (http://bioconductor.org/packages/SC3). We demonstrate that SC3 is capable of identifying subclones from the transcriptomes of neoplastic cells collected from patients.

Bounding the stationary distributions of the chemical master equation via mathematical programming.

The stochastic dynamics of biochemical networks are usually modeled with the chemical master equation (CME). The stationary distributions of CMEs are seldom solvable analytically, and numerical methods typically produce estimates with uncontrolled errors. Here, we introduce mathematical programming approaches that yield approximations of these distributions with computable error bounds which enable the verification of their accuracy. First, we use semidefinite programming to compute increasingly tighter upper and lower bounds on the moments of the stationary distributions for networks with rational propensities. Second, we use these moment bounds to formulate linear programs that yield convergent upper and lower bounds on the stationary distributions themselves, their marginals, and stationary averages. The bounds obtained also provide a computational test for the uniqueness of the distribution. In the unique case, the bounds form an approximation of the stationary distribution with a computable bound on its error. In the nonunique case, our approach yields converging approximations of the ergodic distributions. We illustrate our methodology through several biochemical examples taken from the literature: Schlögl's model for a chemical bifurcation, a two-dimensional toggle switch, a model for bursty gene expression, and a dimerization model with multiple stationary distributions.

GlnK Facilitates the Dynamic Regulation of Bacterial Nitrogen Assimilation.

Ammonium assimilation in Escherichia coli is regulated by two paralogous proteins (GlnB and GlnK), which orchestrate interactions with regulators of gene expression, transport proteins, and metabolic pathways. Yet how they conjointly modulate the activity of glutamine synthetase, the key enzyme for nitrogen assimilation, is poorly understood. We combine experiments and theory to study the dynamic roles of GlnB and GlnK during nitrogen starvation and upshift. We measure time-resolved in vivo concentrations of metabolites, total and posttranslationally modified proteins, and develop a concise biochemical model of GlnB and GlnK that incorporates competition for active and allosteric sites, as well as functional sequestration of GlnK. The model predicts the responses of glutamine synthetase, GlnB, and GlnK under time-varying external ammonium level in the wild-type and two genetic knock-outs. Our results show that GlnK is tightly regulated under nitrogen-rich conditions, yet it is expressed during ammonium run-out and starvation. This suggests a role for GlnK as a buffer of nitrogen shock after starvation, and provides a further functional link between nitrogen and carbon metabolisms.

Flow-Based Network Analysis of the Caenorhabditis elegans Connectome.

We exploit flow propagation on the directed neuronal network of the nematode C. elegans to reveal dynamically relevant features of its connectome. We find flow-based groupings of neurons at different levels of granularity, which we relate to functional and anatomical constituents of its nervous system. A systematic in silico evaluation of the full set of single and double neuron ablations is used to identify deletions that induce the most severe disruptions of the multi-resolution flow structure. Such ablations are linked to functionally relevant neurons, and suggest potential candidates for further in vivo investigation. In addition, we use the directional patterns of incoming and outgoing network flows at all scales to identify flow profiles for the neurons in the connectome, without pre-imposing a priori categories. The four flow roles identified are linked to signal propagation motivated by biological input-response scenarios.

Amplification of small molecule-inducible gene expression via tuning of intracellular receptor densities.

Ligand-responsive transcription factors in prokaryotes found simple small molecule-inducible gene expression systems. These have been extensively used for regulated protein production and associated biosynthesis of fine chemicals. However, the promoter and protein engineering approaches traditionally used often pose significant restrictions to predictably and rapidly tune the expression profiles of inducible expression systems. Here, we present a new unified and rational tuning method to amplify the sensitivity and dynamic ranges of versatile small molecule-inducible expression systems. We employ a systematic variation of the concentration of intracellular receptors for transcriptional control. We show that a low density of the repressor receptor (e.g. TetR and ArsR) in the cell can significantly increase the sensitivity and dynamic range, whereas a high activator receptor (e.g. LuxR) density achieves the same outcome. The intracellular concentration of receptors can be tuned in both discrete and continuous modes by adjusting the strength of their cognate driving promoters. We exemplified this approach in several synthetic receptor-mediated sensing circuits, including a tunable cell-based arsenic sensor. The approach offers a new paradigm to predictably tune and amplify ligand-responsive gene expression with potential applications in synthetic biology and industrial biotechnology.

Emergence of Slow-Switching Assemblies in Structured Neuronal Networks.

Unraveling the interplay between connectivity and spatio-temporal dynamics in neuronal networks is a key step to advance our understanding of neuronal information processing. Here we investigate how particular features of network connectivity underpin the propensity of neural networks to generate slow-switching assembly (SSA) dynamics, i.e., sustained epochs of increased firing within assemblies of neurons which transition slowly between different assemblies throughout the network. We show that the emergence of SSA activity is linked to spectral properties of the asymmetric synaptic weight matrix. In particular, the leading eigenvalues that dictate the slow dynamics exhibit a gap with respect to the bulk of the spectrum, and the associated Schur vectors exhibit a measure of block-localization on groups of neurons, thus resulting in coherent dynamical activity on those groups. Through simple rate models, we gain analytical understanding of the origin and importance of the spectral gap, and use these insights to develop new network topologies with alternative connectivity paradigms which also display SSA activity. Specifically, SSA dynamics involving excitatory and inhibitory neurons can be achieved by modifying the connectivity patterns between both types of neurons. We also show that SSA activity can occur at multiple timescales reflecting a hierarchy in the connectivity, and demonstrate the emergence of SSA in small-world like networks. Our work provides a step towards understanding how network structure (uncovered through advancements in neuroanatomy and connectomics) can impact on spatio-temporal neural activity and constrain the resulting dynamics.

Engineering modular and tunable genetic amplifiers for scaling transcriptional signals in cascaded gene networks.

Synthetic biology aims to control and reprogram signal processing pathways within living cells so as to realize repurposed, beneficial applications. Here we report the design and construction of a set of modular and gain-tunable genetic amplifiers in Escherichia coli capable of amplifying a transcriptional signal with wide tunable-gain control in cascaded gene networks. The devices are engineered using orthogonal genetic components (hrpRS, hrpV and PhrpL) from the hrp (hypersensitive response and pathogenicity) gene regulatory network in Pseudomonas syringae. The amplifiers can linearly scale up to 21-fold the transcriptional input with a large output dynamic range, yet not introducing significant time delay or significant noise during signal amplification. The set of genetic amplifiers achieves different gains and input dynamic ranges by varying the expression levels of the underlying ligand-free activator proteins in the device. As their electronic counterparts, these engineered transcriptional amplifiers can act as fundamental building blocks in the design of biological systems by predictably and dynamically modulating transcriptional signal flows to implement advanced intra- and extra-cellular control functions.

Graph partitions and cluster synchronization in networks of oscillators.

Synchronization over networks depends strongly on the structure of the coupling between the oscillators. When the coupling presents certain regularities, the dynamics can be coarse-grained into clusters by means of External Equitable Partitions of the network graph and their associated quotient graphs. We exploit this graph-theoretical concept to study the phenomenon of cluster synchronization, in which different groups of nodes converge to distinct behaviors. We derive conditions and properties of networks in which such clustered behavior emerges and show that the ensuing dynamics is the result of the localization of the eigenvectors of the associated graph Laplacians linked to the existence of invariant subspaces. The framework is applied to both linear and non-linear models, first for the standard case of networks with positive edges, before being generalized to the case of signed networks with both positive and negative interactions. We illustrate our results with examples of both signed and unsigned graphs for consensus dynamics and for partial synchronization of oscillator networks under the master stability function as well as Kuramoto oscillators.

Hhex and Cer1 mediate the Sox17 pathway for cardiac mesoderm formation in embryonic stem cells.

Cardiac muscle differentiation in vivo is guided by sequential growth factor signals, including endoderm-derived diffusible factors, impinging on cardiogenic genes in the developing mesoderm. Previously, by RNA interference in AB2.2 mouse embryonic stem cells (mESCs), we identified the endodermal transcription factor Sox17 as essential for Mesp1 induction in primitive mesoderm and subsequent cardiac muscle differentiation. However, downstream effectors of Sox17 remained to be proven functionally. In this study, we used genome-wide profiling of Sox17-dependent genes in AB2.2 cells, RNA interference, chromatin immunoprecipitation, and luciferase reporter genes to dissect this pathway. Sox17 was required not only for Hhex (a second endodermal transcription factor) but also for Cer1, a growth factor inhibitor from endoderm that, like Hhex, controls mesoderm patterning in Xenopus toward a cardiac fate. Suppressing Hhex or Cer1 blocked cardiac myogenesis, although at a later stage than induction of Mesp1/2. Hhex was required but not sufficient for Cer1 expression. Over-expression of Sox17 induced endogenous Cer1 and sequence-specific transcription of a Cer1 reporter gene. Forced expression of Cer1 was sufficient to rescue cardiac differentiation in Hhex-deficient cells. Thus, Hhex and Cer1 are indispensable components of the Sox17 pathway for cardiopoiesis in mESCs, acting at a stage downstream from Mesp1/2.

Introduction of a fluorescent probe to amyloid-ß to reveal kinetic insights into its interactions with copper(II).

The kinetics of the interactions between amyloid-ß (Aß) and metal ions are crucial to understanding the physiological and pathological roles of Aß in the normal brain and in Alzheimer's disease. Using the quenching of a fluorescent probe by Cu(2+), the mechanism of Aß/Cu(2+) interactions in physiologically relevant conditions has been elucidated. Cu(2+) binds to Aß at a near diffusion-limited rate, initially forming component I. The switching between component I and II occurs on the second timescale, with a significant energy barrier. Component I is much more reactive towards Cu(2+) ligands and likely responsible for initial Aß dimer formation. Clioquinol (CQ) is shown to sequester Cu(2+) more effectively than other tested ligands. These findings have implications for the potential roles of Aß in regulating neurotransmission, and for the screening of small molecules targeting Aß-metal interactions.