A hierarchically ordered compacted coil scaffold for tissue regeneration

Hierarchically ordered scaffold has a great impact on cell patterning and tissue engineering. The introduction of controllable coils into a scaffold offers an additional unique structural feature compared to conventional linear patterned scaffolds and can greatly increase interior complexity and versatility. In this work, 3D coil compacted scaffolds with hierarchically ordered patterns and tunable coil densities created using speed-programmed melt electrospinning writing (sMEW) successfully led to in vitro cell growth in patterns with tunable cell density. Subcutaneous implantation in mice showed great in vivo biocompatibility, as evidenced by no significant increase in tumor necrosis factor α (TNF-α) and interleukin 6 (IL-6) levels in mouse serum. In addition, a lumbar vertebra was successfully printed for mesenchymal stem cells to grow in the desired pattern. A long-range patterned matrix composed of programmable short-range compacted coils enabled the design of complex structures, e.g., for tailored implants, by readily depositing short-range coil-compacted secondary architectures along with customized primary design. Improvements to a recently developed technique for printing biomaterials could make it easier to fabricate complex scaffolds for tissue regrowth. By shaping biocompatible polymers into thin microfibers and then collecting the strands on a moving platform, researchers can build tissue scaffolds layer-by-layer into 3D structures. Yingchun Su from the Harbin Institute of Technology in Harbin, China, and Aarhus University, Denmark, and colleagues report that materials printed by this ‘speed-programmed melt electrospinning writing’ technique can be composed of coil-shaped microfibers. The team found that the coil-like strands, produced by manipulating platform collection speeds, were effective at introducing pores into 3D scaffolds to better mimic natural tissue. Experiments showed this approach enabled better control over the cell density and growth behavior of stem cells implanted in a scaffold shaped like an artificial lumbar vertebra. Speed-programmed melt electrospinning writing (sMEW) is used to create a hierarchically ordered biomimetic scaffold with long-range patterned and short-range porous architectures for cell growth in patterns with tunable cell density.


Introduction
Cell therapy, especially stem-cell transplantation, has attracted extensive attention as a promising treatment for many diseases since the day it was proposed [1][2][3][4][5][6] . In the typical process of cell transplantation, donor tissue is dissociated into individual cells, and then a proper scaffold is needed to support cell growth, differentiation or functionalization in vitro. Ultimately, the functional cells carried by the scaffold are implanted at the disease site to promote tissue repair or regeneration 7,8 . To date, many three-dimensional (3D) porous scaffolds have been reported as structural supports for cells to promote attachment, proliferation, and differentiation, to yield functionally viable tissues as the ultimate goal [9][10][11] . For example, incorporation of hydroxyapatite, a 3D porous composite scaffold, promoted osteogenic differentiation and provided a favorable microenvironment for good neovascularization and tissue integration 12 . Cell-laden peptide-tethered hydrogel scaffolds with a 3D porous architecture were applied for in vivo spinal cord recovery, which significantly facilitated nerve tissue integration and neuronal regeneration 13 . 3D porous scaffolds with interconnected pore networks provide efficient nutrient transport for cells with great potential in tissue engineering.
Electrospinning has recently become a widespread method to fabricate various fibrous scaffolds because of its easy operation and extracellular matrix mimicry 14,15 . However, traditional solution electrospinning using organic solvent solutions usually forms random, dense, 2D fibrous mats unless a high-speed rotating collector or assisted electrical/ magnetic field is used, and interspacing between adjacent fibers is uncontrollable [16][17][18][19] . With lower conductivity, melt electrospinning generates a more stable jet than solution electrospinning, where solvent evaporation further introduces unstable factors. Recently, inspired by the threedimensional (3D) printing technique, a melt electrospinning writing (MEW) process was proposed [20][21][22] , by which accurate positioning of submicron fibers that are~2 orders of magnitude smaller than 3D printing has been realized. Depending on the collecting speed, these submicron fibers can be either collected as coils at low speed or be stretched into straight lines at high speed 23 .
Poly-ɛ-caprolactone (PCL) has become the most popular scaffold material due to its nontoxicity, low melting point, good mechanical properties, and high biocompatibility 24,25 . Many PCL-based scaffolds have been developed via MEW for tissue engineering applications [26][27][28][29][30][31][32][33][34][35][36][37] : organized microfiber networks contribute to reinforcing hydrogels for chondrocyte differentiation 30 ; scaffolds with box-structured pores are adapted for the simultaneous seeding of multicellular spheroids 31,32 ; stretchable hexagonal microstructures can enhance the beating rate and cardiac maturation-related marker expression as a potential heart patch 33 ; and multiple layers of fibroblasts can be assembled with the aid of MEW scaffolds as a human prostate cancer microtissue model 34 . Many researchers in tissue engineering have paid continuous attention to PCL scaffolds made by MEW.
Here, we present a hierarchically ordered PCL microfibrous architecture via speed-programmed melt electrospinning writing (sMEW). The periodicity and density of the coiling microfibers were tailored by virtue of coiling instability from electrohydrodynamics to achieve a secondary tunable structure along with the primary pattern design. Among many geometrical designs, rhombic geometry has been reported with interesting cellular responses 38,39 . A variety of rhombic patterns were created after optimization of the stepped speed control over coil compactness during the printing process. The programmable speed control provides conventional MEW with an additional tailored secondary coiled structure, introducing tunable scaffold porosity and mechanical properties. Different from the conventional 3D porous scaffolds, which are produced by gas foaming, salt leaching or freezedrying where the interconnected pore networks are dependent on the structure of porogen or phase separation, the porosity of the coil compacted 3D hierarchical architectures can be directly tailored by the stage speed and the layout of the design. The robustness of this technique to form lumbar vertebrae cellular microarrays with tunable density was investigated in vitro. The in vivo biocompatibility was investigated by surgical implantation into female athymic BALB/c nude mice.

Materials and methods
Preparation of PCL scaffolds PCL pellets (Sigma-Aldrich, USA, Mn ∼45 kDa) were loaded into a Spraybase melt electrospinning instrument (CAT000111, Ireland) and heated to 90°C for 1 h to obtain a homogeneous polymer melt. The tip-collector distance, applied voltage, and pressure were 4 mm, 3.5 kV, and 0.25 bar, respectively. The stage was supplied with high voltage, and a 30 G stationary, flat needle was grounded. Aluminum foil was placed on a flat-plate collector mounted on a programmable stage whose motion was controlled by UCCNC software (using G-code). The morphology of the cell-free PCL fibers was characterized by scanning electron microscopy (SEM, HITACHI, Japan, tabletop microscope) at an accelerating voltage of 15 kV. The critical translation speed (CTS) is the transition point between coiled and straight fibers, and after this point, the fibers are stretched by the movement of the collector 40 . The design of a lumbar spine-like scaffold was fabricated by imitating the natural lumbar spine structure. The running G-code with running speeds of 200 mm/min for coils and 5000 mm/min for straight lines are shown in Supplementary note 6.

AFM colloid probe preparation
First, an OLYMPUS microcantilever (Olympus, Japan) was fixed in a micromanipulator. Through microcontrol by moving in the x, y, and z directions, the tip was slightly dipped in UV glue (Norland Optical Adhesives 65, USA) and immediately removed. Then, the dip-removal operation was repeated five times to obtain enough glue on the surface of the tip. After being covered with glue, the tip was moved to touch a single polystyrene particle (50 μm) until the polystyrene particle attached to the tip. Polystyrene particles were monodispersed on the cover glass with absolute ethanol (0.5 mL, VWR, USA) in advance. Finally, the cantilever was removed from the micromanipulator, and the cantilever was irradiated under a UV lamp (365 nm, VILBER-LOURMAT GmbH, Eberhardzell, Germany) to cure for 20 min. Each polystyrene particlemodified cantilever was examined by SEM before use.

Mechanical property measurements
The local mechanical properties of the specific rhombic-patterned scaffold were measured by AFMbased micro indentation under ambient conditions using a Nanoscope Multimode VIII system (Bruck, Santa Barbara, CA, USA). The cantilever sensitivity was calibrated by indenting the mica substrate. The force constant value k of the OLYMPUS microcantilever was experimentally determined before the mechanical experiments. The Derjaguin-Müller-Toporov (DMT) mechanical contact model was used for estimation of the local Young's modulus by Equation 1 from the withdrawal curves.
where F is the force acting on the cantilever, F adh is the adhesion force between the probe and the scaffolds, R is the tip radius, d is the deformation of the scaffold, and E is Young's modulus.

Cell seeding and scaffold construct culture
Before cell seeding, PCL scaffolds (0.6 cm × 0.6 cm) were placed on the bottom of 96-well plates and sterilized with a UV lamp (252 nm, VILBER-LOURMAT GmbH, Eberhardzell, Germany) for 30 min. EGFP-hMSC-TERT cells (an immortalized human telomerase reverse transcriptasetransduced MSC line with an enhanced green fluorescent protein gene) were seeded at a density of 2.5 × 10 5 cells/mL (20 μL). After culturing for 40 min under humid conditions at 37°C with 5% CO 2 , 100 μL of Dulbecco's modified Eagle's medium-high glucose (DMEM with 4.5 g/L D-glucose, L-glutamine and pyruvate) supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin was added, and the cells were cultured under humid conditions at 37°C with 5% CO 2 . After one day of adhesion, the scaffolds were transferred to new 96-well plates to measure the data absorption signal using a CCK-8 kit. New plates were used to ensure that the absorption signal was from the cells attached to the scaffolds rather than the cells on the bottom of plates. Before measuring the cell proliferation value on days 1, 3, 5, and 7, all the seeding scaffolds with cells were transferred to new 96-well plates. For the lumbar spine-like scaffold, 200 μL of cells at a density of 2.5 × 10 5 cells/mL were seeded, and after 40 min of culture, the scaffolds were transferred to a new 24-well plate. The cellladen scaffolds were fixed in 4% formaldehyde in PBS at room temperature for 20 min and treated with 0.2% Triton (Sigma-Aldrich, USA) in PBS for 15 min. Then, the cellladen scaffolds were stained with Hoechst 33258 (Life Technologies, USA) for nuclear visualization. Cells adhered to the fibers were visualized by a Zeiss LSM 700 laser confocal microscope (Carl Zeiss Micro-Imaging GmbH, Germany) and an EVOS FL auto cell imaging system (Life Technologies, Carlsbad, CA, USA). The cytotoxicity test was performed by a lactate dehydrogenase (LDH) assay according to the kit manufacturer's instructions (Roche Diagnostics, Mannheim, Germany), and cell viability was determined by a cell counting kit (CCK-8, Dojindo, Japan). Cytotoxic properties quantified by the LDH assay are a reliable colorimetric reaction measurement at 490 nm to determine the amount of LDH in the media released from damaged cells. LDH activity release analysis was performed from the collected cultured media after 24 h of cell culture. The LDH activity of the cells cultured on TCP without or with 1% Triton X-100 was used as a low toxicity control (0% cell death) and high toxicity control (100% cell death), respectively. The cytotoxicity of the scaffolds was calculated using Equation 2. cytotoxicity % ð Þ ¼ experimental value À low toxicity control high toxicity control À low toxicity control 100% ð2Þ

Surgical implantation
Female athymic BALB/c nude mice aged 4-6 weeks were purchased from Guangdong Medicinal Laboratory Animal Center (Guangzhou, China) and were used in accordance with the regulations of the Animal Ethical and Welfare Committee (AEWC) of Shenzhen University. Analysis was performed with three animals per group. A midline incision was made, and a subcutaneous pocket was created in the abdomen of each mouse in order to implant a 0.36 cm 2 (0.6 cm × 0.6 cm) piece of PCL scaffold. PCL sutures were used to close the incision, and then 0.5 mg/kg Baytril and 0.2 mg/kg Buprenex were administered for 3 days as an antibiotic and an analgesic, respectively. Serum samples were isolated from mouse blood before implantation and on days 3 and 7 after implantation and diluted for analysis. The concentrations of tumor necrosis factor α (TNF-α) and interleukin 6 (IL-6) in the supernatant were detected using the Mouse TNF-α and IL-6 ELISA kits (eBioscience, San Diego, USA) according to manufacturer-recommended procedures.

Results and discussion
Our MEW device, schematically represented in Fig.  1a, consisted of an air pressure pump, a heated extruder, a high voltage supply, and a moving stage collector. The morphologies of the PCL fibers extruded at 90°C, under different stage speeds varying from 200 to 4500 mm/min were investigated (Fig. 1b). Inconsistent with the previously reported coiling instability from electrohydrodynamics 32,35,41 , the fibers showed periodic coils at a low stage speed, lower coil density at an increased speed, and eventually turned into straight lines once the stage speed was increased above a threshold value, namely, the CTS 23,42,43 . Hrynevich et al. provided details about the determination of the CTS value 32 . The CTS of the PCL fiber was approximately 3500 mm/min at 90°C.
The periodic coil density was quantified and marked on the right side of the corresponding SEM images (Fig. 1b), assigned with a color bar. The coil density as a periodic property could be fitted as a function of the collector speed (Fig. 1c). For the coils printed at 90°C, the fitting equations were proposed as Equation 3 and Equation 4 (the fitting computation process is shown in Supplementary note 1 and Fig. S1). From the equations, the theoretical numbers of coils at speeds of 200 mm/min and 400 mm/min (which cannot be manually counted due to experimental conditions) were 66 coils/mm and 31 coils/ mm, respectively.
Coil density n s ¼ 2Q where s is the distance traveled by the collector; Q is the volume flow rate of the melted PCL; c is the course length of one coil; v is the collector speed; r 1 and r 2 are the minor and major axes of cross-section for fibers, respectively; n is the number of coils in distance "s"; and coil density (n/s) is the coil number per unit length. Furthermore, the effects of temperature, distance, and applied voltage on the morphologies, diameters, and coil densities of the PCL fibers were also studied (Supplementary notes 2-3 and Figs. S2-3). With increased temperature, the PCL melt possessed lower viscosity and accordingly showed a higher flow rate, which led to an increase in the diameter. The effect of temperature, and collector speed on the diameter of straight PCL fibers was found to be consistent with previous studies 44,45 . At the same stage speed (800, 1000 or 1500 mm/min), the PCL coil density increased linearly with temperatures from 75°C to 100°C (Fig. S2d). When the temperature was fixed at 90°C and the distance/applied voltage increased from 4 mm, 3.5 kV to 16 mm, 6.5 kV, the coil density at 1000 mm/min increased from 10 coils/mm (Fig. 1c) to 37 coils/mm (Fig. S3f). The coil density along the stage orientation could also be tailored by both temperature and distance/applied voltage. It is noteworthy that Hochleitner et al. used a constant speed below the CTS to print sinusoidal patterns, where the wavelength and peakto-peak distance of sinusoidal shapes were affected by the collector speed 23,46 . The tubular structures of "figure of eight" looping and a sinusoidal fiber pattern were also observed 47 . Various designed soft networks with waving fibers were demonstrated for engineering biomechanically and biologically functional soft tissues 48 . However, those studies applied constant speeds below the CTS, while various speeds (below the CTS) during printing have not been previously systematically explored.
Here, based on the established relationship between stage speed and coil density (Fig. 1b, c), a programmed stepped speed between 200 and 4000 mm/min was applied to print unidirectional (Fig. 1d, e) and bidirectional (Fig. 1f, g) parallel coil gradients. In Fig. 1d, the stage speed was increased from 200 to 4000 mm/min, and the distance between adjacent paths was 300 μm. The number of coils decreased with increasing stage speed and eventually disappeared. For the bidirectional path (Fig. 1f), the stage speed increased in one path from left to right and then inversely increased in the adjacent path. As shown in Fig. 1e, g, after repeating the unidirectional and bidirectional printing 10 times, macroscopic isosceles triangle-like patterns were created along each path. The enlarged images (Fig. 1h, i) showed the fiber morphologies at variable stepped speeds (200, 600, 1000, 2500, and 4000 mm/min) after printing 1 and 10 layers, respectively, where the single layer showed a similar result to that in Fig. 1b. Therefore, using programmed stage speed, triangle-shaped patterns with coil density gradients were successfully fabricated.
Next, by controlling the programmed stepped stage speed along the horizontal printing path, rhombic patterns with tunable geometry and coil density were successfully created. Different structures of the rhombic patterns were characterized by each block having an interstructural gap (g) and a long periodic diagonal length (d) (Fig. 2a). To print the pattern with g = 0.2 mm and d = 1.2 mm, as an example, the stepped stage speed was programmed as shown in Fig. 2b. For each step distance of 0.1 mm, the speed gradually decreased from 4000 to 200 mm/min and then increased symmetrically back to 4000 mm/min to complete one rhombic structure. Various rhombic patterns (p 1.1-1.3) were produced through these blocks, as shown in Fig. 2e-k, with the corresponding stage speed program listed in Fig. S4. The length of the short diagonal for every rhombic structure was 316 μm, which was determined by the intrinsic coiling instability, while g and d can be tailored by programming the stage speed. The faster the stage speed decreases/ increases, the shorter d becomes. A total of 10 layers were used for all patterns; accumulation of layers the corresponding to the pattern in Fig. 2e is shown from the SEM images in Fig. S5, which shows that not only the rhombic shape but also the 3D coil density can be controlled by the tunable layer number. Further combination of both horizontal and vertical printing with rhombus crossing at either their centers (p 2.1-2.3) or gaps (p 3.1-3.3) led to more sophisticated patterns with enhanced durability for practical handling. As shown in Fig. 2l-q, we created a structure with a constant periodic diagonal of 1.2 mm and a tunable inter structure gap from 0.2 to 1.8 mm. These scaffolds were complex systems composed of different patterns and coil compactness, which were arranged in the x and y directions with a highly ordered structure.
Furthermore, three rhombic structures (g = 0.2 mm, d = 1.2 mm) crossing with straight vertical lines, p 1.1, Fig. 2 Stepped speed control over rhombic patterns design. a Schematic of the periodic rhombic pattern design in terms of the parameters of interstructure gap (g) and long diagonal length (d). b A stepped speed curve was used to create the rhombic patterned PCL fibers in (c). c-q SEM images of PCL fibers with rhombic patterns in one direction (c-k, p 1.1-1.3) and in two directions (l-n, p 2.1-2.3, and o-q, p 3.1-3.3) formed by controlled stage speed showing various interstructure gaps and periodic diagonal lengths. The color bar marks the corresponding coil densities.
vertical stepped speed printing crossing at the center, p 2.1, or at the gap, p 3.1, were used to assess cell patterning. SEM images of the three cell-free samples, epifluorescence imaging, and confocal fluorescence z-stack imaging for cell-laden samples are shown in Fig. 3a-c. sMEW technology provides a programmable speedcontrol strategy for the fabrication of fibrous patterned substrates with a proper 3D culture environment for cellular microarrays. The cells grown along the rhombic patterned scaffolds formed rhombic patterned cell constructs with a depth of~50 μm.
Cytotoxicity tests (Fig. 3d) and cell proliferation experiments (Fig. 3e) were also conducted for these rhombic patterned scaffolds. All scaffolds showed toxicity lower than 6% with no significant differences compared to the tissue culture plastic (TCP) low control. It was observed that the initial cell proliferation from days 1 to 5 for cells cultured on the p 1.1, 2.1, and 3.1 samples and TCP were similar, with no significant differences between the groups. After 7 days of culture, cell proliferation in the p 1.1, 2.1, and 3.1 groups was significantly accelerated compared to that of the TCP group, and showed excellent 3D cell culture capacity. No significance was found among the p 1.1, 2.1, and 3.1 groups. The patterned scaffolds yielded high densities of viable cells, and the cell densities were sufficiently high to achieve multicellular architectures. The cell culturing process on grid scaffolds with speeds of 200 mm/min, 600 mm/min, 1000 mm/min, 2500 mm/min, and 4000 mm/min is also presented in Supplementary note 4 and Figs. S6-8. The day 1 adsorption values for all samples at 450 nm from the CCK-8 assay are given to show the initial cell attachment (Fig. S9). Last but not least, due to its flexibility, the scaffolds can also be used as injectable scaffolds (Video S1).
The p 1.1, 2.1, and 3.1 scaffolds were evaluated in vivo by testing the IL-6 response in BALB/c mice for their special patterns. IL-6 is a pro-inflammatory cytokine with multiple functions and is believed to be a pivotal mediator e Cell proliferation on the patterned PCL fabric (n = 6). f, g Relative IL-6 (f) and TNF-α (g) levels in serum from healthy mice before PCL implantation and on days 3 and 7 after implantation. n = 3; *p ≤ 0.05; **p ≤ 0.005; ***p ≤ 0.001. or host responder for injury and infection 49,50 . TNF-α is another cytokine known to play an important role in cellular and inflammatory responses 51 . Overproduction of IL-6 causes neurological disease in the central nervous system 52 . Additionally, expression of TNF-α is related to nuclear factor-κB that involved in tumour initiation and progression 53 . Therefore, the inflammatory response from implantation of the patterned scaffolds was further evaluated by analyzing the relative expression (compared with each mouse group before implantation, considering the individual differences) of TNF-α and IL-6 in mouse serum from blood before implantation and on day 7 after implantation ( Fig. 3f, g). For IL-6 expression, the p 2.1 sample at day 7 showed no significant difference compared with the control group. IL-6 expression was significantly downregulated in groups p 1.1 (from 17.03 to 8.88 pg/mL, a 0.52-fold change on day 7) and p 3.1 (from 32.22 to 14.76 pg/mL, a 0.46-fold change on day 7). The expression levels of TNF-α from mouse serum in groups p 1.1 and p 3.1 before implantation and on day 7 showed no significant difference. For the p 2.1 sample, the level of TNF-α secretion showed a significant reduction from 4.37 to 2.24 pg/mL (a 0.51-fold change). The expression levels, shown in Fig. S10, of TNF-α and IL-6 in mouse serum before implantation and on day 7 after implantation were all in the normal range, confirming the biocompatibility of the PCL scaffolds 54,55 .
The p 2.1 samples exhibited significantly higher IL-6 expression than the p 1.1 and p 3.1 samples, possibly because the p 2.1 scaffold is thicker (~120.8 μm at the thickest point) than the p 1.1 scaffold (~79.4 μm at the thickest point) and p 3.1 scaffold (~91.0 μm at the thickest point). Previous reports have proposed that thicker materials could produce a proportionally higher magnitude of foreign body responses 56,57 . The increasing expression of IL-6 from p 2.1 on day 3 seemed to inhibit TNF-α protein expression, resulting in the downregulation of TNF-α levels on days 3 and 7. This is in agreement with the hypothesis proposed by Starkie et al. 58 that IL-6 can inhibit TNF-α protein production. In turn, the decrease in the IL-6 levels from the p 1.1 and p 3.1 samples also correlated with the upregulation of TNF-α expression with a lower inhibitory effect on day 7, where the level of TNF-α after implantation on the 7th day exhibited a slight increase compared with the control. On day 3, the TNF-α level exhibited the lowest expression as a transitional period.
In order to investigate the microscale mechanical properties of the hierarchical coil compacted scaffolds, a local compression test using atomic force microscopy (AFM) with a modified probe (Fig. 4a) was conducted. The mechanical properties were determined by microindentation using AFM. Young's modulus based on the Derjaguin-Muller-Toporov (DMT) model was calculated from the withdrawal curves using the equation in Fig. 4b, inset. Fig. 4d shows the tip locations on one rhombic pattern, where 10 points were measured. Representative withdrawal curves measured at marked positions 6, 8, and 10 are shown in Fig. 4c. With increasing coil density (Fig.  4e) from tip position 1 to tip position 5, the modulus of the scaffold (Fig. 4f) increased from 0.22 ± 0.03 MPa to 1.10 ± 0.06 MPa and then decreased back to 0.24 ± 0.06 MPa at tip position 10 with decreasing coil density (Fig. 4e). The Young's modulus data at each point were therefore presented (Fig. 4f) in correlation with the theoretical fiber density (Fig. 4e). The mechanical properties of the scaffolds were tunable based on the coil density. The tensile stress-strain curve measurements for the grid scaffolds with different speeds are also shown in Supplementary note 5 and Figs. S11-12. Such coil compacted scaffolds can be designed on top of customized primary shapes with tunable secondary interior structures and mechanical properties 46,48 , which makes sMEW a great contribution to MEW technology.
To demonstrate the tailoring ability, a lumbar vertebra was printed using sMEW, according to the anatomical structure (Fig. 5a), with the design shown in Fig. 5b. The G code is shown in Supplementary note 6. The printed structure with a stage speed of 200 mm/min can be observed by optical microscopy (Fig. 5c) and SEM (Fig. 5e), where short-range compacted coils were also obtained (Fig. 5d, f). Because the running paths are too dense to observe the details in the 10-μm design (Fig. 5b), the overall schematic diagram of G-code with 100 μm spacing is shown in Fig. S13a, which can be helpful to observe the line and piece details clearly. The units of arcs and lines were used to fill in the simplified rectangular, parallelogram, round geometries with different sizes as small pieces. Then, these small pieces were assembled and connected with continuous printing to construct the whole vertebrae structure. The running path and corresponding zoomed-in SEM images of the lumbar spine-like scaffold (Fig. 5e) in different regions are shown in Fig.  S13b-f, and all regions exhibited a similar coiling density. At the speed of 200 mm/min, the pore size was 9.13 ± 4.43 μm, as shown in Fig. S7b, which is similar to the lacunar-canalicular pores surrounding the osteocytes that are in the range of 0.1-10 μm 59 . After culturing with EGFP-hMSC-TERT cells for 7 days, the SEM images (Fig.  5g, h) revealed cell spreading on the scaffold (Fig. 5h). The lumbar spine-like scaffold-patterned EGFP-hMSC-TERT cells monitored by fluorescence microscopy on day 1 (Fig.  5i, j), day 3 (Fig. 5k, l), day 5 (Fig. 5m, n) and day 7 (Fig. 5o, p) showed cell proliferation over the culture time.

Conclusion
In summary, a hierarchically ordered 3D coil compacted scaffold was successfully prepared by using the sMEW method, in which variable long-range patterns and shortrange tunable coil density can be readily achieved. Longrange macroscale shape mimicry is aimed at macroscale support for tissue engineering; moreover, with its abundant surface area and tunable mechanical properties, short-range tailored coil compactness can enhance microscale, cell-cell spatiotemporal modulation, and cellextracellular matrix interactions. In addition to in vivo biocompatibility, hierarchically ordered scaffolds hold great potential for tissue engineering. Such a spatial arrangement of the cell containing coil structural building blocks could potentially facilitate the construction of complex tissues, such as lumbar vertebrae.
Furthermore, cells growing along the patterned biomaterials could generate a cellular microarray, which not only facilitates building tissues with similar or reinforced structural and functional characteristics 60,61 but also offers the opportunity to model disease development and drug screening for precision medicine [62][63][64][65] . Inspired by the cellular microarray reported for high-throughput toxicology assays 66 and anticancer drug screening 67 and assisted by metachip or a microfluidic chip, 3D coil compacted scaffolds with a cellular microarray might be used for early-stage drug discovery screening of drug candidates 68 . It is also believed that spatially patterned scaffolds can provide unique systems to study cell behavior under complex microenvironments [69][70][71] .

Supporting information
Supporting Information is available from the Online Library or from the author.