Materials
To generate FMT and MI organoids, cells extracted from tissue samples were mixed with Matrigel (BD Bioscience, San Jose, CA, USA) and cultured in a stem cell-supporting medium as described previously 19,20,23. FMT samples were portions of mammary tumor tissues from a patient cat that were surgically removed at veterinary clinics (Minamigaoka animal hospital) in Japan, transported immediately to the laboratory in a shipping medium 21,22 and used to generate FMT organoids. The experiments were approved by the Institutional Animal Care and Use Committee of Tokyo University of Agriculture and Technology (approval number: 0020007). Mouse experiments were also performed in accordance with the Institutional Animal Care and Use Committee of the Tokyo University of Agriculture and Technology (approval number: R04-120). Female C57BL/6J mice (6-week-old) were obtained from Japan SLC Inc. (Shizuoka, Japan). The mice were anesthetized by isoflurane and the intestinal tissues were dissected and used for generating MI organoids. The studies were done in agreement with ARRIVE (Animal Research: Reporting of In Vivo Experiments) guidelines.
The anti-cancer drugs used included doxorubicin (Cayman, Ann Arbor, MI, USA) and toceranib sulfate (Sigma-Aldrich, Saint Louis, MO, USA). The antibody source used was HER2 (Novus Biologicals Inc., Centennial CO., USA), whereas the necrosis inhibitor was HS-1371 (MedChem Express Inc., Monmouth, NJ, USA). The fluorescent secondary antibodies used were as follows: Alexa Fluor™ 488 donkey anti-goat IgG; Alexa Fluor 488™ goat anti-rabbit IgG; Alexa Fluor 488™ goat anti-mouse IgG (all Thermo Fisher Scientific Inc., Waltham, MA, USA); Biotinylated goat anti-mouse IgG (Vector Laboratories Inc., Newark, CA, USA); Horseradish peroxidase (HRP)-conjugated anti-rabbit IgG (Cayman); and HRP-conjugated anti-mouse IgG (Millipore, Temecula, CA, USA).
Fabrication of the microfluidic device
The microfluidic device was fabricated using mechanical micromachining. First, a 3D device design was created using computer-aided design (CAD) software and subsequently converted into CNC milling paths using computer-aided manufacturing (CAM) software. The detailed dimensions of the device are provided in Figure S1 (unit: mm). A 2-mm thick acrylic sheet was machined to create microfluidic channels using a CNC milling machine (MiniMiller MM100, Modia Systems Co., Ltd., Japan). Two different drill bits, with diameters of 0.4 and 0.5 mm (MHR230, NS Tool Co., Ltd., Japan), were used to mill various parts of the microfluidic structure, including inlets, outlets, chambers, and channels, to achieve the required dimensions. After milling the microfluidic features, a 1-mm diameter drill bit (MHR230, NS Tool Co., Ltd., Japan) was used to segment the acrylic sheet into individual sections, each measuring 3 cm × 4 cm and containing a complete microfluidic channel network. The bases for the microfluidic device were prepared by machining another 2-mm-thick acrylic sheet into 3 cm × 4 cm blocks.
Thermal bonding was used to assemble the microfluidic device. A hot press (ROMANOFF, USA) was used to apply heat and pressure, which softened the acrylic at the interface and created robust bonds. Before bonding, residual debris within the microchannels was carefully removed using tweezers or adhesive tape to ensure channel clarity. The channel sections were then aligned with the base blocks and secured in position using high-temperature-resistant tape to prevent movement during the bonding process. Glass plates were placed on either side of the assembly to prevent the acrylic material from adhering to the metal plates of the hot press during bonding. The bonding was conducted at 215°C for 21 min.
The connectors for the inlet and outlet of the microfluidic device were fabricated by cutting the tips off the pipettes and inserting 0.75 mm diameter PTFE tubes (F-8007-001, Fron Industry, Japan) into the tips, which were secured with a strong adhesive (High Super 5, Cemedine, Japan). Connectors were then affixed to the inlet and outlet of the microfluidic device. After assembly, the device was left undisturbed overnight to allow the adhesive to fully cure.
Generation and culture of FMT and MI organoids
To prepare the FMT and MI organoids, mammary tumor and intestinal tissues collected from cats with mammary tumors and normal mice were cut and washed three times with PBS (Takara Bio Inc., Shiga, Japan). After washing, the tissues were minced finely, mixed in 0.1 mg/mL LiberaseTH (Sigma-Aldrich), and shaken in a 37 ℃ thermostatic bath for 30 min. They were then centrifuged at 600 × g for 3 min, removing the supernatant. After they were trypsinized in a thermostatic bath at 37 ℃ for 5 min, they were passed through a 70 µm pore size nylon cell strainer, collected in a 15-mL tube with fetal bovine serum (FBS, Sigma-Aldrich), and centrifuged for 3 min. After centrifugation, the supernatant was removed and the cells were washed three times with PBS. The supernatant was then aspirated and the cell pellet was suspended in Matrigel on ice and dropped into wells of a 24-well plate at 40 µL droplet size. After the gel was allowed to solidify in a CO2 incubator at 37°C for 30 min, the organoid culture medium was added (500 µL per well) and the culture plates were maintained at 37 ℃ under 5% CO2. The medium was changed three times per week.
Organoid passage
After culturing the organoids for 7–14 days, they were passaged into new wells using 5 mM EDTA/PBS solution at a ratio of 1:3–4, depending on the size and density of the organoids in the wells. Briefly, to dissolve Matrigel, 500 µL of 5 mM EDTA/PBS was added per well and the culture plate was kept on ice for 90 min. The organoid suspension was collected in a 15 mL tube and centrifuged at 4 ℃ and 600 × g for 3 min. After washing the organoids with PBS, it was treated with 1 mL of pre-warmed (37 ℃) TrypLE Express solution (Life Technologies Co., Grand Island, NY, USA) for 5 min at 37 ℃ and vigorously pipetted, and 100 µL of FBS was added to the tube to neutralize the solution. Cell pellets were collected by centrifugation and mixed with fresh Matrigel on ice as described above. The Matrigel-containing cells were dropped on a 24-well plate at 40 µL/well, solidified in a CO2 incubator at 37 ℃ for 30 min, and then the culture medium was added to each well and replaced three times weekly.
Hematoxylin and eosin (H&E) staining of FMT organoids
H&E staining of organoids was carried out as described previously 19,20. After fixation overnight with 4% paraformaldehyde (4% PFA, FUJIFILM Wako Pure Chemical Corporation, Osaka, Japan), each organoid sample was dehydrated with ethanol and xylene and then embedded in paraffin for microtome slicing into 4-µm-thick sections and mounted onto MS-coated glass slides. Thereafter, the sections were deparaffinized and subjected to H&E staining, according to standard procedures. Images were captured under a light microscope (BX-52; Olympus, Tokyo, Japan).
Immunohistochemical (IHC) staining of FMT organoids
IHC staining of organoids was performed as described previously 19,20. After deparaffinization of organoid sections with xylene and ethanol, the antigen was retrieved in 10 mM citrate buffer with heating at 121 ℃ for 5 min, followed by inactivation of the endogenous peroxidase activity by treating the sections with 1% peroxidase for 30 min. Subsequently, after blocking with 10% normal goat serum (NGS) for 30 min at room temparature, the samples were incubated at 4 ℃ overnight with HER2 (1:200). The sections were then washed three times with PBS, incubated with secondary antibody (EnVision Dual Link System-HRP), and visualized using DAB solution (Nacalai Tesque, Tokyo, Japan). Nuclei were counterstained with Mayer’s hematoxylin. All images were captured using a light microscope (BX-43; Olympus).
Live/dead staining
After perfusion of organoids with anti-cancer drugs, they were transferred from the device to 48-well plates, and 100 µL of LIVE/DEAD™ Cell Imaging Kit (488/570) reagents (Thermo Fisher Scientific) were added per well and incubated for 15 min at room temperature. After incubation, cells were observed under a fluorescence microscope (BZ-9000; Keyence, Osaka, Japan). The cell death rate was calculated by comparing the number of dead cells to the total number of cells.
Quantitative RT-PCR
Total RNA was extracted from the FMT and MI organoids using a NucleoSpin kit (Takara Bio Inc.) following the manufacturer's instructions. First-strand cDNA was synthesized using a QuantiTect Reverse Transcription Kit (QIAGEN). Quantitative real-time PCR was performed using the QuantiTect SYBR I Kit (QIAGEN) and StepOnePlus Real-Time PCR System (Applied Biosystems). The specific primers used for feline Bax, p53, Caspase-9, and Caspase-8 as well as mouse Bax, p53, Caspase-9, Caspase-8, RIPK3, Mlkl, Ulk1, Atg12, and Pick3c3 are listed in Table 1.
Numerical analysis using COMSOL
A COMSOL simulation was conducted to present the transport and diffusion of anti-cancer drugs, providing insights into the time-dependent concentration gradients within gels and organoids. The geometry used in the simulation corresponds to the actual dimensions of the microfluidic device (Supplemental Fig. 1a). The simulation involved a liquid with the same density and viscosity as water, carrying the species entering through an inlet, flowing through four chambers (PAs and PB), and exiting through three outlets. Each chamber contains a hemispherical cap (1.5 mm radius, 1 mm height), embedded with 50 solid microspheres (40–45 µm radius), representing the structural components of gels and organoids, and simulating species diffusion within them (Supplemental Fig. 1b).
The model utilizes Laminar Flow and Transport of Diluted Species interfaces to fully capture the fluid flow and diffusion. The simulations involve solving the fluid flow within the microfluidic device, with an inlet flow rate of approximately 570 µL/min. This results in a low Reynolds number, indicating a laminar flow. The fluid motion is modeled by solving the incompressible Navier-Stokes equations, which govern the behavior of viscous flows under the assumptions of steady state and incompressibility:
$$\:\frac{\partial\:}{\partial\:t}\left(\rho\:\stackrel{⃗}{U}\right)+\nabla\:\cdot\:\left(\rho\:\stackrel{⃗}{U}\stackrel{⃗}{U}\right)=-\nabla\:P+\nabla\:\cdot\:\left[\mu\:\right(\nabla\:\stackrel{⃗}{U}+\nabla\:{\stackrel{⃗}{U}}^{T}\left)\right]$$
$$\:\nabla\:\cdot\:\stackrel{⃗}{U}=0$$
Here ρ denotes density, u is the velocity, µ denotes viscosity, and p equals pressure.
The species are present at relatively low concentrations (10 µM) compared to the solvent, and changes in concentration do not significantly affect the fluid’s density or viscosity. This makes it appropriate to apply Fick’s law to describe the diffusive transport as follows:
$$\:-\nabla\:\cdot\:(-D\nabla\:c+c\stackrel{⃗}{U})=0$$
where D denotes the diffusion coefficient and c represents the concentration. The diffusion coefficient of the species \(\:{D}_{c}\) can be estimated based on the Stokes-Einstein relationship 34:
$$\:{D}_{c}=\frac{kT}{6\pi\:\mu\:R}$$
where k is the Boltzmann constant, T is the Kelvin temperature, R is the Stokes radius, and \(\:\mu\:\) is the dynamic viscosity of media. The diffusion coefficient for the gel \(\:{D}_{g}\) was adjusted by multiplying the media coefficient by a diffusion hindrance factor derived from the area fraction based on the gel fiber diameter in isotropic networks 35. A significantly lower diffusion coefficient \(\:{D}_{o}\) was applied to the organoids to facilitate the establishment of a concentration threshold representing the rate of cell death.
Numerical analysis of liquid flow and drug diffusion
Because the rate of cell death by perfusion of anti-cancer drugs drastically differed before and after branching in the device, we hypothesized that differences in flow velocity and drug infiltration would affect the survival rates of FMT organoids. To explore the influence of flow velocity and drug diffusion on FMT-organoid survival, we conducted a numerical analysis using COMSOL Multiphysics.
Supplemental Fig. 2a shows the velocity distribution within the microfluidic device, revealing significant differences in the flow rate before (PB) and after branching (PA). The flow rate at cross-section B-B (entrance to PB) was 559.8 µL/h, while cross-section A1-A1 (entrance to PA1) and A2-A2 (entrance to PA2) had flow rates of 233.5 and 157.6 µL/h, respectively. These variations resulted in differential extents of drug diffusion across the chambers after 48 h (Supplemental Fig. 2b). PB exhibited the highest degree of species penetration, with minimal concentration differences between the periphery and the center of the gel (Supplemental Fig. 2c). Conversely, PA1 and PA2, which had lower flow rates, exhibited reduced drug diffusion into the gel center. The lowest drug penetration was observed in non-perfused chambers.
To further assess the correlation between drug concentration and organoid mortality, we set concentration thresholds that indicated cellular death when exceeded (Supplemental Fig. 2d). The simulation results were consistent with the experimental findings, where higher drug penetration in PB corresponded to increased cell death, whereas PA regions exhibited reduced mortality owing to lower flow rates.
Statistical analysis
Data are presented as mean ± SD. Statistical evaluation was conducted using SigmaPlot software with one-way analysis of variance (ANOVA) followed by a t-test or Mann-Whitney test. When P values were < 0.05, it was considered statistically significant.