Cell/Colony Motion of Oral Keratinocytes Determined by Non-Invasive and Quantitative Measurement Using Optical Flow Predicts Epithelial Regenerative Capacity

doi:10.21203/rs.3.rs-236508/v1

Download PDF

Research Article

Cell/Colony Motion of Oral Keratinocytes Determined by Non-Invasive and Quantitative Measurement Using Optical Flow Predicts Epithelial Regenerative Capacity

https://doi.org/10.21203/rs.3.rs-236508/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Cell/colony motion determined by non-invasive, quantitative measurements using the optical flow (OF) algorithm can indicate the oral keratinocyte proliferative capacity in early-phase primary cultures. This study aimed to determine a threshold for the cell/colony motion index to detect substandard cell populations in a subsequent subculture before manufacturing a tissue-engineered oral mucosa graft and to investigate the correlation with the epithelial regenerative capacity. The distinctive proliferating pattern of first-passage (passage 1 (p1)) cells reveals the motion of p1 cells/colonies, which can be measured in a non-invasive, quantitative manner using OF with fewer full-screen imaging analyses and cell segmentations. Our results demonstrate that the motion index lower than 40 μm/hour reflects cellular damages by experimental metabolic challenges and can be used as the threshold to determine the quality of cultured cells. Because the p1 cell/colony motion index is correlated with epithelial regenerative capacity, it is a reliable index for quality control of oral keratinocytes.

General Cell Biology & Physiology

Stem Cell & Developmental Cell Biology

quality control

regenerative medicine

oral keratinocytes

image segmentation

proliferative capacity

experimental metabolic challenge

The ability to generate functional equivalents as replacements for damaged or diseased human tissues is rapidly advancing in the fields of tissue engineering and regenerative medicine, which has harnessed the unlimited potential of mono-potent and multi-potent stem cells. By using transplantable tissues/organs manufactured from either autologous or allogeneic stem cells, regenerative technologies have contributed to cure various illnesses and injuries.^1,2 When culturing cells for relatively long periods, both the measurement and the evaluation of cell growth and differentiation potential are crucial to maintaining stable and reliable quality-assured cells for cell-based therapies. To meet the increased demand, quality control (QC) of engineered tissues and organs for clinical applications has become critical to confirm safety and effectiveness for the outcomes of post-grafting.³ However, prior to grafting, conventional cellular evaluation techniques based on molecular biology remain invasive, costly and time consuming. Furthermore, these techniques are incompatible with clinical applications in regenerative medicine because of the inevitable damage or loss of observed cells.⁴ More importantly, intact cells must be used in clinical settings. Thus, a non-invasive, low-cost and rapid cell evaluation method is required to expand the possibilities of regenerative cell-based therapy.⁵ To address such challenges, several types of novel image-based analyses that integrate computational data processing can be applied to a variety of imaging systems and evaluated for efficacy.⁶ These procedures allow investigators to monitor and quantitate specific states of cells during culturing in a non-invasive and high-throughput manner.^7–9

Autologous oral mucosa keratinocytes have been employed in clinical arenas, not only in regenerative dentistry but also in extra-oral regenerative medicine.^10–15 Thus, manufacturing protocols for tissue-engineered constructs using oral keratinocytes are required to confront regulatory challenges pertaining to safety, efficacy, quality etc. The manufacturing protocol comprises two major steps. The initial step establishes a two-dimensional (2D) monolayer cell culture by isolating and culturing cells from living tissue. Next, using a three-dimensional (3D) framework, biomimicking multi-cellular engineered grafts are manufactured as final products for human use. Given that the criteria for the release of final products should be determined immediately before transplantation, real-time and non-invasive examinations of the 3D constructs are crucial. Recently, an advanced imaging system was introduced to evaluate complex 3D cell-based constructs.^16,17 Chen et al. previously reported a non-invasive optical assessment for oral mucosa tissue-engineered constructs using optical molecular microscopy by which metabolic activities were characterised.^18,19 However, because QC needs to be implemented throughout the entire period of cell culture, the state of oral keratinocytes, including proliferation and colony formation before seeding cells onto the scaffold, is critical for manufacturing the quality of cellular-based products. Hence, optimal non-invasive cell measurements are essential for oral keratinocyte regenerative medicine.

Optical flow (OF) is a popular image-based analysis that demonstrates the distribution of the apparent motion of brightness patterns between two consecutive frames, which is caused by the movement of objects.²⁰ Since the OF algorithm is associated with the motility of objects, previous studies have been conducted to assess the ‘motility’ related to specific cell functions, including research on cellular biology^21–24 and cardiology.^25,26 Regarding cardiologic regenerative medicine, OF-based analyses have served as a non-invasive evaluation tool for the QC of cardiomyocytes to monitor contractile function without compromising the condition of the cells.²⁷ OF was also introduced to measure the collective motion speed of keratinocytes to identify skin stem cell colonies that were grown in the culture system with feeder layer cells.²⁸ Similar to that study, we reported that both the OF and normalised cross-correlation (NCC) can be leveraged to evaluate the proliferative capacity of primary oral keratinocytes within a relatively small-sized colony base during an early phase of cell culture.²⁹

To obtain a sufficient number of healthy cells to manufacture 3D cell-based constructs, serial cell cultivation (subculture) from the primary culture, referred to as p0 cells, is indispensable because p0 cells subcultured to the first (p1), second (passage 2) and third (passage 3) passage cultures yield a large number of cells with exponential cell growth.^10,11 It is well-known that the p1 cell/colonies are likely to grow discretely and localise in small, loosely connected colonies, mainly because of a higher plating rate when seeded as p1 cells,^30,31 which results in the rapid proliferation of p1 cells, whereas p0 cells in the feeder layer and serum-free culture take a longer time to reach a higher confluency.^32–36 Although our previous study demonstrated the applicability and feasibility of the OF algorithm to non-invasively measure p0 oral keratinocyte cell/colony motion, the spatio-temporal growth kinetics of p1 oral keratinocyte cells/colonies in a 2D monolayer culture is much different from that of p0 cells. Such an expansive analysis requires an entire microscopic image to target and observe a larger number of cells/colonies, instead of a single-colony base. Since our previous image analyses were not efficient to measure the full extent of the cell area, an improved protocol was needed to implement the time-lapse phase-contrast microphotography to p1 cells.²⁹

This study aimed to improve and enhance our previous protocol for the OF algorithm and to add the capability of determining the threshold of the cell/colony motion speed to differentiate substandard p1 oral keratinocyte populations before manufacturing a tissue-engineered oral mucosa tissue construct. By determining the specific spatio-temporal proliferating pattern of p1 cells, such as rapid growth and narrow space between neighbouring colonies, we could set correlative estimates that reduced the time of microphotography from 24 to 4 h, thereby resulting in fewer image analyses (31 frames compared to 96). Next, the image segmentation under full-screen analysis was used to identify the areas of cells/colonies under the microscopic images before applying the OF algorithm. Hence, 50% confluency was set to start the time-lapse microphotography and assess the overall colony growth and motion speed (MS) of each frame by applying the OF algorithm. We then examined the applicability of the mean motion speed (MMS) to monitor the proliferative capacity of p1 oral keratinocytes by calculating the population doubling times (PDTs) and determining the threshold of MMS that is suitable for transfer to a 3D culture condition non-invasively and quantitatively as a QC tool. Finally, we tested whether the MMS of p1 cells can be used as an index to predict the epithelial regenerative capacity and evaluate histologic features of the tissue-engineered constructs after the p1 oral keratinocytes had received experimental metabolic challenges during a 2D culture condition, thereby demonstrating the threshold of the motion index of the p1 cells to determine whether their proliferation was adequate for manufacturing constructs.

Validation of the time-lapse microphotography protocol for p1 cell/colony using the OF algorithm with image segmentation

Since the spatio-temporal growth (proliferating) kinetics of p1 oral keratinocyte cells/colonies differs from that of p0 cells, mainly due to the plating efficiency, a different protocol of non-invasive measurements was applied in this study. The protocol included time-setting (4 h at an interval of 8 min [31 frames]) for time-lapse microphotography and full-screen analysis before applying the OF algorithm.

To determine what starting condition was appropriate to implement non-invasive measurement for the p1 cells in this protocol, the overall cell/colony growth was investigated. It increased very slowly during 4 h in a similar fashion when time-lapse microphotography started at 30%, 50%, 70% and 90% confluency under manual observation, which indicated that the time-setting of 4 h is appropriate (Fig. 1A). On the other hand, the MMS continued to decline over the period of 4 h when started at 70% and 90% confluency, while it remained unchanged at 50% confluency and showed a slight increase at 30% confluency during the time-lapse imaging (Fig. 1B). Since the MMS is a parameter of the mean of the total frames (31), the increase and decrease of the MMS during microphotography is not appropriate as the representative data of the sample. Thus, as a new protocol, the starting conditions for time-lapse microphotography when the cells reached approximately 50% confluency was set and applied to further measurements used in this study.

Using the same OF algorithm as our previous report and different image segmentation, the motion behaviour of the individual p1 cells/colonies was successfully quantified, which indicated that MS can be applicable and feasible for non-invasive and quantitative measurements of the cell/colony motion (SFig. 1). The sum of the MS of the total 31 frames was divided by the 31 results in MMS (μm/hour), which produced a representative MS for the p1 cells. Validation of the OF algorithm used in this study was completed in the previous study by showing manual tracking of individual cells.²⁹

Correlation of the proliferative capacity of p1 oral keratinocytes with MMS using a standard protocol

When the MMSs of all 32 samples were analysed and their PDTs were plotted on a scatterplot, the diagram showed a negative linear correlation with PDT, which was statistically significant (Fig. 2) (Spearman’s r = –0.6669, 95% confidence interval = −0.8274 to −0.4057, p < 0.0001). Thus, oral keratinocyte p1 cells/colonies possessing a higher proliferative capacity were likely to have higher MMS values, thus indicating a higher locomotive ability, which is more robust correlation than the p0 cell/colony shown in our previous study.

Distinctive distribution of the MMS of p1 oral keratinocytes with their proliferative capacity under metabolic challenge protocols compared with that under the standard protocol

All of the MMS of p1 cells that received two types of metabolic challenges were plotted below the lowest MMS cultured under the standard protocol, which was 40 μm/h, thus demonstrating a distinctive pattern of samples between cells cultured under standard and metabolic challenge protocols (Fig. 3). This differential difference was confirmed by the representative images and movies having vectors (Fig. 4A–D and Supplemental movie 1A–D, compared with Fig. 5D and Supplemental movie 2B) Because the p1 cells cultured under metabolic challenges were damaged during the 2D monolayer culture, the MS of 40 μm/hour is highly likely to be the specific level of MMS and to be used for differentiating the substandard p1 oral keratinocyte populations.

Under the low nutrition challenge, their MMS decreased in the order of 5xPBS and 20xPBS; however, the proliferative capacity of 5xPBS and 20xPBS, which were mostly lower than that of standard protocol, varied and overlapped (Fig. 3; Fig. 4A, B; Supplemental movie 1A, B). Under the no-feeding challenges, the majority of samples showed a similar MMS, which is much lower than the standard protocol (Fig. 3; Figs. 4C, D; Supplemental movie 1C, D). By contrast, although the proliferative capacity of 4dNoF was similar to that of the standard protocol, it was higher than that of 7dNoF.

Histologic examination of the EVPOME received metabolic challenges compared with the standard protocol

We evaluated the histologic findings of EVPOMEs (n = 8) to further examine whether the MMS of the p1 cell/colony has the potential to predict an epithelial regenerative capacity and consequently to confirm whether the MMS of 40 μm/h is a detection threshold for monitoring the quality of the p1 cell before manufacturing an oral mucosa tissue-engineered construct, EVPOME. A well-differentiated, continuous stratified epithelial layer was developed when p1 cells cultured under the standard protocol (Fig. 6A, F) and cells under the 5xPBS and 4dNoF challenge protocols were seeded on AlloDerm^® (Fig. 6B, D, G, I). However, the epithelial layer of cells that received a challenged protocol was thinner and showed less distinct layers than that of the standard protocol, which suggests that their epithelial regenerative capacity was dampened (Fig. 6B, D, G, I). Moreover, although there was a continuous epithelial layer present in some of the EVPOMEs composed of cells cultured under 20xPBS and 7dNoF challenge protocols, some of the EVPOMEs showed more severe metabolic challenges, such as a discontinuous epithelial layer and complete failure to form the epithelium. Such poorly formed epithelial layers were not qualified for clinical use despite a continuous layer (Fig. 6C, E, H, J). The scatterplot of eight samples to examine the epithelial regenerative capacity showed a similar distribution pattern between cells cultured under standard and metabolic challenge protocols (Supplemental Fig. 2).

Proliferative activity of challenged oral keratinocytes in EVPOMEs

Furthermore, to quantify the proliferative activity of EVPOMEs that had received metabolic challenges when cultured in a 2D monolayer, the PI was calculated and compared with that of the standard protocol with Ki-67 immune-reaction (Fig. 7A–E). The result demonstrated that the PIs of EVPOMEs developed under 5xPBS and 4dNoF challenge protocols were significantly lower than that of the standard protocol (Fig. 8 A, B). Additionally, the PIs of EVPOMEs developed under the 20xPBS and 7dNoF challenge protocols were much lower than the 5xPBS and 4dNoF of challenge protocols.

A tissue-engineered cell/product developed ex vivo for use in wound repair, such as EVPOME grafts and various medical therapies, require testing, monitoring and inspection to evaluate its quality from start to completion of the cell culture process, which finally goes through a release testing prior to transplantation into/onto the human body. Current popular technologies usually destroy part of the manufactured cells/tissues for monitoring and testing, although the cells/tissues assessed should be put into the body later on. Thus, non-invasive tools to assess the in vitro viability, metabolic activity and structures in real time are necessary to develop in regenerative medicine. Recent advances in optics and imaging analysis, such as time-lapse microscopic images combined with computational machine learning^4,37,38 and the use of a specific sensor,³⁹ have provided many opportunities to address these issues for the QC of cells to be engineered. Previously, we applied two image analysis algorithms of OF and NCC to primary oral keratinocytes during the early stage of 2D cell culture and confirmed that the two motion indices calculated by both algorithms were applicable and feasible for non-invasive monitoring of cells and can be used as an appropriate tool for QC.²⁹ This study aimed to improve our previous methodology and to expand its application for p1 cells and, hereafter, in a later stage culture as a non-invasive tool for QC. We found that our technique enabled us to evaluate the quality of p1 and subsequent passaged cells non-invasively in a total 2D cell culture period prior to manufacturing the 3D tissue-engineered construct, which allows for thorough monitoring.

The present study applied only the OF algorithm with binary image segmentation to the cell area of p1 cells/colonies for the full-screen analysis as a non-invasive and quantitative measurement, instead of using two different algorithms (OF and NCC). Due to the higher density and narrower spaces between the neighbouring cells/colonies, which are characteristic of p1 cell spatio-temporal growth kinetics compared with p0 cells, we implemented different segmentation imaging techniques that do not require extractions of the cell areas in which targeted colonies are present (conventional approach by using the threshold) before applying the OF algorithm. This pre-processing results in an absolute and immense increase in the background region. In contrast to the decrease in data accuracy of NCC when the microscopic image contains a large background region, the accuracy of OF analysis was not affected in such cases, because the OF algorithm is based on image brightness intensity tracking. Thus, the NCC algorithm was not used in this study in terms of the accuracy of the image analysis. As shown in Fig. 5D and SFig. 1, because displacement vectors to determine MS were successfully visualised and MS was satisfactorily calculated, OF is tolerable and applicable for quantitative and non-invasive measurement of the p1 cell/colony motion in this full-screen analysis.

Although the OF algorithm used in this study was the same as that used in the previous study, the image segmentation as pre-processing prior to applying OF algorithm to calculate MS were different to determine ‘cell area’. In this study, MS recognises the mean value of all the lengths of the displacement vector within the cell area, which were binarised, while it measures only the part where the displacement vector length is 1 pixel or greater than the cell area in the previous study. Technically, in the previous study, in case of the length of the displacement vector smaller than 1 pixel, the targeted area was recognised as the background, not the cell area, although cells moved a little. However, in this study, since the cell area to measure the overall cell/colony growth is directly and entirely visualised and detected by segmentation, MS is accurately analysed despite the image containing a mixture of moving and non-moving cells. Thus, the different approach for image segmentation used in this study is likely to increase the data accuracy of the measurement of cell/colony motion, thereby contributing to a higher correlation coefficient.

Due to the difference in the spatio-temporal proliferating pattern between p0 and p1 cells/colonies, in the present study, time-lapse microscopic imaging started when the cell confluency reached 50%, based on the outcome of overall cell/colony growth and changes in the MS over the 4 h period. Since the reduction of MS was notable when microphotography started at ≥ 70% confluency despite continuous cell proliferation over 4 h, it is highly likely to increase cell–cell contact, which would result in contact inhibition during measurement. By contrast, the MS slightly increased when the imaging started at 30% confluency. This suggested that the cell culture conditions, such as cell confluency at the start of microphotography, affect the MMS as the cell/colony motion index in terms of the definition of MMS in this study.^40,41 Thus, this protocol is appropriate and validated for use in the quantitative measurement of p1 cell/colony motion.

This study revealed that the MMS of p1 cells/colonies had a negative correlation with their proliferative capacity, which indicated that the faster those p1 cells/colonies move, the higher the proliferative capacity they have, which is similar to the findings of our previous study. Notably, their correlation was almost at a ‘strong’ level, which suggests higher accuracy of our outcomes under the protocols used in this study for a non-invasive tool for QC in oral keratinocytes regenerative medicine. In fact, when applying the current protocol to p0 cells/colonies, their MMSs were absolutely lower (data not shown). In the previous study, even though the same objective lens was used, we only analysed zoomed subparts of the image by scaling the limited area where target colonies were present. Thus, the current protocol created for the characteristic growth kinetics of p1 cells/colonies could enhance data accuracy, because the full-screen analysis renders the travelling distance of cells shorter. It also contributed to a decrease in the variation of MSs, since the cell motion ceases when cells divide.

This study was able to determine the threshold of MMS to efficiently differentiate substandard cell populations before manufacturing a tissue-engineered oral mucosa tissue construct by utilising two types of metabolic challenge protocols that are similar to protocols used in previous studies.^18,19 Obviously, as opposed to cells being cultured under the standard protocol, MMSs from cells that received metabolic challenges when grown in a 2D monolayer culture were plotted with an MMS lower than 40 µm/hour. Additionally, the MMSs of cells cultured under more severe metabolic challenge protocols, such as 20xPBS and 7dNoF, tended to be lower than those under milder metabolic challenges, although some overlapped. Accordingly, the distribution of PDT plots was mostly dispersed around 2 days of cell culturing under standard and metabolic challenge protocols, which indicates a correlation with their MMSs, although the PDTs overlapped among cells receiving two types of metabolic challenge protocols. Apart from the primary aim, the nutrition challenge affected MMS more distinctly than the no-feeding challenge did, whereas the no-feeding challenge affected PDT more remarkably. As a result, the MMS of p1 cells/colonies determined by the OF algorithm is reliable as a non-invasive and quantitative index for QC in regenerative medicine.

Most importantly, the histologic evaluation of this study was successful in demonstrating the significant potential of the MMS in predicting the epithelial regenerative capacity of p1 oral keratinocytes. As expected, in contrast to the development of a continuous and well-stratified epithelial layer consisting of cells cultured under the standard protocol, cells cultured under the metabolic challenge protocols generated a thinner epithelial layer in EVPOME or failed to form a continuous epithelial layer, which demonstrated a poor regenerative capacity. Moreover, the significantly lower PI index of EVPOME that consisted of cells receiving 5xPBS and 4dNoF in a 2D monolayer culture condition was in accordance with the histologic findings, which is similar to the correlation with MMS. Considered together, these assessments of the epithelial regenerative capacity will make the MMS of p1 cells/colonies more robust and predictable as a QC tool. These findings confirmed that MMS, which is a non-invasive tool for monitoring the status of the oral keratinocytes culture, can be an informative and quantitative index that ensures that p1 cells are suitable for EVPOME fabrication. Nonetheless, the histologic findings indicated that cells with an MMS lower than 40 µm/hour under the metabolic challenges did not completely lose epithelial regenerative capacity. Thus, an MMS threshold of 40 µm/hour can be used to screen the quality of cultured cells.

However, a molecular mechanism for identifying the correlation between the locomotive ability and proliferative capacity of oral keratinocytes has not been elucidated. Although many factors/signalling pathways should be associated with those cellular characteristics, one that is associated with the epidermal growth factor (EGF) receptor/ligand system can control skin homeostasis and regulation of keratinocyte stem cells.⁴² Since our culture medium contains EGF as a growth supplement, elucidation of the correlation between cell/colony motion and the proliferative capacity of the EGF receptor/ligand system may guarantee more advances in regenerative medicine of oral mucosa.

Procurement of oral mucosa samples and culturing primary oral keratinocytes

The procurement of oral mucosa samples and procedure for oral keratinocyte cultures were described previously.²⁹ All methods were performed in accordance with relevant　guidelines and regulations. The details were described in the Supplementary Methods.

Time-lapse microscopic imaging

Under appropriate conditions, randomly-chosen five locations within the dish of the p1 cell culture were subjected to time-lapse observation. Phase-contrast images were taken at 8-min intervals for 4 h, until a total of 31 images were produced using a ×4 PlanFluor NA0.13 PhL objective lens. The images were converted to video files using a BZ-X analyser (Keyence) (Fig. 5A; Supplemental movie 2A). More detailed information was provided in the Supplementary Methods.

Determination of the cell area using image segmentation

Previously, we utilised a conventional segmentation method of applying the threshold to the intensity value by extracting areas containing targeted cells/colonies within the microscopic image. Briefly, using the MS as an index of cell locomotive ability, we selected areas with an MS of p0 cells/colonies that was equal to or greater than 1.00 pixel/frame, which corresponded to a p0 oral keratinocyte cell/colony area.²⁹ However, it is time consuming and laborious to apply to p1 cell cultures instead of p0 cells/colonies in the current experimental design. Hence, we implemented a different approach to achieve high-precision segmentation.

For image segmentation before determining the MMS of the p1 oral keratinocyte cells/colonies by using the OF algorithm, we used the variational segmentation method of Chan and Vese⁴³ to identify the surface area covered by all cells/colonies. The algorithm developed by Chan and Vese for active contours is a flexible method that enables the segmentation of many types of images. It has been applied to segment biological images such as cells, and its details have been reported.⁴⁴ There are four parameters (λ1, λ2, ν and μ) in this algorithm, and the parameters should be determined by the user. Details of each parameter are explained in a previous study.⁴⁴ Tuning of parameters was conducted by comparing manual and automatic segmentations and was selected to minimise the deviation between manual and automatic segmentations as much as possible. The validation data used was the image data shown in Fig. 5A. In this study, the preferred settings are λ1 = 1.2, λ2 = 1.0, ν = 0.02 and μ = 0.8.

Cell segmentation is conducted using the following three steps:

Step 1. Smoothing: All frames are smoothed using a Gaussian kernel (size is 3 pixel × 3 pixel window).

Step 2. Detect edges: 3 × 3 convolution kernels (Sobel filter) are applied to the image and vertical and horizontal derivatives are generated. Edges are detected by combining the two derivatives using the square root of the sum of the squares (Fig. 5B).

Step 3. Segmentation of p1 cells/colonies: By applying the Chan–Vese algorithm to the obtained image, the image is binarised into the background and cell area (Fig. 5C).

In the Chan–Vese algorithm, the use of small-sized images is recommended, because computation is time consuming. However, since the image size obtained by current time-lapse microscopic imaging is high resolution (1920 pixel × 1440 pixel), the algorithm is programmed in Python language with CuPy⁴⁵ to calculate based on the GPU (graphics processing unit). The runtime of Step 3 is approximately 10 min on a 1920 × 1440 resolution for cell images (31 frames) using a single GPU (GeForce GTX 1080 Ti, 11 GB) core on a common desktop PC (memory, 16 GB).

Measurement of the overall cell/colony growth

Because of the implementation of a different image segmentation, it is necessary to determine what cell growth conditions of p1 oral keratinocytes are appropriate to start time-lapse imaging. The simplest way to evaluate cell growth kinetics is to measure the overall cell/colony growth,⁴⁶ which was obtained via micro-photographing cells grown in the same culture dish at four different cell confluencies of 30%, 50%, 70% and 90% using manual observation. From the binary image after segmentation, the overall cells/colonies growth (%) was calculated, which was obtained by dividing the surface area covered by all cells/colonies shown in the full-screen mode by the entire image size (in pixel) for each frame over 4 h.

Determination of MMS of the p1 cells

The MS was determined in this study by using the identical OF algorithm of our previsou study.²⁹ A total of 31 sequence frames were created, in which vectors are drawn for each full-screen image (Fig. 5D; Supplemental movie 2B). The following procedure was described in the Supplementary Methods.

Sample size estimation

The minimum required sample size was a priori determined by power analysis using G*power software.^47,48 As a result, a total sample size of 26 was required by setting two-tails with an effect size of 0.5, a significance level of 0.05 and a power of 0.8.⁴⁹

Standard protocol for culturing p1 oral keratinocytes and evaluating the correlation between the motion index and proliferative capacity

To evaluate the correlation of the proliferative capacity with the MMS under the standard protocol according to the current experimental design, p1 cells grown at approximately 50% confluence in the 35-mm dish were subject to time-lapse microphotography. The mean MMS of the five locations was represented as the MMS of the sample (n = 32). Regardless of the confluency, 24 h after the completion of time-lapse microphotography, p1 cells were collected, and the number of cells was counted to determine the proliferative capacity. As a parameter of proliferative capacity, population doublings (PDs) and PDT of p1 cells were calculated as described in Supplementary Methods. Subsequently, a total 32 of MMS and PDT were plotted on the scatter plot to evaluate the correlation between MMS and PDT.

Metabolic challenge protocols for oral keratinocyte culture

We hypothesised that there should be a threshold for the MMS of p1 cells that can be differentiated as substandard cell populations. We compared the proliferative capacity of p1 cells under standard protocol and under two types of metabolic challenge protocols that include a lower nutrition challenge created by the dilution of completed EpiLife^® with D-PBS (Wako Chemical, Osaka, Japan) by 1:5 (5xPBS) and 1:20 (20xPBS) and a no-feeding challenge without a fresh culture medium change for four days (4dNoF) and seven days (7dNoF). Different from the standard protocol, which was used as a control for the metabolic challenges, a density of 1.25 × 10⁵ cells was plated into a 35-mm dish for those challenges except for 20x PBS, for which the cells were plated at a density of 1.5 × 10⁵ cells in a 60-mm dish with complete EpiLife^® medium. To secure random sampling, cells were alternately subject to culturing under metabolic challenge protocols using a total of 32 oral keratinocytes (n = 16).

To test a lower nutrition challenge protocol, the cells were fed every 2 days with complete EpiLife^®, and at approximately 50% confluency, they were fed with each medium diluted with D-PBS by the indicated ratio, respectively. The cells of 5xPBS in the 35-mm dish were fed every 2 days with the diluted medium and time-lapse micro-photographed 3 days later. By contrast, the cells of 20xPBS in the 60-mm dish were photographed under the time-lapse microscope 4 h after switching the diluted medium.

For the no-feeding challenge protocol, no medium change was made until they were subjected to time-lapse microphotography, which took three days for 4dNoF and six days for 7dNoF after the medium change, respectively.

Regardless of the confluency rate, 24 h after time-lapse microphotography, the p1 challenged cells were collected using the previously described method, and the number of cells was counted. PD and PDT were calculated using the formula provided in Supplementary Methods in which N = the number of challenged p1 cells collected, N₀ = the cell number inoculated, which is 1.25 × 10⁵ for 5xPBS, 4dNoF and 7dNoF or 1.5 × 10⁵ for 20xPBS) and I = days in the culture of challenged p1 cells. A step-by-step chart of the four metabolic challenge protocols is shown in Supplemental Figs. 3A and 3B.

Manufacturing of ex vivo produced oral mucosa equivalents (EVPOME)

Histologic and immunohistochemical examination of EVPOMEs

Evaluation of the proliferative activity of challenged oral keratinocytes in EVPOME

Manufacturing of EVPOMEs, their histologic and immunohistochemical examinations and evaluation of the proliferative activity of challenged oral keratinocytes in EVPOME were described in the Supplementary Methods.

Statistical analysis

To examine the strength of a linear association between MMS and PDT, Spearman’s correlation coefficient was calculated for non-normal distribution of the PDT, and the coefficient, r, and p values were determined using Prism 7.05 (GraphPad Software, San Diego, CA, USA). The results of PI are presented as the mean ± standard deviation (SD). The comparisons between the cells under the standard protocol and metabolic challenge protocols (5xPBS and 4dNoF) were examined using a paired t-test. A p-value <0.05 was considered statistically significant.

Data availability

The software and datasets generated and analysed during this study can be provided by the corresponding author upon reasonable request due to pending patent application.

In this study, we improved our imaging analysis method to monitor and evaluate the proliferative capacity of oral keratinocytes during the later phase of cell culture, specifically for p1 oral keratinocytes. The present study was successful in measuring p1 cell/colony motion by applying the OF algorithm non-invasively and quantitatively and demonstrated the correlation between the index of cell motion, MMS and the index of proliferative capacity, PDT. This was achieved by a different protocol to determine the cell area of p1 cells/colonies using image segmentation. This includes time-lapse microphotography starting at approximately 50% cell confluency and the fewer frames required for the full-screen analysis before applying the OF algorithm. Combined with the results of the MMS of cells receiving metabolic challenges during the 2D monolayer culture condition, we were able to determine an MMS threshold of approximately 40 µm/hour to differentiate the substandard p1 oral keratinocyte populations before manufacturing a tissue-engineered oral mucosa tissue construct. Additionally, our histologic findings of EVPOME confirmed that the MMS can be used to screen the quality of cells that have sufficient regenerative capacity. Considered together, even as a single parameter, MMS has the potential to predict future epithelial regeneration during a 3D culture and can be used as an informative tool and powerful supportive index for QC in oral keratinocyte regenerative medicine.

Acknowledgements

The authors are grateful to Dr. Hisanori Domon for his assistance of statistical analysis and would also like to thank Mr. Ibuki Takai for his assistance in performing image analysis.

Authorship

Conceived and designed the work: EH, KI.

Performed the experiments: EH, TS, AS, KH, RK, KI.

Analysed the data: EH, TS, KT, KI.

Contributed reagents/materials/analysis tools: EH, TS, AS, KH, RK, KT.

Visualised the data: EH, TS..

Wrote the original paper: EH, TS, KI.

Reviewed and edited the paper: EH, TS, KT, KI.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.

Ethical approval and patient consent

The use of human oral mucosa keratinocytes and the procurement procedure was approved by the Internal Review Board of the Niigata University Medical & Dental Sciences Hospital. Number: 2015-5018, titled ‘Translational research towards advanced regenerative medicine of oral mucosa: From bench to bed side’. Informed consent was obtained from all patients.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship and/or publication of this article: This research was supported by JSPS KAKENHI Grant Number 17H04398G and 20H03870A to K.I. This work was also supported by JSPS KAKENHI Grant Number JP16H06280, Grant-in-Aid for Scientific Research on Innovative Areas–Platforms for Advanced Technologies and Research Resources ‘Advanced Bioimaging Support’.

Behaegel, J., Ní Dhubhghaill, S., Koppen, C. & Zakaria, N. Safety of Cultivated Limbal Epithelial Stem Cell Transplantation for Human Corneal Regeneration. Stem Cells International 2017, (2017).
Mosquera-Perez, R. et al. Stem cells and oral surgery: A systematic review. Journal of Clinical and Experimental Dentistry. 11, 1181–1189 (2019).
Sasaki, T. et al. Evaluation of cell viability and metabolic activity of a 3D cultured human epidermal model using a dynamic autoradiographic technique with a PET radiopharmaceutical. Sci. Rep. 9, 1–10 https://doi.org/10.1038/s41598-019-47153-0 (2019).
Sasaki, K. et al. Non-invasive quality evaluation of confluent cells by image-based orientation heterogeneity analysis. Journal of Bioscience and Bioengineering. 121, 227–234 https://doi.org/10.1016/j.jbiosc.2015.06.012 (2016).
Nishimura, K. et al. Live-cell imaging of subcellular structures for quantitative evaluation of pluripotent stem cells. Sci. Rep. 9, 1–15 https://doi.org/10.1038/s41598-018-37779-x (2019).
Bajcsy, P. et al. Survey statistics of automated segmentations applied to optical imaging of mammalian cells. BMC Bioinformatics. 16, 1–28 https://doi.org/10.1186/s12859-015-0762-2 (2015).
Kato, R. et al. Parametric analysis of colony morphology of non-labelled live human pluripotent stem cells for cell quality control. Sci. Rep. 6, 1–12 https://doi.org/10.1038/srep34009 (2016).
Imai, Y. et al. In-process evaluation of culture errors using morphology-based image analysis. Regenerative Therapy. 9, 15–23 https://doi.org/10.1016/j.reth.2018.06.001 (2018).
Smith, D., Glen, K. & Thomas, R. Automated image analysis with the potential for process quality control applications in stem cell maintenance and differentiation. Biotechnol. Prog. 32, 215–223 (2016).
Izumi, K., Neiva, R. F. & Feinberg, S. E. Intraoral Grafting of Tissue-Engineered Human Oral Mucosa. The International Journal of Oral & Maxillofacial Implants. 28, e295–e303 (2013).
Izumi, K., Feinberg, S. E., Iida, A. & Yoshizawa, M. Intraoral grafting of an ex vivo produced oral mucosa equivalent: A preliminary report. International Journal of Oral and Maxillofacial Surgery. 32, 188–197 (2003).
Burillon, C. et al. Cultured autologous oral mucosal epithelial cell sheet (CAOMECS) transplantation for the treatment of corneal limbal epithelial stem cell deficiency. Investigative Ophthalmology and Visual Science. 53, 1325–1331 (2012).
Ohki, T. et al. Treatment of oesophageal ulcerations using endoscopic transplantation of tissue-engineered autologous oral mucosal epithelial cell sheets in a canine model. Gut. 55, 1704–1710 (2006).
Kim, Y. H. et al. Comparative analysis of substrate-free cultured oral mucosal epithelial cell sheets from cells of subjects with and without Stevens-Johnson syndrome for use in ocular surface reconstruction. PLoS ONE. 11, 1–15 (2016).
Barbagli, G., Heidenreich, A., Zugor, V., Karapanos, L. & Lazzeri, M. Urothelial or oral mucosa cells for tissue-engineered urethroplasty: A critical revision of the clinical outcome. Asian Journal of Urology. 7, 18–23 https://doi.org/10.1016/j.ajur.2018.12.009 (2020).
Teodori, L. et al. Three-dimensional imaging technologies: a priority for the advancement of tissue engineering and a challenge for the imaging community. Journal of Biophotonics. 10, 24–45 (2017).
Zhou, X. et al. Noninvasive real-time monitoring by alamarBlue® during in vitro culture of three-dimensional tissue-engineered bone constructs. Tissue Engineering - Part C: Methods. 19, 720–727 (2013).
Chen, L. C. et al. Optical metric assessed engineered tissues over a range of viability states. Tissue Engineering - Part C: Methods. 25, 305–313 (2019).
Chen, L. C. et al. The potential of label-free nonlinear optical molecular microscopy to non-invasively characterize the viability of engineered human tissue constructs. Biomaterials. 35, 6667–6676 https://doi.org/10.1016/j.biomaterials.2014.04.080 (2014).
Jan Erik Solem. Programming Computer Vision with Python.Program Comput Vis with Python.265–267. (2012)
Huang, Y., Hao, L., Li, H., Liu, Z. & Wang, P. Quantitative Analysis of Intracellular Motility Based on Optical Flow Model. Journal of Healthcare Engineering 2017, (2017).
Zahedi, A. et al. Deep Analysis of Mitochondria and Cell Health Using Machine Learning. Sci. Rep. 8, 1–15 https://doi.org/10.1038/s41598-018-34455-y (2018).
Parrillla, E. et al. Ciliary motility activity measurement using a dense optical flow algorithm. Conference proceedings - IEEE engineering in medicine and biology society 2013, 4446–9(2013).
Tamada, A. & Igarashi, M. Revealing chiral cell motility by 3D Riesz transform-differential interference contrast microscopy and computational kinematic analysis. Nature Communications. 8, https://doi.org/10.1038/s41467-017-02193-w (2017).
Wong, A. O. T. et al. Combinatorial Treatment of Human Cardiac Engineered Tissues With Biomimetic Cues Induces Functional Maturation as Revealed by Optical Mapping of Action Potentials and Calcium Transients. Frontiers in Physiology. 11, 1–11 (2020).
Lee, E. K., Kurokawa, Y. K., Tu, R., George, S. C. & Khine, M. Machine learning plus optical flow: A simple and sensitive method to detect cardioactive drugs. Sci. Rep. 5, 1–12 https://doi.org/10.1038/srep11817 (2015).
Czirok, A. et al. Optical-flow based non-invasive analysis of cardiomyocyte contractility. Sci. Rep. 7, 1–11 https://doi.org/10.1038/s41598-017-10094-7 (2017).
Kinoshita, K. et al. Automated collective motion analysis validates human keratinocyte stem cell cultures. Sci. Rep. 9, 1–12 https://doi.org/10.1038/s41598-019-55279-4 (2019).
Hoshikawa, E. et al. Noninvasive measurement of cell/colony motion using image analysis methods to evaluate the proliferative capacity of oral keratinocytes as a tool for quality control in regenerative medicine. Journal of Tissue Engineering. 10, https://doi.org/10.1177/2041731419881528 (2019).
Lamb, R. & Ambler, C. A. Keratinocytes Propagated in Serum-Free, Feeder-Free Culture Conditions Fail to Form Stratified Epidermis in a Reconstituted Skin Model. PLoS ONE. 8, 1–8 (2013).
Oda, D. & Watson, E. Human oral epithelial cell culture I. Improved conditions for reproducible culture in serum-free medium. In Vitro Cellular & Developmental Biology. 26, 589–595 (1990).
Xu, L., Schantz, S. P., Edelstein, D. & Sacks, P. G. A simplified method for the routine culture of normal oral epithelial (NOE) cells from upper aerodigestive tract mucosa. Methods in Cell Science. 18, 31–39 (1996).
Song, J., Izumi, K., Lanigan, T. & Feinberg, S. E. Development and characterization of a canine oral mucosa equivalent in a serum-free environment. Journal of Biomedical Materials Research - Part A. 71, 143–153 (2004).
Morino, T., Takagi, R., Yamamoto, K., Kojima, H. & Yamato, M. Explant culture of oral mucosal epithelial cells for fabricating transplantable epithelial cell sheet. Regenerative Therapy. 10, 36–45 https://doi.org/10.1016/j.reth.2018.10.006 (2019).
Izumi, K., Takacs, G., Terashi, H. & Feinberg, S. E. Ex vivo development of a composite human oral mucosal equivalent. Journal of Oral and Maxillofacial Surgery. 57, 571–577 (1999).
Izumi, K., Terashi, H., Marcelo, C. L. & Feinberg, S. E. Development and characterization of a tissue-engineered human oral mucosa equivalent produced in a serum-free culture system. Journal of Dental Research. 79, 798–805 (2000).
Matsuoka, F. et al. Morphology-Based Prediction of Osteogenic Differentiation Potential of Human Mesenchymal Stem Cells.PLoS ONE8, (2013).
Moen, E. et al. Deep learning for cellular image analysis. Nature Methods. 16, 1233–1246 https://doi.org/10.1038/s41592-019-0403-1 (2019).
Park, K. et al. Measurement of adherent cell mass and growth. Proceedings of the National Academy of Sciences of the United States of America. 107, 20691–20696 (2010).
Garcia, S. et al. Physics of active jamming during collective cellular motion in a monolayer. Proceedings of the National Academy of Sciences of the United States of America. 112, 15314–15319 (2015).
Tremel, A. et al. Cell migration and proliferation during monolayer formation and wound healing. Chem. Eng. Sci. 64, 247–253 (2009).
Nanba, D., Toki, F., Barrandon, Y. & Higashiyama, S. Recent advances in the epidermal growth factor receptor/ligand system biology on skin homeostasis and keratinocyte stem cell regulation. Journal of Dermatological Science. 72, 81–86 https://doi.org/10.1016/j.jdermsci.2013.05.009 (2013).
Chan, T. F. & Vese, L. A. Active contours without edges. IEEE Transactions on Image Processing. 10, 266–277 (2001).
Möller, M., Burger, M., Dieterich, P. & Schwab, A. A framework for automated cell tracking in phase contrast microscopic videos based on normal velocities. Journal of Visual Communication and Image Representation. 25, 396–409 (2014).
Okuta, R., Unno, Y., Nishino, D., Hido, S. & Loomis, C. CuPy: A NumPy-Compatible Library for NVIDIA GPU Calculations. Workshop on ML Systems at NIPS. 2017, 1–7 (2017).
Scherf, N. et al. Imaging, quantification and visualization of spatio-temporal patterning in mESC colonies under different culture conditions. Bioinformatics. 28, 556–561 (2012).
Faul, F., Erdfelder, E., Lang, A. G. & Buchner, A. G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods. 39, 175–191 (2007).
Beck, T. W. The importance of a priori sample size estimation in strength and conditioning research. The Journal of Strength and Conditioning Research. 27, 2323–2337 (2013).
Erdfelder, E., FAul, F., Buchner, A. & Lang, A. G. Statistical power analyses using G*Power 3.1: Tests for correlation and regression analyses. Behav. Res. Methods. 41, 1149–1160 (2009).

Supplementary Figure 1. MS determined by the current time-lapse microphotography protocol using OF algorithm

Representative change in MS of p1 oral keratinocytes to calculate MMS by 4-hour time-lapse microscopic observation (31 frames) when microphotography started at approximately 50% cell confluency.

Supplementary Figure 2. Distribution of the MMS of p1 oral keratinocytes with PDT under metabolic challenge protocols and the standard protocol used to develop EVPOMEs

Scatterplot showing the MMS and PDT of p1 oral keratinocytes cultured under two types of metabolic challenge protocols (low nutrition and no feeding) and the standard cell culture protocol for the development of EVPOMEs to evaluate the epithelial regenerative capacity after receiving metabolic challenges (n = 8 individuals).

Distinctive pattern of sample distribution is similar to that in Fig. 3 (n = 16 individuals).

Supplementary Figure 3. Step-by-step chart of metabolic challenge protocols for p1 oral keratinocytes

(A) Low nutrition challenge protocol. (B) No-feeding challenge protocol.

Supplementary Movies 1.

Representative video files with displacement vectors for MS calculation using OF algorithm. These movies are converted from phase-contrast images of cultured oral keratinocytes taken at 8 min intervals for 4 h (total 31 frames). Cells were cultured under the metabolic challenge protocols (A) 5xPBS, (B) 20xPBS, (C) 4dNoF and (D) 7dNoF.

Supplementary Movie 2.

These movies are converted from phase-contrast images of cultured oral keratinocytes taken at 8 min intervals for 4 h (total 31 frames).

(A) Representative video file before applying OF algorithm.

(B) Representative video file with displacement vectors for MS calculation using OF algorithm.

No competing interests reported.

Download PDF

Editorial decision: Major revision
30 Mar, 2021
Reviews received at journal
29 Mar, 2021
Reviewers agreed at journal
12 Mar, 2021
Reviews received at journal
09 Mar, 2021
Reviewers agreed at journal
27 Feb, 2021
Reviewers invited by journal
24 Feb, 2021
Editor assigned by journal
24 Feb, 2021
Editor invited by journal
17 Feb, 2021
Submission checks completed at journal
17 Feb, 2021
First submitted to journal
12 Feb, 2021

You are reading this latest preprint version

Cell/Colony Motion of Oral Keratinocytes Determined by Non-Invasive and Quantitative Measurement Using Optical Flow Predicts Epithelial Regenerative Capacity

Status:

Version 1

Abstract

Figures

Introduction

Results

Discussion

Materials And Methods

Conclusions

Declarations

References

Supplementary Materials

Additional Declarations

Supplementary Files

Status:

Version 1