$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
This protocol follows all the requirements set by the Environmental Health and Safety department at the University of Arkansas.
1. Set Up the Tissue Handling Area
- Take a stainless-steel metal tray and cover it with the biohazard bag as shown in Figure 1. Any handling of the biological tissues will be performed within the tray area (i.e., the tissue handling area).
- Prepare laboratory tweezers, tissue wipes, paper towels, filter paper pack, tissue dye bottles, bleach bottle, and ethanol bottle around the tray for easy access when required. Keep any used tissues, wipes, and gloves on the biohazard material surface to dispose of at the end of the protocol.
- Fill a 50 mL centrifuge tube with up to 45 mL of 10% neutral buffered formalin and place it in the centrifuge storage tray near the tissue handling tray.

Figure 1: Setup of tissue handling area. Please click here to view a larger version of this figure.
2. Handling Fresh Breast Cancer Tumor for THz Transmission Spectroscopy
CAUTION: Before handling any live tissues, put on nitrile hand gloves, eye protection goggles, a face mask, and a lab coat. Always use laboratory tweezers to handle tissues and avoid touching them directly with the hands. All work with fresh tissue outside of a sealed container or the scanning stage should be conducted at the tissue handling area established in step 1.1.
NOTE: All tissues handled in this work were shipped in Dulbecco's Modified Eagle's medium (DMEM) and antibiotic solution from the biobank.
- Remove the bulk tumor from the DMEM solution and place it in a Petri dish on the tissue handling area (see Figure 2A).
- From gross inspection, identify distinct tumor regions from which to slice small pieces for transmission characterization. Cut a 0.5 mm thick segment of tumor from the identified points using a stainless steel low profile blade, as shown in Figure 2B. Place this sliced section between two quartz windows with a spacer of 0.1 mm thickness in a liquid sample holder, as shown in Figure 2C.

Figure 2: Tumor sectioning for the THz transmission spectroscopy measurements. (A) Photograph of the bulk tumor. (B) Photograph of the small sections (0.5 mm) of the tumor cut from the bulk tumor. (C) The sliced tumor section placed in the liquid sample holder between the two quartz windows with a 0.1 mm polytetrafluoroethylene spacer for spectroscopy measurement. Figure republished from T. Bowman et al.18 with permission from SPIE. Please click here to view a larger version of this figure.
3. THz Transmission Spectroscopy Measurements
- Set the transmission spectroscopy module inside the THz core chamber by aligning the module handles over the mounting posts in the core system and sliding the stage down into the system. Tighten the two mounting screws in the upper right and lower left corners of the module as shown in Figure 3A.
- Purge the system with dry nitrogen gas at 5 L/min (LPM) during the entire spectroscopy procedure to remove water vapor from the sample space.
- Open the THz transmission spectroscopy measurement software from the desktop connected to the THz system. It will open up the main window.
- Click on the Scan tab on the top of the window. A Spectra Scan Setup window will appear. From the drop-down menu of the Measurement Mode tab on the top right of the window, select Transmission to set up transmission spectroscopy. If the peak is not automatically visible, check the Enable option under the Manual Peak Search tab and manually step the optical delay to bring the peak into view.
- After 30 min of purging, record an air reference signal by following the steps below.
- Under the Scan Settings tab in the spectra scan setup window, input an appropriate Name for the reference file, set Num Scans to 1,800, and set the Start Delay (s) to 0. Leave the other settings as their default values.
- Click on Measure Reference in the scan setup window to take the air reference measurement. Then click on Measure Sample to measure the transmission signal through air as a sample average of 1,800 signals over ~1 min.

Figure 3: THz transmission spectroscopy module setup. (A) THz core chamber with the transmission module mounted on it. (B) A photograph of the liquid sample holder. (C) The sample holder placed inside the core chamber for the measurements. Please click here to view a larger version of this figure.
- Measure the two quartz windows in the liquid sample holder as shown in Figure 3B.
- Place the two quartz windows in the liquid sample holder without a spacer in between.
- Open the THz core chamber. Mount the liquid sample holder on the transmission spectroscopy module, as shown in Figure 3C. Close the chamber.
- Click on the Scan tab on the main window. Repeat steps 3.5.1–3.5.2 for the quartz sample, but update Start Delay (s) to 900. This allows time to purge any water vapor before measurement.
- If the quartz is desired as a reference for additional samples, click on the Clear Reference tab under the Scan Settings. This clears the air reference. Then click on the Measure Reference tab to record the quartz measurements as a new reference.
- Place the sliced tumor section between the two quartz windows inside the liquid sample holder and position the holder inside the chamber for a single point transmission measurement of the tissue. To record the measurement, repeat step 3.6.3.
- Take the liquid sample holder out of the chamber when the measurements are completed and bring it to the area designated for tissue handling. Disassemble the liquid sample holder, wipe the tumor section from the quartz windows with the tissue wipes, and place the used tissue wipes in the same tray to dispose in the biohazard bag along with the other biohazard waste.
- Repeat steps 2.2, 3.7, and 3.8 as necessary to characterize additional tumor slices. When the measurements are completed, go to the main window and click on the File tab to save the measurement data. Close the software window.
4. Handling Fresh Breast Cancer Tumor for THz Reflection Mode Imaging
- Remove the fresh tumor sample from the DMEM and antibiotics solution and place it on a Petri dish. Using gross inspection, select a side of the tumor to be imaged that is sufficiently flat and has little blood and few blood vessels. Avoid imaging tissue with blood or blood vessels if possible.
- Place the tumor with the side to be imaged on grade 1 filter paper to dry the excess DMEM and clear the tissue of fluid or secretions from the tumor, as shown in Figure 4A. Reposition the tumor on the filter paper to a dry spot as the paper saturates. Dry the tumor for ~5 min.

Figure 4: Fresh tumor sample preparation for THz imaging. (A) Tumor placed on filter paper to dry. (B) Tumor placed on polystyrene plate over the imaging window with tissue wipe pads to absorb excess fluids. (C) Tumor viewed from below to track orientation and check for air bubbles. Please click here to view a larger version of this figure.
- Unmount the transmission spectroscopy module and set the reflection imaging module (RIM) mirror base on the THz core system as shown in Figure 5A. Upon setting the mirrors, mount the RIM scanning stage above the mirror base and screw it into the core system (see Figure 5B).
- Purge the system with dry nitrogen gas at 5 LPM for 30 min prior to the imaging procedure to remove water vapor from the sample compartment. After 30 min, reduce the amount of dry nitrogen gas to 3 LPM for the rest of time the system is in use.
- Place a polystyrene plate of thickness ~1.2 mm on the scanning window of diameter ~37 mm. Center the scanning window along with the polystyrene plate on the sample stage.

Figure 5: System setup for reflection imaging. (A) Reflection imaging module mirror base. (B) Scanning stage. Please click here to view a larger version of this figure.
NOTE: Other thicknesses and plate materials are suitable for step 4.5 but should have a uniform thickness and be of low enough absorption to not impede the THz signal.
- Open the THz reflection imaging measurement software from the desktop connected to the THz system. A window will pop up showing several dialog icons for specific functions and two subwindows for THz field plots, (arbitrary units a.u.) against the time and frequency, respectively.
- To set the parameters for the RIM set-up, click on the Image Parameter Dialog icon at the top of the window. An Image Acquisition Parameters window will pop up. Select RIM from the drop-down menu of the Template tab for reflection imaging set up. Hit OK and go back to the main window of the software.
- On the main window, click on the Fixed-point Scan icon. This will activate the THz antennas to start sending the incident THz signal and receiving the reflected THz signal from a single point on the polystyrene plate.
- Click on the Motor Stage Dialog icon on the top of the main window. The motor control window will open up. Adjust the optical delay axis by clicking on the forward/reverse direction arrows to center the reflected pulse from the polystyrene in the main window.
NOTE: After adjusting the optical delay axis, two pulses should appear on the window, as shown in Figure 6: one from the lower interface of the polystyrene plate (primary reflection), and one from the upper interface of the polystyrene plate (secondary reflection).
- Window out the primary reflection from the polystyrene plate and keep the secondary reflection in the window, which will contribute to the reflections from the tissue during the imaging procedure. This is done in two steps.
- First, click on the DAQ Settings button at the top of the main window to open the DAQ settings dialog window. Change the optical delay value from 5 V (default) to 4 V.
- Second, adjust the scanning stage’s vertical position with the micrometer scale on the scanning stage until the minima of the secondary pulse is the strongest. Adjust the optical delay of the axis in the Motor Control Window to put the primary reflection outside of the range of the reflected signal being measured.
NOTE: For a 1.2 mm thick polystyrene plate, the primary reflection is windowed out when the secondary reflection minimum peak is approximately -0.3 mm on the optical delay axis of the time domain window.

Figure 6: THz reflections from the lower and upper interfaces of the polystyrene plate. (A) THz signal incident to and reflected from a 1.2 mm thick polystyrene plate. (B) Measured primary and secondary THz time domain signals from the polystyrene. Please click here to view a larger version of this figure.
- Level the sample stage and record the reference signal.
- Select two points on each axis (A-axis and B-axis) that denote locations on the polystyrene plate near the edge of the sample window. For example, for the A-axis ranging from -15 mm–15 mm, the two position points can be -10 mm and 10 mm; and for the B-axis ranging from -15 mm–15 mm, the two position points can be -10 mm and 10 mm.
- Click on the Motor Control Dialog button to open the motor control window. Reposition the motor control window and the main software window so that the time domain signal is visible while adjusting the motor positions. Set both the A-axis and B-axis to 0 mm.
- Level the A-axis using following steps. A -10 mm–10 mm range is used as an example.
- In the Motor Control Window, change the value of the A-axis from 0 to -10 and hit Enter. The stage moves to the -10 mm position on the A-axis and a shift in the signal position on the main window is observed.
- Use the adjustable micrometer scale on the scanning stage shown in Figure 5B to move the minimum peak of the signal back to the position set in step 4.10.2.
- Change the A-axis value to +10 and hit enter. The stage will now move from the -10 mm position to the +10 mm position on the A-axis and a shift in the signal is observed again. Note the direction and the distance that the signal shifted from its previous position and change the A-axis value again to -10. The signal will go back to the position set in step 4.11.5.
- Rotate the leveling screw on the A-axis of the scanning stage, as shown in Figure 5B and shift the signal to double the distance in the same direction it moved from the original position. Use the micrometer on the scanning stage to shift the signal back to the original position (-0.3 mm for 1.2 mm of polystyrene).
- Repeat steps 4.11.6–4.11.7 until the signal at +10 and -10 are equal and the peak for both positions is focused at the original position (-0.3 mm on the optical axis).
- Once the leveling of the A-axis is achieved, change the A-axis value to 0 and repeat the same procedure for the B-axis. Start by changing the value of the B-axis on the motor control window from 0 to the most positive value (for example +10 mm). Also, while leveling, use the leveling screw on the B-axis of the scanning stage, which is shown in Figure 5B.
- Once both axes are leveled, return both the A-axis and the B-axis to 0 mm. Close the Motor Control Window and verify that the signal is in its original position in case it is shifted a little.
- Record this signal as the reference.
- Go to the set DAQ Properties window. Change the averaging value to 5 and keep all other parameters as default.
- Click on New Reference. The averaging counter in the top right of the window will count from 0–20. Once the counter reaches 20, change the averaging value to 1 and click OK. The reflected signal from the polystyrene will be saved as the reference for any scans taken later.
NOTE: If only the THz imaging procedure has to be performed, then it is best to perform steps 4.3–4.14 before taking the tumor tissue out of the DMEM solution.
- Mount the tumor on the polystyrene plate covering the scanning stage window.
- Remove the imaging window from the scanning stage and bring it to the tissue handling area. Place the tumor on a polystyrene plate, as shown in Figure 4B.
- Ensure that there are no significant air bubbles between the plate and the tumor. If air bubbles are observed, press the tumor with tweezers or lift the tumor and gently roll it onto the polystyrene until the air gaps are minimized.
- Place absorptive spacers at regular intervals around the test sample as shown in Figure 4B. Place another polystyrene plate above the tumor and press gently in order to make the tumor surface as flat as possible. Tape down this polystyrene-tumor-polystyrene arrangement on the sample window.
- Flip the sample window as shown in Figure 4C, and take photos of the tumor to keep a record of its orientation. Return the sample window with the tumor to the scanning stage.
- Click on the Image Parameter Dialog button to open the Image Acquisition Parameters window. Set the values of Axis1min, Axis1max, Axis2min, and Axis2max to fully enclose the position of the tumor in the imaging window
NOTE: By default, Axis1 is the A-axis and Axis2 is the B-axis.
- Set Axis1step and Axis2step to 0.2 mm for the imaging scan.
NOTE: Setting the Axis1step and Axis2step will set the stepper motors' step size to 200 µm increments during the scanning process. The total scan time can be estimated in the Image Acquisition Parameters window.
- Click on the Measure tab on the main window and select the Flyback 2D Scan option. In the window that pops up, indicate the directory and file name under which to save the scan data.
5. Postprocessing the Fresh Tissue in Preparation for Histopathology Procedure
- Upon completion of the scanning process, remove the sample window, polystyrene plates, and sample from the core THz system and move them to the area designated for hazardous waste. Remove the tumor from the polystyrene plate and place it on a flat piece of cardboard of a size comparable to that of the tumor. Make sure the orientation of the tumor is the same as it was on the polystyrene, with the imaging face touching the cardboard.
- Dip a cotton swab in red tissue dye and stain the left side of the tumor down to where the edge of the tumor contacts the cardboard. Similarly, stain the right side of the tumor with blue tissue dye. Stain the exposed surface of the tumor with a line of yellow tissue dye connecting the red stain to the blue stain to denote the back of the sample, as shown in Figure 7A.
NOTE: To prevent the ink from staining the formalin solution, apply only a thin layer to the tissue. This can be accomplished by dabbing the cotton swab on a different surface before staining the tissue or using a clean cotton swab to wipe off any excess dye. Avoid letting the dye contact the skin or clothing. This tumor-staining process is conducted as a reference to provide information about the tumor’s imaging side and its orientation to the pathologist.

Figure 7: Post processing on the tumor after THz imaging. (A) Tumor placed face down on cardboard holder and dyed with tissue marking dye. (B) Filter paper placed over tumor and taped to maintain contact. (C) Stained tumor fixed on the cardboard immersed in 10% neutral buffered formalin solution and sealed with parafilm. Please click here to view a larger version of this figure.
- Let the ink dry for around 3–4 min. Cut a piece of filter paper with the same approximate dimensions as the cardboard. Place it on the tumor and wrap a piece of tape completely around the filter paper and cardboard as shown in Figure 7B. The tape and filter paper should secure the tumor against the cardboard without applying any significant pressure.
- Immerse the stained tissue affixed to the cardboard in 10% neutral buffered formalin solution and seal the centrifuge tube using a paraffin film, as shown in Figure 7C. Designate the sample number, date, tissue type, and tumor number for the sample on the tube label. Send the tumor to the pathologist for further histopathology processing.
6. Hazardous Waste Disposal
- Collect all the waste from the tissue handling tray along with the biohazard bag used to cover the tray and put it in a new biohazard bag, as shown in Figure 8. Bring the bag to the designated biohazardous waste area in the building and set an appointment with the Environmental Health and Safety (EH&S) department for the waste pickup. Clean the tissue handling tray and the surrounding area on the table with 10% bleach solution and ethanol.

Figure 8: Photograph of the biohazardous waste bag. Please click here to view a larger version of this figure.
- Take the liquid sample holder with the spacers and quartz windows, sampling window on which tumor was mounted, polystyrene plates, and laboratory tweezers to the washing area. Rinse all materials with water and then 10% bleach solution, wiping with paper towels as necessary to remove tissue debris. Rinse again with water, scrub with alconox solution, and rinse thoroughly. For glass and plasticware, rinse in 70% isopropyl alcohol and set aside to dry.
NOTE: Once the tumor is in formalin and the sample space is clean, data processing can be handled at the same time as imaging or a later time.
7. Data Processing to Construct THz images
- Export the saved .tvl data files from the THz system. The raw data files obtained from the system are written in Python and are best read in Python before saving as MATLAB data files.
- To construct the THz image of the scanned fresh tissue, convert the raw time domain reflection imaging data into the frequency domain using Fourier transform on the third dimension of the raw data matrix (i.e., the time dimension). Also take the Fourier transform of the reference data.
NOTE: A typical frequency domain spectrum should provide data ranging from 0.1 THz–4 THz.
- Normalize the sample data with the reference data and perform the power spectra based on the integration of the normalized data over the frequency range from f1 = 0.5 THz to f2 = 1.0 THz using the following equation19:

NOTE: Here Esample is the frequency domain reflection imaging data of the tissue sample and Ereference is the frequency domain of a single point reflection data of the reference signal.
- Construct the two-dimensional image by plotting the calculated power spectra data at each point in the matrix defined by the A-axis and B-axis. This is known as the power spectra THz image.
NOTE: The method to obtain a tomographic THz image instead is detailed in steps 7.5–7.7.
- For characterization, calculate the theoretical frequency-dependent reflection for a range of potential tissue properties using the following equation18:

NOTE: Here ρT,ij is the complex Fresnel reflection coefficient between region i and region j;dj is the thickness of region j; and θj is the angle of propagation in region j related to the angle of incidence by Snell’s Law.
is the complex propagation coefficient in region j, where ω is the angular frequency, c is the speed of light in vacuum, nj is the real part of the refractive index, and αabs,j is the absorption coefficient18. Region 1 is air, Region 2 is the polystyrene plate, and Region 3 is the tissue.
- Calculate the reflection in equation (2) for a range of user-defined refractive indexes and absorption coefficients for Region 3 (n3 and αabs,3) and compare with the measured signal at each point to calculate the combined mean squared error for the magnitude and phase.
NOTE: The solution for the refractive index and absorption coefficient is the pair of values that give the lowest error.
- Construct the tomographic THz image from the extracted refractive index and absorption coefficient data (n3 and αabs,3) at each pixel. Analyze the tumor regions by comparing with the pathology slide image obtained from the pathologist. Representative results are shown in Figure 9, with examples of insufficient adherence to the protocol in Figure 10 and Figure 11.
8. Extraction of Electrical Properties of the Tissue Using Transmission Spectroscopy Data
- On the main window of the THz transmission spectroscopy measurement software, go to the File tab and click on the Export option. A window will pop up to select the Data Type and Sample to export. Choose Transmittance and Transmittance Phase data types for the quartz and tissue sample measurements.
- Calculate the theoretical frequency-dependent transmission for a range of potential tissue properties using the following equation15:

NOTE: Here
is the ratio between Fresnel transmission coefficients for the sample and reference setups; γ1 and γ3 are the complex propagation constants of air and tissue, respectively; and d is the thickness of the tissue. The propagation constant in general is defined as
. ñ is the complex refractive index defined as
, where n is the real part of the refractive index; c is the speed of light; ω is the angular frequency; and αabs is the absorption coefficient15.
- Calculate the combined mean squared error between the magnitude and phase of the transmission in equation (3) and the measurement data from the system for a range of user-defined n and αabs values.
NOTE: The solution for the refractive index and absorption coefficient is the pair of values that give the lowest error.
- Plot the extracted refractive index and absorption coefficient data against the frequency range from 0.15–3.5 THz. Representative results are shown in Figure 12.