Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 13 August 2021
Sec. Microbial Physiology and Metabolism
This article is part of the Research Topic Regulation of Prokaryotic Cell Division View all 17 articles

Hyperosmotic Shock Transiently Accelerates Constriction Rate in Escherichia coli

  • 1Department of Bioengineering, Stanford University, Stanford, CA, United States
  • 2Department of Microbiology and Immunology, Stanford University School of Medicine, Stanford, CA, United States
  • 3Chan Zuckerberg Biohub, San Francisco, CA, United States

Bacterial cells in their natural environments encounter rapid and large changes in external osmolality. For instance, enteric bacteria such as Escherichia coli experience a rapid decrease when they exit from host intestines. Changes in osmolality alter the mechanical load on the cell envelope, and previous studies have shown that large osmotic shocks can slow down bacterial growth and impact cytoplasmic diffusion. However, it remains unclear how cells maintain envelope integrity and regulate envelope synthesis in response to osmotic shocks. In this study, we developed an agarose pad-based protocol to assay envelope stiffness by measuring population-averaged cell length before and after a hyperosmotic shock. Pad-based measurements exhibited an apparently larger length change compared with single-cell dynamics in a microfluidic device, which we found was quantitatively explained by a transient increase in division rate after the shock. Inhibiting cell division led to consistent measurements between agarose pad-based and microfluidic measurements. Directly after hyperosmotic shock, FtsZ concentration and Z-ring intensity increased, and the rate of septum constriction increased. These findings establish an agarose pad-based protocol for quantifying cell envelope stiffness, and demonstrate that mechanical perturbations can have profound effects on bacterial physiology.

Introduction

The bacterial cytoplasm contains a dense combination of nucleic acids, proteins, and numerous metabolites that together generate a high internal osmotic pressure relative to the extracellular environment. These turgor pressures of ∼1–10 atm (Cayley et al., 2000; Deng et al., 2011) must be balanced to maintain cellular integrity. To counter this mechanical load, virtually all bacteria have a rigid, highly cross-linked peptidoglycan cell wall (Holtje, 1998) that can bear substantial stress (Yao et al., 1999). Bacterial cells typically encounter a wide range of environmental osmolalities, for example enteric bacteria rapidly transition between high and low osmolalities inside and outside the host, respectively (Gauthier, 2000). Rapid increases in external osmolality cause plasmolysis in which water exits the cytoplasm, leading to reduced stress on the cell envelope [membrane(s) and cell wall] and cytoplasm shrinkage (Cota-Robles, 1963; Pilizota and Shaevitz, 2013). During relatively small hyperosmotic shocks (≤100 mM of osmolyte), cell-wall insertion and elongation continue unaffected (Rojas et al., 2014), while larger hyperosmotic shocks inhibit growth (Walter, 1924; Christian and Scott, 1953; Scott, 1953; Christian, 1955) as well as intracellular diffusion and molecular mobility (van den Bogaart et al., 2007; Boersma et al., 2015). In general, whether intracellular processes are affected by osmolality changes has not been fully investigated, particularly how changes to mechanical load shape cellular structure and physiology.

Microfluidic devices are powerful tools for monitoring growth and division dynamics at single-cell resolution (Taniguchi et al., 2010; Wang et al., 2010; Campos et al., 2014; Taheri-Araghi et al., 2015; Camsund et al., 2020), particularly during acute environmental changes such as osmotic shocks (Pilizota and Shaevitz, 2012; Rojas et al., 2014, 2017; Zhou et al., 2015; Buda et al., 2016; Rojas and Huang, 2018; Yang et al., 2020). A salient example of the utility of microfluidics in bacterial cell biology is the recent discovery that the Gram-negative outer membrane (OM) can bear substantial mechanical stress (Rojas et al., 2018). In a flow cell, Escherichia coli cells were exposed to a large hyperosmotic shock, followed by detergent treatment that induced lysis. The large contraction upon lysis indicated that the stiffness of the OM is comparable to that of the cell wall (Rojas et al., 2018).

Although microfluidic devices can provide dynamical single-cell information and can be used to screen libraries using elaborate designs (Taniguchi et al., 2010; Camsund et al., 2020), cost and throughput is often limiting, with only one strain or species typically tested at a time. Moreover, microfluidic devices limit the movement of cells via rigid physical constraints and hence allow only cells within a particular size range to enter, which imposes additional mechanical stress and makes a single device incompatible with species across a wide range of shapes and sizes. The unintended selection of particular sizes in microfluidic devices also potentially introduces biases (Oliveira et al., 2020). The traditional alternative to microfluidic devices for single-cell imaging is agarose pads, which are versatile platforms that are easy to prepare and applicable for morphologically diverse species. Several recent studies have introduced high-throughput methods for rapidly imaging collections of strains on large-format agarose pads (Kuwada et al., 2015; Shi et al., 2017b), enabling screening of genome-scale libraries. However, it is difficult to track the effects of acute environmental transitions on single cells using agarose pads; instead one must rely on population averages measured pre- and post-transition. Moreover, in the time period required for pad preparation, physiological changes may have taken place that are not captured by snapshots. Thus, it remains unclear whether osmotic shock-related phenomena can be robustly probed in high-throughput on agarose pads.

As a critical part of the bacterial cell cycle, cell division is highly regulated. In bacteria, a ring of the conserved tubulin homolog FtsZ (the “Z-ring”) (Bi and Lutkenhaus, 1991; Dai and Lutkenhaus, 1991) assembles at mid-cell and initiates assembly of the divisome machinery (Goley et al., 2011; Barrows and Goley, 2021). After the Z-ring forms and recruits other division proteins, it progressively constricts the membrane (Osawa and Erickson, 2013) and directs synthesis of the septal cell envelope (Bisson-Filho et al., 2017; Yang et al., 2017) at a constant rate (Reshes et al., 2008b), resulting in growth of new hemispherical endcaps. FtsZ concentration has been linked to the proportion of dividing cells, which changes across nutrient conditions (Ward and Lutkenhaus, 1985; Weart and Levin, 2003; Weart et al., 2007; Hill et al., 2012). Regulation of FtsZ expression affects cell size homeostasis (Si et al., 2019) and FtsZ synthesis and degradation predict the timing of the first division in starved cells supplied with nutrient pulses (Sekar et al., 2018). While it remains unclear whether constriction of the bacterial inner membrane must fight against turgor pressure (Erickson, 2009, 2017), in fission yeast reduction of turgor pressure accelerates cell division (Basu et al., 2014; Chang, 2017). Taken together, bacterial cell division is a natural candidate for processes affected by environmental osmolality.

Here, we developed an agarose pad-based protocol for Gram-negative envelope stiffness measurements, and sought to establish that pad measurements could recapitulate previous findings regarding the stiffness of the OM. To our surprise, the population-averaged length of cells after hyperosmotic shock as observed on agarose pads was significantly smaller than expected based on microfluidic measurements. By tracking single-cell dynamics during an osmotic shock in a microfluidic flow cell, we discovered that the rate of cell division transiently increased after the shock, and the fraction of dividing cells was quantitatively consistent with the apparent decrease in cell length. Treatment with the division inhibitor aztreonam was sufficient to recover the same level of contraction as observed by tracking single cells, thereby enabling quantification of the effects of hyperosmotic shock and detergent treatment in high-throughput on agarose pads. Using a strain expressing a functional msfGFP-tagged FtsZ, we show that FtsZ concentration and Z-ring intensity increased rapidly upon osmotic shock, and that constriction dynamics were significantly accelerated relative to steady-state growth, supporting the hypothesis that the divisome can provide constrictive force that acts against turgor pressure.

Results

Apparent Length Contraction Due to Hyperosmotic Shock Is Much Larger on Agarose Pads Than in a Microfluidic Flow Cell

In a previous study (Rojas et al., 2018), we exposed E. coli cells as well as other Gram-negative bacteria to a large hyperosmotic shock in a microfluidic flow cell to remove turgor pressure (Figure 1A). The plasmolyzed, shrunken cells were then exposed to either EDTA (to remove lipopolysaccharides from the OM) or detergent (to remove the OM altogether); both treatments caused the length of the cell wall to decrease even further (Rojas et al., 2018), indicating that the OM was exerting substantial stretching forces on the cell wall in the plasmolyzed state. We carried out similar experiments and observed similar cell wall mechanical strains (normalized length changes relative to the unextended state) using the wall dye wheat germ agglutinin conjugated to AlexaFluor 488 (WGA-AF488) (Ursell et al., 2014) upon hyperosmotic shock and detergent treatment (Figures 1B,C); unshocked cell walls were extended by 9.6% relative to shocked cells, with a larger contraction upon lysis (Figure 1D).

FIGURE 1
www.frontiersin.org

Figure 1. Population-averaged relative length change after hyperosmotic shock is greater on agarose pads than from single-cell tracking in microfluidic devices. (A) Schematic of microfluidic (top) and agarose pad-based (bottom) measurements of length changes (right) due to hyperosmotic shock (plasmolysis) and detergent-mediated lysis. (B) Relative to normal growth (top), cells exposed to a hyperosmotic shock in a microfluidic flow cell (middle) decreased in length and exhibited plasmolysis (white arrow). Subsequent lysis due to detergent treatment (bottom) resulted in further shrinkage. Shown are phase contrast images (left) and fluorescence from the wall label WGA-AF488 (right). (C) Quantification of cell length dynamics for 5 representative cells for the experiment in panel B. Cell elongation was inhibited after the shock. (D) Quantification of the mechanical strain (extension relative to the smaller state) in cell wall length due to hyperosmotic shock (left) and detergent treatment (right). The additional strain revealed by detergent treatment shows that the stiffness of the outer membrane is comparable to that of the cell wall (Rojas et al., 2018). (E) Population-averaged cell wall length of the same set of cells in a microfluidic device before and after hyperosmotic shock yielded a similar estimate of plasmolysis strain as in panel D. (F) Cells imaged during normal growth in a test tube (top), on agarose pads with 20% sucrose to cause a hyperosmotic shock (middle, white arrows indicate plasmolysis bays), and after detergent treatment to cause lysis (bottom). (G) The relative change in length of cells on an agarose pad before and after a hyperosmotic shock (left) was substantially greater than the relative length change in a microfluidic device shown in panel D, as was the relative length change after lysis (right).

While tracking of cells within the microfluidic device allowed us to measure the relative length change of each cell throughout the hyperosmotic shock and subsequent OM destabilization, we found that the plasmolysis strain could be accurately quantified based on the population-averaged length before and after each treatment (Figure 1E). Moreover, cell growth was inhibited for at least 5 min after the hyperosmotic shock, with most non-dividing cells maintaining a stable length (Figure 1C). Thus, we hypothesized that we could take advantage of traditional agarose pad-based imaging to quantify OM stiffness in high throughput (Figure 1A).

To test whether agarose pad measurements would recapitulate our microfluidics-based measurements of length changes due to hyperosmotic shock, we exposed an exponentially growing E. coli culture in LB to the same hyperosmotic shock as in our microfluidics experiments. To perform the shock as quickly as possible, we directly spotted cells onto agarose pads containing LB+20% sucrose to induce a large (∼0.5 M) osmotic shock, and imaged the cells as quickly as possible (Figures 1A,F). To our surprise, the mean length strain was 32.2% (Figure 1G), far larger than for cells in the microfluidic device (Figures 1D,E).

Next, we diluted exponentially growing cells into a concentrated solution of detergent and sucrose (see section “Materials and Methods”) to induce cell lysis. The lysed cells were spotted onto agarose pads with LB for imaging after 20 min of detergent treatment. The decrease in cell length after lysis (Figure 1G) was also much larger than expected based on our microfluidic measurements (Figure 1D). Thus, we conclude that despite the seemingly simple nature of our length measurements, there was an unidentified factor that was affecting cell length differentially in the microfluidic device and on agarose pads.

Osmotic Shock Transiently Increases the Rate of Cell Division

Despite the slowing of elongation after the hyperosmotic shock (Figure 1C), we noted that some cells continued to divide (Figure 2A). Thus, we hypothesized that one cause of the discrepancy between our microfluidic and agarose pad measurements was the subset of cells that completed constriction in the time interval between the shock and the time of image acquisition, resulting in a decrease in the population-averaged length post-shock. For a population of cells with average length L during exponential growth that is extended by a factor ε relative to its shocked length L0, L = L0(1+ε). If the fraction of cells that divide between the time of the shock and agarose pad imaging is f, the total number of cells increases by the factor (1+f). Therefore, the mean length of shocked cells decreases by (1+f), yielding an apparent extension during exponential growth relative to shocked cells of (1+ε)(1+f) – 1 = ε+ff (Figure 2B). Thus, we hypothesized that the fraction of cells that were dividing after the shock was approximately 20%, based on the difference in the apparent extension on pads and in the microfluidic device.

FIGURE 2
www.frontiersin.org

Figure 2. A transient burst of cell division follows hyperosmotic shock. (A) In the 3 min after exposure to a hyperosmotic shock, a fraction of cells continued to divide (red arrows indicate a few examples) despite growth largely halting (Figure 1C). Shown is a time-lapse in phase contrast (top) and fluorescence of the wall label WGA-AF488 (bottom). (B) If a fraction f of cells divide during the interval between the hyperosmotic shock and imaging on an agarose pad, the average length of cells post-shock will decrease by a factor (1+f), leading to a predicted increase in the plasmolysis strain by (1+f). (C) The fraction of dividing cells (see section “Materials and Methods”) increased from 11.6% pre-shock to 17.3% in the 3 min following the shock. By 10 min after the shock, division had halted completely.

To test this hypothesis, we manually identified cells that divided in a 3-min time interval before and after the hyperosmotic shock using microfluidics. Because it is difficult to pinpoint the time of division initiation due to the diffraction limit of light microscopy, we used WGA-AF488 to label the cell wall and identified cells that formed a clear septum during the 3-min interval (see section “Materials and Methods”). The fraction of dividing cells pre-shock was 11.6% per 3 min (Figure 2C), consistent with our measured doubling time of ∼23 min (Supplementary Figure 1). In the 3 min after the shock, the dividing fraction increased to 17.3% (Figure 2B), in reasonable agreement with our prediction. Division halted completely 10 min after the shock, likely due to the lack of growth (Figure 2C). Thus, we conclude that cell division was a major factor in the discrepancy between contraction estimates.

Inhibition of Cell Division Restores the Same Length Contraction After Hyperosmotic Shock as in Microfluidics

To determine whether cell division was the sole factor distinguishing our microfluidic and agarose pad measurements, we treated cells with the beta-lactam aztreonam, which inhibits the division-specific transpeptidase PBP3 (Spratt, 1975). We exposed cells to 50 μg/mL aztreonam in a microfluidic flow cell to determine the time scale of complete division inhibition, and found that division completely halted within 20 min (Figure 3A). We then exposed cells in a test tube to aztreonam for varying amounts of time before hyperosmotic shock (Figure 3B), and as expected mean length increased monotonically with the duration of treatment (Figure 3C). Moreover, the degree of contraction after hyperosmotic shock relative to initial length decreased with the duration of treatment, plateauing after 20 min (Figure 3C), consistent with the time scale of complete division inhibition (Figure 3A).

FIGURE 3
www.frontiersin.org

Figure 3. Inhibition of cell division allows for accurate measurement of length contraction after hyperosmotic shock. (A) During time-lapse imaging of E. coli cells treated with 50 μg/mL aztreonam, cell division was fully inhibited within 20 min. (B) Schematic of modified protocol for agarose pad-based measurements of cell envelope stiffness, in which cells are treated with 50 μg/mL aztreonam for various amounts of time before imaging on agarose pads. (C) Mean length of cells grown in a test tube with 50 μg/mL aztreonam increased with the duration of treatment, due to the inhibition of cell division. After a hyperosmotic shock on an agarose pad with LB+20% sucrose, the measured length strain gradually decreased with the duration of treatment and plateaued at ∼10%, consistent with microfluidic device measurements (Figure 1D). (D) After detergent treatment, the population-averaged length of aztreonam-treated cells was lower by an amount corresponding to a strain of 37.2%, consistent with microfluidic device measurements (Figure 1D).

We then exposed cells treated for 20 min with aztreonam to a solution of sucrose and detergent. The vast majority of cells lysed (>90%), and the mean length decreased such that unshocked cells were extended by 37.2% relative to the lysed cell wall length (Figure 3D). This value is reasonably consistent with the distribution of extensions of single cells in a microfluidic device (Figure 1D). We also treated cells with detergent in the absence of a hyperosmotic shock, surmising that the removal of the OM should remove both the membranes and turgor pressure, and thereby allow the cell wall to relax to its unstretched state. To our surprise, most cells (>95%) failed to lyse, suggesting that cells are more susceptible to detergent while plasmolyzed.

Thus, we conclude that the transient burst of cell division after hyperosmotic shock is responsible for the large decrease in population-averaged cell length, and that division inhibition is sufficient to enable pad-based measurements of cell envelope mechanical properties.

Septum Formation Time Is Largely Constant Across Mutants With Variable Cell Width

In a previous study, we showed that the concentration of FtsZ is largely constant across a set of MreB mutants with various cell widths and volumes (Shi et al., 2017a). However, the Z-ring in wider mutants was both broader and more intense (Shi et al., 2017a), suggesting that Z-ring assembly is width-dependent, which may lead to different dynamics during septum formation. To further probe the connection between cell width and septum formation, we selected a set of MreB mutants with different cell widths but similar growth rates in exponential phase (Shi et al., 2017a), and performed time-lapse imaging to determine the duration of septum formation in each strain (see section “Materials and Methods”). In these strains, mean cell width varied from 1 μm to 1.6 μm, corresponding to a ∼2.5-fold difference in septum areas. Nonetheless, all strains completed septum formation in ∼19 min, regardless of their width (Figure 4A; r = −0.13, p = 0.68, two-sided Student’s t-test). Thus, the broader and more intense Z-rings in wider cells potentially boost the rate of septum formation and allow cells to synthesize a larger septum within the same duration.

FIGURE 4
www.frontiersin.org

Figure 4. Hyperomostic shock leads to increased FtsZ concentration and ring intensity and faster constriction. (A) Septum formation time (see section “Materials and Methods”) for mreB mutants of various cell widths was approximately constant as a function of mean cell width. (B) Phase-constrast (top) and epifluorescence (middle) images of E. coli cells expressing FtsZ-msfGFP, before (left) and 3 min after (right) a hyperosmotic shock. Fluorescence profiles along the cell contour show that Z-ring intensity increased immediately after the shock (bottom). (C) The concentration of FtsZ-msfGFP (left) and cytoplasmic GFP (right) increased directly after the hyperosmotic shock, consistent with the decrease in cell volume (Figure 1B). (D) Z-ring intensity increased after the shock. Shown are 6 representative cells; the cell colored in green maintained high Z-ring intensity for 10 min after the shock, while Z-ring intensity in the maroon cell increased and then decreased for the ensuing 4 min leading up to the completion of cell division. (E) Z-rings constricted more rapidly after hyperosmotic shock (right) than during normal growth (left). The black curves are representative cells in each category. (F) The time required for Z-ring constriction was substantially shorter after hyperosmotic shock.

FtsZ Concentrates at Midcell After Hyperosmotic Shock

To probe the mechanism underlying the transient increase in division after hyperosmotic shock, we performed time-lapse imaging in a microfluidic flow cell using an E. coli MG1655 strain expressing ftsZ fused to msfGFP as the sole copy of ftsZ (Figure 4B). This strain exhibits approximately normal growth (Moore et al., 2017) although the fusion does affect cell length by altering FtsZ GTPase activity (Yang et al., 2017). We first quantified the total amount and concentration of FtsZ within each cell (see section “Materials and Methods”). FtsZ concentration (total amount normalized by cell volume) increased directly after hyperosmotic shock (Figure 4C), as expected based on cytoplasmic shrinkage. A similar concentration increase was observed for another divisome protein, ZapA (Buss et al., 2013), in a strain expressing GFP-ZapA (Supplementary Figure 2), as well as in a control strain expressing cytoplasmic GFP (Figure 4C). In all cases, the increase was larger than the normalized length decrease of the cell wall, due to plasmolysis collapsing the cytoplasm away from the cell wall (Figure 4B). In addition to FtsZ and ZapA, other divisome proteins likely increase rapidly in concentration given the short time scale of cytoplasmic contraction.

To address whether the increase in FtsZ concentration is associated with accelerated constriction, we quantified the intensity and width of the Z-ring (Shi et al., 2017a). Strikingly, Z-ring intensity increased almost immediately after the shock (Figure 4D). Moreover, Z-ring constriction proceeded at a much faster pace in shocked compared with non-shocked cells (Figures 4E,F), despite the lack of overall volume expansion (Figure 1C). This accelerated Z-ring constriction largely explains the increase in the fraction of cells that complete constriction shortly after osmotic shock (Figure 2C). It is unclear whether the rate of division initiation also increased; regardless, initiation rate is unlikely to be the only cause of increased cell division after shock since septum formation typically requires ∼10–15 min during normal growth. Together, these results support the hypothesis that the enhanced rate of cell division post-hyperosmotic shock is due at least in part to the higher intensity of the Z-ring and faster rate of constriction.

Discussion

In this study, we developed an agarose pad-based measurement protocol to quantify the stiffness of the OM from cell length measurements after hyperosmotic shock and lysis. Our measurements of length strain after plasmolysis were inconsistent between agarose pads and microfluidic devices (Figure 1), which was quantitatively explained by the increase in cell division directly after hyperosmotic shock (Figures 2–4). Treatment with the division inhibitor aztreonam restored the expected degree of length strain (Figure 3C), and thus enables characterization of osmotic shock-induced morphological changes and OM stiffness on agarose pads (Figure 3D). This strategy provides a simpler and more versatile methodology than microfluidics, especially for screening a diverse set of species or mutants. Using single-cell tracking, we discovered that FtsZ (Figure 4C) and ZapA (Supplementary Figure 2) concentration and Z-ring intensity (Figure 4D) increased and Z-ring constriction accelerated (Figures 4E,F) directly after a hyperosmotic shock, providing mechanistic hypotheses for the increased rate of division during plasmolysis. Division may also be affected by the likely increase in the concentration of cell-envelope precursors after the shock.

Our findings suggest that progression of bacterial cell division is coupled to environmental osmolality. The observation that the Z-ring constricts more quickly in the absence of turgor suggests that the divisome exerts constrictive forces on the envelope. During a hyperosmotic shock, both the decrease in turgor and increase in FtsZ concentration may contribute to the acceleration of constriction; it is difficult to decouple the two factors, since the cytoplasmic shrinkage that accompanies the reduced stress on the cell envelope inevitably concentrates FtsZ and presumably all other divisome components. Given the immediate response of cells to osmotic shock, the association between division and osmolarity points to possible ecological implications wherein transitions to a high osmotic environment that limits growth will result in proliferation into a larger population of cells that increases survival odds, reminiscent of the reductive divisions that take place when cells enter stationary phase (Sutterlin et al., 2016). Perhaps turgor decreases along with stationary-phase entry to facilitate division, or in other environmental transitions in which an increase in cell number would be beneficial. In our experiments, division halted almost completely after the initial boost induced by the shock despite the higher concentration of FtsZ relative to normal growth (Figure 2C), potentially due to other limiting factors that prevented growth (Figure 1C). It remains unclear how the activities of FtsZ and other divisome components are generally affected by osmolarity or intracellular crowding.

Our finding that cell division is affected by osmolarity suggests that environmental perturbations such as osmotic shocks may generally induce widespread cellular responses. Indeed, E. coli cells were found to respond in a coupled manner to heat and oxygen shock (Tagkopoulos et al., 2008), likely reflecting the correlation between high temperature and anoxia in human hosts. Certain antibiotics also induce acid stress (Mitosch et al., 2017) or heat-shock response pathways concurrently (Evans et al., 2019). In addition to the divisome, osmolality may affect other cellular structures such as the cytoskeletal filament MreB (Szatmári et al., 2020) and the nucleoid (Finan and Guilak, 2010; Cagliero and Jin, 2013; Wu et al., 2019; Yang et al., 2020), both of which could have widespread downstream physiological consequences. Our findings confirm that agarose pads can be used as a high-throughput platform for bacterial cell mechanics measurements, although microfluidics provides a critical ground truth in which the fate of single cells can be tracked. Thus, rapid screening of large libraries of mutants (Baba et al., 2006; Peters et al., 2016) or species should be straightforward as long as cell division can be inhibited, which should help to uncover the molecular basis of envelope stiffness and understand the role of physical forces in shaping cellular structures and bacterial physiology.

Materials and Methods

Strain Culturing

Escherichia coli cells were grown overnight in LB at 37°C and used to inoculate a test tube with fresh LB. Strains used in this study are listed in Supplementary Table 1.

WGA Staining of the Cell Wall

Cell wall was labeled with wheat germ agglutinin (WGA) conjugated to AlexaFluor-488 (AF488, Invitrogen W11261). WGA-AF488 was added to exponential phase cells at a final concentration of 25 μg/mL and incubated in dark conditions with shaking at 37°C for at least 2.5 h prior to imaging.

Hyperosmotic Shock Application on Agarose Pads

To apply a hyperosmotic shock on an agarose pad, cells were first cultured in liquid LB into steady-state exponential growth using the following serial dilution protocol: starting with a 1:200 dilution from an overnight saturated culture, 1:10 dilutions were performed after 60 and 150 min. After the second dilution, ∼1 μL of liquid culture was spotted onto a 1% agarose pad with LB+20% sucrose, leading to a hyperosmotic shock.

Detergent Treatment for Agarose-Pad Measurements

Exponentially growing cells were exposed to a lysis solution of N-lauroylsarcosine sodium salt (detergent) and sucrose dissolved in 1X PBS. By directly diluting cultures into the lysis solution at a 1:1 ratio, the final concentration of detergent and sucrose was 5% and 10%, respectively. Sucrose was added to promote lysis, recapitulating the perturbation in a microfluidic device in which cells were plasmolyzed by a hyperosmotic shock before being exposed to detergent.

Single-Cell Imaging

One microliter of cells was spotted onto a 1% agarose pad with the appropriate medium, and imaged on a Nikon Eclipse Ti-E inverted fluorescence microscope with a 100X (NA 1.40) oil-immersion objective (Nikon Instruments). For microfluidic experiments, cells were loaded into B04A (CellASIC) microfluidic plates following previous protocols (Rojas et al., 2018). Phase-contrast and epifluorescence images were collected on a DU885 electron-multiplying CCD camera (Andor Technology) or a Neo sCMOS camera (Andor Technology) using μManager v. 1.4 (Edelstein et al., 2010). Cells were maintained at 37°C during imaging with an active-control environmental chamber (Haison Technology).

Identifying Cell Division Events From WGA Signal

In a microfluidic flow cell, cell division events were identified based on the clear formation of a septum in the WGA-AF488 signal. Although septum formation does not represent the end of cell division, the septum can be easily identified visually or from a peak in the fluorescence profile along the long axis of the cell. The fraction of dividing cells was computed based on a combination of visual inspection and fluorescence profile examination, with constriction dynamics in the phase channel serving as further confirmation.

Image Analysis

The MATLAB (MathWorks, Natick, MA, United States) image processing code Morphometrics (Ursell et al., 2017) was used to segment cells and to identify cell outlines from phase-contrast and epifluorescence microscopy images. A local coordinate system was generated for each cell outline using a method adapted from MicrobeTracker (Sliusarenko et al., 2011). Cell length was calculated as the length of the midline from pole to pole. See figure legends for the number of cells analyzed (n) and error bar definitions.

Quantification of FtsZ Fluorescence

Total FtsZ fluorescence was calculated by integrating fluorescence values within the cell contour after background subtraction, and FtsZ concentration was computed as the total fluorescence divided by cross-sectional area. The fluorescence profile along the cell contour had low values at the cell poles and high values at mid-cell, corresponding to the Z-ring. Z-ring intensity was calculated by averaging the fluorescence intensity within the Z-ring (Shi et al., 2017a).

Tracking the Dynamics of Septum Formation

Assuming cells incorporate new septum area at a constant rate k and the shape of the cell poles is hemispherical, the septum area A(t) changes as a function of time according to dA(t)dt=k, and A(t)=4πw021-(w(t)/w0)2, where, w0 is initial septum width (i.e., cell width) and w(t) is the current septum width (Reshes et al., 2008a). Integration yields

t=4πw02k1-(w(t)/w0)2,

where t = 0 corresponds to division onset. Through linear fitting of experimental measurements of 1-(w(t)/w0)2 as a function of t, we obtained the time required for septum formation, T=4πw02k. To reduce the error in estimating septum width, only time points with 0.5w0w(t) ≤ 0.9w0 were used for the fit.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

JS, HS, and KCH designed the research, wrote the manuscript, and analyzedthe data. JS and HS performed the research. All authors contributed to the article and approved the submitted version.

Funding

We acknowledge funding from a James S. McDonnell Postdoctoral Fellowship (to HS) and NIH grant RM1 GM135102 (to KCH). KCH is a Chan Zuckerberg Biohub Investigator.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank members of the Huang lab for helpful discussions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.718600/full#supplementary-material

Supplementary Figure 1 | The distribution of steady-state growth rates is similar in microfluidic devices and on agarose pads.

Supplementary Figure 2 | GFP-ZapA concentration increases directly after a hyperosmotic shock. This increase was consistent with the decrease in cytoplasmic volume (Figure 1B) and the increases in FtsZ-msfGFP and cytoplasmic GFP (Figure 4C).

Supplementary Table 1 | Strains used in this study.

References

Baba, T., Ara, T., Hasegawa, M., Takai, Y., Okumura, Y., Baba, M., et al. (2006). Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol. Syst. Biol. 2:2006.0008.

Google Scholar

Barrows, J. M., and Goley, E. D. (2021). FtsZ dynamics in bacterial division: what, how, and why? Curr. Opin. Cell Biol 68, 163–172. doi: 10.1016/j.ceb.2020.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Basu, R., Munteanu, E. L., and Chang, F. (2014). Role of turgor pressure in endocytosis in fission yeast. Mol. Biol. Cell 25, 679–687. doi: 10.1091/mbc.e13-10-0618

PubMed Abstract | CrossRef Full Text | Google Scholar

Bi, E. F., and Lutkenhaus, J. (1991). FtsZ ring structure associated with division in Escherichia coli. Nature 354, 161–164. doi: 10.1038/354161a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bisson-Filho, A. W., Hsu, Y. P., Squyres, G. R., Kuru, E., Wu, F., Jukes, C., et al. (2017). Treadmilling by FtsZ filaments drives peptidoglycan synthesis and bacterial cell division. Science 355, 739–743. doi: 10.1126/science.aak9973

PubMed Abstract | CrossRef Full Text | Google Scholar

Boersma, A. J., Zuhorn, I. S., and Poolman, B. (2015). A sensor for quantification of macromolecular crowding in living cells. Nat. Methods 12, 227–229. doi: 10.1038/nmeth.3257

PubMed Abstract | CrossRef Full Text | Google Scholar

Buda, R., Liu, Y., Yang, J., Hegde, S., Stevenson, K., Bai, F., et al. (2016). Dynamics of Escherichia coli’s passive response to a sudden decrease in external osmolarity. Proc. Natl. Acad. Sci. U.S.A. 113, E5838–E5846.

Google Scholar

Buss, J., Coltharp, C., Huang, T., Pohlmeyer, C., Wang, S. C., Hatem, C., et al. (2013). In vivo organization of the FtsZ-ring by ZapA and ZapB revealed by quantitative super-resolution microscopy. Mol. Microbiol. 89, 1099–1120. doi: 10.1111/mmi.12331

PubMed Abstract | CrossRef Full Text | Google Scholar

Cagliero, C., and Jin, D. J. (2013). Dissociation and re-association of RNA polymerase with DNA during osmotic stress response in Escherichia coli. Nucleic Acids Res. 41, 315–326. doi: 10.1093/nar/gks988

PubMed Abstract | CrossRef Full Text | Google Scholar

Campos, M., Surovtsev, I. V., Kato, S., Paintdakhi, A., Beltran, B., Ebmeier, S. E., et al. (2014). A constant size extension drives bacterial cell size homeostasis. Cell 159, 1433–1446. doi: 10.1016/j.cell.2014.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Camsund, D., Lawson, M. J., Larsson, J., Jones, D., Zikrin, S., Fange, D., et al. (2020). Time-resolved imaging-based CRISPRi screening. Nat. Methods 17, 86–92. doi: 10.1038/s41592-019-0629-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Cayley, D. S., Guttman, H. J., and Record, M. T. Jr. (2000). Biophysical characterization of changes in amounts and activity of Escherichia coli cell and compartment water and turgor pressure in response to osmotic stress. Biophys. J. 78, 1748–1764. doi: 10.1016/s0006-3495(00)76726-9

CrossRef Full Text | Google Scholar

Chang, F. (2017). Forces that shape fission yeast cells. Mol. Biol. Cell 28, 1819–1824. doi: 10.1091/mbc.e16-09-0671

PubMed Abstract | CrossRef Full Text | Google Scholar

Christian, J. (1955). The influence of nutrition on the water relations of Salmonella oranienburg. Aust. J. Biol. Sci. 8, 75–82. doi: 10.1071/bi9550075

CrossRef Full Text | Google Scholar

Christian, J., and Scott, W. (1953). Water relations of Salmonellae at 30 C. Aust. J. Biol. Sci. 6, 565–573.

Google Scholar

Cota-Robles, E. H. (1963). Electron Microscopy of Plasmolysis in Escherichia Coli. J. Bacteriol. 85, 499–503. doi: 10.1128/jb.85.3.499-503.1963

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, K., and Lutkenhaus, J. (1991). ftsZ is an essential cell division gene in Escherichia coli. J. Bacteriol. 173, 3500–3506. doi: 10.1128/jb.173.11.3500-3506.1991

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Y., Sun, M., and Shaevitz, J. W. (2011). Direct measurement of cell wall stress stiffening and turgor pressure in live bacterial cells. Phys. Rev. Lett. 107:158101.

Google Scholar

Edelstein, A., Amodaj, N., Hoover, K., Vale, R., and Stuurman, N. (2010). Computer control of microscopes using microManager. Curr. Protoc. Mol. Biol. 92, 14.20.1–14.20.17.

Google Scholar

Erickson, H. P. (2009). Modeling the physics of FtsZ assembly and force generation. Proc. Natl. Acad. Sci. U.S.A. 106, 9238–9243. doi: 10.1073/pnas.0902258106

PubMed Abstract | CrossRef Full Text | Google Scholar

Erickson, H. P. (2017). How bacterial cell division might cheat turgor pressure–a unified mechanism of septal division in Gram-positive and Gram-negative bacteria. Bioessays 39:1700045. doi: 10.1002/bies.201700045

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, C. R., Fan, Y., and Ling, J. (2019). Increased mistranslation protects E. coli from protein misfolding stress due to activation of a RpoS-dependent heat shock response. FEBS Lett. 593, 3220–3227. doi: 10.1002/1873-3468.13578

PubMed Abstract | CrossRef Full Text | Google Scholar

Finan, J. D., and Guilak, F. (2010). The effects of osmotic stress on the structure and function of the cell nucleus. J. Cell. Biochem. 109, 460–467.

Google Scholar

Gauthier, M. J. (2000). “Environmental parameters associated with the viable but nonculturable state,” in Nonculturable Microorganisms in the Environment, eds R. R. Colwell and D. J. Grimes (Boston, MA: Springer), 87–112. doi: 10.1007/978-1-4757-0271-2_7

CrossRef Full Text | Google Scholar

Goley, E. D., Yeh, Y. C., Hong, S. H., Fero, M. J., Abeliuk, E., Mcadams, H. H., et al. (2011). Assembly of the Caulobacter cell division machine. Mol. Microbiol. 80, 1680–1698. doi: 10.1111/j.1365-2958.2011.07677.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, N. S., Kadoya, R., Chattoraj, D. K., and Levin, P. A. (2012). Cell size and the initiation of DNA replication in bacteria. PLoS Genet. 8:e1002549. doi: 10.1371/journal.pgen.1002549

PubMed Abstract | CrossRef Full Text | Google Scholar

Holtje, J. V. (1998). Growth of the stress-bearing and shape-maintaining murein sacculus of Escherichia coli. Microbiol. Mol. Biol. Rev. 62, 181–203. doi: 10.1128/mmbr.62.1.181-203.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuwada, N. J., Traxler, B., and Wiggins, P. A. (2015). Genome-scale quantitative characterization of bacterial protein localization dynamics throughout the cell cycle. Mol. Microbiol. 95, 64–79. doi: 10.1111/mmi.12841

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitosch, K., Rieckh, G., and Bollenbach, T. (2017). Noisy response to antibiotic stress predicts subsequent single-cell survival in an acidic environment. Cell Syst. 4:e395.

Google Scholar

Moore, D. A., Whatley, Z. N., Joshi, C. P., Osawa, M., and Erickson, H. P. (2017). Probing for binding regions of the FtsZ protein surface through site-directed insertions: discovery of fully functional FtsZ-fluorescent proteins. J. Bacteriol. 199:e553–16.

Google Scholar

Oliveira, R. A., Ng, K. M., Correia, M. B., Cabral, V., Shi, H., Sonnenburg, J. L., et al. (2020). Klebsiella michiganensis transmission enhances resistance to Enterobacteriaceae gut invasion by nutrition competition. Nat. Microbiol. 5, 630–641. doi: 10.1038/s41564-019-0658-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Osawa, M., and Erickson, H. P. (2013). Liposome division by a simple bacterial division machinery. Proc. Natl. Acad. Sci. U.S.A. 110, 11000–11004. doi: 10.1073/pnas.1222254110

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, J. M., Colavin, A., Shi, H., Czarny, T. L., Larson, M. H., Wong, S., et al. (2016). A comprehensize, CRISPR-based approach to functional analysis of essential genes in bacteria. Cell 165, 1493–1506. doi: 10.1016/j.cell.2016.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Pilizota, T., and Shaevitz, J. W. (2012). Fast, multiphase volume adaptation to hyperosmotic shock by Escherichia coli. PLoS One 7:e35205. doi: 10.1371/journal.pone.0035205

PubMed Abstract | CrossRef Full Text | Google Scholar

Pilizota, T., and Shaevitz, J. W. (2013). Plasmolysis and cell shape depend on solute outer-membrane permeability during hyperosmotic shock in E. coli. Biophys. J. 104, 2733–2742. doi: 10.1016/j.bpj.2013.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Reshes, G., Vanounou, S., Fishov, I., and Feingold, M. (2008a). Cell shape dynamics in Escherichia coli. Biophys. J. 94, 251–264. doi: 10.1529/biophysj.107.104398

PubMed Abstract | CrossRef Full Text | Google Scholar

Reshes, G., Vanounou, S., Fishov, I., and Feingold, M. (2008b). Timing the start of division in E. coli: a single-cell study. Phys. Biol. 5:046001. doi: 10.1088/1478-3975/5/4/046001

CrossRef Full Text | Google Scholar

Rojas, E., Theriot, J. A., and Huang, K. C. (2014). Response of Escherichia coli growth rate to osmotic shock. Proc. Natl. Acad. Sci. U.S.A. 111, 7807–7812. doi: 10.1073/pnas.1402591111

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojas, E. R., Billings, G., Odermatt, P. D., Auer, G. K., Zhu, L., Miguel, A., et al. (2018). The outer membrane is an essential load-bearing element in Gram-negative bacteria. Nature 559, 617–621. doi: 10.1038/s41586-018-0344-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojas, E. R., and Huang, K. C. (2018). Regulation of microbial growth by turgor pressure. Curr. Opin. Microbiol. 42, 62–70. doi: 10.1016/j.mib.2017.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojas, E. R., Huang, K. C., and Theriot, J. A. (2017). Homeostatic cell growth is accomplished mechanically through membrane tension inhibition of cell-wall synthesis. Cell. Syst. 5:e576.

Google Scholar

Scott, W. (1953). Water relations of Staphylococcus aureus at 30 C. Aust. J. Biol. Sci. 6, 549–564. doi: 10.1071/bi9530549

CrossRef Full Text | Google Scholar

Sekar, K., Rusconi, R., Sauls, J. T., Fuhrer, T., Noor, E., Nguyen, J., et al. (2018). Synthesis and degradation of FtsZ quantitatively predict the first cell division in starved bacteria. Mol. Syst. Biol. 14:e8623.

Google Scholar

Shi, H., Colavin, A., Bigos, M., Tropini, C., Monds, R. D., and Huang, K. C. (2017a). Deep phenotypic mapping of bacterial cytoskeletal mutants reveals physiological robustness to cell size. Curr. Biol. 27, 3419–3429. doi: 10.1016/j.cub.2017.09.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, H., Colavin, A., Lee, T. K., and Huang, K. C. (2017b). Strain Library Imaging Protocol: high-throughput, automated single-cell microscopy for large bacterial collections arrayed on multiwell plates. Nat. Protoc. 12, 429–438. doi: 10.1038/nprot.2016.181

PubMed Abstract | CrossRef Full Text | Google Scholar

Si, F., Le Treut, G., Sauls, J. T., Vadia, S., Levin, P. A., and Jun, S. (2019). Mechanistic origin of cell-size control and homeostasis in bacteria. Curr. Biol. 29, 1760.–1770.

Google Scholar

Sliusarenko, O., Heinritz, J., Emonet, T., and Jacobs-Wagner, C. (2011). High-throughput, subpixel precision analysis of bacterial morphogenesis and intracellular spatio-temporal dynamics. Mol. Microbiol. 80, 612–627. doi: 10.1111/j.1365-2958.2011.07579.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Spratt, B. G. (1975). Distinct penicillin binding proteins involved in the division, elongation, and shape of Escherichia coli K12. Proc. Natl. Acad. Sci. U.S.A. 72, 2999–3003. doi: 10.1073/pnas.72.8.2999

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutterlin, H. A., Shi, H., May, K. L., Miguel, A., Khare, S., Huang, K. C., et al. (2016). Disruption of lipid homeostasis in the Gram-negative cell envelope activates a novel cell death pathway. Proc. Natl. Acad. Sci. U.S.A. 113, E1565–E1574.

Google Scholar

Szatmári, D., Sárkány, P., Kocsis, B., Nagy, T., Miseta, A., Barkó, S., et al. (2020). Intracellular ion concentrations and cation-dependent remodelling of bacterial MreB assemblies. Sci. Rep. 10:12002.

Google Scholar

Tagkopoulos, I., Liu, Y. C., and Tavazoie, S. (2008). Predictive behavior within microbial genetic networks. Science 320, 1313–1317. doi: 10.1126/science.1154456

PubMed Abstract | CrossRef Full Text | Google Scholar

Taheri-Araghi, S., Bradde, S., Sauls, J. T., Hill, N. S., Levin, P. A., Paulsson, J., et al. (2015). Cell-size control and homeostasis in bacteria. Curr. Biol. 25, 385–391. doi: 10.1016/j.cub.2014.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Taniguchi, Y., Choi, P. J., Li, G. W., Chen, H., Babu, M., Hearn, J., et al. (2010). Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science 329, 533–538. doi: 10.1126/science.1188308

PubMed Abstract | CrossRef Full Text | Google Scholar

Ursell, T., Lee, T. K., Shiomi, D., Shi, H., Tropini, C., Monds, R. D., et al. (2017). Rapid, precise quantification of bacterial cellular dimensions across a genomic-scale knockout library. BMC Biol. 15:17. doi: 10.1186/s12915-017-0348-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ursell, T. S., Nguyen, J., Monds, R. D., Colavin, A., Billings, G., Ouzounov, N., et al. (2014). Rod-like bacterial shape is maintained by feedback between cell curvature and cytoskeletal localization. Proc. Natl. Acad. Sci. U.S.A. 111, E1025–E1034.

Google Scholar

van den Bogaart, G., Hermans, N., Krasnikov, V., and Poolman, B. (2007). Protein mobility and diffusive barriers in Escherichia coli: consequences of osmotic stress. Mol. Microbiol. 64, 858–871. doi: 10.1111/j.1365-2958.2007.05705.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Walter, H. (1924). Plasmaquellung und wachstum. Z. Bot. 16:1931.

Google Scholar

Wang, P., Robert, L., Pelletier, J., Dang, W. L., Taddei, F., Wright, A., et al. (2010). Robust growth of Escherichia coli. Curr. Biol. 20, 1099–1103.

Google Scholar

Ward, J. E. Jr., and Lutkenhaus, J. (1985). Overproduction of FtsZ induces minicell formation in E. coli. Cell 42, 941–949. doi: 10.1016/0092-8674(85)90290-9

CrossRef Full Text | Google Scholar

Weart, R. B., Lee, A. H., Chien, A. C., Haeusser, D. P., Hill, N. S., and Levin, P. A. (2007). A metabolic sensor governing cell size in bacteria. Cell 130, 335–347. doi: 10.1016/j.cell.2007.05.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Weart, R. B., and Levin, P. A. (2003). Growth rate-dependent regulation of medial FtsZ ring formation. J. Bacteriol. 185, 2826–2834. doi: 10.1128/jb.185.9.2826-2834.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Swain, P., Kuijpers, L., Zheng, X., Felter, K., Guurink, M., et al. (2019). Cell boundary confinement sets the size and position of the E. coli chromosome. Curr. Biol. 29, 2131–2144.e4.

Google Scholar

Yang, D., Männik, J., Retterer, S. T., and Männik, J. (2020). The effects of polydisperse crowders on the compaction of the Escherichia coli nucleoid. Mol. Microbiol. 113, 1022–1037. doi: 10.1111/mmi.14467

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Lyu, Z., Miguel, A., Mcquillen, R., Huang, K. C., and Xiao, J. (2017). GTPase activity-coupled treadmilling of the bacterial tubulin FtsZ organizes septal cell wall synthesis. Science 355, 744–747. doi: 10.1126/science.aak9995

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, X., Jericho, M., Pink, D., and Beveridge, T. (1999). Thickness and elasticity of Gram-negative murein sacculi measured by atomic force microscopy. J. Bacteriol. 181, 6865–6875. doi: 10.1128/jb.181.22.6865-6875.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Halladin, D. K., Rojas, E. R., Koslover, E. F., Lee, T. K., Huang, K. C., et al. (2015). Bacterial division. Mechanical crack propagation drives millisecond daughter cell separation in Staphylococcus aureus. Science 348, 574–578. doi: 10.1126/science.aaa1511

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: FtsZ, divisome, cytokinesis, sucrose, osmolarity, microfluidics

Citation: Sun J, Shi H and Huang KC (2021) Hyperosmotic Shock Transiently Accelerates Constriction Rate in Escherichia coli. Front. Microbiol. 12:718600. doi: 10.3389/fmicb.2021.718600

Received: 01 June 2021; Accepted: 26 July 2021;
Published: 13 August 2021.

Edited by:

Martin Loose, Institute of Science and Technology Austria (IST Austria), Austria

Reviewed by:

Jaan Männik, The University of Tennessee, Knoxville, United States
Ramanujam Srinivasan, National Institute of Science Education and Research (NISER), India

Copyright © 2021 Sun, Shi and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kerwyn Casey Huang, kchuang@stanford.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.