The use of sacrificial ink to print various tubular tissues has been well established by several groups, notably the Lewis laboratory, which has fabricated kidney tubules along with blood vessels67,76. Similar approaches have been used to print alveolar structures77 and vascularized tumors78. Here, these approaches were adapted to generate vascular geometries for investigating endothelial responses under physiological fluid flow (Figure 2). Vessel geometries were converted from a vascular image into a coordinate map and subsequently into G-code to define the printing path (Figure 2A–C). The resulting approach enabled the fabrication of vessel templates for studying how endothelial function and cytoskeletal organization are altered by tortuous vessels and aneurysms (Figure 2D–G). A 3D-printable Rose chamber design was also incorporated to eliminate the need for access to machinist services (Figure 2H). An example G-code file for vessel printing is provided in Supplementary File 1, and the corresponding Rose chamber and stage-mount CAD files are provided in Supplementary File 2 and Supplementary File 3, respectively. STL-based printing can require multiple print-head passes, resulting in braided vessel geometries with separated lumens (Supplemental Figure 1). The resulting platform was designed to be accessible to vascular cell biologists with minimal bioengineering expertise while supporting physiological pressure and flow rates and providing access to high-resolution microscopy (Figure 2I).
Following vessel fabrication and cellularization, the vessel device was incorporated into a closed-loop perfusion system to enable controlled media circulation (Figure 3). The system consisted of a peristaltic pump, a three-way stopcock connected to a priming syringe, a bubble trap, a pulse dampener, the vessel device, and a media reservoir arranged sequentially within the flow circuit (Figure 3A). This configuration allowed media to circulate continuously through the engineered vessel while incorporating components for system priming, bubble removal, and dampening of pulsations generated by the peristaltic pump. The assembled system could be operated with the vessel device positioned within the experimental setup, as shown in Figure 3B, providing a practical configuration for subsequent perfusion of cellularized vessels under controlled flow conditions.
One modification to the previously described sacrificial ink approaches was pasteurization of the ink before printing. This modification was necessary because applying published ink-preparation protocols resulted in mold-contaminated cultures. As with other Pluronic-based inks used for this application, the 25% Pluronic F127/1% PEO ink remains liquid at low temperatures and transitions to a packed micellar gel as the temperature increases. This formulation has a transition temperature of 15°C, which is convenient for templating structures at incubator temperatures and removing the ink at 4°C. Similar to pure Pluronic, this material also displays substantial shear thinning at the printing temperature (Figure 4A). However, even at low temperatures, the ink is too viscous for filter sterilization using standard tissue-culture filters and vacuum. Although the ink can be sterilized by pasteurization at 70°C or higher for 20s, these temperatures irreversibly alter the ink and reduce its viscosity at the printing temperature (data not shown). Lower-temperature, longer-duration pasteurization at 63°C for >1 h reduces contamination by fungi and other microorganisms without altering the ink's complex viscosity at printing-suitable temperatures. Rheological characterization using a rheometer equipped with a 25 mm cone-and-plate geometry showed that raw and pasteurized ink had similar viscosities across a range of deformation speeds at printing temperatures (Figure 4B), demonstrating that pasteurization does not compromise printability. Pasteurized ink also enables longer-term vessel culture without the need for antimycotic drugs, which can alter endothelial function.
Once the vessel is cellularized, confluent endothelial cells should coat the entire vessel, as shown in Figure 5, with the maximum-intensity projection in Figure 5A and the cross-section in Figure 5B. Media flow can then be increased to physiological shear rates, which can be estimated for approximately circular vessels using the Hagen-Poiseuille equation79. In practice, the flow rate required to achieve a desired shear stress can be estimated by accounting for medium viscosity and microfluidic resistance. In a typical 2D culture in vitro, media flow at 12 dyn/cm2induces most endothelial cell types to orient along the flow axis. This response can be visualized by fluorescent phalloidin staining of F-actin (Figure 5D,E) or by analysis of the Golgi-nuclear axis (Figure 5D,F), both of which show anisotropic orientation consistent with endothelial polarization, as previously described18. Although shear stresses exceeding 40 dyn/cm2 occur in vivo, such conditions were not modeled here.
A major consequence of altered vessel geometry is altered fluid flow, to which endothelial cells respond through several signaling pathways80,81. Several methods are available to monitor blood-flow velocity, including 4D Doppler ultrasound82, 4D flow magnetic resonance imaging83, and laser speckle imaging84. Although these techniques offer important in vivo advantages, resolving local flow differences across a single vessel can be challenging. The high numerical aperture imaging capability of this platform enables direct fluorescent particle tracking (Figure 6B) and high-spatial-resolution mapping of flow anomalies induced by vessel curvature using particle image velocimetry (Figure 6C). For example, slower flow is observed along the inner curvature of the vessel than along the outer curvature (Figure 6C), which may contribute to altered endothelial responses85.
Once a confluent vessel is established, patency and barrier function can be assessed by perfusing fluorescently labeled macromolecules, such as dextrans or proteins, through the lumen. Patent vessels lined with HUVECs resist the passage of fluorescently labeled dextrans (Figure 6D,F). In contrast, exposure to barrier-disrupting stimuli such as TNFα increases vessel leakage (Figure 6H, J). Local leakage rates can be visualized using kymograph analysis (Figure 6E,G,I,K). Because the kymograph represents distance on the X-axis and time on the Y-axis, relative leakage rates can be compared by quantifying the slopes of lines fitted to segmented fluorescence intensities (Figure 6E,G,I,K). These analyses enable evaluation of how local vessel geometry influences endothelial barrier function in response to defined stimuli.
A key advantage of this platform is improved microscopic access relative to what is possible in vivo. The system also facilitates immunohistochemical analysis of endothelial molecular organization, particularly the actin cytoskeleton. Endothelial cells are very thin and are surrounded by smooth muscle cells in vivo, which can complicate the localization of cytoskeletal structures specifically within the endothelium. This platform enables straightforward staining and imaging of adhesion molecules and filamentous actin in endothelial cells (Figure 5) and supports multicomponent analyses that are difficult in vivo, including assessment of endothelial polarization and division under physiological flow (Figure 5D–F). Combining these assays on a single platform enables direct examination of relationships among vessel geometry, altered flow, barrier function, endothelial molecular organization, and intercellular adhesion.

Figure 1: Workflow for printing and using 3D blood vessels on a chip. Flowchart of the protocols for generating and using 3D blood vessels on a chip. Black circles indicate the corresponding protocol numbers. Box colors indicate the temperature at which each protocol is performed or the sample is incubated. The workflow includes vessel-shape modeling, ink preparation, extracellular matrix (ECM) preparation, chamber preparation, vessel printing, vessel flushing and coating, cellularization, flow culture, functional assays, and histological analysis. Please click here to view a larger version of this figure.

Figure 2: Vessel printing and chamber design. Vessels can be designed from in vivo vascular images to reproduce patient-specific vascular geometries. (A) Small coronary vessel imaged by contrast-enhanced computed tomography (CT), reproduced with permission from Parekh et al.86 Scale bar = 2 mm. (B) Coordinate map of fiducial points from a portion of the vessel shown in panel A, generated by point recording in FIJI. (C) G-code output generated from the coordinate map using printer-control software. X-, Y-, and Z-coordinates are indicated in the lower right. (D) Schematic of a Rose chamber containing a bioprinted vessel. The design with top thumb nuts is compatible with the 3D-printed polylactic acid (PLA) chamber. (E) Rose chamber containing a vessel and reservoirs printed with sacrificial bioink on an ECM-coated coverslip. Scale bar = 9 mm. (F) Schematic of a Rose chamber configured for the aneurysm model. A wire threaded through the inlet and outlet needles serves as the vessel template, and sacrificial bioink deposited on the wire serves as the aneurysm template. (G) Rose chamber containing the wire and sacrificial bioink template on an ECM-coated coverslip. (H) 3D-printed PLA microscope stage adapter (left) and Rose chamber (right), providing alternatives to machined aluminum components. Metal clamps on the stage adapter are secured to the adapter plate. Scale bar = 30 mm. (I) Printed and flushed vessel in an aluminum Rose chamber mounted on a microscope stage with inlet and outlet needles and tubing attached. Scale bar = 20 mm. Abbreviations: CT = computed tomography; ECM = extracellular matrix; PLA = polylactic acid. Please click here to view a larger version of this figure.

Figure 3: Fluid-flow system for a cellularized vessel in a Rose chamber. (A) Schematic of the tubing configuration used to establish fluid flow through a Rose chamber containing a cellularized vessel. Arrows indicate the direction of flow. Numbered components are (1) peristaltic pump, (2) three-way stopcock with priming syringe, (3) bubble trap, (4) pulse dampener, (5) vessel device, and (6) medium reservoir. All components that contact the culture medium are sterilized before use. (B) Fluid-flow system containing a cellularized vessel in a Rose chamber inside a tissue culture incubator. Numbers correspond to the components shown in panel A. Inlet and outlet needles are secured with tape for stability. Scale bar = 40 mm. Please click here to view a larger version of this figure.

Figure 4: Rheological characterization of the sacrificial bioink. (A) Shear viscosity of unpasteurized 25% Pluronic F127/1% polyethylene oxide (PEO) ink as a function of shear rate at 5 °C, 10 °C, 14.4 °C, 15 °C, and 23 °C. (B) Complex viscosity as a function of angular frequency for unpasteurized and pasteurized Pluronic F127/PEO ink at 25 °C and 37 °C. Abbreviation: PEO = polyethylene oxide. Please click here to view a larger version of this figure.

Figure 5: Imaging and quantitative analysis of fixed and stained endothelialized vessels. (A) Maximum-intensity projection of spinning-disk confocal fluorescence images of a fully cellularized vessel after initial seeding and culture under low shear stress (1 dyn/cm2). The vessel was stained for VE-cadherin (green), F-actin (red), and DNA with DAPI (blue). Scale bar = 50 µm. (B) Reconstructed xz cross-section of a z-stack of spinning-disk confocal fluorescence images of a fully cellularized vessel stained with fluorescent phalloidin to label F-actin (red). Scale bar = 50 µm. (C) Maximum-intensity projection of a z-stack of spinning-disk confocal fluorescence images of a fully cellularized vessel after 24 h of flow at 12 dyn/cm2. F-actin was labeled with fluorescent phalloidin (red). Scale bar = 100 µm. (D) Spinning-disk confocal fluorescence image of a vessel region after 24 h of culture under a shear stress of 16 dyn/cm2. The vessel was stained for VE-cadherin (yellow), F-actin (magenta), the Golgi marker GM130 (green), and DNA with DAPI (blue). The arrow indicates the direction of flow. Scale bar = 30 µm. (E) Representative quantification of mean stress-fiber orientation from the F-actin channel in panel D using FibrilTool74. A grid of 45 regions, each with an area of 1,600 µm2, was analyzed, and mean fiber orientations are displayed as a radial histogram relative to the direction of flow. (F) Representative quantification of the nucleus-Golgi axis from panel D as a measure of endothelial cell polarization. Vectors from the center of each nucleus to the center of the corresponding Golgi apparatus are plotted as a radial histogram relative to the direction of flow; n = 75 vectors. Abbreviation: DAPI = 4′,6-diamidino-2-phenylindole. Please click here to view a larger version of this figure.

Figure 6: Analysis of fluid flow and vascular leakage. (A) Schematic of the vessel geometry used for the analyses in panels B–K. Letters indicate the approximate vessel regions corresponding to the respective panels. (B) Color-encoded temporal projection of a time-lapse spinning-disk confocal fluorescence image sequence showing fluorescent beads moving under flow through a curved vessel region. Color indicates frame number from 1 to 100; the outer and inner curvatures are indicated. (C) Particle image velocimetry (PIV) analysis showing relative flow velocity across the curved vessel region. Colors represent relative velocity, and arrows indicate flow direction and magnitude. Scale bar = 50 µm. (D–G) Representative vascular barrier assay of an untreated vessel perfused with Alexa 647-labeled 10 kDa dextran. (D,D′,D″) Images of a straight vessel region acquired at 0, 70, and 150 min, respectively. The white line indicates the region used for the kymograph in panel E. (E) Kymograph of the straight region in panel D showing a leakage rate of 0.39 µm/min, which is the slope of the dashed line. (F,F′,F″) Images of the curved vessel region acquired at 0, 70, and 150 min. The white line indicates the region used for the kymograph in panel G. Scale bar = 150 µm. (G) Kymograph of the curved region in panel F showing a leakage rate of 1.01 µm/min, which is the slope of the dashed line. (H–K) Representative vascular barrier assay after preincubation with 250 ng/mL of tumor necrosis factor alpha (TNFα). (H,H′,H″) Images of a straight vessel region acquired at 0, 200, and 400 min, respectively. The white line indicates the region used for the kymograph in panel I. Scale bar = 250 µm. (I) Kymograph analysis of the straight vessel region in panel H shows a leakage rate of 2.1 µm/min, which is the slope of the dashed line. (J,J′,J″) Images of the curved vessel region acquired at 0, 200, and 400 min. The white line indicates the region used for the kymograph in panel K. Scale bar = 250 µm. (K) Kymograph analysis of the curved vessel region in panel J shows a leakage rate of 6.7 µm/min, which is the slope of the dashed line. Asterisks in panels D–K identify the corresponding ends of the lines used to generate the kymographs. Abbreviations: PIV = particle image velocimetry; TNFα = tumor necrosis factor alpha. Please click here to view a larger version of this figure.
Supplemental Figure 1: Vessel printing using STL-based printer programming. (A) Differential interference contrast image of a printed vessel generated using STL-based printer programming, which directs the printer head to make two passes, producing a braided vessel structure. Scale bar = 250 µm. (B) 3D rendering of the flushed braided vessel filled with Alexa 647-labeled bovine serum albumin (BSA). A gap is visible between the lumens of the two braided strands. Scale bar = 100 µm. Abbreviations: BSA = bovine serum albumin; STL = standard tessellation language. Please click here to download this file.
Supplementary File 1: G-code for vessel printing. Please click here to download this file.
Supplementary File 2: CAD drawing of the Rose chamber. Please click here to download this file.
Supplementary File 3: CAD drawing of the stage-mount adapter. Please click here to download this file.