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.

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.

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.

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.

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.

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.

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.

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.

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.

Dynamics of Cdc42 network embodies a Turing-type mechanism of yeast cell polarity.

Complex biochemical networks can be understood by identifying their principal regulatory motifs and mode of action. We model the early phase of budding yeast cellular polarization and show that the biochemical processes in the presumptive bud site comprise a Turing-type mechanism. The roles of the prototypical activator and substrate are played by GTPase Cdc42 in its active and inactive states, respectively. We demonstrate that the nucleotide cycling of Cdc42 converts cellular energy into a stable cluster of activated Cdc42. This energy drives a continuous membrane-cytoplasmic exchange of the cluster components to counteract diffusive spread of the cluster. This exchange explains why only one bud forms per cell cycle, because the winner-takes-all competition of candidate sites inevitably selects a single site.

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.

Computational model explains high activity and rapid cycling of Rho GTPases within protein complexes.

Formation of multiprotein complexes on cellular membranes is critically dependent on the cyclic activation of small GTPases. FRAP-based analyses demonstrate that within protein complexes, some small GTPases cycle nearly three orders of magnitude faster than they would spontaneously cycle in vitro. At the same time, experiments report concomitant excess of the activated, GTP-bound form of GTPases over their inactive form. Intuitively, high activity and rapid turnover are contradictory requirements. How the cells manage to maximize both remains poorly understood. Here, using GTPases of the Rab and Rho families as a prototype, we introduce a computational model of the GTPase cycle. We quantitatively investigate several plausible layouts of the cycling control module that consist of GEFs, GAPs, and GTPase effectors. We explain the existing experimental data and predict how the cycling of GTPases is controlled by the regulatory proteins in vivo. Our model explains distinct and separable roles that the activating GEFs and deactivating GAPs play in the GTPase cycling control. While the activity of GTPase is mainly defined by GEF, the turnover rate is a sole function of GAP. Maximization of the GTPase activity and turnover rate places conflicting requirements on the concentration of GAP. Therefore, to achieve a high activity and turnover rate at once, cells must carefully maintain concentrations of GEFs and GAPs within the optimal range. The values of these optimal concentrations indicate that efficient cycling can be achieved only within dense protein complexes typically assembled on the membrane surfaces. We show that the concentration requirement for GEF can be dramatically reduced by a GEF-activating GTPase effector that can also significantly boost the cycling efficiency. Interestingly, we find that the cycling regimes are only weakly dependent on the concentration of GTPase itself.

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.

Spatiotemporal dynamics of contact activation factors of blood coagulation.

A new in vitro model is proposed for studying the spatiotemporal distributions of activated clotting factors, in which clotting is activated in a thin layer of non-stirred plasma supplemented with a fluorogenic substrate and is monitored by fluorescence from its cleavage product. Analysis of the spatiotemporal dynamics of factor XIa and kallikrein in glass-activated human plasma provides evidence that both contact factors remain restricted to the glass surface and possibly a narrow boundary zone (<0.1 mm). The kinetics of factor XIa and kallikrein studied by a new method (in non-stirred plasma) coincided with those studied fluorimetrically with full stirring: their concentrations rapidly rose for the first few minutes after activation and then slowly declined. Factor XI and prekallikrein activation is likely to be restricted by the limited number of sites available for binding to the surface. The maximum concentration of the active factors was estimated at 2 x 10(8) molecules per mm(2) at the glass surface (irrespective of stirring). At the plastic surface, this value was 15--30 times lower.

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.

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.

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.

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.