Histological studies and magnetic resonance imaging were employed to analyze the kidney structure and function of the freshwater mussel, Nodularia douglasiae. The Nodularia kidney consists of proximal, intermediate and distal tubules. The epithelia of the renal tubules were composed of a single layer of cuboidal cells. The proximal and distal tubules run in opposite directions underneath the pericardial cavity. Molecular weight cut-off (MWCO) values for the kidney filtration were detected by MR tracer injections: gadolinium-diethylenetriaminepentaacetic acid (GdDTPA) at 0.55 kDa, an oligomer-based contrast agent (CH3-DTPA-Gd) at 2.2 kDa, as well as Gd-DTPA-polylysine at 10, 22, and 110 kDa. The T1w-MRI intensity and T1 relaxation rate (R1) of the pericardial cavity and renal tubules increased with tracers smaller than 10 kDa. The other tracers showed only minimal or no increase. Thus, we concluded that the MWCO of the kidney is 22 kDa, 50 times larger than that for the Mytilus living in seawater. Since the R1 values of the renal tubules were similar to those of the pericardial cavity, the kidney did not concentrate filtrated tracers. The slow decay of the magnetic resonance (MR) tracers from the renal tubules indicated a low filtration rate, suggesting that the counter-current system reabsorbs useful solutes without reabsorption of water. The higher MWCO may be beneficial to maintain the tubular oncotic pressure and allow excretion of excess water. In conclusion, a main renal function of the freshwater mussel is the excretion of water, opposite to that of the seawater mussel and vertebrates, which preserve water.
The characteristics of the excretion systems of bivalves have not been well analyzed to the present date (Martin and Harrison, 1966; Bayne, 1976). The main function of the kidney in a terrestrial habitat is conserving water, in contrast to the loss of water to the air or urine. Seawater mussels also face the risk of losing body water, as it can be extracted by the hypertonic seawater (Bayne, 1976). We found that a seawater mussel, Mytilus galloprovincialis, can uptake Mn2+ ions from seawater, and can concentrate Mn2+ in urine by more than ten times (Wakashin et al., 2018). Thus, the Mytilus kidney can concentrate waste substances, such as heavy metals, and may reabsorb water to maintain a certain level of body water (Wakashin, et al., 2018, 2019). We also concluded that the molecular weight cut-off (MWCO) for filtration by the auricular wall of the Mytilus is around 0.5 kDa, which is almost 1/100 (Guyton and Hall, 2006) that of the MWCO for filtration by the glomerulus of the kidney in vertebrates, such as humans (Wakashin et al., 2019). Thereafter, we received criticism regarding our data concerning the concentration of Mn2+ in the urine and the small MWCO, compared with results reported for freshwater mussels such as the Anodonta cygnea (Potts, 1954). In freshwater habitats, there is the risk of dilution of body water due to the inflow of water caused by higher osmotic pressure of the body fluid. Therefore, the kidney of the freshwater mollusks has to excrete any excess amount of water to prevent dilution of the body fluid. In addition, the anatomical structure of the exocrine system of the Anodonta varies from that found in the seawater mussel, Mytilus (Martin and Harrison, 1966; Martin, 1983; Andrews and Jennings, 1993). Therefore, this study was conducted in order to analyze the function of the exocrine system of a freshwater mussel, the Nodularia douglasiae, and to compare our results with those reported for the Anodonta and Mytilus.
In order to investigate this hypothesis; (1) we studied the anatomical structure of the exocrine system using histological examinations conducted using a light microscopic method and high-resolution MR imaging, (2) we measured T1w-MRI and R1 values in the urine in the exocrine system of mussels enhanced by five MR tracers with varying molecular weights from 2.2 to 110 kDa: gadolinium-diethylenetriaminepentaacetic acid (GdDTPA) at 0.55 kDa, an oligomer-based contrast agent (CH3-DTPA-Gd) at 2.2 kDa, as well as Gd-DTPA-polylysine (Gd-PL) at 10, 22, and 110 kDa, and (3) we evaluated the relationship between the molecular weight of the MR tracers and the R1 of the urine in the exocrine system. Based on the results of these studies, we were able to determine the MWCO for the filtration conducted by the exocrine system, and, furthermore, to estimate the function of the renal tubules in the freshwater mussel, N. douglasiae.
The structure of the exocrine system
The exocrine system of the N. douglasiae was analyzed using high-resolution three-dimensional (3D) T1w-MR images and light microscopic images of Hematoxylin and Eosin (H&E) staining, and the important anatomical elements of the exocrine system were summarized in a schema of the mid-sagittal section (Fig. 1A).
The excretion system consists of the pericardial gland, auricular wall, pericardial cavity, renopericardial canal, the kidney, and the renal duct. The kidney consists of renal tubules: the proximal renal tubule (PT), the intermediate renal tubule (IT), and the distal renal tubule (DT). A pair of short renopericardial canals connected the near anterior-end of the pericardial cavity and the anterior end of PT, near Keber's valve (Fig. 1B). The PT was convoluted and distributed around the visceral sinus, running to the posterior direction (Fig. 1C). The epithelial cells were cuboidal, and had a shape similar to the brush border cell of the PT seen in vertebrates (Fig. 1F). The luminal space of the PT expanded around the posterior end of the ventricle of the heart, and thereafter, the tubule transferred to the IT. The IT filled the space between the posterior of the pericardium and the anterior of the posterior adductor muscle (Fig. 1D). Then, the IT shrank again and connected to the DT near the end of the pericardium. The DT ran to the anterior direction between the pericardial cavity and the PT (Fig. 1C). At the anterior-end of the pericardium, the DT connected with a pair of renal ducts (Fig. 1B). The epithelial cells were cuboidal, but rather flat compared with cells in the PT (Fig. 1E). Also refer to the high resolution T1w-MRI in the Movie 1.
The pericardial glands and the auricular wall are considered to be filtration sites for the hemolymph (Martin and Harrison, 1966; Bayne, 1976). As shown in Fig. 1E, the auricular wall was thin, and consists of a single flat cell layer. The pericardial glands were found around the anterior end of the pericardium, with epithelial cells similar to those of the auricular wall. Thus, the structure of the exocrine system of the Nodularia freshwater mussel was different from that found in the seawater mussel Mytilus (Bayne, 1976; Wakashin et al., 2018).
Distribution and accumulation of the MR tracers in the kidney
Before injection study of MR tracers, we did some preliminary experiments. The N. douglasiae were incubated in 1 mM GdDTPA containing water for 24 h. The digestive tract was visualized in higher intensity by T1w-MRI, but the image intensity of the kidney remained the same level. Therefore, MR tracers with higher molecular weight were not uptake from the digestive tracts. We also could not detect accumulation of Mn2+ in the kidney of the N. douglasiae after 24 h incubation with 50 µM Mn2+. When the mussel was incubated in water containing 200 µM Mn2+, the T1w-MRI kidney intensities increased to 130% at 12 h, and 200% at 48 h, compared with the control, and were maintained at the same level thereafter for 10 days (Fig. 2). In the digestive tract visualized at the higher signal intensity, Mn2+ might be absorbed slowly from the digestive tract. Since the intensity of the 200 µM Mn2+ water is 180% compared with that of pure water, the Mn2+ concentration in the kidney could be estimated at around 200 µM. In addition, intensities in the foot also increased in a similar fashion, up to 190%. Therefore, the Mn2+ concentration in the kidney might be of the same order as that in the foot, and the kidney did not concentrate the Mn2+. Judging from these results, the function of the exocrine system in a freshwater mussel may be different from that in the seawater mussel, which can concentrate Mn2+ in the kidney (Wakashin et al., 2018).
When MR tracers were injected into the foot, the T1w-MR image intensity of the visceral increased due to the increase in R1 caused by the MR tracers. Therefore, the image intensity of the kidney was depicted as a higher signal intensity at 10 min after the injection of GdDTPA or CH3-DTPA-Gd (P<0.05) (Fig. 3). Then, the image intensity of the kidney remained almost constant for 1 h (P>0.05). Therefore, GdDTPA and CH3-DTPA-Gd were filtrated into the kidney from the hemolymph. Neither the concentration nor excretion of the MR tracers were significant in this time scale.
The T1w-MRI intensity of the pericardial cavity and renal tubules increased with injections of CH3-DTPA-Gd and Gd-PL 10 kDa, while injections of Gd-PL 22 kDa and 110 kDa showed only minimal or no increase. Typical transverse T1w-MR images obtained at 1 mm anterior from the auriculoventricular valve are shown in Fig. 4. As shown in Fig. 4A and B, the T1w-MR image intensities in the area of the hinge were almost the same degree compared with the PC, PT and DT, representing filtration of CH3-DTPA-Gd and PL-Gd 10 kDa in the exocrine system of the mussel. The T1w-MR image intensity of the PC, PT and DT were smaller than those of the hinge. These results indicate that the MR tracers larger than 22 kDa were not filtered in the exocrine system. Coronal images at the level along the PT to the IT are shown in Fig. 4D. These images showed almost the same T1w-MR image intensity from the initial part of the PT to the IT, indicating that the MR tracers were not concentrated along the PT.
R1 relaxation rate of the pericardial cavity and renal tubules
We examined the difference of R1 in the exocrine system after the injection of the MR tracers. These results were summarized in Table 1. Statistically significant differences were detected in the mussel injected with CH3-DTPA-Gd by one-way analysis of variance (ANOVA), and a statistical difference was detected between the PC and the PT by Bonferroni post hoc test (P<0.05). We did not detect any other significant differences for the rest of the MR tracers. A small difference was also detected in the control between the PC and the PT (P<0.05). However, these differences were much smaller than the difference with and without the MR tracers. Thus, thereafter we used average R1 values of the PC, PT and DT for analyzing the effects of the MR tracers.
The means and s.d. of the R1 values after the injection of the MR tracers are summarized in Fig. 5. The R1 of the foot and mantle increased significantly after injection of the MR tracers, compared with the control (P<0.05), but there was no statistical difference shown between the R1 values of the 4 MR tracers. The R1 of the kidney (means of the PC, PT and DT) increased significantly after injection of CH3-DTPA-Gd and 10 kDa Gd-PL, compared with the control (P<0.05), but there were no statistical differences shown after the injection of 22 kDa and 110 kDa Gd-PL (P>0.05). The R1 values of the kidney were significantly smaller than those of the foot and mantle after the injection of 22 kDa and 110 kDa Gd-PL (P<0.05).
The concentrations of the MR tracers were estimated by C=(R1–Rc)/K, where Rc and K are the R1 of the control and the relaxivity value of the MR tracers. Since the increase of R1 was around 0.8 s−1 and the K of the CH3-DTPA-Gd [4.8 (s mM)−1], the concentration of CH3-DTPA-Gd was about 170 µM. From a linear regression of the dose dependency of the CH3-DTPA-Gd injection, a slope of 0.490±0.026 (s mM)−1 (coefficient±s.e.m.) was obtained, which was 1/10 of the relaxivity of the CH3-DTPA-Gd [4.8 (s mM)−1]. Therefore, the injected tracers (0.2 ml) were diluted tenfold in the mussel. Since CH3-DTPA-Gd could not enter cells, the total volume of the hemolymph could be estimated as 2.0 (1.8–2.2) ml (95% confidence limits).
Molecular size dependency for filtration of exocrine system
The MR tracers GdDTPA (0.55 kDa), CH3-DTPA-Gd (2.2 kDa) and 10 kDa Gd-PL (10 kDa) appeared in the PC and renal tubules, but 22 kDa Gd-PL (22 kDa) and 110 kDa Gd-PL (110 kDa) did not appear. From these results we concluded that the MWCO for filtration is around 22 kDa, which means it is possible to filtrate inulin (5.5 kDa). Inulin is the common substance for measurement of the glomerular filtration rate (Gaudino and Levitt, 1949). Therefore, the exocrine system of the N. douglasiae might have a similar or a larger MWCO, compared with the Anodonta Cygnea, which lives in freshwater (Potts, 1954), and the MWCO is 50 times larger than that of the seawater mussel, M. galloprovincialis (Wakashin et al., 2019). Filtration of inulin was reported in the abalone and the octopus, which live in seawater, and the giant African snail, which lives on land (Martin and Harrison, 1966). Therefore, the large MWCO is not likely to be a specific function of the exocrine system of the Mollusca, which lives in freshwater.
Function of the counter-current system of renal tubules
As concerned water balance, the kidney in seawater mussels conserve water against the hyperosmotic seawater. Therefore, the kidney of the Mytilus has a smaller MWCO (0.5 kDa) and can concentrate urine so that reabsorb water to maintain a certain level of body water (Wakashin, et al., 2018, 2019). On the contrary, the freshwater mussel, such as N. douglasiae, has the risk of dilution of body water due to the hypoosmotic freshwater. If filtrated fluid is excreted as it is, the mussel will lose useful substances, such as electrolytes and glucose. Therefore, the exocrine system has to reabsorb useful substances without reabsorbing water, and has to excrete any excess amount of water to prevent dilution of the body fluid. Based on our histological investigation, the renal tubules of N. douglasiae contain a counter-current system. A similar counter-current system was found in the kidney of the Anodonta (Andrews and Jennings, 1993). Andrews commented that “the proximal limb of the kidneys are released for the resorption of ions, excretory products being added in the distal limb” (Andrews, 1988). Indeed, the epithelial cells of the PT had a shape similar to that of the brush border cell of the proximal tubule found in vertebrates, and the epithelial cells of the DT had shape similar to the cell of the collecting duct of vertebrates (Fig. 1). The slow change of the MR tracers in the exocrine system (Fig. 3) suggests a low filtration rate, and this is in good agreement with the low filtration rate of the Anodonta (0.01 ml/hr/g wet weight) estimated by inulin clearance (Potts, 1954). The low filtration rate also allows for the time required to uptake useful solutes along the renal tubules. We have postulated two systems, as follows (Fig. 6).
Model A: (1) the exocrine system filtered out molecules smaller than 20 kDa, (2), so the PT tubule can reabsorb small useful solutes, such as Na+ and amino acids with reabsorption of water via isotonic osmosis, and the concentration of the non-absorbed solutes increased in the tubule, (3) and the DT excretes water due to the oncotic pressure caused by the non-absorbed solutes, and the DT may actively secret waste molecules. This is similar to the reabsorption found in the proximal tubules of the vertebrate (Guyton and Hall, 2006). Since luminal fluid in the PT remains in the isosmotic, reabsorption of solutes works the most efficiently via osmosis. In this system, the apparent concentration of non-absorbed solutes, such as the MR tracers, should be increased in the PT, and should be decreased in the DT (Fig. 6A).
Model B: (1) the exocrine system filters out molecules smaller than 20 kDa, (2) the PT tubule reabsorbs small useful solutes without the reabsorption of water, (3) and the DT tubule reabsorbs electrolytes and may actively secret waste molecules without any transport of water. This is similar to the counter-current system seen in the vertebrate without vasopressin or producing diluted urine (Guyton and Hall, 2006). In this system, the osmotic pressure of luminal fluid in the PT should be decreased by the reabsorption of solutes, while the apparent concentration of the MR tracers should remain constant (Fig. 6B). The long pathway of the IT and PT might be beneficial for additional reabsorption and/or secretion (Andrews, 1988).
As shown in our results for MR tracers smaller than the MWCO in Table 1, the R1 values of the PT and DT were not higher than that of the PC. As explained in the Introduction, the R1 value depends concentration of MR tracers: R1=R0+K·C, where R0, K and C are the intrinsic R1 of the hemolymph, the relaxivity value of the MR tracer (mM−1 s−1), and the concentration of the MR tracer (mM), respectively. Therefore, no increase of R1 values along the renal tubules indicates no increase of concentration of MR tracers in the urine. Thus, there was no water reabsorption in the renal tubules. In the Anodonta kidney, urine was hypotonic by 40% of the salt that was reabsorbed by the kidney (Picken, 1937; Martin and Harrison, 1966). The luminal fluid in the PT might be hypoosmotic at around 20 mosmol kg−1 compared with the hemolymph in Unionoidae (40–50 mosmol kg−1) (Picken, 1937; Burton, 1983; Potts, 1954; Byrne and Dietz, 1997). In addition, even with a small difference, as in the case with CH3-DTPA-Gd (2.2 kDa), the R1 of the PT was smaller than that of the PC (P<0.05). If this is true, the decrease of R1 indicates secretion of water in the PT, because larger molecules such as polypeptides and immunoglobulins are reabsorbed by endocytosis (Boron and Boulpaep, 2004) and it is not likely that the mussel has a specific transporter for CH3-DTPA-Gd. We could not estimate a specific transport system of water. However, the active transport system using ATP can generate a gradient of 100 mosmol kg−1 (Guyton and Hall, 2006), and that could overcome the osmotic gradient caused by even the total osmotic concentration of hemolymph (40–50 mosmol kg−1). The oncotic pressure of the luminal fluid caused by large molecules may also prevent water reabsorption through the epithelia of the renal tubule. Judging by these results, we doubt that model A is the actual case, and feel inclined to believe that model B is actually the true case, or at least close, for transport in the counter-current system of the N. douglasiae. Roles of the IT remains an open question.
In conclusion, we (1) injected five MR tracers in N. douglasiae, and detected the accumulation of the MR tracers in the exocrine system by MRI, and (2) MR tracers smaller than 10 kDa appeared in the exocrine system, but MR tracers larger than 22 kDa did not. Thus, the MWCO is around 22 kDa. (3) MR tracers were not concentrated in the renal tubules, and (4) the renal tubules consist of a counter-current system, which suggested the reabsorption of useful molecules and the excretion of water to prevent any dilution of hemolymph. (5) A main renal function in the freshwater mussel is the excretion of water, which is opposite to that of the seawater mussel and the vertebrate, where one main function is to preserve water.
MATERIALS AND METHODS
Nodularia douglasiae (Gray, 1833) were supplied by Sato Craft (Okayama, Japan). These mussels were collected from a river in Seki, Okayama in June of 2016 and December of 2019. The species of the mussels was identified by sequencing of the mitochondrial cytochrome c oxidase subunit I gene (Folmer et al., 1994; Seo and Seo, 2019). After collection, the mussels were immersed in water, and then transported to the laboratory by a refrigerated transport service, and maintained at 10°C, within a time period of 16 h (Cool Ta-Q-BIN, Yamato Transport Co., Ltd, Tokyo, Japan). At the laboratory, the mussels were housed in aerated water (4 l) with glass beads (0.3–1.0 mm diameter; 5 cm depth) in a 6 l bath at room temperature (20–25°C). The mussels were fed with Chaetoceros calcitrans (WDB Environmental & Biological Research Institute, Tokushima, Japan), and 2.5×106 cells/l water in the bath were applied at intervals of 3–4 days. A total of 21 mussels were used in this study. The length, height and width of the mussels were 38.4±2.2 mm, 20.7±1.2 mm, 15.2±1.4 mm (means and s.d.), respectively. All of the animal experiments in this study were carried out under the rules and regulations of the ‘Guiding Principles for the Care and Use of Animals’ set by the Physiological Society of Japan, and approved by the Animal Research Councils at Dokkyo University School of Medicine.
MR tracers and anesthetic chemicals
Five kinds of MR tracers were employed in the tracer study as follows: GdDTPA (Magnevist, Schering, Berlin, Germany; 0.55 kDa), CH3-DTPA-Gd (NMS-60, Nihon Medi-Physics, Chiba, Japan; 2.2 kDa), and Gd-PL (Gd-DTPA-polylysine 10 kDa, 22 kDa and 110 kDa, BioPAL, Worcester, MA, USA). Two mM Gd-DTPA and CH3-DTPA-Gd dissolved in 20 mM NaCl solution (Burton, 1983), and Gd-PLs solution containing 10 mg polylysine/ml were used for the injections to the mussels.
Tricaine was used as an anesthetic (Lellis et al., 2000). Equivalent volumes of 0.1 g/ml tricaine solution (Tricaine methanesulfonate, FUJIFILM Wako Pure Chemical, Osaka, Japan) and 0.046 g/ml tris buffer solution (Trizma base, Sigma-Aldrich, St Louis, MO, USA) were mixed into water at the final concentration of 500 ppm or 750 ppm, and at a pH of around 7.
The mussels were sedated by 250 ppm and anesthetized by 500 or 750 ppm tricaine (Lellis et al., 2000), followed by injections of 0.2 or 0.1 ml of the MR tracers into the foot using a 30 G needle (Nipro, Osaka, Japan). Thereafter the mussels were incubated in pure water. The MR tracer injections were confirmed by noting a higher signal intensity of the visceral mass in scout T1w-MR images those are a set of images to identify the relative anatomical position of a collection of MR images.
Magnetic resonance imaging
The MRI examination of the N. douglasiae in this study used procedures noted in our previous report (Seo and Seo, 2019). In brief, the mussels were placed in a plastic tube (inner diameter of 22.5 mm), and were immersed in 12 ml of water without aeration, and the temperature was kept at 20°C by variable temperature control units (ECU-20 and BVT-2000, Bruker Biospin, Ettlingen, Germany). Water was exchangeable through another tube set in the bottom of the tube holding the mussel. The 1H MR images were obtained by a 7 T vertical MR microimaging system (AVANCE III, Bruker Biospin, Ettlingen, Germany) with ParaVision operating software (version 5.1) and equipped with an active shielded gradient (micro2.5) and a 25-mm 1H birdcage radiofrequency coil.
R1 relaxation rate was measured by a two-dimensional saturation-recovery imaging (2D SR T1-MRI) method with five relaxation delays from 0.1 s to 4 s. The pixel size was 180×180 µm and the thickness of slice was 1 mm. The total image acquisition time was 9 min 21 s. The R1 image was calculated by the Image Sequence Analysis Tool provided by ParaVision. Images with a large motion artifact were omitted from the following data analysis. In order to include the whole-body structure of the mussel, 3D T1-weighted gradient-echo imaging (T1w-MRI) was used. The MR signal in the kidney was analyzed by 3D T1w-MRI with a voxel size of 240×240×240 µm, and a combination of TR/TE/θ=50 ms/2.5 ms/45°, where TR, TE and θ indicate repetition time, echo-time and flip angle, respectively. The total image acquisition time was 7 min 40 s. Images were Fourier-transformed with a data matrix of 256×128×128 after zero-filling of data, and the final voxel size was 180×180×180 µm. For a higher image resolution, 3D slab T1-weighted gradient-echo imaging was employed with a voxel size 90×90×90 µm with a combination of TR/TE/θ=50 ms/3.2 ms/45°. Motion artifacts of the mussel were reduced by 500 ppm tricaine. The image intensity of the T1w-MR image was quantified by comparison with the image intensity of a reference capillary containing 0.5 mM MnCl2 solution. Statistical analysis of R1 was applied using the unequal variance t-test (Welch's test) or one-way ANOVA employing Excel 2016 (Microsoft, Redmond, WA, USA). P-values less than 0.05 were regarded as significant.
Anatomical information for each 4% paraformaldehyde (PFA) fixed mussel was obtained by 3D T1w-MRI, with a voxel size of 60×60×60 µm, and a combination of TR/TE/θ=50 ms/3.75 ms/45°. For the histological studies, the paraffin sections of the PFA fixed mussel were prepared using a slice thickness of 10 µm. The sections were stained by H&E staining. Images were detected by microscopes (BZ-X710, Keyence, Osaka, Japan; BX63, Olympus, Tokyo, Japan) with an image-stitching mode that is the process of combining multiple images with overlapping fields of view to produce a wide field of image.
Protocols employed in the MR experiments
Experiments on the accumulation of the MR tracers in the kidney started from the low-resolution 3D T1w-MRI (7 min 40 s) before injection of the MR tracers. Then, mussels were anesthetized by tricaine, and thereafter, we waited for 30–60 min until the foot extended. After injection of MR tracers into the foot, the mussel was returned to the pure water and the low-resolution 3D T1w-MRI were measured for 1 h.
Experiments on filtration of the MR tracers started from measurement of the low-resolution 3D T1w-MRI (7 min 40 s) and a pair of 2D SR T1-MRI (18 min 40 s) before the injection of the MR tracers. Then, the mussels were anesthetized by tricaine, and incubated for 30–60 min until the foot extended. After injection of the MR tracers into the foot, the mussel was returned to the pure water for 15 min. The low-resolution 3D T1w-MRI was measured to confirm uniform distribution of the MR tracers. Then, a pair of 2D SR T1-MRI were measured at around 45 min after the injection. When necessary, the high-resolution 3D T1w-MRI was measured thereafter under the anesthesia. The R1 value obtained in separate experiments were included as the R1 for the control condition.
The authors would like to express our thanks to Dr T. Yamaguchi, as well as Ms Y. Imaizumi-Ohashi, Ms M. Yokoi-Hayakawa and Ms Kazuyo Mashiyama (DMUSM) for their technical assistance. We also thank Prof. M. Furuse (National Institute for Physiological Sciences) and the Spectrography and Bioimaging Facility, National Institute for Basic Biology Core Research Facilities for technical support. We must also thank Prof. K. Inoue and Prof. S. Kojima (AORI, UT) for their helpful suggestions and encouragement to E.S.
Conceptualization: E.S., H.W., Y.S.; Methodology: E.S., H.W., Y.S.; Investigation: E.S., H.W., Y.S.; Resources: E.S., H.W., Y.S.; Writing - original draft: E.S., H.W., Y.S.; Writing - review & editing: E.S., Y.S.; Project administration: E.S., Y.S.; Funding acquisition: E.S., Y.S.
Parts of this study were supported by grants from the Japan Society for the Promotion of Science KAKENHI [JP18H02261 to E.S. (AORI, UT) and JP15K08185 to Y.S.].
The authors declare no competing or financial interests.