2024 Volume 21 Issue Supplemental Article ID: e211017
Singularity biology is a scientific field that targets drastic state changes in multicellular systems, aiming to discover the key cells that induce the state change and investigate the mechanisms behind them. To achieve this goal, we developed a trans-scale optical imaging system (trans-scale scope), that is capable of capturing both macroscale changes across the entire system and the micro-scale behavior of individual cells, surpassing the cell observation capabilities of traditional microscopes. We developed two units of the trans-scale scope, named AMATERAS-1 and -2, which demonstrated the ability to observe multicellular systems consisting of over one million cells in a single field of view with sub-cellular resolution. This flagship instrument has been used to observe the dynamics of various cell species, with the advantage of being able to observe a large number of cells, allowing the detection and analysis of rare events and cells such as leader cells in multicellular pattern formation and cells that spontaneously initiate calcium waves. In this paper, we present the design concept of AMATERAS, the optical configuration, and several examples of observations, and demonstrate how the strength-in-numbers works in life sciences.
This article reviews our newly developed trans-scale scope AMATERAS to promote the studies in Singularity Biology. We have realized optical imaging in a centimeter field of view with sub-cellular spatial resolution that allows dynamic observation of more than one million cells. By utilizing image analysis techniques, AMATERAS can accomplish rare cell detection in a massive cell population, their statistical analysis, and statistical analysis of all cells.
Singularity Biology is a novel approach to biological research that primarily focuses on drastic changes in the states of multicellular systems. These changes are triggered by singularities (specific points) involving rare cells or events. In this context, live imaging technology is crucial because it can allow for the direct capture of the moments when specific cells within a population initiate a state transition. Such transitions, which we defined as “singularity phenomena”, occur in various areas of biology, including cell differentiation and tissue development [1,2], macroscale pattern oscillations such as epilepsy and arrhythmia [3,4], disease initiation such as cancer transformation and viral infection [5,6], and brain pathologies such as Alzheimer’s disease [7]. In each case, direct observation of their dynamics by live imaging provides a clear and straightforward method to understand the singularity phenomena. An essential aspect of the observation system is the ability to simultaneously observe a large number of cells. The greater the number of cells observed simultaneously, the greater the likelihood of detecting a cell with a low probability of being present in the population. The performance of an imaging system is determined by its spatial resolution, which is capable of resolving a single cell while maintaining a large field of view (FOV) or observation volume.
Over the past decade, several research groups have reported the development of fluorescence live imaging devices with a large FOV and single cell resolution [8–12]. They are commonly referred to as “mesoscopes”. Most of them target neuronal activity imaging in neuroscience, by combining laser-scanning two-photon excitation fluorescence imaging with a large-aperture lens system to achieve neuronal activity imaging within a FOV of a few millimeters [9–11]. A commercially available two-photon Mesoscope system from Thorlab Inc. [13] makes this technology available to the general research community. Large-aperture lens systems for single-photon excitation fluorescence imaging in the visible wavelength range have also been realized by several groups: the Mesolens system is a pioneering development in this field, with a large FOV of 6 mm and a spatial resolution of sub-micrometers [8]. Now, our goal is to image as many cells as possible to perform singularity biology. To this end, two-photon imaging systems are limited in the number of cells they can image in a single FOV. This limitation stems from their reliance on a laser scanning system, whose pixel throughput cannot exceed the frequency of the excitation pulse train and, typically capped at a few tens of megapixels per second [10,11]. In contrast, CMOS-based wide-field imaging can exceed gigapixels per second thanks to the use of a high-speed interface. In terms of optical design, the Mesolens system achieves intracellular resolution within a 6 mm FOV by choosing a 4x magnification [8]. However, aiming to encompass a greater number of cells, we opted to prioritize a larger FOV, accepting a reduction in the spatial resolution as trade-off. Consequently, we decided to design and develop our own imaging system with a lower magnification.
In this study, we have developed a novel optical imaging technique for the simultaneous observation of massive cell populations [14,15]. We categorized this type of imaging system as “trans-scale scope” because it can provide imaging of spatial information from the micrometer to the centimeter scale, covering a wide range in the scale hierarchy (trans-scale imaging). We named our instrument AMATERAS (A Multiscale/Multimodal Analytical Tool for Every Rare Activity in Singularity) and successfully developed two different systems of AMATERAS with unique configurations. The first configuration is a prototype system, assembled from commercially available components, offering a straightforward approach. The second configuration, however, is a custom-designed unique system with significantly enhanced optical performance. This paper will describe the technological aspects of the imaging systems, explain our design considerations, and provide a comparative analysis with conventional microscopes. We will also discuss the effectiveness of the system and potential contributions to biology.
Design Concept of AMATERASIn order to realize the purpose of this study, i.e., to observe the behavior of individual cells at a single cell resolution during drastic state transitions in a multicellular system, we considered how the optical performance of the imaging system should be set. Obviously, this depends on the cell type to be observed, but for eukaryotic cells we considered a resolution of less than 3 μm as the resolution at which individual cells can be observed separately. This spatial resolution is inferior to the submicron resolution of the mesoscopes of the other groups mentioned above, but it allows for a larger FOV. At this spatial resolution, features of individual cells must be visualized in order to detect characteristic cells. While shape and dynamics obtained by bright-field imaging are important features, it is also essential to visualize the molecular features of the cell (type and amount of proteins, metabolites, messenger molecules, etc.) by using fluorescent probes. For this reason, the wavelength range of observation was fixed in the visible region, where a variety of fluorescent probes (fluorescent dyes and proteins) can be used. As for the configuration of the imaging system, a wide-field illumination type was chosen instead of a laser scanning type in order to achieve high temporal resolution imaging with a large number of pixels. As for the observation modes, we designed an apparatus capable of selecting between fluorescence imaging of specific features using fluorescent labeling and bright-field illumination imaging to observe morphology and dynamics without labeling.
The key components in wide-field visible light imaging are the lens system and the image sensor. The critical factor in increasing the number of cells that can be observed in imaging is the number of sampling points that make up an image, which is equal to the number of pixels in the image sensor used. In contrast to the scientific CMOS sensors with a few megapixels typically used in conventional cell biology, some CMOS sensors used in the machine vision have a significantly higher number of pixels, such as 120 or 250 megapixels. The use of such sensors allows the generation of images with a higher number of sampling points compared to conventional microscopes.
However, because an image sensor with a large number of pixels has a large chip size, it is useless if the lens system cannot cover this size. In addition, if the spread of the point image, which depends on the performance of the lens system, is too large compared to the pixel size, adjacent pixels cannot acquire independent information. That is, the advantage of the number of pixels is reduced. Therefore, it is essential to select an optimal combination of image sensor and lens system.
In conventional biological microscopes, the numerical aperture (NA) of the lens system at the sensor side (sensor-side NA, NAsensor) is at most 0.04, and typically around 0.02 (Fig. 1A). In the green wavelength region (λ~520 nm), the spread of the point image on the imaging plane (d=λ/2 NAsensor) is approximately 13 μm, which is in excellent agreement with the pixel size of commonly used CMOS cameras (a=6.5 μm), since the pixel size is generally considered to be optimal at half the optical resolution (a~d/2). On the other hand, it is larger than the pixel sizes (2.2 μm and 1.5 μm) of 120-megapixel and 250-megapixel CMOS cameras. This means that a lens system with a larger sensor-side NA is required (Fig. 1B), which also means that the object-side NA (NAobj=M×NAsensor; M is the magnification) must be large.
Schematic illustration comparing the scale dimension of the lens system and image sensor between the conventional optical microscope (A) and the two types of trans-scale scopes developed by us (B). SP: sample plane, OL: objective lens, TL: tube lens, CP: conjugate plane, FN: field number, NAsensor: sensor-side numerical aperture, PSF: point-spread function. See main text for the details.
In addition, since the chip size of a CMOS sensor with a large number of pixels is also large, the lens system must have a field number (FN) that can cover the entire chip area, where FN is defined as the diameter of the FOV measured at the image plane. In fact, the size of a 120-megapixel CMOS is 29×20 mm2 (35 mm diagonal), which is considerably larger than the FN of conventional microscope lens systems (FN=22–26 mm), making it impossible to use the periphery of the chip area. To maximize the number of cells, the FN must be larger than the sensor size, i.e., FN≥S (S: diagonal of the image sensor) (Fig. 1B).
Consequently, it is necessary to introduce or newly develop a lens with a large NA and FN. We examined the optical performance required for the lens system by considering the above mentioned two conditions, a~d/2 and FN≥S, as well as the actual pixel size (2.2 μm and 1.5 μm for 120 and 250 megapixel CMOS, respectively) and chip size (S=35 mm). As a result, the required specifications for the lens system are NAsensor>0.059 (for 120 megapixel CMOS) and NAsensor>0.087 (for 250 megapixel CMOS) on the sensor side, and FN>35. Since the required spatial resolution on the object side (d/M, M: magnification) is 3 μm or less, the magnification of the system is M>1.47 and M>1.0 for 120 megapixel and 250 megapixel CMOS sensors, respectively.
Optical Configuration of AMATERASAs a lens system to satisfy the above conditions, we used a telecentric macro lens (LSTL20H-F, Mutron, Japan) with a magnification of 2x and NAobj=0.12 (NAsensor=0.06) in combination with a 120-megapixel CMOS camera (VCC-120CXP1M, CIS, Tokyo, Japan) to develop the first version of AMATERAS (AMATERAS-1) [14]. In addition, to achieve a higher spatial resolution and a three-dimensional resolving capability, we developed an exceptionally giant lens system (Optosigma, Tokyo, Japan) with 2x magnification and NAobj=0.25 (NAsensor=0.125) and combined it with a 250-megapixel CMOS camera (VCC-250CXP1M, CIS, Tokyo, Japan) to build the second version of AMATERAS (AMATERAS-2) [15]. Figures 2A and 2B show the schematic diagrams of the two imaging systems. These configurations provide FOVs of approximately 1.5 cm×1 cm (14.6 mm×10.1 mm for AMATERAS-1 and 14.7 mm×9.4 mm for AMATERAS-2) with a spatial resolution of approximately 2.3 μm and 1.2 μm in the lateral direction, respectively. Both versions are equipped with fluorescence imaging and bright-field imaging capabilities and can be used in either or both modes as required. It is important to note here with regard to the AMATERAS-1 that fluorescence imaging is achieved by dark-field illumination. By illuminating the specimen with LED light from outside the aperture angle (7º) of the camera lens, high-contrast fluorescence imaging can be achieved with preventing excitation light from entering the camera lens. Multiple high-intensity LED light sources can be used for multi-channel fluorescence imaging without the need for a dichroic mirror [14]. On the other hand, it is difficult to achieve dark-field illumination with AMATERAS-2 because the working distance of the objective lens is not large. Therefore, two-color imaging was realized by redesigning a large-size dichroic mirror for epi-illumination. Both the AMATERAS units have motorized stages and mechanical shutters, and are controlled by a self-developed software with an autofocus function. This allows time-lapse observations of cellular dynamics under various measurement conditions (number of wavelength channels, acquisition time, imaging frequency, etc.). Further details can be found in the original papers [14,15].
Optical configuration for fluorescence imaging in AMATERAS. (A) AMATERAS-1. An excitation light consisting of up to five LEDs is directed onto the sample from below with an incident angle of approximately 35 degrees. A fluorescence filter is placed at the entrance of the camera lens. (B) AMATERAS-2w, wide-field imaging configuration of AMATERAS-2. An excitation light consisting of two LEDs or lasers is introduced into the pupil plane of the objective lens through a dual-band dichroic mirror. A fluorescence filter is placed in front of the CMOS camera. (C) AMATERAS-2c, multi-point confocal imaging configuration of AMATERAS-2. A pinhole array disk is placed at the imaging plane of the tube lens. Excitation light is passed through a short-pass dichroic mirror onto the disk. The fluorescence image projected on the disk is transferred to the CMOS camera through another camera lens (1x). (D) An example of fluorescence imaging of MDCK cells stained with NucleoSeeing by AMATERAS-1. The number of cells was analyzed using a machine learning software and found to be approximately 1.19 million cells in the FOV. (E) An example of dual-color fluorescence imaging by AMATERAS-2w. A cardiac cell sheet derived from human-iPS cells was observed with immunostaining for two types of myosin light chains.
In order to apply these imaging tools to a variety of specimens, it is necessary to expand not only the size of the two-dimensional FOV but also the three-dimensional (3D) volume. We have examined how to realize three-dimensional live imaging based on the optical performance of the AMATERAS-2 lens system and proposed two approaches.
The first approach is to image within a 3D volume with a wide-field imaging system and realize 3D imaging computationally; when the 3D volume is imaged, fluorescence from outside the focal plane is superimposed as background light on the fluorescence image at the focal plane. From this image, the out-of-focus fluorescence image component is estimated using a low-pass filter-based algorithm [15] and subtracted from the original image to selectively extract information at the focal plane. This method is effective when the shape of the observed object is simple and uniform, and the fluorescence image from the focal plane is not buried in the background light noise. In fact, it has been applied to the observation of quail embryogenesis to segment the ellipsoidal shape of the nuclei of vascular endothelial cells and track them in three dimensions over time [15].
On the other hand, only a limited number of cell types and organisms can satisfy the above conditions. For example, in thick samples, the fluorescence image of multiple layers overlaps in the depth direction, which can overwhelm the focal plane information. To observe such specimens, optical suppression of the background light is essential. The second approach developed for this purpose is multi-point confocal imaging (Fig. 2C), in which a disk with a pinhole array is rotated in the image plane and this plane is imaged onto the camera by a secondary relay lens system. The pinhole diameter was set to 6 μm, which is comparable to the diameter of an Airy disk at the wavelength of the green and red fluorescence in the disk plane. The performance of this method was verified by imaging chemically fixed and cleared brain slices [15]. and its application to live imaging is being considered for the future.
Number of Observable CellsThe above innovations have made it possible to resolve individual cells in a FOV that is significantly larger than that of conventional microscopes. In addition, the FOV has been extended to three dimensions. We can now discuss how many cells can be captured in a single FOV. In two-dimensional imaging, we have demonstrated the ability to image more than one million eukaryotic cells cultured on a dish substrate in a single image [14]. Figure 2D shows an observation of highly confluent MDCK cells by AMATERAS-1. In this case, the nuclei were stained with a fluorescent dye (NuceloSeeing) and the number of cells was counted by particle analysis to be approximately 1.19×106 in the FOV (14.6×10.1 mm2). The average area occupied by a single cell was estimated to be approximately 11 μm square. This sample had a high confluence, but, of course, the number of cells decreases as the size of individual cells increases. Figure 2E shows an example of the observation of iPS cell-derived cardiac cell sheets [16] by AMATERAS-2. The average cell size was approximately 61 μm square, and the cell number was 3.8×104, assuming that the entire FOV (14.7×9.4 mm2) was filled with the cells (however, the actual cell sheet existed only in an area of approximately 8.4 mm diameter with a cell number of 1.5×104). In this case, intracellular myofibrils can be observed. In contrast, if the target organism is a bacterium such as Escherichia coli (E. coli), the number of observable cells can be much larger because each cell is smaller. Assuming that a single E. coli occupies an area of 1×2 μm2, a simple estimate yields more than 70 million cells in the FOV. Note, however, that it is difficult to spatially resolve individual E. coli with the current AMATERAS system, so it may be necessary to spatially separate individual cells using a microwell array [17] or other means, but even so, it is still possible to observe a large number of cells simultaneously.
In three-dimensional imaging, the potential to observe more cells is even greater because information on multiple cell layers can be obtained. For example, in mouse brain imaging, the combination with tissue clearing techniques [18] would allow for observation at depths of up to several millimeters, making it possible to observe more than 10 million neurons. However, in the case of densely distributed cells, the poor depth resolution (~12 μm) leads to inseparable overlapping of multiple cells along the depth direction. To avoid this problem, the target cells should be sparsely distributed over a distance that allows for 3D separation. A successful example is the three-dimensional time-lapse observation of vascular network formation in quail embryos realized by applying the above-mentioned computational sectioning method to AMATERAS-2 [15]. As shown in Fig. 3A, an entire embryo was observed in a single FOV, allowing the detection of more than 10,000 cells in each time frame and the tracking of cells over 24 hours of time-lapse imaging. This is a much larger number of cells compared to conventional narrow-FOV imaging where approximately 1000 vascular endothelial cells were captured in each time frame in a limited area [19]. Thanks to the ability to view the entire embryo, we can now obtain information that conventional narrow-FOV observations cannot reveal, such as the migration pathways of all cells, correlations between distant sites, and lead to the discovery of specific cells that can trigger drastic morphological changes throughout the entire embryo, which is the target of Singularity Biology.
Applications of AMATERAS. (A) Time-lapse observation of the vascular endothelial cells in tie1:H2B-eYFP transgenic quail embryo [19]. Fluorescence images of the whole embryo and the ventral aortic region at two time points with 24 h difference are shown. Pseudo-colors in individual cell nuclei indicate their vertical positions (z-position) as shown in the inset (right-top). (B) Time-lapse observation of the process of nutrient-depleted Dictyosterium discoideum generating spiral waves of cAMP to eventually form multicellular bodies. Full FOV images of two time points, 5 h 42 m 30 s and 9 h after the nutrient removal, and a magnified image of a local region are shown. Individual particles correspond to individual cells, and cells with high cAMP concentrations are represented in red (indicated by light blue arrows). (Right-bottom) Cells with high cAMP concentrations were automatically detected by an image analysis and their transitions were plotted against time. Intracellular cAMP concentration was visualized by the fluorescent indicator protein, Flamindo2-RFP [21]. (C) Calcium ion dynamics imaging of MDCK cells. Calcium wave propagations were visualized by the calcium indicator O-GECO [23]. The calcium concentration shown in magenta is superimposed on the bright-field image. (Right-top) Stripchart of propagation area of the calcium waves. (Right-bottom) Log-log plot of the complementary cumulative distribution function (CCDF) of propagation area, fitted with a linear function. (D) (Left) Histogram of the fluorescence intensity of the nucleus in the image data of MDCK cells at two different cell densities, 1.19×106 cells/FOV (the image shown in Fig. 2(D)) and 0.23×106 cells/FOV. (Right-top) Fluorescence image of 1/20 FOV at the central region and (Right-bottom) the same image colored according to cell phase.
As described above, AMATERAS increases the possibility of detecting unique cells with a low probability of existence. In fact, in a previous study elucidating spiral pattern formation in the social amoeba Dictyostelium discoideum, AMATERAS was able to find “leader cells” with a less than 0.2% probability of existence [14,20]. By visualizing cAMP, which induces cell chemotaxis, with a fluorescent protein indicator [21], we observed the propagation of cAMP waves over time from the early stages of the spiral formation (Fig. 3B). By retrospectively playing the time-lapse movie and searching for the origin of the spiral pattern, we revealed the mechanism of pattern formation and the cells responsible for it. This result is one of the most symbolic and remarkable achievements in this field and is expected to make a great contribution to the field.
In addition to simply detecting rare cells, we were able to perform statistical analysis of the rare events. Simply detecting one or a few rare cells does not provide much insight into the mechanism of their occurrence. Detecting rare cells in sufficient numbers allows one to speculate about the stochastic process behind the rare cellular activity. Here we present an example of the observation of dynamically changing, spontaneous calcium transients and propagation in MDCK cells [22]. Figure 3C shows the frequency of calcium wave propagations by performing calcium imaging on MDCK cells stably expressing the fluorescent calcium indicator protein [23]. The sample had a high cell density, with approximately 1 million cells in the FOV, and only 2057 cells in 10 minutes were found to cause calcium propagation to surrounding cells. Several parameters related to the wave propagation can be quantified, including the distribution of the wave propagation area. A histogram of the propagation area shows that most calcium waves propagate only over a small area (<5000 μm2), while a very small number of propagation events have an exceptionally large propagation area (Fig. 3C, right-top). We found that it follows the power law distribution (Fig. 3C, right-bottom), which is a statistical distribution that represents self-organized critical phenomena in complex systems [24,25]. Although the detailed analysis is beyond the scope of this article, this result provides a new insight into the understanding of spontaneous calcium wave propagation. Such a result was first realized by the strength-in-number of AMATERAS in this study.
While the above examples introduced statistics for selectively detected minority cells, AMATERAS also has great potential for analyzing the statistics for all cells, including majority cells. Figure 3D shows the results of nuclear DNA content analysis after observing 1 million cultured cells (the image data shown in Fig. 2D). The fluorescence intensity of the nuclei of individual cells was measured. The histogram of fluorescence intensity, representing DNA abundance, was found to be bimodal, indicating the distribution of phases in the cell cycle. The two peaks correspond to the G1 and G2 phase, respectively. A similar measurement on another sample with lower cell density showed that the phase distribution was significantly different. This difference can be explained by the fact that contact inhibition occurs in high-density cells, resulting in a lower percentage of cells transitioning from G1 to G2 phase. Such data are usually obtained by flow cytometry [26]. The reliability of flow cytometry is guaranteed by the large number of cells observed [27], but our method can provide similar data because the number of cells is comparable to that of flow cytometry under imaging. Moreover, information on the morphology and spatial arrangement of individual cells and their dynamics can also be obtained. If a device for selective extraction of single cells [28] or aspiration of cell contents [29] is mounted on AMATERAS, the system would function as an image sorter that can selectively isolate only cells of interest with specific characteristics.
Future Perspectives for AMATERASThere are several perspectives for AMATERAS, both from a technological and application point of view. On the technological side, there is still much room for improvement in imaging performance (sensitivity, resolution, speed, etc.), which is the key to expanding future applications. A worthwhile development is to use this lens system for brain activity imaging by two-photon fluorescence. As mentioned in the introduction, it is difficult to achieve the full FOV of AMATERAS by single-point scanning, which is used by several mesoscopes [9–11]. However, recent innovations in optical technology, such as pulse train multiplexing technology [30] and scanless two-photon fluorescence imaging [31], keep us from losing the possibility of “two-photon AMATERAS”. AMATERAS also has potential as an effective observation tool for whole-brain imaging, which has been actively developed in recent years by combination of tissue clearing and light sheet illumination techniques [18,32,33]. With its capability to encompass an entire mouse brain across-section within its centimeter-scale FOV, AMATERAS achieves comprehensive 3D brain volume mapping through depth scanning alone, eliminating the need for tiling. This method is expected to dramatically enhance the throughput of whole brain imaging. Another possible development is the incorporation of various optical and spectroscopic techniques into the AMATERAS-2 system, such as molecular vibration imaging [34], quantitative phase imaging [35], and light field imaging [36]. In addition to these optical developments, AMATERAS can also be combined with single-cell omics techniques to gain a better understanding of the molecular context of cellular states within a critical transition of multicellular systems. One potential development is the integration with spatial omics technologies [37–39], which have been attracting attention in life sciences in recent years. A possible experimental approach is to observe the dynamics by live imaging in AMATERAS and then perform spatial omics on the same sample. This would allow quantification of the spatial distribution of the cellular state and thus investigate the correlation between the cellular state and the cellular dynamics.
In terms of application, although this paper has presented several examples of applications to systems biology and developmental biology, there are many other possible directions of application. For example, diseases based on genetic mutations, such as cancer, and viral diseases are phenomena that occur extremely rarely in a cell population at the onset of the pathological state and progress beyond the self-repair functions [5,40]. A fundamental understanding of such phenomena can benefit from AMATERAS’ ability to observe a large number of cells simultaneously. In addition, AMATERAS can be a very suitable observation tool for the study of artificial cell sheets in regenerative medicine, where the observation of an entire large area of the sheet with single cell resolution is desirable [41,42]. It also exhibits significant potential as a high-throughput image-based cell screening platform, uniquely suited for the comprehensive study for phenotypic responses to drugs or genetic modifications. These are just a few of the possible applications and our group has received many more research proposals. We expect that new potential applications will be discovered in the future as this technology becomes more widely known to the public.
This paper has reviewed the development and application of a trans scale scope AMATERAS, which is one of the key technologies for promoting Singularity Biology. The development strategy to realize imaging observation of a much larger number of cells than conventional microscopy techniques was described, and the effect of “strength in numbers” in practice was demonstrated with several examples. The development of AMATERAS was conceived driven by a clear need in the life sciences and was finally realized through the technological development described in this paper. We started this research with a specific purpose in Singularity Biology, but over time we found various directions of application, and finally it has become a multi-purpose tool. We believe that this technology will be applied to various research fields in the future and will become one of the standard tools in the life sciences.
T.I and T.N have patent applications related to AMATERAS. The authors declare no additional conflict of interests.
TI wrote the manuscript. TI, TK, and TN reviewed and edited it. TI, TK prepared figures. TK, YS, YF prepared samples and acquired data. YS, YF, YO and KH provided comments to improve the manuscript.
The evidence data generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
This work was supported by a Grant-in-Aid for Scientific Research on Innovative Areas “Singularity Biology (No. 8007)” Grant No. 21H00431 (YS), 21H00413 (YF), 19H05411 (YO), 18H05415 (KH), 18H05410 and 18H05408 (TN) from The Ministry of Education, Culture, Sports, Science, and Technology (MEXT), Grant-in-Aid for Transformative Research Areas (A) “Seeing through Scattering Media (No. 20A207)” Grant No. 21H05590 and 23H04135 (TI), the Research Program of “Five-star Alliance” in “NJRC Mater. & Dev.” (TN), The Uehara Memorial Foundation (TN), and Takeda Science Foundation (TN). We acknowledge the use of DeepL Write for enhancing the quality of the English language in this manuscript.