<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>A portable pen-sized instrumentation to measure stiffness of soft tissues in vivo</title></titleStmt>
			<publicationStmt>
				<publisher></publisher>
				<date>12/01/2021</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10356603</idno>
					<idno type="doi">10.1038/s41598-020-79735-8</idno>
					<title level='j'>Scientific Reports</title>
<idno>2045-2322</idno>
<biblScope unit="volume">11</biblScope>
<biblScope unit="issue">1</biblScope>					

					<author>Zhengwei Li</author><author>Alireza Tofangchi</author><author>Robert A. Stavins</author><author>Bashar Emon</author><author>Ronald D. McKinney</author><author>Paul J. Grippo</author><author>M. Taher Saif</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[Abstract            Quantitative assessment of soft tissue elasticity is crucial to a broad range of applications, such as biomechanical modeling, physiological monitoring, and tissue diseases diagnosing. However, the modulus measurement of soft tissues, particularly in vivo, has proved challenging since the instrument has to reach the site of soft tissue and be able to measure in a very short time. Here, we present a simple method to measure the elastic modulus of soft tissues on site by exploiting buckling of a long slender bar to quantify the applied force and a spherical indentation to extract the elastic modulus. The method is realized by developing a portable pen-sized instrument (EPen: Elastic modulus pen). The measurement accuracies are verified by independent modulus measures using commercial nanoindenter. Quantitative measurements of the elastic modulus of mouse pancreas, healthy and cancerous, surgically exposed but attached to the body further confirm the potential clinical utility of the EPen.]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><p>Considerable progress has been made over the past two decades in measuring elastic modulus of hard materials at macro <ref type="bibr">1,</ref><ref type="bibr">2</ref> , micro <ref type="bibr">3,</ref><ref type="bibr">4</ref> and nano <ref type="bibr">5,</ref><ref type="bibr">6</ref> scales. However, there is a long-term unmet need for measuring the local elastic modulus of soft materials and tissues on site, including cell culture substrates for mechanobiological studies <ref type="bibr">7</ref> and drug screening <ref type="bibr">8</ref> , engineered scaffolds for tissue engineering <ref type="bibr">9</ref> , bone regeneration <ref type="bibr">10</ref> and biohybrid machines <ref type="bibr">11,</ref><ref type="bibr">12</ref> , and in situ healthy and pathological tissues for disease diagnosis and prognosis <ref type="bibr">13,</ref><ref type="bibr">14</ref> . In contrast to nonliving hard materials, most soft biomaterials and tissues, engineered or in situ, do not lend themselves to fabrication of test specimens with standard geometry used by current material testing instruments. In addition, there are many challenges associated with the elastic modulus measurements. For example, the methods to measure their modulus need to satisfy the following conditions in regards to the instrument: (1) has to reach the site of the materials or tissues (not the other way around), therefore, the instrument has to be portable and preferably hand-held, and (2) must be able to measure the modulus in a very short time to preserve the change of mechanical properties, particularly for the soft tissues after their surgical resection from the body or after opening the body during a surgical procedure, and (3) should be versatile for a wide range of stiffness measurements for different types of soft tissues or different stages of pathological tissues (e.g. malignant tumors) <ref type="bibr">15</ref> . Currently, there are no method available to meet these criteria.</p><p>Indentation <ref type="bibr">16,</ref><ref type="bibr">17</ref> has been widely used to measure the mechanical properties of biological tissue systems. Much of the reported literature work focused on the skin and subcutaneous tissues, such as human forehead skin <ref type="bibr">18</ref> , residual limbs <ref type="bibr">19</ref> , fibrosis neck tissues <ref type="bibr">20</ref> and plantar foot tissues <ref type="bibr">21</ref> . Several generations of indentation instruments have been developed which usually employ a load cell to measure the loading force and a displacement transducer to record the tissue deformation <ref type="bibr">[22]</ref><ref type="bibr">[23]</ref><ref type="bibr">[24]</ref><ref type="bibr">[25]</ref> . However, the structures of their mechanical indentation apparatuses are either complicated or large which are not convenient for in vivo tests, such as that of soft pancreatic organ. Some portable, hand-held indentation instruments <ref type="bibr">[26]</ref><ref type="bibr">[27]</ref><ref type="bibr">[28]</ref> have also been reported. For example, an ultrasound indentation system <ref type="bibr">26</ref> has been developed for measuring biomechanical properties of plantar-foot tissues and breast tissue in vivo. A portable pen sized indentation probe <ref type="bibr">27</ref> has been reported by tracking the ultrasound echo that is reflected from the substrate bone to measure the tissue deformation. However, these ultrasound indentation systems are very expensive, sensitive to the noise, and not suitable for the measurements of thick tissues due to the attenuation of ultrasound <ref type="bibr">28</ref> . Furthermore, almost all current indentation systems have a fixed stiffness of the indenter which greatly hinder their usage for diverse modulus measurements.</p><p>Here, we propose a simple and novel method for measuring elastic modulus of soft materials or tissues on site. Such method is based on two well-established theories of mechanics: (1) buckling of a long slender bar, where the critical threshold force is determined by its geometry and elastic properties, and (2) Hertz contact theory between a hard sphere and soft materials, which offers the relationship between the applied force and deformation (i.e., indentation) of the soft materials or tissues. Our method exploits buckling to quantify the applied force on the hard sphere, and spherical indentation on the soft materials to extract the elastic modulus. The method is realized by developing a hand-held portable instrument (EPen: Elastic modulus pen) with a wide measurement range, and validated by measuring elastic modulus of Polyacrylamide (PA) gel with known modulus values (0.5, 2.5 and 5 kPa). The measured elastic modulus by the EPen matches well with those by both commercial nanoindenter and standard fabrication protocol. The potential clinical application of our EPen is further demonstrated by measuring the modulus in situ in living tissues: healthy and cancerous pancreatic tissues of healthy and KPC mice respectively.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Results</head><p>Theory. The method exploits two well-established theories of mechanics to measure the elasticity of the soft materials or tissues (Fig. <ref type="figure">1</ref>). (1) Buckling of long slender bar: a bar subjected to an axial compressive force P at both ends buckles when the force exceeds a critical threshold value P cr . When P &gt; P cr , the bar continues to buckle with very small increase of force &#948;P(&#948;P &lt; &lt; P cr ) due to the ignorable axial stiffness <ref type="bibr">29,</ref><ref type="bibr">30</ref> . Therefore, P cr can be considered as the force imposed by the slender bar at any post buckled state. This introduces a small error in the measurement of elastic modulus, yet can be quantified (See Supplementary Note S1). ( <ref type="formula">2</ref>) Hertz contact theory: a hard sphere in contact with an elastic material subjected to a compressive force deforms the material or indents it. If the deformation is small, the soft material recovers to its original configuration upon removal of the force, and the amount of indentation depends on the contact load, geometry of sphere and elasticity of the indented material. Consider a long slender bar of length l and dimeter d with one end fixed to a solid glass bead of radius r and the other end held by a movable piston, as shown in Fig. <ref type="figure">1a</ref>. Firstly, the sphere is brought in contact with the surface of a soft sample. Then, the piston is used to push the bar towards the sample until the sphere indents the sample by &#948; i when the bar begins to buckle (Fig. <ref type="figure">1b</ref>). The critical buckling force P cr is determined by the geometry and elasticity of the bar and given by, where the moment of inertia, I = &#960;d 4 64 and E g is the elastic modulus of glass fiber. When P &gt; P cr , the bar continues to buckle with a buckling angle &#952; with a small increase of force &#948;P (Fig. <ref type="figure">1c</ref>). Hence, P cr could be consid- ered simply as the force supported by the bar at any post-buckling state. The total translation of the piston from the initial contact to buckling is &#948; m = &#948; i + &#948; b , where &#948; b is the part of the translation needed to accommodate the buckling (Fig. <ref type="figure">1b</ref>) which can be determined from the deformed shape of the bar after buckling. Once the length l is given, &#948; b can be obtained analytically in the function of buckling angle &#952; , i.e., &#948; b = 0.16 &#8226; l &#8226; sin 2 (&#952;) (See Supplementary Note S2). Therefore, indentation &#948; i on the soft materials or tissues is obtained from &#948; m and &#948; b by, Although the bar shortens elastically by a small amount due to the axial force, the axial shortening is negligible compared to the displacement &#948; b due to the buckling (see Supplementary Note S3).</p><p>Hertz contact theory is then used to extract the elastic modulus of soft material based on the indentation on the soft materials or tissues. A hard sphere of radius r contacting an elastic substrate with a contact force P , deforms the material elastically by &#948; i , as long as the strains induced are within elastic limit of the indented mate- rial (Fig. <ref type="figure">1d</ref>). The force-indentation relation is given by where E s and &#957; are the elastic modulus and Poisson's ratio of the soft material. It is known that most soft materials and tissues remain linearly elastic up to a strain of 10% <ref type="bibr">31,</ref><ref type="bibr">32</ref> , and have Poisson's ratio between 0.4 and 0.5 <ref type="bibr">33,</ref><ref type="bibr">34</ref> . We use &#957; = 0.5 for all materials to be tested similar to other standard nanoindenter testing machines (e.g. Piuma nanoindenter). It should be noted that this assumption introduces a small error in the measurement of elastic modulus E . Alternatively, one can consider E (1-&#957; 2 ) as an effective elastic modulus of the material without introducing any value for &#957; . With the known values of applied force and indentation, we can evaluate the elastic modulus of soft tissues by It should be noted that our method exploits buckling as a threshold force readout and avoids using any force sensor.</p><p>Implementation. The method is realized by developing a portable pen-sized instrument that is EPen: Elastic modulus pen (Fig. <ref type="figure">2a,</ref><ref type="figure">b</ref>). It consists of a long slender glass fiber with 80 &#956;m diameter encased in a glass/acrylic hollow cylinder (i.e. outer tubing) which serves as the body of the EPen (Fig. <ref type="figure">2a</ref>). One end of the glass fiber is fixed to a glass bead/sphere with 2 mm diameter, and the other end is connected to a spring-loaded piston, the position of which can be precisely controlled by a Vernier micrometer. Two concentric tubes (i.e. glass tube and metal tube) are used to precisely align the glass fiber along the axial direction, as well as to selectively create "effective length" for buckling and thus an adjustable stiffness. The glass tube is fixed at its upper end to the plunger section, while the metal tube is attached at its lower end to a stainless steel piston enclosed in the outer acrylic cylinder tubing. The motion of steel piston (Fig. <ref type="figure">2d</ref>) is controlled by a small magnet outside the enclosure adjacent to it, which can move the metal tube inside the concentric glass tube. This unique feature enables to create variable effective glass fiber length for the buckling and thus adjustable stiffness, with one end clamped near the piston at desired position and the other end probing the substrate. Longer the glass fiber length, lower is the critical buckling load. We made the length as an adjustable parameter so that the same EPen can be used to measure elastic modulus of samples with diverse stiffness. A pair of arc-shaped spacer/aligner (Fig. <ref type="figure">2f</ref>) are further designed and attached to the bottom end of EPen body to guide and retain the plane of buckling for the sake of consistency and accurate optical imaging of the glass fiber.</p><p>To measure the elastic modulus, the EPen is held vertically by a clamp (Fig. <ref type="figure">2c</ref>). The sample is brought in proximity of the glass sphere by adjusting the clamp height. Once approached close enough, the Vernier micrometer of the EPen is then used to bring the glass sphere in contact with the sample. And a camera is used to visualize the fine motion of the sphere and monitor the contact. The micrometer is then used to move the glass fiber further by &#948; m (Fig. <ref type="figure">1c</ref>) until the bar buckles with a buckling angle &#952; that can be obtained directly based on the image of buckled glass fiber (Fig. <ref type="figure">2e</ref>). The total effective length of glass fiber l can readily be determined based on a scale ruler in mm range (Fig. <ref type="figure">2b</ref>) attached to the EPen body. The total translation &#948; m of the glass fiber is directly deter- mined by taking two measurements of the Vernier micrometer by the buckling bar before and after application of force on the sample. The buckling rotation angle &#952; is obtained by two images of the glass fiber taken before and after the buckling, which gives &#948; b . Thus, the indentation &#948; i is obtained by &#948; m -&#948; b , which provides the elastic (1)</p><p>i modulous of soft materials or tissues. Therefore, the entire testing process involves taking two measurements and two images, which takes less than 1 min.</p><p>Calibration. The critical buckling load P cr of glass fiber in EPen is given by Eq. (1). However, the load is for an ideal bar without any imperfection, i.e., there is no initial curvature in the fiber (perfectly straight) and the contact force between the glass sphere and glass fiber is aligned with its longitudinal axis. In reality, these conditions might be violated, and the glass fiber may show transverse deformation before the force is applied. Hence, the glass fiber with a sphere need to be calibrated experimentally for its "real" force-deformation relation. We calibrate the glass fiber with a prescribed length by a force sensor, namely a micro balance (Fig. <ref type="figure">3a</ref>). We use the EPen to apply force on the balance until the bar buckles along the transverse direction. The critical buckling force P cr can be directly read from the micro balance. By calibrating the buckling load for different glass fiber length l , we get a relation between l and P cr (Fig. <ref type="figure">3b</ref>). Corresponding theoretical relation is also shown in Fig. <ref type="figure">3b</ref>. As expected, the theoretical prediction is slightly higher (8.2-15.5% within the fiber lengths considered) than the experimental values, since the theory does not account for (a) any imperfection in the buckling bar and (b) any misalignment between the bar and the glass sphere. In reality, these conditions are violated. Therefore, we calibrated the buckling load for different glass fiber lengths. When P &gt; P cr , the bar continues to buckle with very small increase of force &#948;P(&#948;P &#8810; P cr ). Simply consider- ing P cr as the applied force will introduce a small error in the measurement of elastic modulus. Therefore, we further investigate the effect of buckling angle &#952; on the applied force both experimentally and theoretically. The force scale is again used to measure the applied force with the change of buckling angles. Figure <ref type="figure">3c</ref> shows the schematic illustration of the measuring process. It should be noted that the indentation &#948; i is zero for the solid surface, so &#948; m equals to &#948; b . The accuracy of our experimental measurements, shown by the black dots in Fig. <ref type="figure">3e</ref>, is further validated by both finite element analysis (FEA) (Fig. <ref type="figure">3d</ref> and Supplementary Fig. <ref type="figure">S4a</ref>) and analytical buckling analysis (see Supplementary Note S1). The results show that the estimations of P = P cr results in up to 5% error as the buckling angle is below 5&#176;, yet the error decreases as small as 2% when the buckling angle exceeds 5&#176; (Fig. <ref type="figure">3e</ref>). Therefore, we allow the glass fiber to deflect up to 5&#176; before the read out for sake of accuracy and ease of optical recording.</p><p>Validation. We use the EPen to measure the elastic modulus of PA gel with different known moduli to validate the accuracy (Fig. <ref type="figure">4a</ref>). The gel samples are prepared by curing liquid PA gel in 3D printed circular wells of 25 mm in diameter and 5 mm in depth, following the standard fabrication protocol <ref type="bibr">35</ref> . We prepare the PA gels with nominal elastic moduli of 0.5, 2.5, and 5 kPa, and Poisson's ratio is assumed to be 0.5. Nanoindenter (Optics Piuma, 2015) with spherical tip of 25 &#181;m radius is firstly used to measure the prepared 5 kPa gel samples. The measured modulus by nanoindenter is 5.16 kPa verifying the accuracy of the fabrication protocol. EPen is then used to measure the elastic moduli of each gel sample at 5 to 8 different positions. The measured modulus by EPen for 0.5, 2.5, and 5 kPa PA gel samples are all very close (within 10%) to the values obtained by standard fabrication protocols (Fig. <ref type="figure">4b</ref>).</p><p>FEA is further used to study the indentation mechanics and simulate the spherical indentation on PA gel substrates with different elastic modulus. Figure <ref type="figure">4c</ref> gives the net indentation force,P versus the indentation depth,&#948; i , by both FEA and the EPen. FEA results agree very well with the experiment, verifying the accuracy of the EPen measurements. We also find that the boundaries of stress field induced by the indentation stay far away from the bottom surface of PA gel (Fig. <ref type="figure">4d</ref>) and maximum strain level remains under 10% (Fig. <ref type="figure">4e</ref>) when the spherical indentation is 134 &#181;m, similar to the EPen measurements for 5 kPa PA gel in the experiment. However, as the indentation increases to 500 &#181;m, the maximum strain in the substrate increases to 17.7% due to the boundary effect (Fig. <ref type="figure">4f</ref>). In the EPen measures, the indentation depths &#948; i are restricted to 250 &#181;m to ensure the elastic strains below 15% to avoid the effect of bottom boundary (See Supplementary Note and Supplementary Fig. <ref type="figure">S5</ref>). Furthermore, we find that thinner substrate requires higher applied net force (Fig. <ref type="figure">4g</ref>), causing a higher effective modulus deviating from the actual values (Fig. <ref type="figure">4h</ref>). FEA results also show that the friction or adhesion between the glass sphere and PA gel substrate does not affect the net indentation force (Supplementary Fig. <ref type="figure">S6</ref>). Here, the effect of adhesion between the sphere and the substrate on the indentation load was simulated by introducing no-slip condition. The results show that the adhesion does not noticeably affect the indentation force, especially when the indentation depth is small (SI Fig. <ref type="figure">6</ref>). Finally, we performed a 3D finite element analysis to investigate the effect of adhesion on the indentation force as the sphere rotates after buckling. Here, we applied an adhesive contact condition between the sphere and soft substrate. Thus, there is no separation or sliding between them. We used a 2 mm diameter glass sphere to indent 5 kPa substrate by 150 &#181;m. We then rotate the sphere by 5&#176; and quantify the net force needed to keep the sphere indented at a constant 150 &#181;m. The results show that adhesion between sphere and the substrate does not noticeably affect the indentation force. It increases by only about 0.35% (Supplementary Fig. <ref type="figure">S7</ref>). This suggests that the effect of adhesion or friction between the sphere and the soft substrate has negligible effect on the measurement of elastic modulus by E-Pen.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Application: measuring elastic modulus of soft tissue in vivo. We demonstrate the application of</head><p>EPen to measure the elastic modulus of pancreatic tissues of healthy and KPC mice in vivo, where the green dotted area indicates the pancreatic tissue (Fig. <ref type="figure">5a,</ref><ref type="figure">b</ref>). The detailed description for each mouse could be found in the Methods section. To perform testing, the EPen is held vertically by a clamp with a camera fixed on the side view to monitor the contact of glass sphere with the tissue. The Vernier micrometer is used to induce the buckling of glass fiber and indenting the tissue surface to extract the elastic modulus of the pancreatic tissue in vivo (Fig. <ref type="figure">5c</ref>). The measured elastic modulus of healthy and cancerous mouse pancreas by EPen are given in Fig. <ref type="figure">5d</ref>. Mouse 1 is a control and has a healthy looking pancreas with the measured stiffness of 0.43 kPa. Mouse 2 and 3 show no statistically significant difference with the measured stiffness of 0.41 and 0.50 kPa respectively, which match with the fact that they are in early stages of disease. Interestingly, we notice that the mouse 4, which is reported "clearly diseased" by the clinic specialist, shows a perceptible increase in the stiffness of 0.78 kPa. Mouse 5 exhibits the smallest stiffness of 0.33 kPa, which match the fact that it has the youngest age (2 month old) and still in the early stage of disease without symptoms. Quantitative assessments of elastic modulus of mouse pancreas, healthy and cancerous, all of which surgically exposed but stayed intact in the live body confirm the potential clinical utility of the EPen.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Discussion</head><p>A unique feature of our EPen, in contrast to existing linear spring-based method of measuring material behavior, is that the stiffness of EPen's nonlinear spring (i.e., the buckling bar) is tunable. The underlying mechanics of nonlinear buckling and indentation embedded in the nature of the EPen device makes it capable for a wide range of in-situ stiffness measurement of soft materials or tissues. To bring this notion into a quantitative perspective, Eq. ( <ref type="formula">4</ref>) gives where the glass fiber length l and indentation amount &#948; i can be considered as operational variables, sphere radius r and fiber cross-sectional diameter d as structural variables. All of these parameters influence the EPen's measurement range of elastic modulus for soft materials or tissues. In the current design, the effective length l for the buckling of the glass fiber can be easily tuned from 20 to 110 mm by means of an external magnet based on the needs. Glass fiber diameter d and the glass bead radius r can be chosen from 80 to 200 um, and from 1 to 3 mm respectively. And the indentation amount &#948; i can vary within 0.05 to 0.3 mm in the actual tests. Therefore, the same EPen can measure elastic modulus of materials that span 3 orders of magnitude. For example, an 80 &#956;m diameter glass fiber can be used to test materials with elastic modulus varying over three orders of magnitude (from 0.1 to 100 kPa). Larger glass fibers are used to measure the stiffer soft samples (Supplementary Fig. <ref type="figure">S8</ref>). Soft materials such as polydimethylsiloxane (PDMS) <ref type="bibr">36,</ref><ref type="bibr">37</ref> with elastic modulus as high as 4.2 MPa can also be tested by the same device by replacing the glass fiber of 200 &#956;m diameter.</p><p>Measurement time for a modulus measurement is less than 1 min. It involves taking two measurements of the Vernier scale and two images by the buckling bar, before and after application of force on the sample. The measurement time can be significantly reduced to few seconds if the Vernier micrometer readings are done electronically. In that case no manual measurements will be necessary. Furthermore, the EPen is inexpensive, and it does not require any specialized facility for its assembly. It only requires a straight bar with known geometry and material properties so that the critical buckling load can be estimated, a Vernier micrometer to measure its end motion as it applies indentation on a sample. Owing to these beneficial characteristics, our EPen is suitable for measuring elastic modulus of soft tissues, particularly diseased ones for prognosis. The tissues can be in situ, intact with the body. However, the tissue needs to be accessible to the probe of the EPen. Hence, EPen is appropriate for use when the tissue is exposed during surgery or after it is resected (e.g., human colon tumor <ref type="bibr">38</ref> ).</p><p>To summarize, we have developed a hand-held pen type instrument (EPen: Elastic modulus Pen) to measure the local elastic modulus of soft materials and tissues on site, where the tip of the pen serves as the modulus probe. To investigate the accuracy and reliability of the EPen measurements, we provide a systematic experimental, computational and analytical studies of the method and demonstrate that it can be used as an accurate reliable tool for rapidly measuring the mechanical properties of soft materials and tissues. Assessments of the elastic modulus of mouse pancreas, healthy and cancerous, surgically exposed but attached to the body further confirm the clinical utility of the EPen. Our findings establish routes for the design and fabrication of stiffness measurement instruments that meet requirements for clinical use.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Methods</head><p>Assembly of glass bead-fiber. In the EPen, the glass bead is glued to the glass fiber. A glass rod with 80 &#181;m diameter is firstly obtained by stripping off plastic coating from a fiber optic cable (Supplementary Fig. <ref type="figure">S1a</ref>). And a 10 cm long segment of the glass fiber is then cut out transversely using fiber cleaver (Supplementary Fig. <ref type="figure">S1b</ref>). A 3D printed aluminum funnel fixture with a hollow glass needle is used to align the fiber and the bead (Supplementary Fig. <ref type="figure">S1c</ref>). A glass sphere with 2 mm in diameter is placed on the funnel hole and glued to the glass fiber at the cut end, such that the fixture holds the sphere in a conical cavity. Finally, the glass sphere and glass fiber are assembled and glued together, guided by a 3D-printed fixture (Supplementary Fig. <ref type="figure">S1c</ref>). The fabricated product of glass bead-glass fiber assembly is shown in Supplementary Fig. <ref type="figure">S1e</ref>.</p><p>Fabrication of the EPen. The end of glass bead-glass fiber assembly is glued to a spring-based piston enclosed into an aluminum tube that is secured to the body of the EPen (Supplementary Fig. <ref type="figure">S2</ref>). The forward motion of the piston is driven and controlled by a Vernier Micrometer (SM-13LH) attached to the EPen body, while the back spring (Supplementary Fig. <ref type="figure">S2c</ref>) allows the piston to return to initial state after the micrometer motion is reversed. Two concentric glass tube and metal tube are used to ensure the alignment and translation of glass fiber along the axial direction, while the inner tube also anchored to a stainless steel piston through its central hole. The translation of this piston (Supplementary Fig. <ref type="figure">S2d,</ref><ref type="figure">e</ref>) is controlled by an external magnet that can freely move from outside the EPen glass enclosure, which thus can move the metal tube inside the concentric glass tube, and hence creating an adjustable length of the glass fiber for the buckling. A pair of arc-shaped wires are attached to the end side of EPen body to guide the direction and plane of buckling of the glass fiber.</p><p>The operation of EPen. In the measurements, the pen is held vertically by using a clamp. The samples are placed on a movable z-stage which brings the sample close to the glass bead. The Vernier micrometer is used to bring the glass bead in contact with the sample, and then move the glass fiber farther until the fiber buckles. Two measurements of the Vernier micrometer by the buckling bar before and after application of force on the sample give the total translation &#948; m of the glass fiber. A camera aligned orthogonally to the plane of arc-shaped aligner is used to visualize the fine motion of the glass bead and monitor the contact. The buckling angle &#952; is obtained based on two optical images of glass fiber in proximity of samples surface taken by the camera before and after the buckling.</p><p>PA gel preparation. Polyacrylamide (PA) gels are synthesized following a protocol by Tse and Engler <ref type="bibr">35</ref> .</p><p>Acrylamide (40% sol., Sigma-Aldrich), N, N&#8242;-Methylenebis-acrylamide (2% sol., Sigma-Aldrich) and Dul-</p><p>becco's Phosphate Buffered Saline (DPBS, Corning) are mixed at different ratios according to the protocol for various Young's moduli (0.5, 2.5 and 5 kPa). 1% (v/v) of ammonium persulfate (APS) solution (10% w/v APS, Bio-Rad) and 0.1% (v/v) of Tetramethylethylenediamine (TEMED) solution (Bio-Rad) is used to catalyze the polymerization reaction in the gel precursor solution. The gel precursor is then poured into a mold and a large coverless (Fisher Scientific) is placed on top to form a plane test surface. After polymerization for about 20 min, the coverglass is removed and the gel was rinsed 3 times with Dulbecco's phosphate-buffered saline (DPBS).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Finite element analysis.</head><p>Three-dimensional finite element model (FEM) is used to perform the buckling and postbuckling analysis of long slender glass fiber and investigate how the buckling force changes at any post buckled state (Fig. <ref type="figure">3e</ref>). Buckling analysis is firstly performed to determine the critical buckling loads and corresponding buckling modes of the long slender glass fiber with one end is pinned and the other end fully fixed.</p><p>The first buckling mode is then used as initial small geometrical imperfection to trigger the buckling. The beam elements (B21 in Abaqus) are used to discretize the geometry of glass fiber with the length of 50 mm and diameter 80 &#956;m. A total of 21,000 elements are used to ensure mesh convergence. FEA results quantitively give the buckling forces at any post buckled states with different buckling angles (Supplementary Fig. <ref type="figure">S4</ref>). A 2D axisymmetric finite element model is further established to simulate the spherical indentation on PA gel substrate and study the indentation mechanics. For simulation, four-node bilinear axisymmetric quadrilateral (CAX4) elements are used to discretize the geometry of PA gel substrate with the thickness of 3 and 5 mm. Refined meshes are adopted to ensure the accuracy. The Young's modulus and Poisson's ratio used for PA gel are E = 0.5, 2.5 or 5 kPa and &#957; = 0.49. The glass bead with the diameter of 2 mm is considered to be rigid. "Hard contact" is implemented to simulate the interface between the glass bead and the PA gel substrate. Displacement boundary conditions are assigned to the rigid glass bead to apply different levels of compression. In addition, the effect of interaction property, adhesive or frictionless, between the glass sphere and PA gel substrate is also studied. FEA results show that friction does not affect the indentation force (Supplementary Fig. <ref type="figure">S6</ref>).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Description of mice.</head><p>All animal procedures were approved by Institutional Animal Care and Use Committee (IACUC) at the University of Illinois and were conducted following the Guide for the Care and Use of Laboratory Animals. Mouse 1: 9 month old female HK2 fl/+ (no modification of allele in the absence of Cre; wild type control) with a body weight of 30.4 g. Having previously collected these control mice, their pancreas is histologically indistinct from wild type mice. Mouse 2: 10 month old p48-Cre/LSL-Kras G12D /Akt2 fl/+ female with a body weight of 20.6 g. p48Cre and LSL-Kras G12D positive (KC mice <ref type="bibr">39</ref> ) and heterozygous for Akt2 as determined by PCR genotyping. This mouse has a mild morphological phenotype in regards to its pancreas where mutant KRAS alters tissue architecture, though there could be some protective effect from reduced AKT2 to suppress disease development. Mouse 3: 9 month old p48-Cre/LSL-Kras/Akt2 fl/+ male with a body weight of 31.5 g. It has the same genotype as the mouse 2, with mild/slight pancreas stiffness. Mouse 4: 5.5 month old p48-Cre/ LSL-Kras G12D /LSL-p53 R172H/+ (KPC <ref type="bibr">40</ref> ) female with a body weight of 25.5 g. This mouse is PCR positive for both p48Cre and LSL-Kras and heterozygous for one mutant LSL-p53 R172H allele. KPC mice are the gold standard of mouse models of pancreatic cancer with a distinct cancer phenotype in the pancreas that resembles human pancreatic cancer on both morphological and molecular levels <ref type="bibr">40</ref> . Mouse 5: 2 month old KPC male (no body weight available) that is Pdx1-Cre and LSL-Kras G12D positive and heterozygous for LSL-p53 G12D . This mouse has a more mild phenotype but will likely have worse diseased pancreas within 2-3 months. These mice were selected to include the control mouse 1, which has a floxed allele that is not recombined. Ideally, this would have been Akt2, but those mice were not generated, and HK2 served as its substitute. The next two mice (2 and 3) are KC mice with loss of one allele of Akt2, which does have a histotype/phenotype similar to KC mice, and we remain interested in how loss of this allele may impact tissue stiffness. The final two mice (4 and 5) are KPC mice at older and younger ages to determine the relative stiffness during the progression of disease.</p></div><note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_0"><p>Scientific Reports | (2021) 11:378 | https://doi.org/10.1038/s41598-020-79735-8</p></note>
			<note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_1"><p>Scientific Reports | (2021) 11:378 | https://doi.org/10.1038/s41598-020-79735-8</p></note>
		</body>
		</text>
</TEI>
