Full text
Interlaboratory Reproducibility in Growth and Reporter Expression in the Cyanobacterium Synechocystis sp. PCC 6803 Maurice Mager, ¶ Hugo Pineda Hernandez, ¶ Fabian Brandenburg, Luis López-Maury, Alistair J. McCormick, Dennis J. Nurnberg, Tim Orthwein, David A. Russo, Angelo Joshua Victoria, Xiaoran Wang, Julie A. Z. Zedler, Filipe Branco dos Santos,*and Nicolas M. Schmelling* Cite This: ACS Synth. Biol. 2023, 12, 1823−1835 Read Online ACCESS Metrics & More Article Recommendations * sı Supporting Information ABSTRACT: In recent years, a plethora of new synthetic biology tools for use in cyanobacteria have been published; however, their reported characterizations often cannot be reproduced, greatly limiting the comparability of results and hindering their applicability. In this interlaboratory study, the reproducibility of a standard microbiological experiment for the cyanobacterial model organism Synechocystis sp. PCC 6803 was assessed. Participants from eight different laboratories quantified the fluorescence intensity of mVENUS as a proxy for the transcription activity of the three promoters PJ23100, PrhaBAD, and PpetE over time. In addition, growth rates were measured to compare growth conditions between laboratories. By establishing strict and standardized laboratory protocols, reflecting frequently reported methods, we aimed to identify issues with state-of-the-art procedures and assess their effect on reproducibility. Significant differences in spectrophotometer measurements across laboratories from identical samples were found, suggesting that commonly used reporting practices of optical density values need to be supplemented by cell count or biomass measurements. Further, despite standardized light intensity in the incubators, significantly different growth rates between incubators used in this study were observed, highlighting the need for additional reporting requirements of growth conditions for phototrophic organisms beyond the light intensity and CO2supply. Despite the use of a regulatory system orthogonal to Synechocystis sp. PCC 6803, PrhaBAD, and a high level of protocol standardization, ∼32% variation in promoter activity under induced conditions was found across laboratories, suggesting that the reproducibility of other data in the field of cyanobacteria might be affected similarly. KEYWORDS: promoter, cyanobacteria, reproducibility, interlab ■INTRODUCTION As oxygenic phototrophs, cyanobacteria can potentially mitigate climate change when utilized by a carbon-neutral biotechnological industry. However, robust molecular biology tools are needed to study their physiology, enable their manipulation, and tap into their potential as green biotechnology platforms. Over the past decades, a plethora of new molecular tools like molecular manipulation tools and CRISPR/Cas systems have been developed for various model species. 1−3 Those tools have been adapted and further extended for cyanobacteria at an unprecedented rate in recent years. RSF1010-based replicative plasmid for various cyanobacterial species 4 and chromosomal integration for introducing transgenes have been used for decades to introduce heterologous genes into cyanobacteria. Many constitutive and inducible promoters as well as novel terminators have been characterized for transgenic expression in cyanobacteria. 5−7 Inducible promoters have been successfully implemented in more sophisticated applications such as CRISPR and CRISPRi systems as well as for metabolic engineering. 3,8 However, to ensure a seamless transfer of inducible promoter tools to other investigators, their performance needs to be optimally characterized under different experimental conditions to give users a sense of their robustness and limitations. Usually, those promoters are described and tested in cyanobacterial model strains like Synechocystis sp. PCC 6803 (hereafter PCC 6803) or Synechococcus elongatus PCC 7942. However, Cyanobacteria are a diverse phylum, represented by multiple model strains, which show differences in their optimal cultivation conditions, their ability for genetic manipulation, and the behavior of genetic parts. While lacI-based inducible expression systems have been commonly used in Synechococ- Received: March 11, 2023 Published: May 29, 2023 Research Articlepubs.acs.org/synthbio © 2023 The Authors. Published by American Chemical Society 1823 https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 Downloaded via UNIV DE SEVILLA on July 7, 2023 at 14:34:55 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.
cus 9 and Anabaena, 10 attempts to design similar systems for PCC 6803 showed very little tunability and high leaky expression. 11 Further, different investigators have reported extreme variability in promoter activity among PCC 6803 substrains. Reported fold changes for those inducible promoters range from 3- to 32-fold for the PpetE promoter and 30- to 55-fold for variations on the PrhaBAD promoter in PCC 6803 3,6,12 (see Table S1). For the PpetE promoter, varying basal activity levels have been reported 13−15 (see Table S2). Due to the lack of a standardized measurement unit for promoter activity, comparisons among different reports in the literature are nearly impossible. Thus, despite the efforts to create robust genetic tools for cyanobacteria, the implementation of inducible promoters often requires readaptation of the promoter sequence and sometimes also the protocols, even when used in the same cyanobacterial model strain under seemingly identical cultivation conditions. Reproducibility of research results has been noted as a problem in various fields and dubbed as the Reproducibility Crisis, including improper documentation of equipment and detailed methods. 16 This problem is much more extensive for cyanobacterial research as culture conditions for phototrophs are very complex. Variations in light intensity and quality, as well as CO2availability in different incubators, are additional factors that influence growth rates and cellular metabolism. Furthermore, it has been shown that optical density (OD), a standard reporting unit for the growth of bacteria, is not comparable across devices, strains, and, thus, laboratories. 17 This fact is especially an issue when using dose-responsive inducible promoters, as the ratio of inducer molecules to the number of cells will determine the actual induction rather than a defined inducer concentration. 6 However, OD is still the most dominant form used to report bacterial growth and define the starting point for experiments. Further, the respective mRNA sequence can influence the observed promoter activity, 18 as well as the type of assay used to measure promoter activity, such as fluorescence or luminescence-based reporters, protein quantification, or Northern blot. All these factors lead to different reports on promoter strengths and fold changes of inducible promoters. Interlaboratory comparisons (ILCs) are, according to the European Commission Science Hub, organized to either “check the ability of laboratories to deliver accurate testing results” or to “find out whether a certain analytical method performs well and is fit for its intended purposes”. 19 Within an ILC, the same analytical method is performed by multiple laboratories on the same samples. Subsequently, the results of each independent laboratory are compared in terms of conformity and deviation. As such, ILCs are used to assess each laboratory’s accuracy and the method’s accuracy in general. The term “interlab study” was recently popularized within the synthetic biology community by the international Genetically Engineered Machines (iGEM) competition as an ILC study to assess the reproducibility of properties from single genetic elements inside a defined context. 20−23 Such studies are a valuable tool to investigate the robustness of methods and expression platforms and identify possible sources of variation. However, those studies are seldom reported and, to our knowledge, have not been performed with cyanobacteria to test the reproducibility of widely used techniques such as promoter activity or growth. We, therefore, set up an interlab study in eight different laboratories to investigate the reproducibility of a relatively simple growth and promoter activity quantification experiment in the cyanobacterial model strain PCC 6803. We designed a protocol reflecting state-of-the-art methods reported in the literature and performed the experiment four times in each laboratory. Three commonly used promoters for heterologous expression in PCC 6803, PJ23100, PpetE, and PrhaBAD were chosen, and their transcriptional activity was quantified using the fluorescent reporter protein mVENUS as a proxy. While actual transcript levels and to a lesser extent protein levels may fluctuate in response to minute changes in environmental conditions, we chose mVENUS expression as a simplified form to approximate transcript and resulting protein levels, which we further refer to as “mVENUS expression”. Experimental conditions were standardized, leaving equipment and investigators as the primary sources of variability. OD measurements were highly reproducible within replications in a single laboratoy. However, we noticed that dilutions of the initially identical cultures to the starting OD resulted in vastly different cell biomass concentrations across laboratories. This effect results from taking OD values as a proxy for cell concentration instead of cell count. In the following, we observed significant differences in growth rate across laboratories that could not be found to correlate with the initial difference in cell biomass concentration. Even when expression of the promoters of interest was induced with concentrations above saturation level, high variability of promoter strengths was observed across different laboratories. With this information, we aim to formulate best practices for reporting these parameters to ensure better reproducibility and robustness of research results in the future. ■RESULTS AND DISCUSSION The aim of this interlab study was to assess the reproducibility of routinely performed microbiological experiments for the cyanobacterial model organism PCC 6803. The experiment of choice was a time series of the transcription strength of three promoters in PCC 6803 using the fluorescence reporter mVENUS, representing common procedures used frequently in molecular biology laboratories. Those experiments include the growth of PCC 6803 under “standard” cultivation conditions and fluorescence measurements in a plate reader and should not require overly sophisticated equipment. Furthermore, the flask type as well as the flask cap, light intensity, and shaking speed were defined to reduce the number of confounding factors. Thus, we performed the same predefined experiments (see Methods) under conditions as identical as possible across all eight independent laboratories routinely working with cyanobacteria. This study is hence aimed to give a broad, unbiased picture of the current state of reproducibility of some of the published methods in cyanobacterial research. Selection and Design of Genetic Parts and Constructs. Two widely used inducible promoters have been chosen as representative candidates. The rhamnose-dependent rhaBAD system is based on the E. coli native rhaBAD operon regulated by the AraC-like positive transcription regulator RhaS. Upon addition of rhamnose, RhaS dimers bind to the RhaS regulon and recruit RNA polymerases by interaction with E. coli sigma 70 factor RpoD. 24 Based on sequence similarity, it has been hypothesized that in PCC 6803, RhaS recruits the main sigma factor SigA instead. This is further supported by conserved sites between RpoD and SigA, which were previously described as relevant for RhaS-RpoD interaction. 13 ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1824
RhaS-based regulation of the rhaBAD promoter was utilized in PCC 6803 to reach a 55×induction when RhaS was expressed from a strong constitutive promoter J23119. 6 The copper-dependent petE system is based on the PCC 6803 native promoter of the petE gene, which is responsible for plastocyanin expression. The petE promoter is regulated by an interplay of PetR and PetP: PetR represses transcription of petE by binding to the petE promoter, while PetP is a protease that degrades PetR in the presence of copper. 14 Englund et al. observed a 5×induction of the petE promoter in PCC 6803 in the presence of copper 15 using the same experimental approach as the one used in this manuscript. The J23100 promoter is a constitutive promoter from the Anderson promoter collection. 25 This collection is a small combinatorial library of J23119 derivatives and covers multiple orders of magnitude in transcription strength in E. coli. The entire collection has been characterized in PCC 6803 by Vasudevan et al. 5 Each individual promoter was cloned upstream of the bicistronic design (BCD2) ribosomal binding site (RBS). 25 This module consists of a ribosomal binding site followed by a short reading frame, a stop codon, and a secondary ribosomal binding site. Through translational coupling, this sequence has been reported to minimize coding sequence-based bias to translation activity. Additionally, the insulation effect of the BCD2 ribosomal binding site reduces background activity from upstream genes, improving the fold changes of inducible promoter systems. We chose this RBS to reduce deviation through genetic design further. In addition, each individual promoter was cloned upstream of the mVENUS coding sequence on the RSF1010 origin of replication, a commonly used shuttle vector for PCC 6803. As an empty vector control (hereafter: EVC), the RSF1010 background only harboring the chloramphenicol resistance gene was used. A single designated laboratory was chosen to clone the desired constructs. Further, a defined PCC 6803 background strain was selected to mitigate any effects that different PCC 6803 background strains might have on the results. The designated laboratory transformed the PCC 6803 background strain and later shipped cryopreserved stock cultures on dry ice to each participating laboratory to ensure that each laboratory had genetically identical strains for the experiments. Experimental Setup of the Interlab Study. We created a set of detailed protocols to standardize the experimental conditions and data collection across all participants. The protocols contained the required information to handle the (frozen) stocks of the strains, prepare the growth medium, set up the incubator conditions, handle the cultures during the experiment, and perform the measurements. In short, glycerol stocks of the four strains, EVC, PJ23100, PpetE, and PrhaBAD (Table 1), were inoculated in liquid culture and grown for 36 to 48 h. The day before the assay, each preculture was diluted to OD730 of 0.3 to ensure cells grew exponentially at the onset of the assay. On the next day, each culture was divided into two flasks (induced and uninduced), the OD730 was adjusted to 0.5 when necessary, and inducers were added accordingly. From this moment, samples were taken during the first seven hours and after 24 h (see Figure S1 for a graphical representation of the protocol). In these samples, growth was recorded by measuring OD730 in a benchtop spectrophotometer and promoter activity by measuring fluorescence and OD730 in a plate reader. Additionally, full absorption spectra were measured as a control for the vitality of the cultures and as an internal quality control. However, these data were not included in this analysis because no notable changes in absorption spectra were observed. The raw data can be found at 10.6084/m9.figshare.21525747.v5. OD730 Measurements Are Highly Reproducible. The measurements performed by all participating laboratories, spectrophotometer OD730 and the normalized relative fluorescence units (nRFU) were used as proxies for biomass concentration and mVENUS expression, respectively (see the section Fluorescence Analysis and Normalization in Methods for a detailed explanation of our normalization strategy). We estimated the coefficient of variation (CV) to assess the reproducibility of our protocols by calculating the ratio of the standard deviation to the mean for each time point, strain, and induction regime in the first seven hours of the assay, either for each participant (intralab) or over all locations (interlab) (Figure 1). The CV, inversely proportional to the precision of replicate measurements, was lower at both the intra- and interlab level for the spectrophotometer OD730. Measurements from this data set presented a median CV of 6.1% and 11.5% (Figure 1), with 95.4% and 100% of measurements with a CV lower than 20% for intra- and interlab, respectively (Figure S2). On the other hand, the nRFU showed a median CV of 23.9% and 60.6% at the intra- and interlab levels, respectively (Figure 1). In this data set, 43% of intralab replicates had a CV lower than 20%, but at the interlab level, all replicates showed a CV higher than 20% (Figure S2). Even though the nRFU values were less reproducible than the spectrophotometer measurements, this normalization strategy clearly improved the comparability of results across laboratories. This can be seen from the much higher interlab CV values calculated from not normalized values as the background-corrected fluorescence units (FUbc) or the relative fluorescence units (RFU) (Figure S3). Interestingly, including the OD730 in the normalization procedure did not lead to lower CV at the interlab level (Figure S3). This analysis shows that the measurements performed in the spectrophotometer were the most reproducible in our protocol, which can be partially explained by the fact that the starting conditions of the assay were based on measurements from these devices. Discrepancies in Growth Rates Are Not Explained by Different Initial Biomass Concentrations. Even though the spectrophotometer results showed the highest reproducibility, we further evaluated if this resulted from highly reproducible culture conditions. To assess this, we used the spectrophotometer OD730 to estimate growth rates for each biological replicate over 24 h (Figure 2). This metric is not dependent on the actual OD730 values per se but rather on their relationship. The growth rate is known to be influenced by environmental conditions. We performed ANOVA to test if significant differences were found at the strain, induction regime, or laboratory level. The latter was the only factor Table 1. Strains Used in This Study strain name fluorescent reporter reporter promoter inducer EVC − − − PJ23100 mVENUS PJ23100 constitutive PpetE mVENUS PpetE CuSO4 PrhaBAD mVENUS PrhaBAD rhamnose ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1825
showing a significant influence (P-value < 0.0001) on the measured growth rates. Therefore, we compared the average growth rate between laboratories using a post hoc Tukey’s test. The results of this analysis showed significant differences (P-value < 0.05) between most of the participants. However, the most striking and biologically relevant differences were found between the laboratories, with the highest (Seville, Tuebingen, Jena, and Amsterdam) and lowest growth rates (Duesseldorf and Leipzig). In the first group, growth rates varied between 0.056 and 0.058 h−1, while in the second, it ranged between 0.035 and 0.039 h−1. This translates into a 36% reduction in the average growth rate. Our protocols precisely defined the growth conditions (temperature, shaking speed, light intensity) and media composition. However, since we used the spectrophotometer OD730 as a reference value to set the initial biomass concentration in the assay, the starting amount of photons per cell could have differed among laboratories. For example, this could have been the case if spectrophotometers used across laboratories had different relationships between OD730 and the number of cells. To test this hypothesis, we performed additional measurements in each spectrophotometer by preparing a dilution of the four glycerol stocks and directly measuring it without allowing cells to grow. The results (Figure S4A) showed that, indeed, there were significant differences (ANOVA, P-value < 0.05) across spectrophotometers when measuring the same amount of cells. With this information, we could determine if the differences in growth rates we observed during the assays resulted from the initial biomass concentration. Therefore, we normalized the initial OD730 values of the assays by the measurements obtained from the glycerol stocks (Figure S4B). Next, we tested whether growth rates correlated with this relative OD730 but found no significant correlation (Pearson’s coefficient: 0.089, P-value: 0.834, Figure S4C). Thus, we can conclude that even if the starting amount of photons per cell was not the same across laboratories, this could not explain the observed differences in growth rates. PpetE Expression Across Laboratories Is Less Reproducible than PrhaBAD When No Inducer Is Added. The uninduced PJ23100 (PJ23100−, see the Promoter Quantification Assay section) RFU was used as an internal standard to compare the expression of mVENUS across laboratories (Figure 3). As expected, we observed an average increase in fluorescence during the first seven hours after induction in PpetE and PrhaBAD strains. However, these promoters’ expression patterns differed in reproducibility across laboratories, leakiness, or magnitude of induction. The reproducibility of the nRFU within each laboratory was comparable for both strains and induction regimes (Figure 4). The median CV of the cultures without induction was 19.3% and 29.3% for PrhaBAD and PpetE, respectively. Interestingly, the reproducibility across laboratories was higher when an inducer was added (median CV of 16.4% and 17.1% for PrhaBAD and PpetE, respectively). Since we used the PJ23100−RFU to normalize the RFU values of the other strains, it was impossible to use the nRFU CV to compare the reproducibility of fluorescence measurements across all strains. Thus, we additionally calculated the CV of the not normalized RFU at the intralab level (Figure S5). PJ23100 measurements were the most reproducible among all strains and induction regimes. When no inducer was added, EVC, PpetE, and PrhaBAD showed similar reproducibility. However, when the inducer was added, the PpetE and PrhaBAD measurements were more reproducible than the EVC RFU. As expected, the observed variation of nRFU at the interlab level was higher than within each laboratory (Figure 4). When cultures were induced, increased reproducibility was observed compared to uninduced. A pattern that was already observed on the intralab level comparison. However, while CV values were comparable between both induced strains (median CV of 32% and 31.2% for PrhaBAD and PpetE, respectively), the Figure 1. Coefficient of variation (%) of nRFU in EVC, PrhaBAD, and PpetE (left) and spectrophotometer OD730 (right) data sets. The coefficient of variation was calculated for either all the replicates within a laboratory (intralab, in green) or for all the replicates across all laboratories (interlab, orange). In each panel, a boxplot summarizes all data points by showing the median as a horizontal line, the 25th and 75th percentiles as the bottom and top of the box, respectively, and whiskers extending 1.5 ×IQR from the box margins. In addition, all data points are shown and distributed over the area of a violin plot. Figure 2. Measured growth rates in each laboratory. Bars represent the mean growth rate and error bars 95% CI (n= 32 or 24). Text boxes on the right side of the bars show the results from Tukey’s test. Laboratories with significantly different growth rates (P-value < 0.05) are labeled with different letters. ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1826
difference between induced and uninduced was much more pronounced. Further, the interlab variability observed in the uninduced PpetE cultures was much higher than for the uninduced PrhaBAD cultures. We observed nonoverlapping distribution of CV values between the two strains, where PpetE’s median CV value was 67% compared to the 46.8% PrhaBAD’s median CV. Inducible promoters can be characterized by measuring a change in expression strength between induced “ON” and uninduced “OFF” states. We, therefore, investigated how these parameters varied between the PrhaBAD and PpetE strains and how they varied between laboratories. We calculated the changes in transcription activity upon induction by dividing the nRFU at seven hours by the value at the beginning of the assay for both induced and uninduced cultures (Figure 5A). When the inducer was absent, the transcriptional activity did not change over the experiment as expected. Changes for all strains in the absence of an inducer were close to 1, with a mean change of 0.86, 95% CI [0.77, 0.95] for PrhaBAD and 0.93, 95% CI [0.73, 1.13] for PpetE. When an inducer was added, PpetE cultures showed the largest average change (2.87, 95% CI [1.46, 4.28]) while also showing high variability between laboratories (Figure 5A). On the other hand, PrhaBAD cultures showed an average change of 2.35, 95% CI [1.84, 3.76] with less variability across laboratories. Regarding expression strength, we observed that after seven hours of induction and with the addition of 1 μM of CuSO4, the PpetE promoter led to similar mVENUS expression levels as the PJ23100 promoter (average 106% of PJ23100, 95% CI [89%, 123%]) (Figure 5B). Under the same conditions and with 10 mM rhamnose, the average expression of the PrhaBAD was 32% of PJ23100 (95% CI [27%, 37%]). The leakiness of an inducible promoter refers to its basal expression level in the absence of the inducer. To assess the PpetE and PrhaBAD promoters’ leakiness, we calculated the ratio of RFU OD730−1at seven hours in the uninduced cultures to the EVC (Figure 5C). This ratio is equal to 1 if no signal from the fluorescent reporter is measured in the uninduced condition and will increase proportionally to the leakiness of the promoter. The uninduced PrhaBAD cultures showed an average ratio of 1.6, 95% CI [1.3, 1.9]. Conversely, the mean ratio in the PpetE cultures was 8.4, 95% CI [0.9, 15.9] when no inducer was added, suggesting that either the promoter was very leaky or the environment contained too many residual copper ions (Figure 5C). The high variability observed across the different laboratories (Figure 5D), with an average ratio within laboratories ranging from 1.2 to 38.5 and a median value of 4.8, suggests that residual copper present in the medium of some of the participating laboratories is responsible. Across eight different laboratories, the growth of PCC 6803 was quantified via OD730 measurements, and the transcriptional activity of three promoters (PJ23100, PrhaBAD, and PpetE) was quantified using fluorescence intensity from mVENUS as a Figure 3. Time series of promoter assay. The nRFU over time is shown for the four strains, each depicted in an individual panel. Cultures, where an inducer was added, are shown in yellow, and those without are shown in black. Smaller data points represent the average values for a single laboratory (n= 4 or 3). Larger data points show the overall average and error bars of the overall 95% CI (n= 7). The induced and uninduced conditions for EVC and PJ23100 refer to specific aspects of the assay preparation. See the Promoter Quantification Assay section in Methods for further details. ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1827
proxy. Reproducibility was evaluated by using inter- and intralaboratory coefficients of variation. We discovered that OD730 values showed high reproducibility, albeit the absolute OD730 values were not normalized between different laboratories. The fluorescence intensity values of reporters were less reproducible than OD730. Generally, fluorescence intensity values of induced promoters were more reproducible than uninduced promoters. Reproducibility Is Not Homogeneous Across the Different Measurement Methods. The spectrophotometer measurements in the growth assay showed the highest reproducibility (Figure 1). However, this was expected since the initial OD730 of the assay was determined from this device. Nevertheless, the growth rates estimated from these measurements showed significant differences across the participating laboratories (Figure 2). It is well-known that OD measurements are only comparable across devices with additional calibrations. 17 However, as this study is intended to investigate experimental variation according to commonly used experimental practices, we were interested in the effect of this variation on the experimental results. Thus, we chose not to calibrate the starting OD730 by cell counts across laboratories. Instead, calibration measurements were initially performed to later infer the influence of differences in OD730 measurements on the overall reproducibility of the study. In these calibration measurements, the reported OD730 varied more than 5-fold between participants, confirming the high variability of bulk OD730 measurements between spectrophotometer devices. Although this shows that internally, OD730 measurements were very reproducible (Figure 1), they were practically incomparable across laboratories without normalization. Therefore, the starting OD730 target at 0 h does not represent the same biomass concentration across different laboratories. Initially, we expected to find a correlation between the growth rate and the normalized initial biomass concentration. However, our data did not reveal such a correlation (Figure S4C), and it remains a matter of speculation which factor caused the differences in the growth rates. It is possible that these differences arose from variations in the light spectra since these are known to affect this physiological parameter. 26,27 In our protocol, the light color was specified as “white”, which is how it is commonly defined in scientific literature. White light, however, is a combination of different wavelengths, and light spectra across laboratories could have been different. Thus, we compared the light spectra across laboratories and incubators. However, most of the laboratories did not have the instruments to measure the light spectra accurately; thus, we relied, in most cases, on manufacturer specifications. Unfortunately, those specifications were, at best, if at all available, a picture of the light spectrum. After aligning all available spectra and comparing the peaks, we did not find any conclusive trend that explains the observed differences in growth rate. The plate-reader measurements resulted in less reproducible data at the intra- and interlab level than the spectrophotometer measurements (Figure 1). This observation can be partially explained by the fact that the initial conditions of the assay were set based on the latter device. In addition, several other factors could have influenced the reproducibility of the estimated RFU. First, two different measurements were performed to obtain these values and, thus, two different sensors (OD and fluorescence), increasing the chances of Figure 4. Coefficient of variation (%) of the nRFU measured for all time points in PpetE and PrhaBAD cultures without induction (−, left column) and with induction (+, right column) at the interlab (top row) and intralab level (bottom row). The coefficient of variation was calculated for either all the replicates within a laboratory (intralab, in green) or for all the replicates across all laboratories (interlab, orange). In each panel, a boxplot summarizes all data points by showing the median as a horizontal line, the 25th and 75th percentiles as the bottom and top of the box, respectively, and whiskers extending 1.5 ×IQR from the box margins. In addition, all data points are shown and distributed over the area of a violin plot. Notice that the Y-axis range is different in the top and bottom rows. ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1828
measurement errors. Next, it has been shown that fluorescence measurements of weak promoters (such as PrhaBAD) are less reproducible due to the poor sensitivity of plate readers at low signals. 21,22 Lastly, the inclusion of two different induction regimes, which certainly affect the fluorescent output but not the growth behavior, might have contributed to the lower reproducibility of the estimated RFU ×OD730−1compared to the spectrophotometer measurements. Uninduced Reporter Systems Show a Higher Level of Variation Compared to Systems Induced above Saturation Levels. Interestingly, the promoter activity of inducible promoters was more reproducible in the presence than in the absence of the respective inducer. Introducing an additional experimental step (addition of the inducer) did not lead to a higher level of variation. Further, the observed significant differences between measured growth rates did not seem to translate into high variation in promoter activity, highlighting the robustness of those promoters across different physiological stages of cells. Lastly, the differences in initial cell biomass, inferred by differences in OD730 values of the identical stock between laboratories, resulted in differences in the inducer/cell count ratio. However, as induction was performed above saturation concentration for the inducer based on previous experiments, 6 the investigated promoters both seem to perform relatively robustly toward variations in inducer/cell count ratio. It is important to note that the expression level of the PrhaBAD promoter in the OFF state is extremely low and barely above background signals. Thus, much of the intralab CV of the PrhaBAD can be attributed to the background noise and not variations in expression level (compare Supplemental Figure S5). The PpetE system shows higher location-dependent variation. Whereas residual copper could explain the higher variance in the PpetE promoter system, this can be excluded for PrhaBAD as rhamnose is unlikely to be naturally present in trace concentrations. Varying copper residuals on glassware or in filtered water are expected to be the main culprit of this issue. Copper is a ubiquitously present metal ion that is difficult to remove altogether. Acid-washing glassware prior to experiments 28 could help here. Additionally, the use of chelators such as bathocuproinedisulfonic acid disodium salt 14 in the medium is a way to reduce copper availability efficiently during the experiment. For future studies employing the PpetE promoter in real applications, we recommend supporting the data with a fluorescence activity-based assay, as done in this study. This additional experiment should be done in the same strain used in the application and under similar experimental conditions. This process control not only helps to approximate the activity of the PpetE promoter in the actual experiment but also serves as a process control that can be used by investigators trying to replicate the results. As our results not only show high activity and variance in the OFF state of the PpetE promoter, we recommend carefully observing experimental results of a strain Figure 5. Characterization of PpetE and PrhaBAD promoters. (A) nRFU ratio of time 7 h over time 0 h. (B) nRFU values at time 7 h of induced cultures. (C) nRFU ratio of uninduced PpetE and PrhaBAD to EVC culture at 7 h. (D) nRFU ratio of uninduced PpetE to EVC culture at 7 h in each laboratory. From A to C, points represent the mean value per laboratory (n= 3 or 4), bars show the average across all laboratories, and error bars show 95% CI (n= 7). In D, points show the mean value of technical replicates from a single experiment (n= 3), and bars show the average from each laboratory (n= 3 or 4). The color of bars and points indicate the induction regime, black for uninduced and yellow for induced cultures. In C and D, the horizontal and vertical black lines, respectively, represent a ratio of 1. ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1829
with a seemingly uninduced PpetE promoter. Furthermore, these native problems of the PpetE promoter might be mitigated in the future by further engineering the promoter system: heterologous coexpression of the petR/petP regulators may positively impact the sensitivity to residual copper. Strategies to Improve Comparability and Reproducibility. With the increasing throughput of experiments, improving reporting and reproducibility standards becomes increasingly important. None of the following recommendations are definite answers to this question. Instead, they are intended to provoke discussions and reflections about the reproducibility of cyanobacterial research. A good strain-handling policy is important for keeping research results reproducible over time. Especially for bacterial strains, this is essential since genetic changes occur rapidly. Those changes resulting from improper strain maintenance lead to genotypic diversity and consequently to conflicts in research findings. The effect of the genetic diversity of PCC 6803 and its impact on phenotype has been analyzed in the past. 29,30 Cryopreservation techniques have been established for long-term microbial evolution experiments to reduce genotypic variability between experiments with the same strain. 31 To our knowledge, such techniques are not standard practice in the cyanobacterial field. In this study, we aimed to minimize strain variability not only between participants but also between experimental runs within a single location. Therefore, we adapted a cryopreservation and inoculation protocol based on glycerol conservation from Price et al. 31 to the PCC 6803 strains used in this interlab study. With this protocol, we hope to contribute to standardizing strain handling for cyanobacteria, especially for PCC 6803. Next to varying practices in strain-keeping, the preparation of BG11 medium is a source for variations. During the establishment of the experimental procedures, we encountered that the preparation and final mass concentration of almost all solutes in BG11 media differed across participating laboratories. In addition, we obtained similar outcomes when comparing BG11 media protocols published in the scientific literature. Regarding standardization, we agreed on the formula developed by van Alphen et al., 32 except for removing copper from the media for our particular purposes. Additionally, as outlined in the methods, a weak HEPES buffer was added to standardize starting conditions after inoculation. As cyanobacterial cultures rapidly increase the pH of their environment during cultivation, this buffer is assumed to be nonsignificant in prolonged cultivation. However, we assumed it would be better to buffer initial pH disruptions during inoculation. It remains open if those changes in BG11 media would substantially affect the experiments. However, we would encourage everyone also to standardize the preparation of BG11 as most differences in the preparation have no fundamental reason other than long-lasting traditions in respective laboratories. In addition, when working with phototrophs, reporting the exact specification of light spectra used would be an important addition to reporting growth conditions. However, measuring light spectra requires a light meter, which is not standard laboratory equipment. Furthermore, manufacturers of light bulbs do not disclose those specifications in most cases or in an inappropriate format. Thus, we encourage manufacturers of light bulbs to report the light spectra as a CSV file with raw values. Further, we would encourage scientists to ask for the light spectra when buying a new incubator or contact the manufacturer to request the light spectra of your current incubators. The iGEM interlab studies have gathered a considerable amount of data and protocols that support the use of external calibrants for both optical density and fluorescence for promoter characterization. 20−23 There are indeed clear benefits in defining such protocols. By calibrating plate readers (or flow cytometers) with external standards for fluorescence and cell concentration, promoter activity can be reported in absolute units, allowing direct comparisons across laboratories and facilitating the detection of biological and technical errors. 20−23 However, these protocols are not universally applicable. Regarding fluorescence reporters, GFP is not the only fluorescent protein used in cyanobacterial research, and YFP derivatives such as mVENUS have become popular in recent years. 33 Therefore, this requires that new calibrants are found for each new fluorescent protein, matching its excitation and emission wavelengths. It has been known that OD measurements differ between instruments and, thus, laboratories. 17,34 Hence, only reporting OD values as a measure for cell biomass is problematic but still common practice in microbiology. The Stevenson et al. and later the iGEM interlab studies proposed using materials that match the refractive index of E. coli to calibrate OD. 20−23,34 In addition to being species-specific, this approach presents the disadvantage of assuming that cells of a given species are always the same size and shape. An option to overcome this problem would be to report cell counts and size with the starting OD values so that individual researchers can calibrate their OD measurements to the respective cell count. 35 We are aware that cell counters are not always available in all laboratories, as was the case among the participants of this study. However, we encourage reporting cell counts in combination with OD for better reproducibility in the future. Our alternative to overcome these issues was using an internal biological standard. In our experimental design, we employed PJ23100 as the normalizing promoter within each location, therefore accounting for differences between instruments. In addition, this approach also excludes factors such as the intracellular abundance of RNA polymerases and sigma factors, ribosomes, and other components of protein production impacting the mVENUS expression and maturation. Increased mVENUS expression from one of the inducible promoters related to any of these factors was expected to affect the mVENUS expression in the J23100 strain equally. As all our RFU values are reported as RFUstrain/RFUJ23100, these factors are presumably either buffered or completely negated. Considering that no significant differences in growth rates between the strains were found within each laboratory and identical absorption spectra across all the used strains in all laboratories, highly similar growth conditions can be assumed for all data sets. By correcting these factors, this normalization approach enables enhanced comparability across laboratories. This is demonstrated by the reduced interlab CV in nRFU, compared to the non-normalized metrics (Figure S3). It is worth noting that we observed a slightly lower coefficient of variation (CV) between laboratories when employing FUbc instead of RFU for our normalization method (Figure S3). Essentially, including OD730 in the procedure led to an increased variation due to the intrinsic error of these measurements propagating into nRFU. Although this small increase in variation is present, we consider it justified as it addresses potential variation in fluorescence due to cultures ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1830
with different optical densities. Therefore, we recommend such a normalization strategy to report more reproducible data. Lastly, our experimental setup used the native copper repressible promoter PpetE and the heterologous PrhaBAD system as a somewhat orthologous regulation system. Even though the BG11 medium used in this study was prepared without adding copper salts, copper is hard to remove from glassware, and even a slight amount of residual copper could have caused the differences in basal activity observed in our experiments. Additionally, as a native promoter, PpetE may be subjected to more levels of regulation than so far known, which could lead to greater variability across different laboratories. Thus, we encourage the use of orthogonal systems like the PrhaBAD promoter, which overall performed better in this study in terms of reproducibility. Considering the orthogonal nature of the PrhaBAD promoter in PCC 6803 and assuming no residual rhamnose in uninduced media, we hypothesize that the interlaboratory reproducibility of this experimental procedure and the multitude of measures taken to ensure maximum comparability of results across this study, we propose that the coefficients of variation of promoter activity for this system could indicate where the baseline for maximum reproducibility of quantitative data in PCC 6803 lies. Essentially, the coefficient of variation for the PrhaBAD promoter serves as a call for caution when it comes to the overall reproducibility of fluorescence activity�and potentially other�data in the field of cyanobacterial research and a call for action to aid in improving comparability by following the discussed strategies: standardized and rigid strain keeping policies, standardized media composition, thorough reporting of cultivation conditions and the correct use of normalization agents. ■CONCLUSION In conclusion, the reproducibility of promoter activity in cyanobacteria within this study across multiple laboratories was better than expected. However, we also identified some causes for errors that could be improved in future studies. The remaining problem when using cyanobacteria is that most options to correct for those errors are understudied or completely missing, and those tools would first need further development to increase reproducibility. ■METHODS The strains for this study were generated in a single laboratory (Duesseldorf) and distributed to all participants. The interlab experiment was performed following a detailed set of protocols, which can be found at 10.17504/protocols.io. 3byl4j69rlo5/v1. See Figure S1 for a visual description of the protocol. Plasmid and Strain Construction. The plasmids were constructed via Golden Gate Cloning. 1 Complete and annotated sequences of all plasmids are available in supplementary data on Figshare (10.6084/m9.figshare. 21525747.v5). According to standard procedures, Escherichia coli DH5αcells were transformed with Golden Gate mixes according to the MoCloFlex protocol. 36 Sequences were verified via Sanger sequencing. PCC 6803 was conjugated via triparental mating, 37 and conjugants were confirmed via colony PCR. The PCC 6803 strain (glucose tolerant, nonmotile) was obtained from D. Bhaya (Carnegie Institution for Science, Stanford, USA). Media. All participants prepared a BG11 medium without CuSO4supplemented with 10 mM NaHCO3and 5 mM HEPES-NaOH (pH 8). The exact composition of BG11 and the detailed recipe for the stock’s preparation can be found in the supplements (Table S3 and S4) and was extracted from the Supporting Data S3 document of van Alphen et al. 32 It is hereafter referred to as BG11. For each repetition of the assay, 1 L of BG11 was prepared fresh for each round of experiments: First, approximately 500 mL of ultrapure H2O were autoclaved in a 1 L volumetric bottle. Then stock solutions were acclimated to room temperature and added in the following order: 5 mL of 1 M HEPES-NaOH, 2.5 mL BG11 S1, 2.5 mL BG11 S2 (without CuSO4), 2.5 mL BG11 S3, and 10.5 mL 0.95 M NaHCO3. After adding the solutions, sterile ultrapure H2O was added up to the 1 L mark of the bottle. Cultivation for Strain Conservation. For the preparation of cryoconserved cultures, PCC 6803 was grown under the same conditions outlined in Culture Conditions in Interlab Experiment, except for increasing the light intensity to 80 μ mol photons ×m−2×s−1. When the OD measured at 730 nm wavelength (hereafter OD730) reached a value of 3, cells were centrifuged at 12,000g(fixed angle rotor) for 15 min, washed with fresh BG11 medium, centrifuged once again, and resuspended in 10% of the initial volume with BG11 and 15% (v/v) final concentration glycerol. Cells were stored at −80 °C and shipped to participating laboratories on dry ice. The incubator and photometer listed in Table S5 under “Duesseldorf I” were used to prepare cryoconserved cultures. Participants. Eight different laboratories participated in the study and are described throughout the study by the geographical location of the research facility. In alphabetical order, the participants were from Amsterdam (University of Amsterdam), Berlin (Freie Universitat Berlin), Duesseldorf (Heinrich Heine University), Edinburgh (University of Edinburgh), Jena (Friedrich Schiller University), Leipzig (Helmholtz Centre for Environmental Research), Seville (University of Seville), and Tuebingen (University of Tubingen), as included in the author list. In Duesseldorf, the complete experiment was independently performed by two researchers, indicated in the text as Duesseldorf (I) and Duesseldorf (II). Strains were prepared and grown in Duesseldorf and shipped to participants on dry ice. Culture Conditions in Interlab Experiment. Cultures were grown in 100 mL Erlenmeyer nonbaffled glass flasks with cotton plugs in shaking incubators set at 100 rpm (see Table S5) under 50 μmol photons ×m−2×s−1constant white light illumination, ambient CO2, and 30 °C (except for Tuebingen where the temperature was 28 °C). Cultures were grown in BG11 medium, as outlined above. Cultures were always grown with 10 μg mL−1chloramphenicol, which was added individually to each flask before inoculation. Promoter Quantification Assay. PCC 6803 strains harboring plasmids containing the investigated promoterreporter cassette were inoculated from cryoconserved cultures by adding 330 μL of inoculum to 10 mL of copper-free BG11 medium supplied with chloramphenicol and cultivated for 48 h. Cultures were then diluted to an OD730 of 0.3 in 35 mL and grown overnight to OD730 of 0.5−0.6. The following day, the complete volume of each preculture was transferred to a sterile tube and diluted to OD730 0.5 if necessary. The flasks where ACS Synthetic Biology pubs.acs.org/synthbio Research Article https://doi.org/10.1021/acssynbio.3c00150 ACS Synth. Biol. 2023, 12, 1823−1835 1831