The σB alternative sigma factor circuit modulates noise to generate different types of pulsing dynamics.

Single-cell approaches are revealing a high degree of heterogeneity, or noise, in gene expression in isogenic bacteria. How gene circuits modulate this noise in gene expression to generate robust output dynamics is unclear. Here we use the Bacillus subtilis alternative sigma factor σB as a model system for understanding the role of noise in generating circuit output dynamics. σB controls the general stress response in B. subtilis and is activated by a range of energy and environmental stresses. Recent single-cell studies have revealed that the circuit can generate two distinct outputs, stochastic pulsing and a single pulse response, but the conditions under which each response is generated are under debate. We implement a stochastic mathematical model of the σB circuit to investigate this and find that the system's core circuit can generate both response types. This is despite one response (stochastic pulsing) being stochastic in nature, and the other (single response pulse) being deterministic. We demonstrate that the main determinant for whichever response is generated is the degree with which the input pathway activates the core circuit, although the noise properties of the input pathway also biases the system towards one or the other type of output. Thus, our work shows how stochastic modelling can reveal the mechanisms behind non-intuitive gene circuit output dynamics.

Prediction of photoperiodic regulators from quantitative gene circuit models.

Photoperiod sensors allow physiological adaptation to the changing seasons. The prevalent hypothesis is that day length perception is mediated through coupling of an endogenous rhythm with an external light signal. Sufficient molecular data are available to test this quantitatively in plants, though not yet in mammals. In Arabidopsis, the clock-regulated genes CONSTANS (CO) and FLAVIN, KELCH, F-BOX (FKF1) and their light-sensitive proteins are thought to form an external coincidence sensor. Here, we model the integration of light and timing information by CO, its target gene FLOWERING LOCUS T (FT), and the circadian clock. Among other predictions, our models show that FKF1 activates FT. We demonstrate experimentally that this effect is independent of the known activation of CO by FKF1, thus we locate a major, novel controller of photoperiodism. External coincidence is part of a complex photoperiod sensor: modeling makes this complexity explicit and may thus contribute to crop improvement.

A spatial model of the plant circadian clock reveals design principles for coordinated timing.

Individual plant cells possess a genetic network, the circadian clock, that times internal processes to the day-night cycle. Mathematical models of the clock are typically either "whole-plant" that ignore tissue or cell type-specific clock behavior, or "phase-only" that do not include molecular components. To address the complex spatial coordination observed in experiments, here we implemented a clock network model on a template of a seedling. In our model, the sensitivity to light varies across the plant, and cells communicate their timing via local or long-distance sharing of clock components, causing their rhythms to couple. We found that both varied light sensitivity and long-distance coupling could generate period differences between organs, while local coupling was required to generate the spatial waves of clock gene expression observed experimentally. We then examined our model under noisy light-dark cycles and found that local coupling minimized timing errors caused by the noise while allowing each plant region to maintain a different clock phase. Thus, local sensitivity to environmental inputs combined with local coupling enables flexible yet robust circadian timing.

Spatially specific mechanisms and functions of the plant circadian clock.

Like many organisms, plants have evolved a genetic network, the circadian clock, to coordinate processes with day/night cycles. In plants, the clock is a pervasive regulator of development and modulates many aspects of physiology. Clock-regulated processes range from the correct timing of growth and cell division to interactions with the root microbiome. Recently developed techniques, such as single-cell time-lapse microscopy and single-cell RNA-seq, are beginning to revolutionize our understanding of this clock regulation, revealing a surprising degree of organ, tissue, and cell-type specificity. In this review, we highlight recent advances in our spatial view of the clock across the plant, both in terms of how it is regulated and how it regulates a diversity of output processes. We outline how understanding these spatially specific functions will help reveal the range of ways that the clock provides a fitness benefit for the plant.

Tunable phenotypic variability through an autoregulatory alternative sigma factor circuit.

Genetically identical individuals in bacterial populations can display significant phenotypic variability. This variability can be functional, for example by allowing a fraction of stress prepared cells to survive an otherwise lethal stress. The optimal fraction of stress prepared cells depends on environmental conditions. However, how bacterial populations modulate their level of phenotypic variability remains unclear. Here we show that the alternative sigma factor σV circuit in Bacillus subtilis generates functional phenotypic variability that can be tuned by stress level, environmental history and genetic perturbations. Using single-cell time-lapse microscopy and microfluidics, we find the fraction of cells that immediately activate σV under lysozyme stress depends on stress level and on a transcriptional memory of previous stress. Iteration between model and experiment reveals that this tunability can be explained by the autoregulatory feedback structure of the sigV operon. As predicted by the model, genetic perturbations to the operon also modulate the response variability. The conserved sigma-anti-sigma autoregulation motif is thus a simple mechanism for bacterial populations to modulate their heterogeneity based on their environment.

Coordinated circadian timing through the integration of local inputs in Arabidopsis thaliana.

Individual plant cells have a genetic circuit, the circadian clock, that times key processes to the day-night cycle. These clocks are aligned to the day-night cycle by multiple environmental signals that vary across the plant. How does the plant integrate clock rhythms, both within and between organs, to ensure coordinated timing? To address this question, we examined the clock at the sub-tissue level across Arabidopsis thaliana seedlings under multiple environmental conditions and genetic backgrounds. Our results show that the clock runs at different speeds (periods) in each organ, which causes the clock to peak at different times across the plant in both constant environmental conditions and light-dark (LD) cycles. Closer examination reveals that spatial waves of clock gene expression propagate both within and between organs. Using a combination of modeling and experiment, we reveal that these spatial waves are the result of the period differences between organs and local coupling, rather than long-distance signaling. With further experiments we show that the endogenous period differences, and thus the spatial waves, can be generated by the organ specificity of inputs into the clock. We demonstrate this by modulating periods using light and metabolic signals, as well as with genetic perturbations. Our results reveal that plant clocks can be set locally by organ-specific inputs but coordinated globally via spatial waves of clock gene expression.

Cell size control driven by the circadian clock and environment in cyanobacteria.

How cells maintain their size has been extensively studied under constant conditions. In the wild, however, cells rarely experience constant environments. Here, we examine how the 24-h circadian clock and environmental cycles modulate cell size control and division timings in the cyanobacterium Synechococcus elongatus using single-cell time-lapse microscopy. Under constant light, wild-type cells follow an apparent sizer-like principle. Closer inspection reveals that the clock generates two subpopulations, with cells born in the subjective day following different division rules from cells born in subjective night. A stochastic model explains how this behavior emerges from the interaction of cell size control with the clock. We demonstrate that the clock continuously modulates the probability of cell division throughout day and night, rather than solely applying an on-off gate to division, as previously proposed. Iterating between modeling and experiments, we go on to identify an effective coupling of the division rate to time of day through the combined effects of the environment and the clock on cell division. Under naturally graded light-dark cycles, this coupling narrows the time window of cell divisions and shifts divisions away from when light levels are low and cell growth is reduced. Our analysis allows us to disentangle, and predict the effects of, the complex interactions between the environment, clock, and cell size control.

Nitrate modulates stem cell dynamics in Arabidopsis shoot meristems through cytokinins.

The shoot apical meristem (SAM) is responsible for the generation of all the aerial parts of plants. Given its critical role, dynamical changes in SAM activity should play a central role in the adaptation of plant architecture to the environment. Using quantitative microscopy, grafting experiments, and genetic perturbations, we connect the plant environment to the SAM by describing the molecular mechanism by which cytokinins signal the level of nutrient availability to the SAM. We show that a systemic signal of cytokinin precursors mediates the adaptation of SAM size and organogenesis rate to the availability of mineral nutrients by modulating the expression of WUSCHEL, a key regulator of stem cell homeostasis. In time-lapse experiments, we further show that this mechanism allows meristems to adapt to rapid changes in nitrate concentration, and thereby modulate their rate of organ production to the availability of mineral nutrients within a few days. Our work sheds light on the role of the stem cell regulatory network by showing that it not only maintains meristem homeostasis but also allows plants to adapt to rapid changes in the environment.

Widespread inter-individual gene expression variability in Arabidopsis thaliana.

A fundamental question in biology is how gene expression is regulated to give rise to a phenotype. However, transcriptional variability is rarely considered although it could influence the relationship between genotype and phenotype. It is known in unicellular organisms that gene expression is often noisy rather than uniform, and this has been proposed to be beneficial when environmental conditions are unpredictable. However, little is known about inter-individual transcriptional variability in multicellular organisms. Using transcriptomic approaches, we analysed gene expression variability between individual Arabidopsis thaliana plants growing in identical conditions over a 24-h time course. We identified hundreds of genes that exhibit high inter-individual variability and found that many are involved in environmental responses, with different classes of genes variable between the day and night. We also identified factors that might facilitate gene expression variability, such as gene length, the number of transcription factors regulating the genes and the chromatin environment. These results shed new light on the impact of transcriptional variability in gene expression regulation in plants.

Serum CXCL13 levels are associated with lymphoma risk and lymphoma occurrence in primary Sjögren's syndrome.

Primary Sjögren's syndrome (pSS) is an autoimmune disease characterised by an increased risk for non-Hodgkin lymphoma (NHL) development. Ectopic germinal centre (GC) in the salivary gland is associated with increased NHL risk in pSS, and the chemokine CXCL13 is implicated in B-cell migration and GC formation. Serum CXCL13 concentrations were quantified by ELISA in 48 healthy individuals, 273 pSS patients without NHL (pSS-nonL), and 38 pSS patients with NHL (pSS-NHL+) from the United Kingdom Primary Sjögren's Syndrome Registry cohort. PSS-nonL patients were stratified into low risk (LR), moderate risk (MR) and high risk (HR) groups according to the lymphoma risk score proposed by Fragkioudaki et al. Differences in serum CXCL13 levels among groups were analysed using the Wilcoxon method. Also, changes in serum CXCL13 over a time period of at least 1 year and a median 4 years were assessed for 200 pSS-nonL and 8 pSS-NHL+ patients. In addition, associations of serum CXCL13 with B-cell and inflammatory markers were investigated by correlation analyses and logistic regression. Serum CXCL13 levels were higher in all pSS groups compared to controls (p < 0.0001), and in pSS-NHL+ compared to pSS-nonL patients (p = 0.0204). LR patients had lower CXCL13 levels than MR patients (p < 0.0001) and pSS-NHL+ patients (p = 0.0008). CXCL13 levels remained stable over the study period for all pSS groups. CXCL13 was associated (p < 0.0005) with Immunoglobulin G (IgG), B-cell activating factor, ß2 microglobulin, combined free light chains, κ and λ light chains, anti-Ro/SSA, anti-La/SSB, and erythrocyte sedimentation rate. IgG and C3 controlled for age and gender were significantly associated with NHL risk in pSS. Serum CXCL13 levels were elevated in pSS-NHL+ and MR patients compared to LR patients and remained stable over time. Further study is required to investigate the role of CXCL13 in pSS-associated NHL risk.

Astronaut Cardiovascular Health and Risk Modification (Astro-CHARM) Coronary Calcium Atherosclerotic Cardiovascular Disease Risk Calculator.

BACKGROUND: Coronary artery calcium (CAC) is a powerful novel risk indicator for atherosclerotic cardiovascular disease (ASCVD). Currently, there is no available ASCVD risk prediction tool that integrates traditional risk factors and CAC. METHODS: To develop a CAC ASCVD risk tool for younger individuals in the general population, subjects aged 40 to 65 without prior cardiovascular disease from 3 population-based cohorts were included. Cox proportional hazards models were developed incorporating age, sex, systolic blood pressure, total and high-density lipoprotein cholesterol, smoking, diabetes mellitus, hypertension treatment, family history of myocardial infarction, high-sensitivity C-reactive protein, and CAC scores (Astro-CHARM model [Astronaut Cardiovascular Health and Risk Modification]) as dependent variables and ASCVD (nonfatal/fatal myocardial infarction or stroke) as the outcome. Model performance was assessed internally, and validated externally in a fourth cohort. RESULTS: The derivation study comprised 7382 individuals with a mean age 51 years, 45% women, and 55% nonwhite. The median CAC was 0 (25th, 75th [0,9]), and 304 ASCVD events occurred in a median 10.9 years of follow-up. The c-statistic was 0.784 for the risk factor model, and 0.817 for Astro-CHARM ( P<0.0001). In comparison with the risk factor model, the Astro-CHARM model resulted in integrated discrimination improvement (0.0252), and net reclassification improvement (0.121; P<0.0001), as well. The Astro-CHARM model demonstrated good discrimination (c=0.78) and calibration (Nam-D'Agostino χ2, 13.2; P=0.16) in the validation cohort (n=2057; 55 events). A mobile application and web-based tool were developed to facilitate clinical application of this tool ( www.AstroCHARM.org ). CONCLUSION: The Astro-CHARM tool is the first integrated ASCVD risk calculator to incorporate risk factors, including high-sensitivity C-reactive protein and family history, and CAC data. It improves risk prediction in comparison with traditional risk factor equations and could be useful in risk-based decision making for cardiovascular disease prevention in the middle-aged general population.

A new MHC-linked susceptibility locus for primary Sjögren's syndrome: MICA.

The association of primary Sjögren's syndrome (pSS) with Major Histocompatibility Complex (MHC) alleles is quintessential of MHC-disease associations. Indeed, although disease associations with classical HLA class I and II alleles/haplotypes are amply documented, further dissection is often prevented by the strong linkage disequilibrium across the entire MHC complex. Here we study the association of pSS, not with HLA genes, but with the non-conventional MHC encoded class I gene, MICA (MHC class I chain-related gene A). MICA is selectively expressed within epithelia, and is the major ligand for the activatory receptor, NKG2D, both attributes relevant to pSS' etiology. MICA-pSS association was studied in two independent (French and UK) cohorts representing a total of 959 cases and 1,043 controls. MICA*008 allele was shown to be significantly associated with pSS (pcor=2.61 × 10-35). A multivariate logistic regression showed that this association was independent of all major known MHC-linked risk loci/alleles, as well as other relevant candidate loci that are in linkage disequilibrium with MICA*008 i.e. HLA-B*08:01, rs3131619 (T), MICB*008, TNF308A, HLA-DRB1*03:01 and HLA-DRB1*15:01 (P = 1.84 × 10-04). Furthermore, independently of the MICA*008 allele, higher levels of soluble MICA proteins were detected in sera of pSS patients compared to healthy controls. This study hence defines MICA as a new, MHC-linked, yet HLA-independent, pSS risk locus and opens a new front in our understanding of the still enigmatic pathophysiology of this disease. The fact that the soluble MICA protein is further amplified in MICA*008 carrying individuals, might also be relevant in other auto-immune diseases and cancer.

Reducing the complexity of mathematical models for the plant circadian clock by distributed delays.

A major bottleneck in the modelling of biological networks is the parameter explosion problem - the exponential increase in the number of parameters that need to be optimised to data as the size of the model increases. Here, we address this problem in the context of the plant circadian clock by applying the method of distributed delays. We show that using this approach, the system architecture can be simplified efficiently - reducing the number of parameters - whilst still preserving the core mechanistic dynamics of the gene regulatory network. Compared to models with discrete time-delays, which are governed by functional differential equations, the distributed delay models can be converted into sets of equivalent ordinary differential equations, enabling the use of standard methods for numerical integration, and for stability and bifurcation analyses. We demonstrate the efficiency of our modelling approach by applying it to three exemplar mathematical models of the Arabidopsis circadian clock of varying complexity, obtaining significant reductions in complexity in each case. Moreover, we revise one of the most up-to-date Arabidopsis models, updating the regulation of the PRR9 and PRR7 genes by LHY in accordance with recent experimental data. The revised model more accurately reproduces the LHY-induction experiments of core clock genes, compared with the original model. Our work thus shows that the method of distributed delays facilitates the optimisation and reformulation of genetic network models.

Reduction of Glucocorticoid Receptor Function in Chronic Fatigue Syndrome.

Glucocorticoid receptor (GR) function may have aetiopathogenic significance in chronic fatigue syndrome (CFS), via its essential role in mediating inflammatory responses as well as in hypothalamic-pituitary-adrenal axis regulation. GR function can be estimated ex vivo by measuring dexamethasone (dex) modulation of cytokine response to lipopolysaccharide (LPS), and in vivo using the impact of dex on cortisol levels. This study aimed to compare the GR function between CFS (n = 48), primary Sjögren's syndrome (a disease group control) (n = 27), and sedentary healthy controls (HCs) (n = 20), and to investigate its relationship with clinical measures. In the GR ex vivo response assay, whole blood was diluted and incubated with LPS (to stimulate cytokine production), with or without 10 or 100 nanomolar concentrations of dex. Cytometric bead array (CBA) and flow cytometry enabled quantification of cytokine levels (TNFα, interleukin- (IL-) 6, and IL-10) in the supernatants. In the in vivo response assay, five plasma samples were taken for determination of total cortisol concentration using ELISA at half-hourly intervals on two consecutive mornings separated by ingestion of 0.5 mg of dex at 11 pm. The association of the data from the in vivo and ex vivo analyses with reported childhood adversity was also examined. CFS patients had reduced LPS-induced IL-6 and TNFα production compared to both control groups and reduced suppression of TNFα by the higher dose of dex compared to HCs. Cortisol levels, before or after dex, did not differ between CFS and HCs. Cortisol levels were more variable in CFS than HCs. In the combined group (CFS plus HC), cortisol concentrations positively and ex vivo GR function (determined by dex-mediated suppression of IL-10) negatively correlated with childhood adversity score. The results do not support the hypothesis that GR dysregulation is aetiopathogenic in CFS and suggest that current and future endocrine cross-sectional studies in CFS may be vulnerable to the confounding influence of childhood trauma which is likely increased by comorbid depression.

Frequency doubling in the cyanobacterial circadian clock.

Organisms use circadian clocks to generate 24-h rhythms in gene expression. However, the clock can interact with other pathways to generate shorter period oscillations. It remains unclear how these different frequencies are generated. Here, we examine this problem by studying the coupling of the clock to the alternative sigma factor sigC in the cyanobacterium Synechococcus elongatus Using single-cell microscopy, we find that psbAI, a key photosynthesis gene regulated by both sigC and the clock, is activated with two peaks of gene expression every circadian cycle under constant low light. This two-peak oscillation is dependent on sigC, without which psbAI rhythms revert to one oscillatory peak per day. We also observe two circadian peaks of elongation rate, which are dependent on sigC, suggesting a role for the frequency doubling in modulating growth. We propose that the two-peak rhythm in psbAI expression is generated by an incoherent feedforward loop between the clock, sigC and psbAI Modelling and experiments suggest that this could be a general network motif to allow frequency doubling of outputs.

Germline variation of TNFAIP3 in primary Sjögren's syndrome-associated lymphoma.

BACKGROUND AND OBJECTIVE: A germline and coding polymorphism (rs2230926) of TNFAIP3 (A20), a central gatekeeper of nuclear factor-kappa B (NF-kB) activation, was recently found associated with primary Sjögren's syndrome (pSS)-associated lymphoma in a French cohort. We aimed to replicate this association. PATIENTS AND METHODS: The rs2230926 polymorphism was genotyped in cases and controls of European ancestry from two independent cohorts from UK and France. Case control association tests were performed (Fisher's test) in the two cohorts, followed by a meta-analysis of the two cohorts. RESULTS: The UK cohort included 308 controls and 590 patients with pSS including 31 with a history of lymphoma. The French cohort consisted of 448 controls and 589 patients with pSS including 47 with lymphoma. In both cohorts, the rs2230926 missense polymorphism was not associated with pSS. However, in the UK cohort, the rs2230926G variant was significantly associated with pSS-associated lymphoma (OR=2.74, 95% CI (1.07 to 7.03), p=0.0423, compared with patients with pSS without lymphoma, and OR=3.12, 95% CI (1.16 to 8.41), p=0.0314, compared with healthy controls) as observed in the French cohort. The meta-analysis of the two cohorts confirmed these results (OR=2.48, 95% CI (1.87 to 3.28) p=0.0037 and OR=2.60, 95% CI (1.91 to 3.53) p=0.0031, respectively). CONCLUSIONS: This study confirms the role of A20 impairment in pSS-associated lymphoma. Subtle germline abnormalities of genes leading to impaired control of NF-kB activation in B cells continuously stimulated by autoimmunity enhance the risk of lymphoma.

Developmental mechanisms underlying variable, invariant and plastic phenotypes.

BACKGROUND: Discussions of phenotypic robustness often consider scenarios where invariant phenotypes are optimal and assume that developmental mechanisms have evolved to buffer the phenotypes of specific traits against stochastic and environmental perturbations. However, plastic plant phenotypes that vary between environments or variable phenotypes that vary stochastically within an environment may also be advantageous in some scenarios. SCOPE: Here the conditions under which invariant, plastic and variable phenotypes of specific traits may confer a selective advantage in plants are examined. Drawing on work from microbes and multicellular organisms, the mechanisms that may give rise to each type of phenotype are discussed. CONCLUSION: In contrast to the view of robustness as being the ability of a genotype to produce a single, invariant phenotype, changes in a phenotype in response to the environment, or phenotypic variability within an environment, may also be delivered consistently (i.e. robustly). Thus, for some plant traits, mechanisms have probably evolved to produce plasticity or variability in a reliable manner.

Rate of environmental change determines stress response specificity.

Cells use general stress response pathways to activate diverse target genes in response to a variety of stresses. However, general stress responses coexist with more specific pathways that are activated by individual stresses, provoking the fundamental question of whether and how cells control the generality or specificity of their response to a particular stress. Here we address this issue using quantitative time-lapse microscopy of the Bacillus subtilis environmental stress response, mediated by σ(B). We analyzed σ(B) activation in response to stresses such as salt and ethanol imposed at varying rates of increase. Dynamically, σ(B) responded to these stresses with a single adaptive activity pulse, whose amplitude depended on the rate at which the stress increased. This rate-responsive behavior can be understood from mathematical modeling of a key negative feedback loop in the underlying regulatory circuit. Using RNAseq we analyzed the effects of both rapid and gradual increases of ethanol and salt stress across the genome. Because of the rate responsiveness of σ(B) activation, salt and ethanol regulons overlap under rapid, but not gradual, increases in stress. Thus, the cell responds specifically to individual stresses that appear gradually, while using σ(B) to broaden the cellular response under more rapidly deteriorating conditions. Such dynamic control of specificity could be a critical function of other general stress response pathways.

Standards for plant synthetic biology: a common syntax for exchange of DNA parts.

Inventors in the field of mechanical and electronic engineering can access multitudes of components and, thanks to standardization, parts from different manufacturers can be used in combination with each other. The introduction of BioBrick standards for the assembly of characterized DNA sequences was a landmark in microbial engineering, shaping the field of synthetic biology. Here, we describe a standard for Type IIS restriction endonuclease-mediated assembly, defining a common syntax of 12 fusion sites to enable the facile assembly of eukaryotic transcriptional units. This standard has been developed and agreed by representatives and leaders of the international plant science and synthetic biology communities, including inventors, developers and adopters of Type IIS cloning methods. Our vision is of an extensive catalogue of standardized, characterized DNA parts that will accelerate plant bioengineering.