Microﬂuidic-Enabled Multi-Cell-Densities-Patterning and Culture Device for Characterization of Yeast Strains’ Growth Rates under Mating Pheromone

: Yeast studies usually focus on exploring diversity in terms of a speciﬁc trait (such as growth rate, antibiotic resistance, or fertility) among extensive strains. Microﬂuidic chips improve these biological studies in a manner of high throughput and high efﬁciency. For a population study of yeast, it is of great signiﬁcance to set a proper initial cell density for every strain under speciﬁc circumstances. Herein, we introduced a novel design of chip, which enables users to load cells in a gradient order (six alternatives) of initial cell density within one channel. We discussed several guidelines to choose the appropriate chamber to ensure successful data recording. With this chip, we successfully studied the growth rate of yeast strains under a mating response, which is crucial for yeasts to control growth behaviors for prosperous mating. We investigated the growth rate of eight different yeast strains under three different mating pheromone levels (0.3 µ M, 1 µ M, and 10 µ M). Strains with, even, a six-fold in growth rate can be recorded, with the available data produced simultaneously. This work has provided an efﬁcient and time-saving microﬂuidic platform, which enables loading cells in a pattern of multi-cell densities for a yeast population experiment, especially for a high-throughput study. Besides, a quantitatively analyzed growth rate of different yeast strains shall reveal inspiring perspectives for studies concerning yeast population behavior with a stimulated mating pheromone.


Introduction
The budding yeast S. cerevisiae has long been studied as a model organism. In contrast with gathering information from a single cell [1][2][3][4][5][6][7], interest in cell behaviors at the population scale have been rising recently [8][9][10][11][12][13][14][15], in which precise measurement of the growth rate has become a considerable and fundamental task. Microfluidic chips are often used to study questions about yeasts or other organisms [16][17][18][19][20][21][22][23][24][25]. However, measuring the growth rate of various strains with different genes knocked-out meets several problems (or have several requirements): (1) The difference in growth rate, as counted by cell number between the fastest-growth strain and the slowest one under diverse conditions, is extremely dramatic, when both size and shape of observation chamber is uniform for each strain. If initial cell density is too high, yeast cells will grow out of the visual field for each of the selected measuring points at the end of the experiment. Thus, initial cell density of every strain is of great importance when given a fixed measuring time under all conditions; (2) Morphology of yeast cells may change a lot when facing stimulus, which means the number of cells or optical density of cell suspension are not sufficient to give information about growth status of yeast cells; (3) Manually added inducer (mostly small molecules) may be degradable by specific enzymes, such as protease secreted by yeast cells, and, thus, is hard to influence yeast cells, peculiarly when the concentration of this factor is low and cell density is high, which causes an unstable biochemical environment. To solve similar problems, as mentioned above, a series of microfluidic chips have been developed. Unfortunately, many of them [17][18][19]24] are designed for fetching data from a single cell, such that they may only have one observation for each strain under a specific condition (thus, they cannot adjust loaded cell density and do not fulfill requirement (1), or exclusively satisfied the requirements of (2) and (3)). In this paper, we are eager to introduce a newly designed microfluidic chip, which conforms to all these requirements, especially by providing users multiple choices for proper initial cell density, for specific experiment requirements.
Mating is a vital biological pathway, which regulates yeasts behaviors [26][27][28][29], so dissecting the growth rate of yeasts with extensive response intensities to a mating pheromone is quite captivating. Haploid yeasts include two mating types: MATa and MATα. Mating pheromone a-factor and α-factor are small peptides naturally secreted by MATa yeasts and MATα yeasts for mediating mating behaviors, respectively [26,30]. Yeast study focusing on mating usually applies MATa haploid yeast and α-factor. Since MATa cell secretes Bar1 protein, a protease that can cleave and inactivate α-factor, requirement (3) is a pivotal condition to be satisfied. What is more, yeasts in response to mating pheromone would encounter cell cycle arrest. During cell cycle arrest, yeasts exhibit morphological changes mainly in two forms: elongation and shmoo [31]. As a previous study [32][33][34][35] pointed out, elongation appears where the mating pheromone level is low, when downstream proteins Fus3 and Kss1 take part, while shmoo happens only when the mating pheromone is high, when only Fus3 participates. Elongated yeast is like a line, while shmoo is like a polygon with every side having curvature. Although yeasts stop cell division during cell cycle arrest, the size of the individual continues increasing, undetectable by traditional methods that merely count cell number. This indicates the significance of requirement (2). Both of these two requirements have been successfully achieved by our microfluidic device in previous work [18], combined with a microscope to observe and determine the mass growth rate of GFP-labeled strains (which are all wild type) under two mating pheromone concentrations.
To solve the challenge of requirement (1), additionally, we introduced a newly designed microfluidic chip specialized for yeast population study in this paper. Chip design fulfills self-distribution of cell density in a gradient order, which results in multiple choices for specific initial cell density, especially suitable for high-throughput study. On the other hand, we established a set of quantitative methods for accurate analyzation of yeast growth rate with large diversity. In this study, we simultaneously set different α-factor levels (0.3 µM, 1 µM and 10 µM) besides the control group (0 µM) and compared the growth rate of the wild-type strain and eight different knock-out strains, whose genes were related with the cell cycle, cytoskeleton, protein synthesis, and mating response pathway. Based on our analyzation, we discussed a series of yeast's growth status under the environments containing α-factor. We expect that this platform might inspire microfluidic design for cell population study and provide guidance to quantitative and high-throughput exploration of the yeast growth rate in the future.

System Setup and Microfluidic Chip Design
Our experimental platform is shown in Figure 1a, which consists of four parts: (1) a novel designed microfluidic chip for this study; (2) syringe pumps for controlling culturing medium with 4 different concentrations of yeast mating pheromone α-factor; (3) inverted microscope imaging system (Nikon Instruments Inc., Melville, NY, USA) for capturing phase contrast images of yeasts; and (4) a computer in charge of controlling syringe pumps Chemosensors 2022, 10, 141 3 of 16 and microscope system. When we were designing the chip in this study, the very first purpose was that researchers can choose specific chambers with an expected initial cell density from one set of observation chambers. Based on such reason, we ingeniously designed the structure of microfluidic chip for yeast experiments (Figure 1b). The overview of our chip displayed loading wells (black dotted box) and medium inlets (brown dotted box). The medium outlet is the same as the loading well, since we load yeasts through these wells into the observation chambers, while during the experiment, the medium was injected through the medium inlets into chambers in a contrary direction, so that the waste came out from the same wells as loading wells. In observation region (blue dotted box) of this chip, there are four identical areas for setting different medium input, each of which includes six parallel observation regions. Current design realizes observation of 6 different strains in 4 different culture conditions, and it has the capability to extend to higher throughput.

System Setup and Microfluidic Chip Design
Our experimental platform is shown in Figure 1a, which consists of four parts: (1) a novel designed microfluidic chip for this study; (2) syringe pumps for controlling culturing medium with 4 different concentrations of yeast mating pheromone α-factor; (3) inverted microscope imaging system (Nikon Instruments Inc., Melville, NY, USA) for capturing phase contrast images of yeasts; and (4) a computer in charge of controlling syringe pumps and microscope system. When we were designing the chip in this study, the very first purpose was that researchers can choose specific chambers with an expected initial cell density from one set of observation chambers. Based on such reason, we ingeniously designed the structure of microfluidic chip for yeast experiments (Figure 1b). The overview of our chip displayed loading wells (black dotted box) and medium inlets (brown dotted box). The medium outlet is the same as the loading well, since we load yeasts through these wells into the observation chambers, while during the experiment, the medium was injected through the medium inlets into chambers in a contrary direction, so that the waste came out from the same wells as loading wells. In observation region (blue dotted box) of this chip, there are four identical areas for setting different medium input, each of which includes six parallel observation regions. Current design realizes observation of 6 different strains in 4 different culture conditions, and it has the capability to extend to higher throughput. Overview of the chip. Different colored dotted boxes present different functional parts of the chip. Current design allows 6 different strains growing under same environment. (c) Schematic diagram of 1 set of observation chambers, which includes 6 observation chambers, fence, main channel, and by-pass channel. During preparation, the loaded cells were flowing with the medium from the right side (as the arrow on the right shows). For each strain, there are six chambers for researchers to choose for an appropriate initial cell density. Circled number 1-6 marks a gradient distribution of initial cell density. During experiment, the culture medium would be injected by the syringe pump from the left side as the 7 arrows on the left shows. The various heights of each region were denoted by different colors in the diagram. (d) The microfluidic chip with 2 different inks load into, placed with a ruler. (e) Inside channels and observation chambers of this microfluidic chip observed by a 4× objective.
Each observation channel comprises of 6 chambers, which were presented in Figure  1c. A single observation chamber is a 120 μm-wide square with a height of 4 μm that is Each observation channel comprises of 6 chambers, which were presented in Figure 1c. A single observation chamber is a 120 µm-wide square with a height of 4 µm that is slightly smaller than the diameter of yeast cells, so that yeast cells in loading medium can be trapped when entering the observation chamber form main channel. We appended the by-pass channel within every observation region, which helps provide both a stable and relatively slow flow velocity of the culture medium. The height of by-pass channel is no difference from the main channel. At the distal end of the observation area, there are fences (h = 1.7 µm), which prevent the cells from flowing into the other side and causing cross-contamination among different strains. Compared with previous chips [18,24], the most significant advantage of this chip is that it has a multi-chamber observation region, which enables gradient distribution of the number of one strain within 6 chambers. It extremely increased the possibility of loading most suitable initial cell density required for experiments.
To get the microfluidic chip for this study, we first made 3 chrome plate masks, which corresponded to three layers with different heights (fence, observation chamber, main channel) of the chip. We then used mask aligners to exposure each layer with mercury lamp. The photoresist we used to make the mold is the SU-8 3000 series. After exposure of each layer, we washed the unexposed area and began next turn of exposure. When we finished the chip mold, we poured well-mixed PDMS (Momentive RTV615, A: B = 8:1) on the mold. After degassed, the PDMS along with mold were placed in air oven at 70 • C overnight. The second day, we removed the bonded PDMS from the mold and punched all the inlets or outlets on the chip. At last, we used the plasma cleaner to attach the PDMS part of the chip with the cover glass part. Figure 2d shows the final product of the chip was placed along with a ruler to make acknowledge of the chip size. The 2 different colors of ink were injected into 4 separate areas of the chip, which enables 4 different culture conditions during experiment. Moreover, Figure 2e exhibits the detail of main channels and observation chambers of the chip observed by a 4× objective.
be trapped when entering the observation chamber form main channel. We appended the by-pass channel within every observation region, which helps provide both a stable and relatively slow flow velocity of the culture medium. The height of by-pass channel is no difference from the main channel. At the distal end of the observation area, there are fences (h = 1.7 μm), which prevent the cells from flowing into the other side and causing cross-contamination among different strains. Compared with previous chips [18,24], the most significant advantage of this chip is that it has a multi-chamber observation region, which enables gradient distribution of the number of one strain within 6 chambers. It extremely increased the possibility of loading most suitable initial cell density required for experiments.
To get the microfluidic chip for this study, we first made 3 chrome plate masks, which corresponded to three layers with different heights (fence, observation chamber, main channel) of the chip. We then used mask aligners to exposure each layer with mercury lamp. The photoresist we used to make the mold is the SU-8 3000 series. After exposure of each layer, we washed the unexposed area and began next turn of exposure. When we finished the chip mold, we poured well-mixed PDMS (Momentive RTV615, A: B = 8:1) on the mold. After degassed, the PDMS along with mold were placed in air oven at 70 °C overnight. The second day, we removed the bonded PDMS from the mold and punched all the inlets or outlets on the chip. At last, we used the plasma cleaner to attach the PDMS part of the chip with the cover glass part. Figure 2d shows the final product of the chip was placed along with a ruler to make acknowledge of the chip size. The 2 different colors of ink were injected into 4 separate areas of the chip, which enables 4 different culture conditions during experiment. Moreover, Figure 2e exhibits the detail of main channels and observation chambers of the chip observed by a 4× objective.  The principle of gradient density distribution is illustrated as follows. Each observation chamber has various distance from the entrance of the cell loading wells, so that it allows decreasing flow velocity within the observation chambers with increasing distance from loading entrance when loading cells. To confirm this tentativeness, we employed COMSOL to simulate several cases when using the chip ( Figure 2). Figure 2a lists several pieces of information about our simulation process, including software detail and fluid velocity in different cases. Figure 2b illustrates the status when loading cells. As we consider a situation where there are relatively few cells, we made this simplification that we only simulate the condition where there is only the medium inside this chip. In this figure, the medium-loading direction is from the right. As shown in Figure 2b, medium velocity increases gradually from chamber No. 1 (most blue, low velocity) to chamber No. 6 (most yellow, high velocity). This is a good proof, as expected.
Then we determined to simulate status of culturing yeasts with by-pass channel during experiment, as Figure 2c exhibits. The medium's incoming direction is from the left. In Figure 2c, it is clear that the medium velocity in almost all 6 chambers were quite low (deep blue). Moreover, velocity within different chambers is alike, which means there's little bias on medium velocity among different chamber during experiment, which enhanced reliability to experiment result. To verify the contribution of the by-pass channel for decreasing velocity in the chamber, we continued simulating the condition of culturing yeasts without the by-pass channel ( Figure 2d). We then found that the velocity within the chambers was higher than the situation with the by-pass channel. Compared with Figure 2c, we can also distinguish that the velocity among each chamber is no longer uniform. The by-pass channel brings a more uniform and slower velocity, theoretically. Thus, such a design enables us to select chambers with the expected initial cell density as observation points, according to specific experiment requirements.

Selection of Yeast Strains
The mating pheromone response pathway [36][37][38][39] is activated by α-factor, and MATa haploid cells undergo cell cycle arrest. Figure 3a shows how this pathway works: when attached by the α-factor, the Ste2 protein is activated. Ste4 protein (Gβ subunit), accompanying with Ste18 (Gγ subunit), releases from Gpa1 (Gα subunit) and starts the following activation steps of mating pathway until the activation of Fus3. Ste5 protein plays as scaffold which provides a platform for MAPK kinases and other proteins. After activated by Cdc42, Ste20 starts sequential phosphorylation of Ste11, Ste7, and Fus3 until turn-on of downstream proteins such as Ste12. Therefore, we chose genes STE4 and STE5, which encode 2 important components, proteins Ste4 and Ste5, located in the mating pathway. We then utilized the corresponding single gene knock-out strains to observe their behaviors under different α-factor levels.
Yeast growth is affected by absence of certain genes [17]. Without these genes, yeasts undergo abnormal growth. For example, SWE1, YKU70 is related with cell cycle, LOC1 is related with translation, and RVS161, CIN8 is related with actin cytoskeleton, which is involved with morphological changes resulted from cell cycle arrest [40][41][42]. Thus, we tried to figure out how different yeasts without these genes behaved from the wild type strain, under a different mating pheromone level. We finalized the genes listed above and chose corresponding single-gene knock-out strains from the yeast strain library.
It has also been reported that removing some key components of such mating response pathway results in desensitization of yeasts to mating pheromone [43], contributing to growth changes. On the other hand, lack of LOC1 made yeast vegetate at lower growth rate. It is fascinating to find out how absence of genes directly controlling growth such as LOC1 interacted with deletion of genes located at mating pathway such as STE4 under various α-factor level. Therefore, we innovatively constructed a double-gene-knocked-out strain ste4-∆/loc1-∆ to investigate how growth status influenced by these two pathways under different intensities of mating pheromone.  Yeast growth is affected by absence of certain genes [17]. Without these genes, yeasts undergo abnormal growth. For example, SWE1, YKU70 is related with cell cycle, LOC1 is related with translation, and RVS161, CIN8 is related with actin cytoskeleton, which is involved with morphological changes resulted from cell cycle arrest [40][41][42]. Thus, we tried to figure out how different yeasts without these genes behaved from the wild type strain, under a different mating pheromone level. We finalized the genes listed above and chose corresponding single-gene knock-out strains from the yeast strain library.
It has also been reported that removing some key components of such mating response pathway results in desensitization of yeasts to mating pheromone [43], contributing to growth changes. On the other hand, lack of LOC1 made yeast vegetate at lower growth rate. It is fascinating to find out how absence of genes directly controlling growth such as LOC1 interacted with deletion of genes located at mating pathway such as STE4 under various α-factor level. Therefore, we innovatively constructed a double-geneknocked-out strain ste4-Δ/loc1-Δ to investigate how growth status influenced by these two pathways under different intensities of mating pheromone.
These thoughts promoted us to choose strains listed in Figure 3b. The haploid yeast cell of the single-gene knock-out used in this study is from the deletion yeast strain library [44] provided by Jef D. Boeke at Johns Hopkins University, the wild type strain came form GFP-tagged-protein yeast library [45] provided by Dr. Erin O' Shea and Dr. Jonathan Weissman at University of California, San Francisco. Besides, we constructed dual-gene knock-out strain (ste4-Δ/loc1-Δ) by LiAc transformation method [46], with a knock-out plasmid containing a His-auxotroph selection site, based on the single-gene knocked-out ste4-Δ strain. These thoughts promoted us to choose strains listed in Figure 3b. The haploid yeast cell of the single-gene knock-out used in this study is from the deletion yeast strain library [44] provided by Jef D. Boeke at Johns Hopkins University, the wild type strain came form GFP-tagged-protein yeast library [45] provided by Dr. Erin O' Shea and Dr. Jonathan Weissman at University of California, San Francisco. Besides, we constructed dual-gene knock-out strain (ste4-∆/loc1-∆) by LiAc transformation method [46], with a knock-out plasmid containing a His-auxotroph selection site, based on the single-gene knocked-out ste4-∆ strain.

Experiment Workflow and Image Acquisition
Before executing the experiment, we inoculated the strains in synthetic complete (SC) medium and cultured them at 30 • C overnight. On next day, the cell suspensions were diluted into fresh SC medium at a ratio of 5%. After culturing about 5-6 h, we injected cell suspension into chip. Before loading cells, we used a vacuum to degas the microfluidic chip for at least 15 min. Then, we used a 10-µL pipette to load the yeast cells directly from the loading wells into the chip. To make yeasts adapt to new environment, we continued culturing all strains with medium containing no α-factor about 3 h before imaging.
When culture medium switched to experimental condition (0.3 µM, 1 µM, 10 µM), we used a computer system to automatically control Nikon Ti-E inverted microscope and syringe pump. By setting up customized acquisition condition, we were able to capture images of all the selected positions with a fixed time interval (10 min) in an entire loop (10 h). Meanwhile, the temperature around the chip was maintained at 30 • C, by a live cell station. A phase contrast objective and a corresponding phase condenser annulus 3 (Ph3) were employed for yeast imaging, since the phase difference images provide unambiguous periphery of a single cell, specifically the intracellular and extracellular region that appears dark, whereas the margin region appears bright.

Validation on Function of Loading Yeasts with a Multi-Cell-Density Pattern
To verify the novel functions of our chip, we conducted several tests on loading yeasts, to observe whether the initial cell distribution was loaded with a pattern of multicell-densities, as expected. Figure 4a provided here is one ideal example. There were 6 numbered chambers in Figure 4a. It is obvious that chamber No. 1 had the lowest cell density, while chamber No. 6 had the highest cell density, which was consistent with the simulation result in tendency.  , which corresponds to the increasing distance for these chambers to loading well. (b) An ideal example for selecting chambers for successfully recording available data. In this situation, yeasts will not grow out of the chamber at the end of the experiment. (c) One of the situations that will not be selected where initial cell density is too high, so that yeasts grow out of the chamber when the experiment ends. (d) Another situation that will not be picked, where yeasts stay at the edge of the chamber after being loaded. Yeasts will grow directly into the main channel during the experiment. The left half of (b-d) is the schematic diagram, and the other half is from loading trials.

Growth Rate Analyzation on the Wild Type Strain
Cell number and area are crucial parameters in growth rate calculation. In this study, our procedure of analyzing the cellular growth rate included three steps: (1) Select the suitable cell populations within the microfluidic observation chamber, which should have a clear boundary and maintain a visual field during the whole experiment period (and also follow the 3 rules discussed above); (2) Utilize ImageJ build-in tools to measure the area and cell number of all selected populations in terms of 4 time points (0 h, 1.5 h, 3 h, and 6 h); and (3) Perform curve fitting by MATLAB to illustrate the cell growth dynamic in terms of the population's area and number. The relationship between parameters is calculated as follows:  6, which corresponds to the increasing distance for these chambers to loading well. (b) An ideal example for selecting chambers for successfully recording available data. In this situation, yeasts will not grow out of the chamber at the end of the experiment. (c) One of the situations that will not be selected where initial cell density is too high, so that yeasts grow out of the chamber when the experiment ends. (d) Another situation that will not be picked, where yeasts stay at the edge of the chamber after being loaded. Yeasts will grow directly into the main channel during the experiment. The left half of (b-d) is the schematic diagram, and the other half is from loading trials.
We were trying to summarize several guidelines for choosing a chamber with suitable cell density, to finally obtain the available data. Figure 4b-d exhibited our points about this topic. Figure 4b shows Rule No. 1: choose chambers with cells located in the center or near the center. Additionally, one may choose chambers with relatively low cell density, if the strain loaded in the chamber has a known large growth rate. Figure 4c exhibits the situation that we should avoid, when cells are too much located (Rule No. 2). Chambers with lower cell density had a higher possibility to avoid failure by holding all the populations, so that it shall increase the possibility of obtaining perfect raw data. Figure 4d illustrates another circumstance that we should be aware of (Rule No. 3), which is that one should be careful to choose chambers with the yeast cell staying at the edge of the chamber, especially at the boundary between the inner edge and the main channel. In this case, it is very possible that yeasts grow into the main channel, so that the new-born yeasts will be washed away through the main channel.

Growth Rate Analyzation on the Wild Type Strain
Cell number and area are crucial parameters in growth rate calculation. In this study, our procedure of analyzing the cellular growth rate included three steps: (1) Select the suitable cell populations within the microfluidic observation chamber, which should have a clear boundary and maintain a visual field during the whole experiment period (and also follow the 3 rules discussed above); (2) Utilize ImageJ build-in tools to measure the area and cell number of all selected populations in terms of 4 time points (0 h, 1.5 h, 3 h, and 6 h); and (3) Perform curve fitting by MATLAB to illustrate the cell growth dynamic in terms of the population's area and number. The relationship between parameters is calculated as follows: Here, GRa is growth rate, resulted from fitting population area. GRn is growth rate calculated from fitting population number. t is time, A represents area of population, and N is on behalf of population number. A 0 and N 0 are the initial values of the fitting process. The larger GRa or GRn is, the faster the area or number of the population changes. In the following analyzation and comparison, we mainly used GRa and GRn to estimate the growth status of strains.
For the wild type strain (Figure 5a), it is obvious that yeasts were growing well when there's no α-factor in the culture medium. At 1 µM, the yeasts had noticeable morphological changes (elongation) [32][33][34][35] at 6 h, which had no significant contradictions with other works. At 0.3 µM, cells had nearly no morphological change; while at 10 µM, yeasts exhibited the shape of shmoo [33], and the morphological change was very severe, which is a typical trait when a yeast is living in an extremely high level of mating pheromone ( Figure S1). Figure  10 μM, 6.090 × 10 −8 ). At 1 μM, yeast populations suffered from growth arrest, which can be reflected by the lower GRn and GRa (GRn = 0.146, GRa = 0.259). The reason why the decreasing of GRa is not as obvious as GRn is that yeasts in the status of cell cycle arrest continue increasing individuals' area. At 0.3 μM, we obtained GRn (0.280) intermediate between the results of 1 μM (0.146) and 0 μM (0.365); at 10 μM, yeast populations suffered from the most severe cell cycle arrest, with the lowest GRn (0.0964) and GRa (0.192). We noticed that the GRa of the four conditions was larger than the GRn (Figure 5c). To investigate the difference between the GRa and GRn change rates, we used formula (2) to calculate the change of GRa (or GRn), as the experimental concentration increased per 1 μM named "drop rate per μM (DRPU)". In this formula, the subscript of every GR or We noticed that the GRa of the four conditions was larger than the GRn (Figure 5c). To investigate the difference between the GRa and GRn change rates, we used formula (2) to calculate the change of GRa (or GRn), as the experimental concentration increased per 1 µM named "drop rate per µM (DRPU)". In this formula, the subscript of every GR or [α-factor] ("con") represents a specific experimental condition (1 as 0 µM, 2 as 0.3 µM, 3 as 1 µM, and 4 as 10 µM). Thus, we got DRPU of GRa (DRPUa) and GRn (DRPUn), respectively.
By calculating DRPU, we were able to evaluate the change rate of GRa and GRn between every neighboring concentration ("0 µM to 0.3 µM", "0.3 µM to 1 µM", and "1 µM to 10 µM") ( Figure 5d). As a result, DRPUa is smaller than DRPUn, except for "1 µM to 10 µM". Such variation between DRPUa and DRPUn indicated that the change rate of population number is larger than the rate of population area. The reason may be that area change mostly results from the accumulation of biomass, no matter if the cell is divided or not, while number change is mainly influenced by division. When undergoing cell cycle arrest, the size of single yeasts, rather than the number, increases what it may spend in energy on the synthesis of proteins relating to the cytoskeleton rather than proliferation. Thus, we saw a larger value of DRPUa than DRPUn at "0 µM to 0.3 µM" and "0.3 µM to 1 µM", which means more sensitivity to number than area, in response to the α-factor. Analyzation of the wild type strain demonstrates that our analyzation method is suitable to estimate the growth status of other yeast cells with gene knock-out.

Growth Rate Analyzation on Strains with Growth-Related Genes Absent
Besides related genes, cellular growth is also affected by a mating pheromone. To find out the effect of a mating pheromone on the growth status of growth-gene knocked-out strains, we tested several strains, including swe1-∆ and loc1-∆ strains. Figure 6a is images captured for these 2 strains, under stimulated 0 µM and 1 µM α-factor. By observing theses images, we found swe1-∆ and loc1-∆ grew just like the wild type, at 0 µM. At 1 µM, both strains had morphological change of elongation (6 h), like the wild type. Figure 6b shows the GRa and GRn fitting results, corresponding to the microscope images. Moreover, at 0.3 µM, both strains had elongation change, more distinctly than the wild type, while these two strains encountered shmoo at 10 µM ( Figure S2). energy on the synthesis of proteins relating to the cytoskeleton rather than proliferation Thus, we saw a larger value of DRPUa than DRPUn at "0 μM to 0.3 μM" and "0.3 μM t 1 μM", which means more sensitivity to number than area, in response to the α-factor Analyzation of the wild type strain demonstrates that our analyzation method is suitabl to estimate the growth status of other yeast cells with gene knock-out.

Growth Rate Analyzation on Strains with Growth-Related Genes Absent
Besides related genes, cellular growth is also affected by a mating pheromone. To find out the effect of a mating pheromone on the growth status of growth-gene knocked out strains, we tested several strains, including swe1-Δ and loc1-Δ strains. Figure 6a is im ages captured for these 2 strains, under stimulated 0 μM and 1 μM α-factor. By observing theses images, we found swe1-Δ and loc1-Δ grew just like the wild type, at 0 μM. At 1 μM both strains had morphological change of elongation (6 h), like the wild type. Figure 6b shows the GRa and GRn fitting results, corresponding to the microscope images. Moreo ver, at 0.3 μM, both strains had elongation change, more distinctly than the wild type while these two strains encountered shmoo at 10 μM ( Figure S2). Then, we compared the growing status of both strains with the wild type under four experiment conditions, based on GRn and GRa (Figure 6c-f). For swe1-∆, GRa (0.421, Figure 6c) is bigger than the wild type (0.342, Figure 6c) at 0 µM (p-value: 1.068 × 10 −2 ). At this concentration, the GRn of swe1-∆ had no clear difference from the wild type (Figure 6d). At 1 µM, the GRa of swe1-∆ is like the wild type (Figure 6c, p-value: 7.501 × 10 −1 ), while the GRn is visibly smaller than the wild type at 0.3 µM and 10 µM (Figure 6d; p-value, 0 µM vs. 0.3 µM: 6.850 × 10 −3 , 0 µM vs. 10 µM: 4.118 × 10 −2 ). We also noticed that even though both swe1-∆ and the wild type are undergoing cell cycle arrest, the difference in GRn between 0 µM and 0.3 µM is larger than that of the wild type (Figure 6d). This may imply that the swe1-∆ strain had a higher sensitivity to the α-factor than the wild type. Herein, the growth rate of swe1-∆ under 0 µM (GRn = 0.375, Figure 6d) is six-fold the rate of swe1-∆ under 10 µM (0.055, Figure 6d). We also analyzed the loc1-∆ strain. GRa of loc1-∆ are smaller than ones of the wild type at four conditions (Figure 6e; p-value of GRa, 0 µM: 4.654 × 10 −6 , 0.3 µM: 5.786 × 10 −5 , 1 µM: 1.282 × 10 −7 , 10 µM: 1.144 × 10 −4 ). Except for 10 µM, the GRn of loc1-∆ are smaller than the ones of the wild type at three conditions (Figure 6f; p-value of GRn, 0 µM: 2.172 × 10 −4 , 0.3 µM: 9.764 × 10 −5 , 1 µM: 2.086 × 10 −4 ). This means that strains of loc1-∆ grew more slowly than the wild type strain, even when there was mating pheromone around. Like the wild type, the GRn of both swe1-∆ and loc1-∆ were influenced more than the GRa.
Absence of SWE1 resulted in a faster growth pattern than the wild type strain, because the Swe1 protein functions in a cell cycle point to prevent early terminating of the cell cycle [47]. We assumed that the disfunction of the cell cycle checkpoint aggravates the influences of the higher mating pheromone environment. This also explains why the strain of swe1-∆ was more sensitive to the mating pheromone, compared with the wild type strain. For the loc1-∆ strain, it suffered from cell cycle arrest like the wild type, while it always grew slower than the wild type. We inferred that the gene LOC1 is related with protein translation [48], so that the absence of LOC1 may slacken the generation rate of almost all the proteins. This may be a global influence, which impacted mating-pheromone-response proteins, so that the mating pheromone response under all three experimental conditions was weakened globally, in strains with LOC1 absent.

Growth Rate Analyzation on Strains with Defective Mating Pheromone Response Pathway
Since cell cycle arrest is guaranteed by a complete mating response pathway, we wondered how gene knock-out, of genes located on this pathway, would influence the growth status of yeasts. Thus, we observed ste4-∆ strain and ste4-∆/loc1-∆ strains at 0 µM and 1 µM (Figure 7a), along with fitting results (Figure 7b). Images of these two strains are notable, since we did not find any morphological changes in any of them (Figures 7a and S3). This means these strains can still survive without STE4. It also indicated that the absence of both STE4 and LOC1 at the same time is not lethal to yeasts. Figure 6b showed fitting results of all conditions corresponding to microscope images (Figure 7a).
We then decided to compare growth rate of both strains at all four conditions. By comparing the p-value between 0 µM with the other three concentrations, we found that the ste4-∆ strain had a similar GRa and GRn, no matter if there is a mating pheromone or not (p-value of GRa, 0 µM vs. 0.3 µM: 2.419 × 10 −1 , 0 µM vs. 1 µM: 8.398 × 10 −1 , 0 µM vs. 10 µM: 2.313 × 10 −1 ; p-value of GRn, 0 µM vs. 0.3 µM: 2.647 × 10 −1 , 0 µM vs. 1 µM: 6.001 × 10 −1 , 0 µM vs. 10 µM: 4.096 × 10 −1 ), which also proved that this strain is not affected by the α-factor in terms of growth rate (Figure 7c,d). Although the ste4-∆ strain seems to grow slightly faster than the wild type at 0 µM, as reflected by the GRa of ste4-∆ (0.398) and the GRa of the wild type (0.342), this difference is not significant (p-value: 3.612 × 10 −1 ). Similarly, the difference between the GRn of the ste4-∆ (0.427) and the GRn of the wild type (0.365) at 0 µM is not significant (Figure 7d, p-value: 2.986 × 10 −1 ). Immunity to the α-factor described above may facilitate the ste4-∆ strains' proliferation, rather than responding to a mating signal. Next, we studied the ste4-∆/loc1-∆ strain. This strain also was not influenced by the α-factor and had a similar growth rate at all three conditions, except 0 µM (Figure 7e . This result is riveting, since this strain had traits like the ste4-∆ strain and the loc1-∆ strain, kind of an overlayed trait (slower growth and insensitive to mating pheromone). Such a trait, exhibited by the ste4-∆/loc1-∆ strain, implied that the functions of STE4 and LOC1 are probably independent. One thing, which should be pointed out, is that the cause of slower growth of the ste4-∆/loc1-∆ strain is due to the absence of the gene LOC1, since the absence of gene STE4 prevents this strain from being affected by the α-factor. , which also proved that this strain is not affected by the α factor in terms of growth rate (Figure 7c,d). Although the ste4-Δ strain seems to grow slightly faster than the wild type at 0 μM, as reflected by the GRa of ste4-Δ (0.398) and the GRa of the wild type (0.342), this difference is not significant (p-value: 3.612 × 10 −1 ). Simi larly, the difference between the GRn of the ste4-Δ (0.427) and the GRn of the wild type (0.365) at 0 μM is not significant (Figure 7d, p-value: 2.986 × 10 −1 ). Immunity to the α-facto described above may facilitate the ste4-Δ strains' proliferation, rather than responding to a mating signal. Next, we studied the ste4-Δ/loc1-Δ strain. This strain also was not influ

Landscape View of Eight Gene Knock-out Strains under Four Levels of Mating Pheromone A-Factor
Besides the wild type strains and the four gene knock-out strains mentioned above, we measured the growth rate of the other four gene knock-out strains under different α-factor conditions. They were rvs161-∆, yku70-∆, the cin8-∆ strain (related to abnormal growth), and the ste5-∆ strain (defect on the mating pathway). In order to make a more detailed comparison in a manner of the landscape view, we calculated ∆GRa and ∆GRn, the difference between each strain and the wild type strain, as formula (3) below: In formula (3), the GR strain is the GRa or GRn of a gene knock-out strain, while the GRwt is the GRa or GRn of the wild type strain. Then, we obtained the ∆GR of area or number, respectively. Figure 7 is built by considering the ∆GR ratio of area (∆GRa ratio) as the x-axis and the ∆GRn ratio as the y-axis, so that the point of (∆GRa ratio, ∆GRn ratio) determined a specific location for a specific strain on a 2D plane, as exhibited by Figure 8. In this figure, strains located in the blue regions have a negative ∆GRa ratio or ∆GRn ratio, indicating slower growth than the wild type; while in the yellow regions, strains have a positive ∆GRa ratio or ∆GRn ratio, indicating faster growth than the wild type. growth), and the ste5-∆ strain (defect on the mating pathway). In order to make a mor detailed comparison in a manner of the landscape view, we calculated ∆GRa and ∆GRn the difference between each strain and the wild type strain, as formula (3) below: ∆GR ratio = GR strain − GR wt GR wt (3 In formula (3), the GR strain is the GRa or GRn of a gene knock-out strain, while th GRwt is the GRa or GRn of the wild type strain. Then, we obtained the ∆GR of area o number, respectively. Figure 7 is built by considering the ∆GR ratio of area (∆GRa ratio as the x-axis and the ∆GRn ratio as the y-axis, so that the point of (∆GRa ratio, ∆GRn ratio determined a specific location for a specific strain on a 2D plane, as exhibited by Figure 8 In this figure, strains located in the blue regions have a negative ∆GRa ratio or ∆GRn ratio indicating slower growth than the wild type; while in the yellow regions, strains have positive ∆GRa ratio or ∆GRn ratio, indicating faster growth than the wild type. Figure 8. The growth difference ΔGRa ratio and ΔGRn ratio of all strains, as compared with the wil type strains under four α-factor concentrations. The x axis represents ΔGRa ratio and y axis is ΔGR ratio. Each point represented a specific strain and is obtained by calculating the average da Figure 8. The growth difference ∆GRa ratio and ∆GRn ratio of all strains, as compared with the wild type strains under four α-factor concentrations. The x axis represents ∆GRa ratio and y axis is ∆GRn ratio. Each point represented a specific strain and is obtained by calculating the average data calculated from at least three parallel observation points from the same strain. The color blocks in each figure represent the area where strains have a faster (yellow area) or smaller (blue area) growth rate than the wild type strains. The hollow circle inside each figure represents (0, 0), where the strain has the same growth rate as the wild type strains in terms of both area and number. The line of y = x represents the set of points where the ∆GRa ratio is equal to the ∆GRn ratio.
At 0 µM, we found that roughly all strains, except swe1-∆, are located near the line of y = x. Especially, ste4-∆ and ste5-∆ strains grow faster than the wild type strain, even when there is no α-factor. At 0.3 µM, we found that nearly all strains, except ste4-∆ or ste5-∆, stay in blue, which means these strains did not grow faster than the wild type strains, in terms of both area and number. The strain of ste4-∆/loc1-∆ stays in the blue regions, apart from the ste4-∆ and ste5-∆ strains. This is because this strain has a lower basal growth other than the cell cycle arrest. At 0.3 µM, the swe1-∆ and loc1-∆ strains had similar growth status (p-value: ∆GRa, 4.017 × 10 −1 ; ∆GRn, 6.610 × 10 −1 ), even they had a huge difference at 0 µM (p-value: ∆GRa, 5.294 × 10 −8 ; ∆GRn, 3.549 × 10 −4 ). At 1 µM, locations of all strains, except ste4-∆ and ste5-∆, are closer to each other in the blue region. For strains with a mating-related gene knock-out, only the ste4-∆/loc1-∆ strain is an exception, located outside the yellow region, although the other two strains stay in the yellow region. At 10 µM, all strains were almost the same as in 1 µM, except the ste4-∆/loc1-∆ strain that moved into the yellow region. Here, we were convinced that this dual-gene knock-out strain had the potential to serve as a footstone to build a synthetic eukaryotic cell competition system. Unlike the ste4-∆ or ste5-∆ strains though, the growth rate of the ste4-∆/loc1-∆ strain had a progressive pattern. In other words, at the very beginning it grew more weakly than many other strains, until it gradually catches up with the other strains at higher mating pheromone levels. The results above shall enlighten future studies concerning interactions between strains with different adaptivity to a mating pheromone environment.

Conclusions
In this work, we introduced a newly designed microfluidic device to realize a gradient distribution of cell density, with only one trial of loading operation. Users can easily choose from six observation chambers with gradient initial cell density, for an expected cell density based on a simple pre-experiment for their own experiment design. Based on this microfluidic platform, strains with even a six-fold variation in growth rate can provide available data for the analyzation of one chip. We explored the growth nature of eight gene knock-out yeast strains under three different mating pheromone levels, ranging from mild (0.3 µM) to extremely high (10 µM). Then, we analyzed yeast growth by a quantitative method, which allowed for the estimation of yeast population growth in terms of population area and number together, which may bring a deeper understanding of yeast growth. Besides, to explore the potential interaction between growth and the mating response pathway, we constructed and studied the double-gene knocked-out ste4-∆/loc1-∆ strain. This shall be helpful for biological research focusing on behaviors of cooperation or competition. Our method shall contribute to a more throughput and comprehensive way to analyze a yeast population and provide insights about the growth status of relating genes or pathways.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/chemosensors10040141/s1. Figure S1: Growth rate analyzation of the wild type strain, in terms of population number and area; Figure S2: Growth rate analyzation of the swe1-∆ and loc1-∆ strains in terms of population area (GRa) and number (GRn); Figure S3: Growth rate analyzation of the ste4-∆ and ste4-∆/loc1-∆ of strains in terms of population area (GRa) and number (GRn).