Targeted single-cell RNA sequencing of transcription factors enhances the identification of cell types and trajectories.

Single-cell RNA sequencing (scRNA-seq) is a widely used method for identifying cell types and trajectories in biologically heterogeneous samples, but it is limited in its detection and quantification of lowly expressed genes. This results in missing important biological signals, such as the expression of key transcription factors (TFs) driving cellular differentiation. We show that targeted sequencing of ∼1000 TFs (scCapture-seq) in iPSC-derived neuronal cultures greatly improves the biological information garnered from scRNA-seq. Increased TF resolution enhanced cell type identification, developmental trajectories, and gene regulatory networks. This allowed us to resolve differences among neuronal populations, which were generated in two different laboratories using the same differentiation protocol. ScCapture-seq improved TF-gene regulatory network inference and thus identified divergent patterns of neurogenesis into either excitatory cortical neurons or inhibitory interneurons. Furthermore, scCapture-seq revealed a role for of retinoic acid signaling in the developmental divergence between these different neuronal populations. Our results show that TF targeting improves the characterization of human cellular models and allows identification of the essential differences between cellular populations, which would otherwise be missed in traditional scRNA-seq. scCapture-seq TF targeting represents a cost-effective enhancement of scRNA-seq, which could be broadly applied to improve scRNA-seq resolution.

Global proteomic analysis of extracellular matrix in mouse and human brain highlights relevance to cerebrovascular disease.

The extracellular matrix (ECM) is a key interface between the cerebrovasculature and adjacent brain tissues. Deregulation of the ECM contributes to a broad range of neurological disorders. However, despite this importance, our understanding of the ECM composition remains very limited mainly due to difficulties in its isolation. To address this, we developed an approach to extract the cerebrovascular ECM from mouse and human post-mortem normal brain tissues. We then used mass spectrometry with off-line high-pH reversed-phase fractionation to increase the protein detection. This identified more than 1000 proteins in the ECM-enriched fraction, with > 66% of the proteins being common between the species. We report 147 core ECM proteins of the human brain vascular matrisome, including collagens, laminins, fibronectin and nidogens. We next used network analysis to identify the connection between the brain ECM proteins and cerebrovascular diseases. We found that genes related to cerebrovascular diseases, such as COL4A1, COL4A2, VCAN and APOE were significantly enriched in the cerebrovascular ECM network. This provides unique mechanistic insight into cerebrovascular disease and potential drug targets. Overall, we provide a powerful resource to study the functions of brain ECM and highlight a specific role for brain vascular ECM in cerebral vascular disease.

Common transcriptional signatures of neuropathic pain.

ABSTARCT: The dorsal root ganglia (DRG) are key structures in nociception and chronic pain disorders. Several gene expression studies of DRG in preclinical pain models have been performed, but it is unclear if consistent gene changes are identifiable. We, therefore, compared several recent RNA-Seq data sets on the whole DRG in rodent models of nerve injury. Contrary to previous findings, we show hundreds of common differentially expressed genes and high positive correlation between studies, despite model and species differences. We also find, in contrast to previous studies, that 60% of the common rodent gene response after injury is likely to occur in nociceptors of the DRG. Substantial expression changes are observed at a 1-week time-point, with smaller changes in the same genes at a later 3- to 4-week time-point. However, a subset of genes shows a similar magnitude of changes at both early and late time-points, suggesting their potential involvement in the maintenance of chronic pain. These genes are centred around suppression of endogenous opioid signalling. Reversal of this suppression could allow endogenous and exogenous opioids to exert their analgesic functions and may be an important strategy for treating chronic pain disorders. Currently used drugs, such as amitriptyline and duloxetine, do not seem to appropriately modulate many of the critical pain genes and indeed may transcriptionally suppress endogenous opioid signalling further.

A single-input binary counting module based on serine integrase site-specific recombination.

A device that counts and records the number of events experienced by an individual cell could have many uses in experimental biology and biotechnology. Here, we report a DNA-based 'latch' that switches between two states upon each exposure to a repeated stimulus. The key component of the latch is a DNA segment whose orientation is inverted by the actions of ϕC31 integrase and its recombination directionality factor (RDF). Integrase expression is regulated by an external input, while RDF expression is controlled by the state of the latch, such that the orientation of the invertible segment switches efficiently each time the device receives an input pulse. Recombination occurs over a time scale of minutes after initiation of integrase expression. The latch requires a delay circuit, implemented with a transcriptional repressor expressed in only one state, to ensure that each input pulse results in only one inversion of the DNA segment. Development and optimization of the latch in living cells was driven by mathematical modelling of the recombination reactions and gene expression regulated by the switch. We discuss how N latches built with orthogonal site-specific recombination systems could be chained together to form a binary ripple counter that could count to 2N - 1.

Mathematical model of a serine integrase-controlled toggle switch with a single input.

Dual-state genetic switches that can change their state in response to input signals can be used in synthetic biology to encode memory and control gene expression. A transcriptional toggle switch (TTS), with two mutually repressing transcription regulators, was previously used for switching between two expression states. In other studies, serine integrases have been used to control DNA inversion switches that can alternate between two different states. Both of these switches use two different inputs to switch ON or OFF. Here, we use mathematical modelling to design a robust one-input binary switch, which combines a TTS with a DNA inversion switch. This combined circuit switches between the two states every time it receives a pulse of a single-input signal. The robustness of the switch is based on the bistability of its TTS, while integrase recombination allows single-input control. Unidirectional integrase-RDF-mediated recombination is provided by a recently developed integrase-RDF fusion protein. We show that the switch is stable against parameter variations and molecular noise, making it a promising candidate for further use as a basic element of binary counting devices.

Monitoring of nutrient limitation in growing E. coli: a mathematical model of a ppGpp-based biosensor.

BACKGROUND: E. coli can be used as bacterial cell factories for production of biofuels and other useful compounds. The efficient production of the desired products requires careful monitoring of growth conditions and the optimization of metabolic fluxes. To avoid nutrient depletion and maximize product yields we suggest using a natural mechanism for sensing nutrient limitation, related to biosynthesis of an intracellular messenger - guanosine tetraphosphate (ppGpp). RESULTS: We propose a design for a biosensor, which monitors changes in the intracellular concentration of ppGpp by coupling it to a fluorescent output. We used mathematical modelling to analyse the intracellular dynamics of ppGpp, its fluorescent reporter, and cell growth in normal and fatty acid-producing E. coli lines. The model integrates existing mechanisms of ppGpp regulation and predicts the biosensor response to changes in nutrient state. In particular, the model predicts that excessive stimulation of fatty acid production depletes fatty acid intermediates, downregulates growth and increases the levels of ppGpp-related fluorescence. CONCLUSIONS: Our analysis demonstrates that the ppGpp sensor can be used for early detection of nutrient limitation during cell growth and for testing productivity of engineered lines.

A simplified mathematical model of directional DNA site-specific recombination by serine integrases.

Serine integrases catalyse site-specific recombination to integrate and excise bacteriophage genomes into and out of their host's genome. These enzymes exhibit remarkable directionality; in the presence of the integrase alone, recombination between attP and attB DNA sites is efficient and irreversible, giving attL and attR products which do not recombine further. However, in the presence of the bacteriophage-encoded recombination directionality factor (RDF), integrase efficiently promotes recombination between attL and attR to re-form attP and attB The DNA substrates and products of both reactions are approximately isoenergetic, and no cofactors (such as adenosine triphosphate) are required for recombination. The thermodynamic driving force for directionality of these reactions is thus enigmatic. Here, we present a minimal mathematical model which can explain the directionality and regulation of both 'forward' and 'reverse' reactions. In this model, the substrates of the 'forbidden' reactions (between attL and attR in the absence of RDF, attP and attB in the presence of RDF) are trapped as inactive protein-DNA complexes, ensuring that these 'forbidden' reactions are extremely slow. The model is in good agreement with the observed in vitro kinetics of recombination by ϕC31 integrase, and defines core features of the system necessary and sufficient for directionality.

The mechanism of ϕC31 integrase directionality: experimental analysis and computational modelling.

Serine integrases, DNA site-specific recombinases used by bacteriophages for integration and excision of their DNA to and from their host genomes, are increasingly being used as tools for programmed rearrangements of DNA molecules for biotechnology and synthetic biology. A useful feature of serine integrases is the simple regulation and unidirectionality of their reactions. Recombination between the phage attP and host attB sites is promoted by the serine integrase alone, giving recombinant attL and attR sites, whereas the 'reverse' reaction (between attL and attR) requires an additional protein, the recombination directionality factor (RDF). Here, we present new experimental data on the kinetics and regulation of recombination reactions mediated by ϕC31 integrase and its RDF, and use these data as the basis for a mathematical model of the reactions. The model accounts for the unidirectionality of the attP × attB and attL × attR reactions by hypothesizing the formation of structurally distinct, kinetically stable integrase-DNA product complexes, dependent on the presence or absence of RDF. The model accounts for all the available experimental data, and predicts how mutations of the proteins or alterations of reaction conditions might increase the conversion efficiency of recombination.

Defining the robust behaviour of the plant clock gene circuit with absolute RNA timeseries and open infrastructure.

Our understanding of the complex, transcriptional feedback loops in the circadian clock mechanism has depended upon quantitative, timeseries data from disparate sources. We measure clock gene RNA profiles in Arabidopsis thaliana seedlings, grown with or without exogenous sucrose, or in soil-grown plants and in wild-type and mutant backgrounds. The RNA profiles were strikingly robust across the experimental conditions, so current mathematical models are likely to be broadly applicable in leaf tissue. In addition to providing reference data, unexpected behaviours included co-expression of PRR9 and ELF4, and regulation of PRR5 by GI. Absolute RNA quantification revealed low levels of PRR9 transcripts (peak approx. 50 copies cell(-1)) compared with other clock genes, and threefold higher levels of LHY RNA (more than 1500 copies cell(-1)) than of its close relative CCA1. The data are disseminated from BioDare, an online repository for focused timeseries data, which is expected to benefit mechanistic modelling. One data subset successfully constrained clock gene expression in a complex model, using publicly available software on parallel computers, without expert tuning or programming. We outline the empirical and mathematical justification for data aggregation in understanding highly interconnected, dynamic networks such as the clock, and the observed design constraints on the resources required to make this approach widely accessible.

Mathematical modelling of the diurnal regulation of the MEP pathway in Arabidopsis.

Isoprenoid molecules are essential elements of plant metabolism. Many important plant isoprenoids, such as chlorophylls, carotenoids, tocopherols, prenylated quinones and hormones are synthesised in chloroplasts via the 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway. Here we develop a mathematical model of diurnal regulation of the MEP pathway in Arabidopsis thaliana. We used both experimental and theoretical approaches to integrate mechanisms potentially involved in the diurnal control of the pathway. Our data show that flux through the MEP pathway is accelerated in light due to the photosynthesis-dependent supply of metabolic substrates of the pathway and the transcriptional regulation of key biosynthetic genes by the circadian clock. We also demonstrate that feedback regulation of both the activity and the abundance of the first enzyme of the MEP pathway (1-deoxy-D-xylulose 5-phosphate synthase, DXS) by pathway products stabilizes the flux against changes in substrate supply and adjusts the flux according to product demand under normal growth conditions. These data illustrate the central relevance of photosynthesis, the circadian clock and feedback control of DXS for the diurnal regulation of the MEP pathway.

Mathematical modelling of diurnal regulation of carbohydrate allocation by osmo-related processes in plants.

Plants synthesize sucrose in source tissues (mainly mature leafs) and supply it for growth of sink tissues (young leafs). Sucrose is derived from photosynthesis during daytime and from starch at night. Because the diurnal regulation of sucrose fluxes is not completely understood, we built a mathematical model designed to reproduce all key experimental observations. For this, assumptions were made about the molecular mechanisms underlying the regulations, which are all motivated by experimental facts. The key regulators in our model are two kinases (SnRK1 and osmo-sensitive kinase OsmK) under the control of the circadian clock. SnRK1 is activated in the night to prepare for regularly occurring carbon-limiting conditions, whereas OsmK is activated during the day to prepare for water deficit, which often occurs in the afternoon. Decrease of SnRK1 and increase of OsmK result in partitioning of carbon towards sucrose to supply growing sink tissues. Concomitantly, increasing levels of the growth regulator trehalose-6-phosphate stimulates the development of new sink tissues and thus sink demand, which further activates sucrose supply in a positive feedback loop. We propose that OsmK acts as a timer to measure the length of the photoperiod and suggest experiments how this hypothesis can be validated.

Adjustment of carbon fluxes to light conditions regulates the daily turnover of starch in plants: a computational model.

In the light, photosynthesis provides carbon for metabolism and growth. In the dark, plant growth depends on carbon reserves that were accumulated during previous light periods. Many plants accumulate part of their newly-fixed carbon as starch in their leaves in the day and remobilise it to support metabolism and growth at night. The daily rhythms of starch accumulation and degradation are dynamically adjusted to the changing light conditions such that starch is almost but not totally exhausted at dawn. This requires the allocation of a larger proportion of the newly fixed carbon to starch under low carbon conditions, and the use of information about the carbon status at the end of the light period and the length of the night to pace the rate of starch degradation. This regulation occurs in a circadian clock-dependent manner, through unknown mechanisms. We use mathematical modelling to explore possible diurnal mechanisms regulating the starch level. Our model combines the main reactions of carbon fixation, starch and sucrose synthesis, starch degradation and consumption of carbon by sink tissues. To describe the dynamic adjustment of starch to daily conditions, we introduce diurnal regulators of carbon fluxes, which modulate the activities of the key steps of starch metabolism. The sensing of the diurnal conditions is mediated in our model by the timer α and the "dark sensor"ß, which integrate daily information about the light conditions and time of the day through the circadian clock. Our data identify the ß subunit of SnRK1 kinase as a good candidate for the role of the dark-accumulated component ß of our model. The developed novel approach for understanding starch kinetics through diurnal metabolic and circadian sensors allowed us to explain starch time-courses in plants and predict the kinetics of the proposed diurnal regulators under various genetic and environmental perturbations.

Regulatory principles and experimental approaches to the circadian control of starch turnover.

In many plants, starch is synthesized during the day and degraded during the night to avoid carbohydrate starvation in darkness. The circadian clock participates in a dynamic adjustment of starch turnover to changing environmental condition through unknown mechanisms. We used mathematical modelling to explore the possible scenarios for the control of starch turnover by the molecular components of the plant circadian clock. Several classes of plausible models were capable of describing the starch dynamics observed in a range of clock mutant plants and light conditions, including discriminating circadian protocols. Three example models of these classes are studied in detail, differing in several important ways. First, the clock components directly responsible for regulating starch degradation are different in each model. Second, the intermediate species in the pathway may play either an activating or inhibiting role on starch degradation. Third, the system may include a light-dependent interaction between the clock and downstream processes. Finally, the clock may be involved in the regulation of starch synthesis. We discuss the differences among the models' predictions for diel starch profiles and the properties of the circadian regulators. These suggest additional experiments to elucidate the pathway structure, avoid confounding results and identify the molecular components involved.

Modelling the widespread effects of TOC1 signalling on the plant circadian clock and its outputs.

BACKGROUND: 24-hour biological clocks are intimately connected to the cellular signalling network, which complicates the analysis of clock mechanisms. The transcriptional regulator TOC1 (TIMING OF CAB EXPRESSION 1) is a founding component of the gene circuit in the plant circadian clock. Recent results show that TOC1 suppresses transcription of multiple target genes within the clock circuit, far beyond its previously-described regulation of the morning transcription factors LHY (LATE ELONGATED HYPOCOTYL) and CCA1 (CIRCADIAN CLOCK ASSOCIATED 1). It is unclear how this pervasive effect of TOC1 affects the dynamics of the clock and its outputs. TOC1 also appears to function in a nested feedback loop that includes signalling by the plant hormone Abscisic Acid (ABA), which is upregulated by abiotic stresses, such as drought. ABA treatments both alter TOC1 levels and affect the clock's timing behaviour. Conversely, the clock rhythmically modulates physiological processes induced by ABA, such as the closing of stomata in the leaf epidermis. In order to understand the dynamics of the clock and its outputs under changing environmental conditions, the reciprocal interactions between the clock and other signalling pathways must be integrated. RESULTS: We extended the mathematical model of the plant clock gene circuit by incorporating the repression of multiple clock genes by TOC1, observed experimentally. The revised model more accurately matches the data on the clock's molecular profiles and timing behaviour, explaining the clock's responses in TOC1 over-expression and toc1 mutant plants. A simplified representation of ABA signalling allowed us to investigate the interactions of ABA and circadian pathways. Increased ABA levels lengthen the free-running period of the clock, consistent with the experimental data. Adding stomatal closure to the model, as a key ABA- and clock-regulated downstream process allowed to describe TOC1 effects on the rhythmic gating of stomatal closure. CONCLUSIONS: The integrated model of the circadian clock circuit and ABA-regulated environmental sensing allowed us to explain multiple experimental observations on the timing and stomatal responses to genetic and environmental perturbations. These results crystallise a new role of TOC1 as an environmental sensor, which both affects the pace of the central oscillator and modulates the kinetics of downstream processes.

The clock gene circuit in Arabidopsis includes a repressilator with additional feedback loops.

Circadian clocks synchronise biological processes with the day/night cycle, using molecular mechanisms that include interlocked, transcriptional feedback loops. Recent experiments identified the evening complex (EC) as a repressor that can be essential for gene expression rhythms in plants. Integrating the EC components in this role significantly alters our mechanistic, mathematical model of the clock gene circuit. Negative autoregulation of the EC genes constitutes the clock's evening loop, replacing the hypothetical component Y. The EC explains our earlier conjecture that the morning gene Pseudo-Response Regulator 9 was repressed by an evening gene, previously identified with Timing Of CAB Expression1 (TOC1). Our computational analysis suggests that TOC1 is a repressor of the morning genes Late Elongated Hypocotyl and Circadian Clock Associated1 rather than an activator as first conceived. This removes the necessity for the unknown component X (or TOC1mod) from previous clock models. As well as matching timeseries and phase-response data, the model provides a new conceptual framework for the plant clock that includes a three-component repressilator circuit in its complex structure.

Stochastic properties of the plant circadian clock.

Circadian clocks are gene regulatory networks whose role is to help the organisms to cope with variations in environmental conditions such as the day/night cycle. In this work, we explored the effects of molecular noise in single cells on the behaviour of the circadian clock in the plant model species Arabidopsis thaliana. The computational modelling language Bio-PEPA enabled us to give a stochastic interpretation of an existing deterministic model of the clock, and to easily compare the results obtained via stochastic simulation and via numerical solution of the deterministic model. First, the introduction of stochasticity in the model allowed us to estimate the unknown size of the system. Moreover, stochasticity improved the description of the available experimental data in several light conditions: noise-induced fluctuations yield a faster entrainment of the plant clock under certain photoperiods and are able to explain the experimentally observed dampening of the oscillations in plants under constant light conditions. The model predicts that the desynchronization between noisy oscillations in single cells contributes to the observed damped oscillations at the level of the cell population. Analysis of the phase, period and amplitude distributions under various light conditions demonstrated robust entrainment of the plant clock to light/dark cycles which closely matched the available experimental data.

Temporal repression of core circadian genes is mediated through EARLY FLOWERING 3 in Arabidopsis.

The circadian clock provides robust, ∼24 hr biological rhythms throughout the eukaryotes. The clock gene circuit in plants comprises interlocking transcriptional feedback loops, reviewed in [1], whereby the morning-expressed transcription factors CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY) repress the expression of evening genes, notably TIMING OF CAB EXPRESSION 1 (TOC1). EARLY FLOWERING 3 (ELF3) has been implicated as a repressor of light signaling to the clock [2, 3] and, paradoxically, as an activator of the light-induced genes CCA1 and LHY [4, 5]. We use cca1-11 lhy-21 elf3-4 plants to separate the repressive function of ELF3 from its downstream targets CCA1 and LHY. We further demonstrate that ELF3 associates physically with the promoter of PSEUDO-RESPONSE REGULATOR 9 (PRR9), a repressor of CCA1 and LHY expression, in a time-dependent fashion. The repressive function of ELF3 is thus consistent with indirect activation of LHY and CCA1, in a double-negative connection via a direct ELF3 target, PRR9. This mechanism reconciles the functions of ELF3 in the clock network during the night and points to further effects of ELF3 during the day.

Ubiquitin ligase switch in plant photomorphogenesis: A hypothesis.

The E3 ubiquitin ligase COP1 (CONSTITUTIVE PHOTOMORPHOGENIC1) plays a key role in the repression of the plant photomorphogenic development in darkness. In the presence of light, COP1 is inactivated by a mechanism which is not completely understood. This leads to accumulation of COP1's target transcription factors, which initiates photomorphogenesis, resulting in dramatic changes of the seedling's physiology. Here we use a mathematical model to explore the possible mechanism of COP1 modulation upon dark/light transition in Arabidopsis thaliana based upon data for two COP1 target proteins: HY5 and HFR1, which play critical roles in photomorphogenesis. The main reactions in our model are the inactivation of COP1 by a proposed photoreceptor-related inhibitor I and interactions between COP1 and a CUL4 (CULLIN4)-based ligase. For building and verification of the model, we used the available published and our new data on the kinetics of HY5 and HFR1 together with the data on COP1 abundance. HY5 has been shown to accumulate at a slower rate than HFR1. To describe the observed differences in the timecourses of the "slow" target HY5 and the "fast" target HFR1, we hypothesize a switch between the activities of COP1 and CUL4 ligases upon dark/light transition, with COP1 being active mostly in darkness and CUL4 in light. The model predicts a bi-phasic kinetics of COP1 activity upon the exposure of plants to light, with its restoration after the initial decline and the following slow depletion of the total COP1 content. CUL4 activity is predicted to increase in the presence of light. We propose that the ubiquitin ligase switch is important for the complex regulation of multiple transcription factors during plants development. In addition, this provides a new mechanism for sensing the duration of light period, which is important for seasonal changes in plant development.

Quantitative analysis of regulatory flexibility under changing environmental conditions.

The circadian clock controls 24-h rhythms in many biological processes, allowing appropriate timing of biological rhythms relative to dawn and dusk. Known clock circuits include multiple, interlocked feedback loops. Theory suggested that multiple loops contribute the flexibility for molecular rhythms to track multiple phases of the external cycle. Clear dawn- and dusk-tracking rhythms illustrate the flexibility of timing in Ipomoea nil. Molecular clock components in Arabidopsis thaliana showed complex, photoperiod-dependent regulation, which was analysed by comparison with three contrasting models. A simple, quantitative measure, Dusk Sensitivity, was introduced to compare the behaviour of clock models with varying loop complexity. Evening-expressed clock genes showed photoperiod-dependent dusk sensitivity, as predicted by the three-loop model, whereas the one- and two-loop models tracked dawn and dusk, respectively. Output genes for starch degradation achieved dusk-tracking expression through light regulation, rather than a dusk-tracking rhythm. Model analysis predicted which biochemical processes could be manipulated to extend dusk tracking. Our results reveal how an operating principle of biological regulators applies specifically to the plant circadian clock.

Data assimilation constrains new connections and components in a complex, eukaryotic circadian clock model.

Circadian clocks generate 24-h rhythms that are entrained by the day/night cycle. Clock circuits include several light inputs and interlocked feedback loops, with complex dynamics. Multiple biological components can contribute to each part of the circuit in higher organisms. Mechanistic models with morning, evening and central feedback loops have provided a heuristic framework for the clock in plants, but were based on transcriptional control. Here, we model observed, post-transcriptional and post-translational regulation and constrain many parameter values based on experimental data. The model's feedback circuit is revised and now includes PSEUDO-RESPONSE REGULATOR 7 (PRR7) and ZEITLUPE. The revised model matches data in varying environments and mutants, and gains robustness to parameter variation. Our results suggest that the activation of important morning-expressed genes follows their release from a night inhibitor (NI). Experiments inspired by the new model support the predicted NI function and show that the PRR5 gene contributes to the NI. The multiple PRR genes of Arabidopsis uncouple events in the late night from light-driven responses in the day, increasing the flexibility of rhythmic regulation.