Natural changes in light interact with circadian regulation at promoters to control gene expression in cyanobacteria

  1. Joseph Robert Piechura
  2. Kapil Amarnath
  3. Erin K O'Shea  Is a corresponding author
  1. Harvard University, United States
8 figures, 6 tables and 1 additional file

Figures

The circadian and light response pathways in cyanobacteria.

(A) Schematic of gene expression output of the circadian clock under Constant Light conditions. Under Constant Light intensity (dashed navy blue line), dawn gene expression (dashed maroon line) and …

https://doi.org/10.7554/eLife.32032.003
Figure 2 with 4 supplements
Natural clear day conditions sharpen the expression of dusk genes to peak just before expected darkness.

(A) Experimental setup for testing the effects of Clear Day conditions on circadian gene expression. The upper plot shows the light intensity profiles of Low Light (black) and Clear Day (magenta) …

https://doi.org/10.7554/eLife.32032.004
Figure 2—source data 1

Normalized gene expression in Low Light and Clear Day conditions.

https://doi.org/10.7554/eLife.32032.009
Figure 2—figure supplement 1
Pigment levels of cyanobacteria grown under Low Light or Clear Day conditions reveal adjustments in the photosynthetic apparatus to optimize growth in different light conditions.

(A) Estimation of phycocyanin and chlorophyll levels in cells grown under Low Light (black) or Clear Day (magenta) conditions for two days, measured at midday of the third light period. Phycocyanin …

https://doi.org/10.7554/eLife.32032.005
Figure 2—figure supplement 2
Gene expression dynamics of dusk and dawn circadian genes under Constant Light conditions (data from Markson et al., 2013).

(A) Gene expression dynamics of circadian genes over 24 hr in Constant Light conditions in wildtype cells (left heat map) and over 12 hr in OX-D53E cells (rpaA-, kaiBC-, Ptrc::rpaA(D53E)) (middle …

https://doi.org/10.7554/eLife.32032.006
Figure 2—figure supplement 3
Dawn gene expression increases during the early part of Clear Day relative to Low Light conditions.

(A) Gene expression dynamics of dawn genes (n=169) under Low Light (top) and Clear Day (bottom) conditions. Gene expression is quantified as the log2 fold change from the average expression of the …

https://doi.org/10.7554/eLife.32032.007
Figure 2—figure supplement 4
The gene expression dynamics of glycogen production and breakdown enzymes change in Clear Day conditions relative to Low Light conditions.

(A) Gene expression dynamics of the dusk gene glgP, encoding a key enzyme in glycogen breakdown, under Low Light (black) and Clear Day (magenta) conditions as measured by RNA sequencing (left …

https://doi.org/10.7554/eLife.32032.008
Figure 3 with 2 supplements
Rapid changes in light intensity modulate the recruitment of RNA polymerase to dusk genes to control dusk gene expression.

(A) Light intensity profiles of Low Light (black) and High Light pulse (orange) conditions, in units of μmol photons m-2 s-1 (see Materials and methods - Calibrating light conditions for more …

https://doi.org/10.7554/eLife.32032.010
Figure 3—source data 1

Normalized gene expression in High Light pulse and Shade pulse conditions.

https://doi.org/10.7554/eLife.32032.013
Figure 3—source data 2

List of RNAP peaks, gene targets, and quantification of enrichment under High Light pulse and Shade pulse conditions.

https://doi.org/10.7554/eLife.32032.014
Figure 3—figure supplement 1
Rapid changes in light intensity affect dawn gene expression in an opposite direction compared to dusk gene expression.

(A) Gene expression dynamics of dawn genes (n=169) under High Light pulse conditions. Gene expression is quantified as the log2 fold change from the average expression of the gene over all time …

https://doi.org/10.7554/eLife.32032.011
Figure 3—figure supplement 2
Changes in RNAP enrichment and downstream dusk gene expression after rapid changes in light intensity.

(A) Changes in enrichment of RNAP upstream of dusk genes during High Light pulse conditions (left heat map) and corresponding changes in target dusk gene expression (right heat map) for the 82 dusk …

https://doi.org/10.7554/eLife.32032.012
Figure 4 with 4 supplements
Changes in environmental light intensity regulate RpaAP DNA binding activity and RNAP recruitment to control dusk gene expression downstream of clock regulation of RpaA.

(A) Phosphorylation dynamics of RpaA under Low Light vs High Light pulse. Relative levels of phosphorylated RpaA were measured using Phos-tag Western blotting (left y-axis) in cells grown under Low …

https://doi.org/10.7554/eLife.32032.015
Figure 4—source data 1

Quantification of relative RpaA∼P levels.

https://doi.org/10.7554/eLife.32032.020
Figure 4—source data 2

List of RpaA peaks, gene targets, and quantification of enrichment under High Light pulse and Shade pulse conditions.

https://doi.org/10.7554/eLife.32032.021
Figure 4—figure supplement 1
Representative Western blots used to quantify relative levels of RpaAP under dynamic light conditions.

(A) Representative Western Blot used to quantify levels of RpaAP under Low Light and Clear Day conditions. Lysates were prepared from cells harvested from either Low Light (L) or Clear Day (C) …

https://doi.org/10.7554/eLife.32032.016
Figure 4—figure supplement 2
Changes in RpaA enrichment and downstream dusk gene expression after rapid changes in light intensity.

(A) Changes in enrichment of RpaA upstream of dusk genes during High Light pulse conditions (left heat map) and corresponding changes in target dusk gene expression (right heat map) for the 56 dusk …

https://doi.org/10.7554/eLife.32032.017
Figure 4—figure supplement 3
Changes in RpaA and RNA polymerase enrichment upstream of dusk genes after rapid changes in light intensity.

(A) Changes in enrichment of RpaA upstream of dusk genes during High Light pulse conditions (left heat map) and corresponding changes in RNAP enrichment upstream of the same gene (right heat map) …

https://doi.org/10.7554/eLife.32032.018
Figure 4—figure supplement 4
Multifactorial behavior of RpaAP at select promoters under changes in light intensity.

(A)-(C) Normalized ChIP-seq signal of RpaA (red), RpaB (blue), RNAP (green) and mock IP (black) upstream of the (A) the representative dusk gene Synpcc7942_2267, (B) the kaiBC operon and (C) another …

https://doi.org/10.7554/eLife.32032.019
Figure 5 with 3 supplements
Light-induced changes in RpaBP levels modulate RpaB and RNAP binding upstream of dusk genes to directly regulate dusk gene expression in response to light.

(A) Phosphorylation dynamics of RpaB under Low Light vs High Light pulse. Relative levels of phosphorylated RpaB were measured using Phos-tag Western blotting (left y-axis) in cells grown under Low …

https://doi.org/10.7554/eLife.32032.022
Figure 5—source data 1

Quantification of relative RpaB∼P levels.

https://doi.org/10.7554/eLife.32032.026
Figure 5—source data 2

List of RpaB peaks, gene targets, and quantification of enrichment under High Light pulse and Shade pulse conditions.

https://doi.org/10.7554/eLife.32032.027
Figure 5—figure supplement 1
Representative Western blots used to quantify relative levels of RpaBP under dynamic light conditions.

(A) Representative Western Blot used to quantify levels of RpaBP under Low Light and Clear Day conditions. Lysates were prepared from cells harvested from either Low Light (L) or Clear Day (C) …

https://doi.org/10.7554/eLife.32032.023
Figure 5—figure supplement 2
Changes in RpaB enrichment and downstream dusk gene expression after rapid changes in light intensity.

(A) Changes in enrichment of RpaB upstream of dusk genes during High Light pulse conditions (left heat map) and corresponding changes in target dusk gene expression (right heat map) for the 42 dusk …

https://doi.org/10.7554/eLife.32032.024
Figure 5—figure supplement 3
Changes in RpaB and RNA polymerase enrichment upstream of dusk genes after rapid changes in light intensity.

(A) Changes in enrichment of RpaB upstream of dusk genes during High Light pulse conditions (left heat map) and corresponding changes in RNAP enrichment upstream of the same gene (right heat map) …

https://doi.org/10.7554/eLife.32032.025
Figure 6 with 1 supplement
Global regulation of dusk gene expression in response to light changes.

(A) Number of dusk gene targets of RpaA only (red), RpaB only (blue), RpaA and RpaB (yellow), or neither (black). Target genes of binding sites of RpaA and RpaB were determined using chromatin …

https://doi.org/10.7554/eLife.32032.028
Figure 6—figure supplement 1
Regulation of dusk sigma factor gene expression by RpaA and RpaB.

(A)-(C) Normalized ChIP-seq signal of RpaA (red), RpaB (blue), RNAP (green) and mock IP (black) upstream of the sigma factor genes (A) rpoD6, (B) rpoD5, and (C) sigF2. The location of the gene is …

https://doi.org/10.7554/eLife.32032.029
Figure 7 with 1 supplement
Dusk genes group into three major clusters that show distinct and coordinated responses to changes in light intensity.

(A) Average expression profiles of genes belonging to the Early dusk gene cluster under Clear Day (magenta) and Shade pulse (gray) conditions (left y-axis). Dusk genes were grouped using k-means …

https://doi.org/10.7554/eLife.32032.030
Figure 7—source data 1

Lists of genes belonging to the Early, Middle, and Late dusk clusters, and scaled gene expression values.

https://doi.org/10.7554/eLife.32032.032
Figure 7—figure supplement 1
Average expression profiles of the major dusk gene clusters under various conditions.

(A) Average expression profiles of the Early (left plot), Middle (middle plot), and Late (right plot) dusk gene clusters under Low Light (black) and High Light pulse (orange) conditions (left …

https://doi.org/10.7554/eLife.32032.031
Figure 8 with 4 supplements
Phenomenological modeling of the activation of clusters of light-responsive dusk genes.

(A) Normalized RpaAP levels (left plot) and RpaBP levels (right plot) under Clear Day (magenta) and Shade pulse (gray) conditions used as input for mathematical models of dusk gene expression. …

https://doi.org/10.7554/eLife.32032.033
Figure 8—figure supplement 1
Best fit simulations of ‘RpaA-only’ and ‘RpaB-only’ models in which RpaAP or RpaBP solely activates the expression of the dusk gene clusters.

(A) Normalized RpaAP levels under Clear Day (magenta) and Shade pulse (gray) conditions used as input for mathematical models of dusk gene expression. RpaAP levels from all four light conditions …

https://doi.org/10.7554/eLife.32032.034
Figure 8—figure supplement 2
Models in which either the Middle or Late cluster feeds back to influence Early cluster expression.

(A) Feedback model in which the expression of the Early dusk cluster is an activation Hill function of Middle gene expression and an activation Hill function of both RpaAP and RpaBP. The left …

https://doi.org/10.7554/eLife.32032.035
Figure 8—figure supplement 3
Models in which either the Early or Late cluster feeds back to influence Middle cluster expression.

(A) Feedback model in which the expression of the Middle dusk cluster is an activation Hill function of Early gene expression and an activation Hill function of both RpaAP and RpaBP. The left …

https://doi.org/10.7554/eLife.32032.036
Figure 8—figure supplement 4
Models in which either the Early or Middle cluster feeds back to influence Late cluster expression.

(A) Feedback model in which the expression of the Late dusk cluster is an activation Hill function of Early gene expression and an activation Hill function of both RpaAP and RpaBP. The left plot …

https://doi.org/10.7554/eLife.32032.037

Tables

Key resources table
Reagent type or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Synechococcus elongatus)PCC 7942 (wild-type)ATCCCat. Num. 33912
Strain, strain background (Escherichia coli)Tuner (DE3)EMD MilliporeCat. Num. 70263
Gene (S. elongatus)RNA polymerase Beta’ subunitN/ACyanobase: Synpcc7942_1524
Gene (S. elongatus)rpaBN/ACyanobase: Synpcc7942_1453
Recombinant DNA reagentRNA polymerase beta prime subunit FLAGThis paperAddgene: 102337Plasmid encoding C-terminal FLAG tag RNA polymerase Beta’ subunit (Synpcc7942_1524) with Kan selection marker, targeted to integrate at native gene locus
Recombinant DNA reagentpET-48b(+)EMD MilliporeCat. Num. 71462
Renetic reagent (S. elongatus)EOC 398 and EOC 399This paperS. elongatus PCC7942 transformed with RNA polymerase beta prime subunit FLAG plasmid. Confirmed by PCR and Western blot.
Antibodyanti-RpaBThis paperAnti-RpaB serum was produced by Cocalico Biologicals. Anti-RpaB was affinity purified as described in this work.
Antibodyanti-RpaAThis paperAnti-RpaA serum was produced by Cocalico Biologicals as described in Markson et al., 2013. Anti-RpaA was affinity purified as described in this work.
AntibodyFLAG M2 mouse monoclonal antibodySigma AldrichCat. Num. F3165
Software, algorithmImagequantGE Healthcare
Software, algorithmBowtiePMID: 19261174
Software, algorithmPeak-SeqPMID: 19122651
Software, algorithmMATLABMathWorks
Commercial assay or kitRNeasy Mini kitQiagenCat. Num. 74104
Commercial assay or kitRibo-Zero bacteria rRNA removal kitIlluminaCat. Num. MRZMB126
Commercial assay or kitTruseq Stranded mRNA sample prep kitIlluminaCat. Num. 20020594
Commercial assay or kitNEBNext Ultra II DNA library prep kitNew England BiolabsCat. Num. E7645S
Chemical compound, drugPhos-tagAcrylamide AAL-107Wako Pure Chemical IndustriesCat. Num. 304–93521
Table 1
Fitting bounds.

Bounds used for fitting the variables in our simple model of gene expression. H is the Hill coefficient, β is the max transcription rate, α is the decay/dilution rate, B is the background …

https://doi.org/10.7554/eLife.32032.038
VariableLower boundUpper bound
H07
β080
α080
B010
K01
Table 2
Fitting results.

The definitions of the variables are given in Equations 1-3, p. 1–3. The error is defined as the square root of the sum of the squared deviations between simulation and data.

https://doi.org/10.7554/eLife.32032.039
ModelClusterFigureBXβXαXKAXHAXKBXHBXKYXHYXError
RpaA-onlyEarly7D0.7137.5472.710.716.76----0.85
RpaB-onlyEarly7-Fig. Supp. 2C0.3724.0378.62--0.370.78--1.01
RpaA and RpaBEarly7G0.3551.2837.760.354.190.82.5--0.41
Feedback, M act.Early7-Fig. Supp. 3A0.0155.8530.010.010.30.872.380.062.470.37
Feedback, M rep.Early7-Fig. Supp. 3B0.6758.6938.890.676.960.622.470.9670.24
Feedback, L act.Early7-Fig. Supp. 3C0.235.8719.030.24.430.983.350.056.150.38
Feedback, L rep.Early7I, 7-Fig. Supp. 3D0.7569.3442.680.756.220.593.530.712.390.21
RpaA-onlyMiddle7D0.7937.95630.796.76----0.86
RpaB-onlyMiddle7-Fig. Supp. 2C0.260.03---0.265.6--0.85
RpaA and RpaBMiddle7G157.4625.9714.960.524.12--0.29
Feedback, E act.Middle7-Fig. Supp. 4A0.823.7322.190.86.960.494.530.216.350.32
Feedback, E rep.Middle7-Fig. Supp. 4B0.7371.0839.240.735.140.536.580.740.880.35
Feedback, L act.Middle7I, 7-Fig. Supp. 4C0.1878.6376.50.186.090.332.640.161.550.16
Feedback, L rep.Middle7-Fig. Supp. 4D0.6831.0217.980.683.340.576.79100.44
RpaA-onlyLate7D0.9639.8264.370.966.7----0.78
RpaB-onlyLate7-Fig. Supp. 2C0.0500--0.050.68--0.79
RpaA and RpaBLate7G0.9577.6567.10.9570.485.9--0.5
Feedback, E act.Late7-Fig. Supp. 5A0.9923.9320.010.995.80.46.950.186.770.53
Feedback, E rep.Late7-Fig. Supp. 5B0.7659.8118.430.766.220.696.130.473.120.29
Feedback, M act.Late7I, 7-Fig. Supp. 5C0.3727.316.090.373.720.013.460.916.230.22
Feedback, M rep.Late7-Fig. Supp. 5D0.8625.114.460.866.920.487100.52
Table 3
Parts for controllable light source.

The table includes the parts chosen for their specific properties. The remaining parts, such as wires, heat shrink tubing, thermal paste for mounting the LEDs on the heat sinks, proto-boards, and …

https://doi.org/10.7554/eLife.32032.040
Part nameDigikey part numberCurrent price ($)Quantity
PWR SUP MEDICAL 18V 8.3A 150WEPS439-ND73.711
CONN RCPT 8CONT DIN SLD PNL MNTSC2007-ND5.641
LEDDynamics Flexblock BUCK BOOST 48V, 700 mA788–1038-ND19.994
AD7376 digital potentiometerAD7376ARWZ10-ND8.664
AC to DC power supply, 10VDC, 275 mA993–1233-ND4.682
BXRA-30E1200-B-03, Bridgelux, Warm white, LEDNot sold at Digikey.
Need to order from:10.474
AMBIT ELECTRONICS, INC.
Aavid thermalloy Spotlight 47W heat sink1061–1092-ND9.504
Arduino Uno Board Rev31050–1024-ND21.491
Table 4
Wiring the FlexBlock LED driver.

The FlexBlock LED driver needs to be connected in a ’boost only’ configuration (see spec sheet for more details), with connections as shown.

https://doi.org/10.7554/eLife.32032.041
LineConnection
DIM GNDGND of 10 V power supply/Arduino
DIMWipe of AD7376 potentiometer (Pin 16)
Vin++of 18V power supply AND + of LED array
Vin-GND of 18V power supply
LED+NC (not connected)
LED-- of LED array
Table 5
Wiring the AD7376 potentiometer.

We used the SOIC-16 housing for the AD7376 potentiometer for ease of soldering to wires. The table indicates how each pin was connected. The length of the GND wire from the Arduino board to the …

https://doi.org/10.7554/eLife.32032.042
PinConnection
1+of 10 V power supply
2GND (shared GND between that of 10V power supply and Arduino
3GND
4GND
5pin 10 on Arduino (or any other pin designated as a Slave Select, such as 5, 6, or 9
6+5V of Arduino
7pin 13 on Arduino (SCLK)
8NC (not connected)
9NC
10NC
11pin 11 on Arduino (MOSI)
12+5V of Arduino
13NC
14+of 10V power supply
15NC
16DIM line of FlexBlock

Additional files

Download links