Samantha M Fix

Bhanu P Koppolu

Anthony Novell

Jared Hopkins

Thomas M Kierski

David A Zaharoff

Paul A Dayton

Virginie Papadopoulou
email: papadopoulou@unc.edu

Ultrasound-Stimulated Phase-Change Contrast Agents for Transepithelial Delivery of Macromolecules, Toward Gastrointestinal Drug Delivery

Keywords: Ultrasound, Transepithelial Delivery, Acoustic Droplet Vaporization, Phase-change contrast agents, Nanodroplet

come   L'archive ouverte pluridisciplinaire

INTRODUCTION

The gastrointestinal (GI) tract has evolved to efficiently breakdown food products, allow nutrient uptake, and simultaneously provide a robust defense against ingested toxins and pathogens [START_REF] Barrett | Epithelial biology in the gastrointestinal system: insights into normal physiology and disease pathogenesis[END_REF], Schoellhammer, et al. 2016). This presents a complex barrier for both local and systemic drug delivery through the GI tract, particularly for biologics (e.g., therapeutic peptides, proteins, and nucleic acids) [START_REF] Goldberg | Challenges for the oral delivery of macromolecules[END_REF]. The first hurdle is maintaining drug stability in the degradative environment. This often requires specialized formulations to protect the active pharmaceutical ingredient from (1) the harsh and fluctuating acidic pH and (2) an array of degradative species and enzymes [START_REF] Moss | The biological challenges and pharmacological opportunities of orally administered nanomedicine delivery[END_REF]. If drug integrity is maintained, the second hurdle is permeability. The GI tract is blanketed with mucus that traps and protects the underlying epithelium from foreign particles [START_REF] Cone | Barrier properties of mucus[END_REF].

Furthermore, GI epithelial cells are linked by tight junctions that restrict paracellular transport of large molecules (> 200 Da), and the hydrophilic nature of most biologics limits passive transcellular diffusion [START_REF] Goldberg | Challenges for the oral delivery of macromolecules[END_REF]. Due to these challenges, generally only small molecule drugs can be delivered through the GI tract.

Therapeutic ultrasound has recently been proposed as a universal method to physically enhance GI permeability to macromolecular drugs (Schoellhammer, et al. 2017, Schoellhammer, et al. 2016, Schoellhammer, et al. 2017[START_REF] Schoellhammer | Ultrasound-mediated gastrointestinal drug delivery[END_REF], Schoellhammer and Traverso 2016). This technique involves colonic insertion of an ultrasound probe with co-administration of a medicated enema (Schoellhammer and Traverso 2016). Low-frequency (≤ 100 kHz) ultrasound exposure for one minute or less has been shown to significantly enhance macromolecule delivery into surrounding GI tissue (Schoellhammer, et al. 2017, Schoellhammer, et al. 2017). This technology has the potential to transform the treatment landscape of several GI diseases including ulcerative colitis, cancers of the GI tract, infection, and Crohn's disease, and offers the following unique benefits (Schoellhammer, et al. 2017, Schoellhammer, et al. 2016[START_REF] Schoellhammer | Ultrasound-mediated gastrointestinal drug delivery[END_REF]. First, rapid delivery limits concerns of drug degradation by intestinal enzymes, eliminating the need for complex protective formulations. Second, this technique may allow for local GI delivery of biologic drugs previously restricted to intravenous administration (including peptides, proteins, monoclonal antibodies, nucleic acids), which would expand therapeutic options for GI diseases. Finally, concentrated local delivery has the potential to enhance drug efficacy and reduce systemic toxicities.

The mechanism of ultrasound-mediated drug delivery involves in situ generation and acoustic excitation of vapor cavities near a biological barrier through the process of cavitation [START_REF] Lentacker | Understanding ultrasound induced sonoporation: definitions and underlying mechanisms[END_REF][START_REF] Miller | Sonoporation: mechanical DNA delivery by ultrasonic cavitation[END_REF], Schoellhammer and Traverso 2016). When a fluid is subjected to ultrasound of sufficient energy, vapor cavities are created (cavitation). These bubbles oscillate with the passing acoustic wave (stable cavitation) and cause flow disturbances in the surrounding fluid, and this 'microstreaming' exerts shear stress on nearby structures (Kooiman, et al. 2014[START_REF] Pereno | Layered acoustofluidic resonators for the simultaneous optical and acoustic characterisation of cavitation dynamics, microstreaming, and biological effects[END_REF]. With sufficient acoustic amplitude, the bubble becomes unstable, violently collapses (inertial cavitation), and creates a fluid microjet towards nearby boundaries. Stresses exerted by these physical phenomena transiently enhance the permeability of biological barriers [START_REF] Coussios | Cavitation-mediated extravasation and transtumoral drug delivery by microstreaming: What role do the gas nuclei and the physical properties of the therapeutic play?[END_REF][START_REF] Escoffre | Mini-review -Biophysical mechanisms of cell membrane sonopermeabilization: Knowns and unknowns[END_REF], Kooiman, et al. 2014[START_REF] Miller | Sonoporation: mechanical DNA delivery by ultrasonic cavitation[END_REF][START_REF] Postema | High-speed photography during ultrasound illustrates potential therapeutic applications of microbubbles[END_REF]. In addition to enhanced permeability, acoustic cavitation and radiation force enhance the transport of molecules through newly permeabilized membranes [START_REF] Bhatnagar | Exploitation of acoustic cavitation-induced microstreaming to enhance molecular transport[END_REF].

The likelihood of ultrasound to induce cavitation is related to the wave's mechanical index (MI) which is defined as the peak rarefactional pressure (in MPa) divided by the square root of the center frequency (in MHz) [START_REF] Apfel | Gauging the likelihood of cavitation from short-pulse, low-duty cycle diagnostic ultrasound[END_REF]Holland 1991, Bhatnagar, et al. 2014). To achieve ultrasound-mediated drug delivery, low frequencies (≤ 100 kHz) and consequently high MIs (>2) are often needed (Schoellhammer, et al. 2017, Schoellhammer, et al. 2017[START_REF] Schoellhammer | Ultrasound-mediated gastrointestinal drug delivery[END_REF]. However, introducing exogenous cavitation nuclei into the acoustic field substantially lowers the energy required to induce cavitation-mediated bioeffects as initial bubble formation is no longer required [START_REF] Bhatnagar | Exploitation of acoustic cavitation-induced microstreaming to enhance molecular transport[END_REF]). This can be achieved through the use of microbubble-based ultrasound contrast agents.

Microbubbles are 1-10 µm gas-filled particles stabilized by a lipid, protein or polymer shell [START_REF] Sirsi | Microbubble Compositions, Properties and Biomedical Applications[END_REF]. When subjected to ultrasound, they oscillate and/or collapse, providing biological permeabilization similar to that described above, but with acoustic parameters in line with what is employed for clinical diagnostic imaging (i.e. with MIs less than the U.S. Food and Drug Administration's limit of 1.9) [START_REF] Delalande | Sonoporation at a low mechanical index[END_REF][START_REF] Lentacker | Understanding ultrasound induced sonoporation: definitions and underlying mechanisms[END_REF]. Indeed, drug delivery by this technique can be achieved with diagnostic ultrasound imaging equipment, which provides the added benefit of simultaneous ultrasound image guidance [START_REF] Kotopoulis | Treatment of human pancreatic cancer using combined ultrasound, microbubbles, and gemcitabine: a clinical case study[END_REF].

The use of ultrasound-stimulated microbubbles for drug delivery has been studied extensively for the disruption of vascular barriers (Martin andDayton 2013, Sennoga, et al. 2017). For applications such as opening the blood brain barrier [START_REF] Burgess | Drug delivery across the blood-brain barrier using focused ultrasound[END_REF][START_REF] Wu | Characterization of Different Microbubbles in Assisting Focused Ultrasound-Induced Blood-Brain Barrier Opening[END_REF][START_REF] Zhao | Ultrasound-mediated strategies in opening brain barriers for drug brain delivery[END_REF] or enhancing vascular permeability within solid tumors [START_REF] Chowdhury | Ultrasound-guided drug delivery in cancer[END_REF][START_REF] Qin | Sonoporation: Applications for Cancer Therapy[END_REF], microbubbles are administered intravenously and are stimulated with an external ultrasound source. We envision an adaptation of this technology for GI drug delivery, which would involve co-administration of contrast agents with a medicated enema in the colon followed by ultrasound stimulation using an endoscopic probe. With the use of contrast agents, we believe that higher frequencies and lower pressures may be used compared to the previously described, contrast-free methods of ultrasound-mediated GI drug delivery. This approach may theoretically localize bioeffects more precisely, as permeabilization would only occur where contrast agents are present (e.g., at the GI wall). This would reduce concerns of unsuppressed cavitation activity outside the desired treatment area (e.g., surrounding tissues).

As a first step towards this goal, we studied the dynamics of epithelial monolayer disruption and recovery in vitro using phase-change ultrasound contrast agents (PCCAs) stimulated with 1.0 MHz ultrasound pulses. PCCAs are liquid perfluorocarbon-filled particles that can be vaporized into acoustically active microbubbles with the application of an ultrasound pulse of sufficient amplitude. In the liquid state, PCCAs are characterized by nanometer-range size distributions (100-300 nm) that may allow more thorough permeation through the GI mucus mesh compared to microbubbles (0.5-10 μm [START_REF] Sirsi | Microbubble Compositions, Properties and Biomedical Applications[END_REF]). Once vaporized, PCCAs form microbubbles several microns in diameter that can be utilized to enhance drug delivery.

Typically, PCCAs are created using perfluorocarbons with bulk boiling points near body temperature (e.g. dodecafluoropentane, b.p. = 29 °C). These agents require acoustic pressures in the range of 3-6 MPa to initiate vaporization (Burgess andPorter 2015, Williams, et al. 2013), which may induce undesired tissue bioeffects (e.g., cell lysis or tissue heating) for applications where transient drug delivery is desired. We have developed a low boiling point formulation of PCCAs filled with octofluoropropane (OFP), which has a bulk boiling point of -36.7°C. These agents can be vaporized with peak negative pressures as low as 300 kPa (at 1.0 MHz) [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF]) and with clinically available hardware. We believe that this low boiling point PCCA formulation is ideal for GI drug delivery applications where transient permeabilization is desired without permanent tissue damage.

The primary objective of this study is to demonstrate the potential of using ultrasound-stimulated PCCAs to cause transient disruption of confluent colorectal adenocarcinoma (Caco-2) epithelial monolayers and enhance the permeation of a model macromolecular drug. This represents a first step towards the goal of in vivo GI drug delivery applications.

MATERIALS AND METHODS

Phase-change ultrasound contrast agent (PCCA) fabrication and characterization

PCCAs were generated as previously described [START_REF] Sheeran | Formulation and acoustic studies of a new phase-shift agent for diagnostic and therapeutic ultrasound[END_REF], Sheeran, et al. 2012). Briefly, 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and 1,2-distearoyl-sn-glycero-3phosphoethanolamine-N-methoxy(polyethylene-glycol)-2000 (DSPE-PEG2000) (both lipids >99% purity, Avanti Polar Lipids, Alabaster, AL, USA) were combined in a 9:1 molar ratio and solubilized in a phosphate-buffered saline (PBS)-based solution (1× PBS: 11.9 mM phosphates, 137 mM solium chloride, 2.7 mM potassium chloride) containing 15% propylene glycol (v/v) (≥99.5% purity, Fisher Chemical, Waltham, MA, USA) and 5% glycerol (v/v) (≥99.5% purity, Fisher Chemical, Waltham, MA, USA) for a final combined concentration of 1.0 mg/mL. Aliquots (1.5 mL) of this solution were dispensed into 3.0 mL glass vials, the headspace air was replaced with octofluoropropane (OFP) gas (99.8% purity, FluoroMed, Round Rock, TX, USA), and microbubbles were generated via shaking for 45 seconds in a VialMix device, which operates at a frequency of 4530 ± 100 oscillations per minute (Bristol-Myers-Squibb, New York, NY, USA). Finally, PCCAs were formed via microbubble condensation through simultaneous cooling and pressurization of the microbubble suspension (Sheeran, et al. 2012). The size distribution of resultant PCCAs was characterized using a NanoSight NS500 (Malvern Instruments, Westborough, MA, USA) (detection capability 50-2000 nm) (n=3 representative vials).

Do note that during in vitro experiments performed at 37°C, a portion of the PCCA sample is lost to spontaneous vaporization due to the very low boiling point of OFP. However, we previously showed that stimulating cells with ultrasound and OFP-filled PCCAs below the vaporization threshold of the PCCAs does not cause barrier permeabilization and does not result in substantial cavitation activity [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF]. This indicates that spontaneous vaporization alone likely does not contribute a substantial enough bubble population to influence drug delivery outcome. Furthermore, in vivo imaging studies from our group have shown that while some spontaneous vaporization occurs after intravenous injection of OFP PCCAs, substantial vaporization can be achieved in a controllable manner via ultrasound-induced vaporization (Sheeran, et al. 2015).

PCCA vaporization captured via high-speed optical microscopy

PCCA vaporization events were captured at a range of acoustic settings using high-speed optical microscopy, using a previously described experimental set-up [START_REF] Sheeran | Phase-transition thresholds and vaporization phenomena for ultrasound phase-change nanoemulsions assessed via high-speed optical microscopy[END_REF]. Briefly, a temperature-controlled water bath (37°C) was mounted on an inverted microscope with a 100× water immersion objective (Olympus IX71, Center Valley, PA, USA), and the microscope was interfaced with a high-speed camera to capture vaporization events (FastCam SA1.1, Photron USA, Inc., San Diego CA, USA).

First, a calibrated needle hydrophone (HNA-0400, Onda Corp., Sunnyvale, CA, USA) was aligned to the optical focus of the microscope and subsequently used to characterize the peak negative pressure (PNP) of a 1.0 MHz unfocused piston transducer (IP0102HP, Valpey Fisher Corp., Hopkinton, MA, USA). The transducer was excited with sinusoidal pulses generated by an arbitrary function generator (AFG3021, Tektronix, Inc., Beaverton, OR, USA) and amplified approximately 55 dB with a power amplifier (A500, ENI, Rochester, NY, USA). After the transducer was calibrated, the hydrophone was replaced with a 200-μm inner diameter microcellulose tube (Spectrum Labs, Inc., Rancho Dominguez, CA, USA), which was also aligned to the optical focus of the microscope. This process ensured that contents of the microcellulose tube that were visible in the microscope focus would experience the acoustic pressure calibrated at that location. The microscope field of view was 130 µm × 130 µm.

PCCAs were diluted 1:4 in PBS and perfused through the microcellulose tube. Injection of the PCCAs was followed by a brief waiting period to allow particles in the field of view to become nearly stationary, which minimized blurring in the captured images. Previous work has shown that PCCA vaporization is a pressure-threshold dependent phenomenon, independent of pulse length when submillisecond pulses are used [START_REF] Fabiilli | The role of inertial cavitation in acoustic droplet vaporization[END_REF][START_REF] Lo | Acoustic droplet vaporization threshold: effects of pulse duration and contrast agent[END_REF]. Therefore, to determine an approximate pressure threshold for vaporization, PCCAs were exposed to single 20-cycle acoustic pulses with PNPs of 100, 200, 300, 400, 500, or 600 kPa.

A synchronization pulse from the function generator was used to trigger the high-speed camera and initiate video recording at 500 frames per second. Videos were set to begin recording just before the triggered ultrasound pulse, which allowed complete PCCA vaporization events within the optical field of view to be captured. Videos were reviewed offline using ImageJ to determine whether or not vaporization was achieved at the various acoustic excitation parameters (National Institutes of Health, Bethesda, MD, USA).

Cell culture

Human colorectal adenocarcinoma (Caco-2) cells were chosen as an in vitro model of intestinal absorption. When Caco-2 monolayers are grown to confluence, they differentiate, developing polarity and tight junctions between cells. This accurately mimics the human intestinal epithelial barrier, making Caco-2 permeability testing a reliable tool for predicting in vivo intestinal absorption of various compounds (van Breemen and Li 2005).

Mycoplasma-free caco-2 cells were purchased from the American Type Culture Collection (ATCC) and cultured using minimum essential medium without L-glutamine (Corning Inc., Corning, NY, USA) supplemented with 100 U/mL penicillin/streptomycin, 1× non-essential amino acid supplement (NEAA), 1× sodium pyruvate, 1× L-glutamate, and 20% heat-inactivated fetal bovine serum (FBS) (VWR, Radnor, PA, USA) in a humidified, 37°C incubator with 5% CO2 atmosphere. Low passage number (passage 15) Caco-2 cells were seeded at 150,000 cells per well onto Transwell supports with 0.4 µm pores (12-well plate, 12 mm membrane diameter, 1.12 cm 2 cell growth area, Corning Inc., Corning, NY, USA) and allowed to grow for approximately 10 days to form confluent monolayers. Monolayer integrity was regularly monitored via the measurement of transepithelial electrical resistance (TEER) (World Precision Instruments, Sarasota, FL, USA). Monolayers were used for permeability experiments after TEER values had reached at least 500 Ohm (560 Ohm•cm 2 ), which is similar to TEER values reported in the literature for this model [START_REF] Kudsiova | A comparison of the effect of chitosan and chitosan-coated vesicles on monolayer integrity and permeability across Caco-2 and 16HBE14o-cells[END_REF]Lawrence 2008, Srinivasan, et al. 2015).

Furthermore, above this threshold, macromolecular tracers (e.g., the 70 kDa dextran used here) are unable to permeate through the Caco-2 monolayers.

Ultrasound alignment strategy

Caco-2 monolayers were sonicated from the apical side using a 1.0 MHz unfocused piston transducer (diameter = 8 mm, -6 dB 'focal spot' = ~6 mm diameter, IP0102HP, Valpey Fisher Corp., Hopkinton, MA, USA). A custom water bath was designed to (1) facilitate consistent alignment of the transducer above the center of each well and (2) reduce acoustic reflections from the base of the well plate (Figure 1). The 14 cm×10 cm×8.5 cm (L×W×H) water bath was created using 3.2 mm acrylic sheets. For each experiment, a 12-well culture plate was placed on an internal ledge within the bath, providing coupling between the bottom of the well plate and the pre-heated 37°C water below. Water was maintained at 37°C for the duration of the experiment by placing the apparatus on a heated plate. Coupling between the well plate and water was designed to minimize acoustic reflections that would have otherwise occurred at the air-plastic interface, although we do note that it is unlikely that this strategy entirely eliminated reflections and standing waves. A lid for the water bath was created with 1.0 cm-diameter circular holes centered above each well. This served as a method to hold the ultrasound transducer centered within each well at a consistent height above the cells (6 mm). The full design of this water bath can be found in the supplemental materials, which includes complete instructions for fabrication (available in Mendeley Data repository online).

Ultrasound-mediated FITC-dextran delivery through epithelial monolayers

Twelve-well plates containing confluent Caco-2 monolayers on Transwell ® supports were positioned in the custom water bath as described above. Complete cell media (400 µl) containing PCCAs (0.5 μl; 1.75×10 9 particles/mL; roughly matching the circulating concentration of OFP-PCCAs used as vascular ultrasound contrast agents (Sheeran, et al. 2015)) and 70 kDa FITC-dextran (50 µg total; 0.125 mg/mL; similar in concentration to prior reports of colonic delivery of fluorescent dextran in mice (Schoellhammer, et al. 2017)) was added to the donor chamber (apical side) of a well to simulate absorption from the intestinal lumen. The 70 kDa dextran was chosen to model macromolecular drugs typically unable to pass the GI epithelial barrier. Subsequently, the 1.0 MHz transducer was placed in center of the well and excited with amplified sinusoidal signals from an arbitrary function generator (AFG3021C, Tektronics, Inc., Beaverton, OR, USA and 3100LA Power Amplifier, ENI, Rochester, NY, USA) (Figure 1). A pulse repetition frequency of 5000 Hz and exposure time of 30 seconds was used for all conditions. Pulse lengths of 20 cycles (10% duty cycle) or 40 cycles (20% duty cycle) were employed, and peak negative pressures were varied from 300 to 600 kPa. Control samples were exposed to (1) sham treatment (without PCCAs or ultrasound), (2) PCCAs only, and (3) ultrasound only (using the highest energy condition of 600 kPa and 40 cycle pulse length). A summary of the conditions tested can be found in Table 1, including experimental replicates for each group.

Samples (100 µl) were collected from the receiving chamber (basolateral side) at the following time points after treatment and replaced with cell media: ~5 min, 3 hr, 8 hr, 24 hr, 48 hr, and 72 hr. Dextran concentration was determined via fluorescence intensity (excitation: 485 nm, emission: 528 nm) using a Cytation 5 Plate Reader (BioTek, Winooski, VT, USA). Percent dextran delivery was calculated based on the mass initially added to the donor compartment, correcting for the mass removed through sampling at each time point and considering 100% to be the receiving well concentration if equilibrium had been established between both compartments.

TEER values were recorded at the same time points and used to monitor monolayer disruption and recovery dynamics. For all experiments, monolayers with a TEER value of at least 560 Ohms⋅cm 2 were considered intact, and monolayers below this threshold were considered disrupted. "Time to monolayer recovery" was calculated as a summary metric to quantify the degree of monolayer disruption caused by sonication. For monolayers that recovered by the last measurement time point (72 hours), time to recovery was calculated as the time that TEER values remained below 560 Ohms⋅cm 2 , using a linear interpolation between measured time points. Monolayers that never dropped below 560 Ohms⋅cm 2 therefore had a time to monolayer recovery of 0 hours.

Detection of cavitation signals and their persistence during 30 second ultrasound exposures

Cavitation is thought to be the main mechanical driver of biological barrier permeabilization with ultrasound. As such, we were interested in characterizing the generation and persistence of stable and inertial cavitation energy generated by ultrasound-stimulated PCCAs during 30 second exposures.

We detected cavitation signals using a method adapted from our previous work [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF].

PCCAs (1.75×10 9 particles/mL) were suspended in 800 µl of 37°C PBS in a disposable polystyrene cuvette (One-way transmission factor at 1.0 MHz = 93±2% (n=10); thus, attenuation ≅ 7%) (FisherBrand, Thermo Fisher Scientific, Waltham, MA, USA). The unfocused 1.0 MHz transducer was shallowly submerged in the PCCA suspension and set to transmit the same ultrasound pulses that were used for transepithelial delivery (Table 1) with the only difference being that the pulse repetition frequency was reduced to 2000 Hz. Unfortunately, hardware limitations prevented us from saving data with a pulse repetition frequency of 5000 Hz to match the cell experiments. We do not believe that this discrepancy significantly alters the interpretation of our results, as discussed in subsequent sections.

A separate, single element, spherically focused transducer (7. We saved signals over a total of 30 seconds to monitor the dynamics of cavitation energy generated from the acoustically-stimulated PCCAs. Ultrasound was transmitted continuously throughout the 30second period; however, signals were saved intermittently to reduce data to a manageable size. Each second, 20 individual signals (i.e., 10 ms of data) were saved followed by a 990 ms saving delay. This was repeated every second for a total of 30 seconds.

All saved data was post-processed using an adaptation of our previously described post-processing analysis script in MATLAB (MathWorks Inc., Natick, MA, USA) to quantify stable and inertial cavitation doses [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF]. Briefly, for each individual radiofrequency line, the PCCA signal was selected by applying a 30 μs window from the point corresponding to the beginning of the acoustic pulse. Fast Fourier transform was used to convert individual time-domain signals to the frequency domain in order to estimate the level of stable and inertial cavitation. For stable cavitation, the area under curve (AUC) of the second harmonic component was calculated considering a spectral window from 1.9 to 2.1 MHz. For inertial cavitation, the broadband emission was quantified by calculating the AUC of frequency content ranging from 5.25 to 7.75 MHz (while excluding the 6 and 7 MHz harmonic components). Stable and inertial cavitation doses were then calculated by normalizing AUC values obtained for PCCAs with those calculated for a PBS-filled cuvette (without PCCAs) exposed to the same acoustic parameters. This was repeated for three independent vials of PCCAs; average cavitation doses (relative AUC) with inter-vial standard deviation are reported. Finally, average of relative AUC for stable and inertial cavitation are plotted over the 30 second period, and total cavitation dose for each was calculated by taking the area under these curves.

Statistical analyses

All statistical analyses were performed in GraphPad Prism 6 (GraphPad Software, Inc., La Jolla, CA, USA). P-values of <0.05 were considered statistically significant, and data are presented as average ± standard deviation throughout this work.

To evaluate which ultrasound parameters resulted in statistically significant dextran delivery, we Linear regressions were performed to assess the relationships between pressure and (1) total stable cavitation dose, (2) total inertial cavitation dose, and (3) the sum of total stable and inertial cavitation.

This was done separately for 20-cycle and 40-cycle conditions. Again, goodness-of-fit (r 2 ) and p-values are reported.

RESULTS AND DISCUSSION

PCCA characterization

The PCCA formulation was characterized by a polydisperse size distribution with an average concentration of (1.4±0.2) × 10 12 particles/mL and mean particle size of 170±20 nm (n=3 independent vials) (Figure 2A). Figure 2B shows representative optical microscopy images captured before (t = -2 ms)

and after (t = 8 ms) a single 20-cycle ultrasound pulse at t = 0 ms. We found that PCCAs consistently vaporized into microbubbles at and above a peak negative pressure of 300 kPa, but showed little to no activation at lower pressures (at 1.0 MHz), which is consistent with our previous findings [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF]). Qualitatively, we observed an increase in the number of bubbles generated with a single ultrasound pulse as pressure was increased from 300 kPa to 600 kPa (Figure 2A). From these experiments, we concluded that the pressure threshold for PCCA vaporization is between 100 -300 kPa at 1.0 MHz, and therefore pressures at or above 300 kPa were used for all subsequent experiments.

Ultrasound-mediated FITC-dextran delivery through epithelial monolayers

As expected, 70 kDa dextran was not able to penetrate sham-treated Caco-2 monolayers over the 72 hr incubation period, confirming that effective barrier function had been established in our model (Figure 3A). Similarly, we did not observe appreciable delivery through control monolayers treated with

(1) PCCAs only (no ultrasound) or (2) ultrasound only (no PCCAs, ultrasound stimulation at the highest setting of 600 kPa and 40 cycles) (Figure 3A). This indicated that neither our agents nor ultrasound alone altered monolayer integrity with the parameters used in our study. Correspondingly, TEER values for all control groups remained well above the previously chosen threshold for an intact monolayer (560 Ohm•cm 2 ) throughout the 72-hour observation period, further demonstrating that control treatments did not significantly alter membrane integrity (Figure 4A). Neither PCCAs alone nor ultrasound alone resulted in statistically significant changes in TEER values at any time point compared to the time-point respective sham controls (Table 2).

When ultrasound was combined with PCCAs we were able to achieve significant delivery (vs.

time-point respective sham controls) of both dextran sizes with select parameter combinations (Figure 3B and 3C). When substantial delivery was achieved, it most often reached statistically significant levels between 8-24 hours post-treatment with cumulative percent delivery increasing slightly thereafter. We note that we do not see appreciable dextran accumulation before the 8 hr time point in any case, and we speculate that the lag time may be due to the diffusional characteristics of the dextran molecules. Future experiments will look into the origin of this lag time with greater detail. In general, the dextran delivery results correspond with TEER values indicating significant monolayer disruption (Figure 4 and Table 2).

For all ultrasound-stimulated PCCA treatments, we observed an immediate and statistically significant reduction in TEER, with all groups characterized by average TEER values in the range of 200-280

Ohm⋅cm 2 approximately 5 minutes post-treatment (Figure 4 and Table 2). TEER values tended to recover within 24-48 hours post-treatment, regaining values ≥560 Ohm•cm 2 . The 300-kPa treatment groups were the only conditions that did not show 100% recovery of all monolayers by the 72-hour time point. For the 300 kPa treatment groups, we achieved 33% and 50% recovery for the 20 and 40 cycle cases, respectively.

Treatment groups diverged with respect to monolayer recovery time. To investigate this further, we interpolated "time to monolayer recovery" for each treatment (i.e., time for a monolayer to regain a TEER value ≥560 Ohm•cm 2 ). This summary value, which captures the monolayer recovery dynamics, correlated well with the cumulative percent dextran delivered (r = 0.8475, p < 0.0001) (Figure 5).

Biologically, the main factor dictating dextran delivery efficiency appears to be how long the monolayers remain leaky.

In terms of acoustic parameters, rarefactional pressure has a greater influence on dextran delivery outcome than pulse length. To test this observation statistically, we performed a two-way ANOVA to quantify the relative contributions of acoustic pressure and pulse length to overall variation in cumulative dextran delivery (72 hr time point). We found that pressure (p = 0.0002) but not pulse length (p = 0.1476) was a significant source of variability in percent dextran delivery, with pressure accounting for 53.5% of the total variation compared to 3.5% attributed to pulse length. This can be assessed visually in Figure 6, where the cumulative percent delivery varies substantially with respect to pressure but remains fairly consistent between pulse lengths. Indeed, multiple comparison testing showed that differences in cumulative dextran delivery between 20 and 40 cycle conditions were statistically insignificant at each rarefactional pressure.

The most interesting and unexpected finding was that the 300 kPa exposure conditions consistently provided significantly greater dextran delivery compared to higher-pressure groups, while the higherpressure conditions (400, 500, and 600 kPa) provided similar results (statistically insignificantly different) (Figure 6). These findings are counter to our initial hypothesis that as acoustic energy was increased, a greater proportion of PCCAs would vaporize into microbubbles and cavitation dose would be more pronounced, together leading to more dramatic biological effects and ultimately greater dextran delivery efficiency. We hypothesize that these unexpected findings may be explained by differences in cavitation persistence between ultrasound exposure conditions, as is explored in the subsequent section.

Persistence of cavitation over 30 second exposures

In our experiments, as would be the case for GI drug delivery applications, the same sample of PCCAs were stimulated for the duration of the ultrasound exposure time (30 seconds in this study).

Therefore, at high pressure, PCCAs may be rapidly converted to microbubbles but also subsequently destroyed via inertial cavitation. This would limit the effective cavitation duration to a fraction of the ultrasound on time (30 s). Conversely, if PCCAs were being vaporized more slowly and were able to persist for longer at lower pressures, the effects of microbubble cavitation would have been felt by the cells for a longer period of time, enhancing the biological effect.

To test this hypothesis, we quantified stable and inertial cavitation dose generated from acoustically-stimulated PCCAs in vitro over the course of 30 seconds. All parameters matched those of our dextran delivery experiments except for pulse repetition frequency, which was reduced to 2000 Hz.

For all acoustic exposure conditions (300-600 kPa with 20 and 40 cycle pulse lengths), we observed substantial stable and inertial cavitation immediately after the ultrasound transmission was initiated (Figure 7A-D). For 20 cycle pulse lengths, there is a clear separation in cavitation dynamics with respect to rarefactional pressures, with lower pressures resulting in more sustained stable and inertial cavitation compared to the higher pressures. The lower pressures (especially 300 kPa) provide sustained stable and inertial cavitation activity over the entire 30 second exposure, while cavitation dose decreases dramatically within the first 15 s for the higher pressures (especially 600 kPa) (Figure 7A & 7C). This trend is also observed for stable cavitation dose with respect to pressure at 40 cycles (Figure 7B). This is less pronounced for inertial cavitation dose at 40 cycles, as cavitation dose decreases more rapidly for all pressures (Figure 7D). In general, these data demonstrate that both peak rarefactional pressure and pulse length influence cavitation persistence for PCCAs, as has previously been reported for microbubbles [START_REF] Pouliopoulos | Rapid short-pulse sequences enhance the spatiotemporal uniformity of acoustically driven microbubble activity during flow conditions[END_REF]. The lower pulse repetition frequency used for these experiments compared to the dextran delivery experiments may have biased results slightly. At a higher pulse repetition frequency, cavitation dose would be expected to decrease more rapidly for all conditions, thus this effect is likely more pronounced in our cell experiments at higher pulse repetition frequency.

In general, we hypothesize that two main factors contribute to the inverse trend between total cavitation dose and pressure. First, we qualitatively observed that PCCA vaporization efficiency increases with increasing pressure. This can be observed in the representative images shown in Figure 2, where more bubbles are generated after a single ultrasound pulse at 600 kPa vs. 300 kPa. This qualitative finding is supported by previous work showing that relative PCCA vaporization efficiency increases with increasing pressure (Reznik, et al. 2013[START_REF] Wu | Focused ultrasound-facilitated brain drug delivery using optimized nanodroplets: vaporization efficiency dictates large molecular delivery[END_REF]). As a result, it will take longer to vaporize all PCCAs in the sample with lower pressures, contributing to sustained cavitation activity over the 30-second exposure time. Secondly, as pressure increases, we believe that generated microbubbles will be destroyed more rapidly, ultimately limiting microbubble lifetime and cavitation persistence. A similar inverse trend between generated microbubble stability/survival and pressure has previously been described by Reznik et al. (Reznik, et al. 2013).

To quantitatively analyze the trends in acoustic cavitation with respect to pressure, we first calculated the cumulative stable cavitation dose (SCD) and cumulative inertial cavitation dose (ICD) by calculating the area under each of the cavitation curves over the 30 s ultrasound exposure. This provided single values summarizing the total SCD and ICD generated by each condition. We also summed these cumulative SCD and ICD values to provide an estimate of the overall cavitation dose generated per condition (SCD+ICD). Linear regressions were performed between each of these summary values and acoustic pressure (Figure 7E-F). For the 20 cycle cases, we find statistically significant inverse trends between pressure and total SCD+ICD and total SCD (p = 0.0007 -0.0100). For the 40 cycle cases, all three summary cavitation metrics were found to significantly and inversely trend with pressure (p = 0.0049 -0.0100). For both pulse lengths, 300 kPa qualitatively provided the greatest cumulative SCD, ICD and SCD+ICD compared to all higher-pressure groups, which corresponds with our finding that 300 kPa results in the greatest dextran delivery efficiency.

Limitations

Total cavitation dose alone cannot predict dextran delivery efficiency as evidenced by the discrepancy we find between 20 and 40 cycle conditions: dextran delivery outcomes are comparable across pulse lengths, while total cavitation doses were in general lower for 40 cycle conditions compared to respective 20 cycle conditions. For the 20 cycle conditions, we find a significant correlation between cumulative dextran delivery and (1) total SCD (Pearson r = 0.9624, p = 0.0376) and (2) total SCD+ICD (Pearson r = 0.9638, p = 0.0362). However, for the 40 cycle cases we do not find significant correlations between dextran delivery and any of the total cavitation metrics. This discrepancy may be due to various confounding variables that were not captured by the total cavitation dose metric. For instance, we hypothesize that peak cavitation dose and the balance between stable and inertial cavitation may be important variables contributing to biological outcome.

Furthermore, acoustic pressure and pulse length are known to influence the stability and size distribution of microbubbles generated by PCCA vaporization (Reznik, et al. 2013[START_REF] Sheeran | Phase-transition thresholds and vaporization phenomena for ultrasound phase-change nanoemulsions assessed via high-speed optical microscopy[END_REF]), and differences in microbubble populations may have influenced dextran delivery efficiency. The proportion of small microbubbles generated by PCCA vaporization has previously been show to increase with increasing acoustic pressure [START_REF] Fix | An evaluation of the sonoporation potential of lowboiling point phase-change ultrasound contrast agents in vitro[END_REF][START_REF] Sheeran | Phase-transition thresholds and vaporization phenomena for ultrasound phase-change nanoemulsions assessed via high-speed optical microscopy[END_REF], which is attributed to the inverse relationship between vaporization pressure threshold and PCCA size (Sheeran andDayton 2012, Sheeran, et al. 2013). Additionally, the likelihood of microbubble fusion (i.e., generation of larger microbubbles) increases with increasing pulse length [START_REF] Sheeran | Phase-transition thresholds and vaporization phenomena for ultrasound phase-change nanoemulsions assessed via high-speed optical microscopy[END_REF]. Furthermore, differences in bubble aggregation and coalescence due to secondary Bjerknes forces may alter the number and size of bubbles between treatment groups. These expected differences in microbubble size distribution as a function of acoustic excitation parameters may influence monolayer disruption efficiency.

After ultrasound-mediated monolayer disruption, TEER value analysis indicated that it typically took 24-48 hours for the epithelial cells to regain their integrity. While recovery mechanism was not explicitly studied in this work, this rather long recovery time may indicate that recovery is due to cell repopulation and growth rather than resealing of transiently opened tight junctions between cells. A limitation of this study is that we did not quantify cell death after treatment. Furthermore, not all monolayers treated with the 300 kPa conditions recovered within the 72-hour observation window. We do not claim that we have found optimal acoustic parameters for transepithelial delivery in vivo. Rather, this study provides an encouraging proof of principle demonstration that transepithelial delivery is feasible with ultrasound stimulated PCCAs and supports further optimization of acoustic parameters to maximize transient disruption while minimizing recovery time and cell death.

It is important to note that cell-based assays of intestinal permeabilization are known to overestimate damage compared to the what would be observed in viable intestinal tissue (in the context of chemically induced permeability enhancement) [START_REF] Mccartney | Safety concerns over the use of intestinal permeation enhancers: A mini-review[END_REF][START_REF] Petersen | Colonic absorption of salmon calcitonin using tetradecyl maltoside (TDM) as a permeation enhancer[END_REF]. This is attributed to the lack of complete intestinal repair mechanisms and protective mucus in simple cell culture assays. As a result, any future optimization of acoustic parameters should be carefully validated (e.g., through the in situ intestinal perfusion [START_REF] Escribano | Single-pass intestinal perfusion to establish the intestinal permeability of model drugs in mouse[END_REF])) to ensure that the chosen parameters are effective and safe when working with viable intestinal tissue. One of the greatest safety concern for intestinal barrier disruption in vivo is the risk of facilitating absorption of harmful bystanders such as bacteria, viruses, and toxins [START_REF] Mccartney | Safety concerns over the use of intestinal permeation enhancers: A mini-review[END_REF]. This risk should be carefully evaluated and minimized upon in vivo translation.

Future directions -potential for ultrasound-mediated GI drug delivery with PCCAs

The data presented herein represent a first step towards the ultimate goal of ultrasound-mediated GI drug delivery with PCCAs, and the finding that lower pressure results in greater delivery efficiency is important for successful in vivo translation. In the in vivo setting, we envision co-administering a medicated enema with a solution of PCCAs. Subsequently, an endoscopic ultrasound probe would be inserted and used to stimulate PCCAs for GI permeabilization. Without natural replenishment of the PCCAs (e.g., through blood flow), achieving persistent cavitation by limiting the destruction of generated microbubbles will be of utmost importance. This is in contrast to an analogous application where ultrasound-stimulated contrast agents are used to enhance drug delivery through vascular barriers. For vascular disruption, microbubbles are administered intravenously, and blood flow provides continuous replenishment of intact microbubbles. Indeed, delivery efficiency has been shown to increase with increasing pressure for this application [START_REF] Chen | The size of blood-brain barrier opening induced by focused ultrasound is dictated by the acoustic pressure[END_REF][START_REF] Liu | Low-pressure pulsed focused ultrasound with microbubbles promotes an anticancer immunological response[END_REF][START_REF] Wang | Ultrasound-guided delivery of microRNA loaded nanoparticles into cancer[END_REF], and microbubble reperfusion into the focal zone can be achieved by alternating sonication with rest time [START_REF] Chowdhury | Ultrasound-guided therapeutic modulation of hepatocellular carcinoma using complementary microRNAs[END_REF][START_REF] Wang | Ultrasound-guided delivery of microRNA loaded nanoparticles into cancer[END_REF] or theoretically by using short pulses and long pulse repetition periods [START_REF] Shamout | Enhancement of noninvasive trans-membrane drug delivery using ultrasound and microbubbles during physiologically relevant flow[END_REF]. However, using high pressures for vascular permeabilization is associated with an increased risk of undesirable bioeffects [START_REF] Baseri | Multi-modality safety assessment of blood-brain barrier opening using focused ultrasound and definity microbubbles: a short-term study[END_REF], pointing to another motivation for using reduced pressures for any drug delivery application. To this end, Pouliopoulos and colleagues have demonstrated the value of redefining the acoustic pulses typically used for vascular permeabilization to provide sustained, controllable, and safe cavitation activity with low pressures and short pulse lengths [START_REF] Pouliopoulos | Exploiting flow to control the in vitro spatiotemporal distribution of microbubble-seeded acoustic cavitation activity in ultrasound therapy[END_REF][START_REF] Pouliopoulos | Rapid short-pulse sequences enhance the spatiotemporal uniformity of acoustically driven microbubble activity during flow conditions[END_REF].

The likelihood that ultrasound will induce cavitation is related to the mechanical index (MI = peak rarefactional pressure divided by center frequency) of the transmitted pulse. The United States Food and Drug Administration stipulates that MI cannot exceed 1.9 for clinical diagnostic ultrasound, which is intended to prevent cavitation-mediated biological effects in the absence of exogenous contrast agents.

Therefore, when ultrasound-mediated GI permeabilization is performed in the absence of contrast agents, MIs above this threshold (MI~2.2 -3.3) are often employed to generate the cavitation activity necessary for permeabilization (Schoellhammer, et al. 2017, Schoellhammer, et al. 2017[START_REF] Schoellhammer | Ultrasound-mediated gastrointestinal drug delivery[END_REF]. (Note: these MIs were estimated based on the reported acoustic intensities for ex vivo experiments, assuming plane wave transmission, an acoustic impedance of water of 1.48×10 6 kg⋅m -2 ⋅s -1 , negligible tissue attenuation, and that the reported intensity represented instantaneous acoustic intensity.) While these protocols have generally been described as safe in preclinical studies (Schoellhammer, et al. 2017[START_REF] Schoellhammer | Ultrasound-mediated gastrointestinal drug delivery[END_REF], cavitation-mediated biological effects could be induced in tissues beyond the GI epithelium. Conversely, we were able to achieve efficient drug delivery at a mechanical index of 0.3 (300 kPa at 1.0 MHz), which is well below the FDA limit of MI = 1.9, indicating that biological effects will likely only be induced in areas with direct contact to the contrast agents. This claim is supported by the absence of monolayer disruption (as evidenced by stably high TEER values and impermeability to dextran) when cells were stimulated with ultrasound alone in our experiments. This offers an opportunity for improved treatment localization and the potential for enhanced safety compared to ultrasoundmediated GI drug delivery in the absence of contrast agents.

In the GI space, we believe low boiling point PCCAs will provide greater success than standard microbubbles. Unlike microbubble contrast agents, the small size of PCCAs may allow permeation through the GI mucus mesh and permeabilization of the underlying GI epithelium. In future studies, we will test the penetration of PCCAs through mucus and optimize acoustic parameters for GI drug delivery in a model that more accurately captures the in vivo setting (e.g., in situ intestinal perfusion).

The use of ultrasound for drug delivery offers an opportunity for simultaneous image guidance. This will be particularly important for GI drug delivery applications where site-specific (e.g., immunotherapy delivery to a colon tumor) rather than global GI drug delivery is desired. Ultrasound image guidance could be achieved in one of two ways. First, it may be possible to achieve ultrasound imaging and subsequent drug delivery with clinically available endoscopic ultrasound probes (e.g. with 5 MHz center frequency).

A similar approach has been published for externally applied ultrasound for drug delivery to pancreatic tumors with a clinically available machine [START_REF] Dimcevski | A human clinical trial using ultrasound and microbubbles to enhance gemcitabine treatment of inoperable pancreatic cancer[END_REF][START_REF] Kotopoulis | Treatment of human pancreatic cancer using combined ultrasound, microbubbles, and gemcitabine: a clinical case study[END_REF]. Perhaps more exciting is the prospect of developing dual frequency endoscopic ultrasound probes that incorporate aligned low frequency elements for optimal drug delivery and high frequency elements for high-resolution ultrasound imaging. The prospect of high-frequency and therefore high-resolution image guidance may be particularly advantageous for small animal studies where anatomical features are difficult to resolve with relatively low-frequency imaging. Dual frequency transducers are currently under development for contrast-enhanced intravascular and intracavity imaging [START_REF] Kim | Phantom evaluation of stacked-type dual-frequency 1-3 composite transducers: A feasibility study on intracavitary acoustic angiography[END_REF][START_REF] Lindsey | Dual-Frequency Piezoelectric Endoscopic Transducer for Imaging Vascular Invasion in Pancreatic Cancer[END_REF][START_REF] Wang | Contrast Enhanced Superharmonic Imaging for Acoustic Angiography Using Reduced Form-factor Lateral Mode Transmitters for Intravascular and Intracavity Applications[END_REF]), which could be modified for image-guided therapy purposes. A dual frequency approach can also be employed to initiate PCCA phase change at high frequency and detect the unique acoustic signature of PCCA vaporization at low frequency (Arena, et al. 2015). This would offer a method to image the efficiency and duration of PCCA vaporization events during the treatment.

CONCLUSIONS

We have demonstrated efficient delivery of a macromolecular drug mimic (70 kDa FITC-dextran) through otherwise impermeable Caco-2 epithelial monolayers using ultrasound-stimulated PCCAs. We found that the lowest pressure conditions (300 kPa) consistently provided the greatest dextran delivery efficiency (vs. 400-600 kPa), which is explained in part by the observation that cavitation is more persistent during ultrasound exposure at lower pressures. While we believe that persistent cavitation activity is important for achieving efficient epithelial disruption, we have not fully explored other acoustic metrics that may also play a role, such as acoustic radiation force, peak cavitation dose, and variable size distributions of acoustically generated microbubbles. Further experimentation should be conducted to evaluate the relative importance of these variables and their contribution to dextran delivery outcome.

Insight gleaned from these experiments will allow for rational and thorough optimization of acoustic parameters for in vivo drug delivery through the GI epithelial barrier. (SCD) and inertial cavitation (ICD) for all acoustic conditions tested. The persistence of this cavitation over 30 seconds varied between conditions. In general, SCD and ICD was most persistent at the lowpressure conditions (300 or 400 kPa) with short pulse length (20 cycles). Total SCD and ICD dose was calculated as the area under the cavitation curves over 30 seconds. (E) For the 20 cycle conditions, we found significant inverse trend between rarefactional pressure and total SCD+ICD dose and SCD dose.

(F) This trend held true for SCD+ICD, SCD and ICD doses for the 40 cycle conditions. All data presented as mean±SD with n=3 replicates.

  performed a two-way repeated-measures ANOVA with respect to time and acoustic treatment. At each timepoint, the mean of each treatment group was compared to the mean of the sham control at the same time point via Dunnett's multiple comparisons testing. The same statistical procedure was applied to analyze changes in TEER values over time in various ultrasound treatment groups vs. time pointrespective sham controls.Cumulative dextran delivery achieved at 72 hours was analyzed via an additional two-way ANOVA with respect to acoustic pressure and pulse length. The purpose of this test was to evaluate the extent to which dextran delivery is affected by these factors (i.e. the percent of variability in dextran delivery that can be attributed to pressure vs. number of cycles). Following this two-way ANOVA, two multiple comparisons tests were performed. First, a Sidak multiple comparison test was used to compare mean delivery efficiency between 20 and 40 cycle cases at each pressure. Second, a Tukey multiple comparison test was used to compare mean delivery efficiency between all rarefactional pressures separately for the 20 and 40 cycle cases.The relationship between monolayer disruption and percent dextran delivery was assessed by calculating the correlation between time to monolayer recovery and cumulative percent dextran delivery at the 72-hour time point. Spearman correlation coefficients (r) and p-values are reported.

Figure 1 :

 1 Figure 1: (Color online) Sonication strategy. Colorectal adenocarcinoma (Caco-2) cells were cultured

Figure 2 :Figure 3 :Figure 4 :

 234 Figure 2: (A) PCCA characterization. The size distribution and concentration of PCCAs were

Figure 5 :

 5 Figure 5: Correlation between cumulative dextran delivery and time to monolayer recovery. A

Figure 6 :

 6 Figure 6: Influence of pulse length and rarefactional pressure on overall dextran delivery outcome.

  * **** **** **** **** ****

  5 MHz nominal frequency, measured center frequency = 8.740 MHz, diameter = 19 mm, focal distance = 50 mm, -6 dB bandwidth = 1.270

	MHz) (IL0706HP, Valpey Fisher Corp., Hopkinton, MA, USA) was arranged perpendicular to the
	transmit transducer with its focus positioned within the PCCA solution. Acoustic signals received by this
	transducer were passed through a receive amplifier (BR-640A, RITEC, Inc., Warwick, RI, USA) with 27
	dB gain and a 500 kHz high pass filter. Signals were subsequently digitized using a 12-bit analog-to-
	digital conversion card with a sampling rate of 200 MHz (GaGe model #CSE1222, DynamicSignals LLC,

Lockport, IL, USA) installed in a computer (Dell, Round Rock, TX, USA) running a custom acquisition program (LabVIEW, National Instruments Corp., Austin, TX, USA).

Table 1 :

 1 Summary of conditions tested. Summary of the acoustic conditions tested for ultrasoundmediated Caco-2 permeabilization with PCCAs, including the number of experimental replicates performed.

	Abbreviations: PCCA -phase change contrast agent; kPa -kilopascal; kDa -kilodalton; US -
	ultrasound; Y -yes; N -no; NA -not applicable

Table 2 : Statistics describing transepithelial resistance (TEER) values before and after monolayer sonication.

 2 Stars indicate the results of Dunnett's multiple comparison testing at each timepoint comparing mean TEER value of each group to that of the sham control for that timepoint after two-way ANOVA. Key: * p ≤ 0.05; ** p ≤ 0.01; *** p ≤0.001; **** p ≤ 0.0001 Abbreviations: PCCA -phase change contrast agent; kPa -kilopascal; US -ultrasound; cyc -cycles.

		Significantly different than respective sham
		control?
		Time point (hours)
	Group	pre 0	3	8	24	48	72
	US only	ns	ns	ns	ns	ns	ns	ns

ACKNOWLEDGEMENTS

The authors would like to thank Eric Markley for assistance adapting the LabVIEW code used in this study, as well as Catherine Feng and Dr. Nicholas Shaheen for useful discussions on the gastrointestinal drug delivery background relating to this work. Financial support was provided by an NRSA Individual Predoctoral Fellowship (F31CA220970, S.M.F) and the University of North Carolina General Assembly Research Opportunities Initiative: Pharmacoengineering (V.P & B.P.K). P.A.D declares that he is a co-inventor on a patent describing the formulation of low-boiling point perfluorocarbon agents and a co-founder of Triangle Biotechnology, a company that has licensed this patent. All the other authors declare that they have no competing interests.


Influence of Water Content on Spectral Reflectance

of Leaves in the 3-15-µm Domain Sophie Fabre, Audrey Lesaignoux, Albert Olioso, and Xavier Briottet

Abstract-This letter describes a laboratory experiment where reflectance signatures of three plant species are measured at a leaf level in the 3-15-µm spectral domain. The leaf samples are progressively dried in order to analyze the behavior of their spectral signature according to the variations in their water content. Our first objective aims to underline leaf water content (LWC) impact on the spectral signatures. This work is a necessary step toward further studies dealing with interpretation of multispectral remote sensing data or estimation of water stress and energy budget. The drying process and measurement method are detailed. This letter deals with dry and fresh leaves (as found in literature) and considers intermediate water content levels as well. For intermediate LWC levels, our analysis outlines some important results: The spectral domain may be divided into two parts, namely, 3-5.5 and 5.5-15 µm, each corresponding to different impacts of LWC variation; both sides of a cherry tree leaf have not the same behavior according to the water content amount; and, in the 8-15-µm, the drying process impacts when the LWC becomes lower than a threshold value (around 30%). Index Terms-Leaf, optical properties, spectral reflectance, thermal infrared domain, water content.

I. INTRODUCTION

T HE IMPACT of leaf water content (LWC) on leaf re- flectance has been widely documented in the 0.4-2.5-µm solar spectrum for a long time [START_REF] Gausman | Leaf optical properties in visible and near-infrared light[END_REF]. It has been shown that the shortwave infrared (1.4-2.5 µm) is strongly influenced by the water in plant tissue, particularly at 1.45 and 1.94 µm. Radiometric measurements in this domain have been suggested for estimating vegetation water content [START_REF] Zarco-Tejada | Water content estimation in vegetation with MODIS reflectance data and model inversion methods[END_REF].

Conversely, little information is available for larger wavelengths situated in the 3-20 µm range. However, this information can be helpful for remote sensing interpretation or land surface and atmospheric studies (for example, in the 7-15-µm range [START_REF] Olioso | Evidences of low land surface thermal infrared emissivity in presence of dry vegetation[END_REF] or around 4.7 µm [START_REF] Ho | Estimates of 4.7 µm surface emissivity and their impact on the retrieval of tropospheric carbon monoxide by measurements of pollutions in the troposphere (MOPITT)[END_REF]). Only few measurements of leaf spectral reflectance and emissivity were presented in the past S. Fabre and X. Briottet are with ONERA, 31055 Toulouse, France (e-mail: Sophie.Fabre@onera.fr; Xavier.Briottet@onera.fr).

A. Lesaignoux is with ONERA, 31055 Toulouse, France, and also with the Institut Supérieur de l'Aéronautique et de l'Espace, Université de Toulouse, 31055 Toulouse, France (e-mail: Audrey.Lesaignoux@onera.fr).

A. Olioso is with the Environnement Méditerranéen et Modélisation des Agro-Hydrosystèmes, Unité Mixte de Recherche 1114, Institut National de la Recherche Agronomique, 84914 Avignon, France (e-mail: olioso@avignon. inra.fr). Fig. 1. Specific absorption coefficient of water. Redrawn from [START_REF] Wieliczka | Wedge shaped cell for highly absorbent liquids: Infrared optical constants of water[END_REF].

[5]- [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF]. Moreover, these measurements considered fresh or fully dried leaves only. Intermediate moisture levels (that means variation of leaf reflectance as a function of water content) were not considered. These measurements showed that the differences between wet and dry leaf reflectances in the 7-15-µm spectral range are limited in comparison to those in the 1.4-2.5-µm range. Similar differences between wet and dry leaf reflectances are expected in the 3-6-µm window [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF]. These low differences may be explained by the high absorption characteristics of the other leaf materials which limit the effect of water. As shown in Fig. 1, the absorption coefficient of water is usually very high, which is above 2.5 µm with fundamental molecular vibrations of water located at 2.86, 3.05, and 6.08 µm [START_REF] Wieliczka | Wedge shaped cell for highly absorbent liquids: Infrared optical constants of water[END_REF].

Furthermore, active mid-infrared backscatter characteristics of various agricultural crop leaves were measured at different wavelengths in the 9-11-µm spectral domain. Both copolarized and cross-polarized reflectance dependences on leaf moisture were analyzed by Narayanan et al. [START_REF] Narayanan | Mid infrared laser reflectance of crop leaves subjected to water stress[END_REF]. They pointed out the interest of using this range to study the impact of leaf moisture. However, these works were based on directional active measurements, and their results cannot be transposed to the analysis of remote sensing data which are currently limited to passive acquisitions and which are not obtained in the same geometrical conditions.

In order to analyze the impact of LWC on spectral reflectance, we conducted a laboratory experiment. Spectral reflectances were measured over the entire thermal infrared spectral domain between 3 and 15 µm for leaves which were progressively dried. The analysis of these laboratory measurements is provided for the middle infrared (3-5 µm) and the long-wave infrared [(LWIR); 5.5-15 µm] domains.

First, a laboratory measurement method and the acquired spectral reflectance at several water contents are described in Section II. Results are then analyzed in Section III. Discussions will follow in Section IV. Finally, conclusions and perspectives are given in Section V.

II. METHOD AND MEASUREMENTS

A. Overview of the Measurement Sequence and the Drying Process

The effect of water content on leaf optical properties was investigated for three plant species: one monocotyledon, namely, sorghum (Sorghum bicolor), and two dicotyledons, namely, sunflower (Helianthus annuus) and cherry tree (Prunus avium). A preliminary visual analysis showed that their leaves differed by their thickness, color, surface composition (waxes of abaxial leaves of Prunus was more visually apparent), and roughness. All samples were collected in agricultural fields near Toulouse (France) during summer 2006. Branch and plants selected for their good health and physiologic conditions were collected and stored in a fresh environment for an efficient conservation. Because of instrumental constraints, a part of leaf was sampled for each species before proceeding to the spectral measurements.

In order to measure the reflectance spectra at different moisture levels, successive LWCs were obtained by progressively drying the leaf samples in a laboratory oven (Sapratin brand). In order to avoid leaf rolling, the oven temperature was set to 60 • C, the drying time was limited to 30 min, and a rigid perforated sheet was put on the sample during the drying (experimental conditions inferred by preliminary tests). Thus, spectral reflectances of each leaf sample were measured at a minimum of five LWC levels. In order to obtain fully dried leaves, the samples were put in the oven for a period of 24 h between the two last spectral measurements. At each drying step, the reflectance spectrum of the leaf sample was measured using a spectrometer operating in the 3-15-µm range. In order to control the leaf moisture variations during the spectral measurements and to provide an accurate moisture assessment, the sample was weighed before and after each acquisition of reflectance data. The spectra were always measured using the same leaf on both sides for each species, as well as the same location on the leaf. The spatial variability within the leaf was not investigated. The dissimilar spectral behaviors are presented in the next paragraph.

B. Measurement of LWC

Leaf weights were measured using a laboratory scale with 100-µg accuracy. The LWC (in percent) was deduced from dry (m D ) and wet (m W ) weights according to the following relation: LW C = 100 (m W -m D )/m W . The weight m D of the fully dried sample was obtained after 24 h of drying (assuming that residual LWC was 2% [START_REF] Peters | On-farm moisture testing of corn silage[END_REF]).

C. Measurement of Spectral Reflectances

In the 3-15-µm domain, spectral reflectances were directly measured with a Fourier transform infrared instrument: Bruker Equinox 55 spectrometer (http://www.brukeroptics.com/). Owing to several mercury-cadmium-telluride detectors, which were cooled with liquid nitrogen (77 K) in order to reduce thermal noise, it made it possible to scan with a waveband of 7000 to 370 cm -1 , i.e., from 1.4 to 27 µm, with a standard spectral resolution better than 0.5 cm -1 . One measure resulted from the average of 256 acquisitions (scans) done with an 8-cm -1 resolution (providing a noise level less than 2% [START_REF] Tanguy | Intercomparaison, évaluation et amélioration de l'étalonnage d'instruments de mesure de réflectance spectrale en infrarouge[END_REF]). In order to keep a good signal-to-noise ratio, only the 3-15-µm range was used. The instrument was set using an Infragold integration sphere, with a 127-mm diameter, and measured directional-hemispherical reflectances. The sphere was purged to remove carbon dioxide and water vapor. The angle of the incident beam was 13 • . Under this sphere, the sample size had a minimum diameter of 36 mm. Due to the long duration of the measurement sequence, we limited our measures to only one sample per species. The reflectance was obtained by comparing the measurement to an Infragold diffuse reference (annually certified by Labsphere).

III. RESULTS:MEASURED SPECTRA AT SEVERAL WATER CONTENTS

A. Fresh Leaf Samples Fig. 2 shows the spectral reflectances acquired in the 3-15-µm domain for different LWCs. For fresh leaves, the spectral reflectance levels were low and ranged between 0.02 and 0.07. The sunflower leaf sample exhibited weaker spectral reflectance levels than those related to cherry tree and sorghum leaf samples. The behaviors of both leaf sides were dissimilar for the cherry tree leaf sample (see spectral ranges 3-5 and 8-10 µminFig.2).

In the middle infrared domain, i.e., 3-5.5 µm, our measurements were similar to spectra of typical green foliage measured by Salisbury and D'Aria [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], with low reflectance due to intense leaf water absorption [START_REF] Elvidge | Thermal infrared reflectance of dry plant materials: 2.5-20 µm[END_REF] (Fig. 1). In this spectral domain, leaf samples had similar global behaviors with two apparent maxima around 4 and 5.2 µm, due to cellulose absorption for dry leaves and water absorption at 4.65 µm for fresh leaves (Fig. 1), and the Reststrahlen bands near 3.43 and 3.51 µm due to the H-C vibration bands [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF]. Furthermore, for all the spectra, the strong absorption at 3.05 µm linked to water (Fig. 1) is well noticed.

In the 5.5-15-µm range, the absorption peak linked to water at 6 µm (Fig. 1) is observed for the abaxial cherry, adaxial sunflower, and adaxial sorghum samples.

In the entire spectral domain 3-15 µm, a more detailed analysis showed that the spectral features varied from one species to another. Furthermore, in the case of the cherry leaf sample, both sides of the leaf had dissimilar spectral signatures: In the 3-5.5-µm range, the reflectance level difference between both sides was around 0.02, and double reflectance peaks became less apparent for the adaxial side of the leaf; in the 8-15-µm range, the reflectance peak around 8.7 µmw a sl e s s pronounced than for the abaxial side of the leaf [Fig. 2(a) → 1]. This was in agreement with previously published data [START_REF] Da | Spectral reflectance and emissivity features of broad leaf plants: Prospects for remote sensing in the thermal infrared (8.0-14.0 µm)[END_REF]- [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF] and with the data available in the MODIS (http://www. icess.ucsb.edu/modis/EMIS/html/em.html) spectral library. Indeed, the spectral signatures of abaxial and adaxial side leaves of green cherry (Prunus serotina) presented in [START_REF] Salisbury | Thermal infrared (2.5 to 13.5 µm) directional and hemispherical of leaves[END_REF] and those of Prunus serotina measured by Ribeiro da Luz [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF] were in good agreement with our measurements: the maximum of reflectance at 5.94 µm [Fig. 2 The sunflower leaf sample had similar spectral signature as the one of Red oak leaf in [START_REF] Salisbury | Emissivity of terrestrial materials in the 8-14 µm atmospheric window[END_REF] and was close to those of Heliantus annuus and Cucurbita pepo in [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF]. For the Heliantus annuus variety, the only difference was related to the absorption peak near 9 µm [Fig. 2(c

) → 1].
The sorghum spectral signature looked like Indiangrass blades in [START_REF] Salisbury | Emissivity of terrestrial materials in the 8-14 µm atmospheric window[END_REF] and Zea mays in [START_REF] Da | Spectral reflectance and emissivity features of broad leaf plants: Prospects for remote sensing in the thermal infrared (8.0-14.0 µm)[END_REF], [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF]: absorption peaks near 6 µm [Fig. 2(d

B. Drying Leaf Samples

For fully dried leaves (Fig. 2; a residual LWC of 2%), the spectral reflectance levels increased and ranged between 0.02 and 0.18. The spectral reflectance signature of cherry tree sample was in good agreement with dried cherry leaf in [START_REF] Salisbury | Thermal infrared (2.5 to 13.5 µm) directional and hemispherical of leaves[END_REF]. Dried sorghum and dried sunflower spectra were similar to senescent foliage of rye grass in [START_REF] Salisbury | Emissivity of terrestrial materials in the 8-14 µm atmospheric window[END_REF] and [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF] with weaker reflectance levels, particularly in the 3-5-µm domain. Furthermore, our measurements tallied with the spectra of dry plant materials in [START_REF] Gausman | Leaf optical properties in visible and near-infrared light[END_REF]. For example, the signature due to the mixture of holocellulose and lignin was well observed with a maximum around 4-5.5 µm.

In complement to the spectra in Fig. 2, the impact of leaf drying is shown for two specific wavelengths in Fig. 3. For our data set, the water had more impact in the 3-5 µm domain. Indeed, the spectral reflectance increased when the water content decreased, and the levels between fresh and dry spectra differed from 0.1 to 0.14 [Figs. 2 and3(a)]. Furthermore, the reflectance levels were more influenced when LWC varied below 30% to 40% (in relation to the high absorption by water). In Fig. 3(a), the reflectance at 4.1 µm increased with LWC, as an exponential function. There were few spectral variations in the 5.5-10-µm spectral domain due to the strong absorption of cellulose.

The impact of water content increased beyond 10 µm. This increase is underlined in particular for low LWC (LWC < 20%), with variations of reflectances between 0.03 and 0.01 at 11 µm [Fig. 3(b)]. This behavior was not observed for the abaxial side of the cherry tree, and the drying impact was different on both sides of the cherry tree leaf. This might be due to dissimilar leaf compounds and textures on the two sides of the leaf.

IV. DISCUSSION

Works conducted in [START_REF] Elvidge | Thermal infrared reflectance of dry plant materials: 2.5-20 µm[END_REF], [START_REF] Da | Spectral reflectance and emissivity features of broad leaf plants: Prospects for remote sensing in the thermal infrared (8.0-14.0 µm)[END_REF]- [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], and [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF] provide information for the interpretation of fresh and dry leaf spectra and allow one to validate our measurements by comparison. The leaf spectra behavior might be then defined, owing to its specific components (see, for example, the spectral impact of wax on cherry tree sample or of hemicellulose on sorghum sample). For fresh and dry conditions, some differences in terms of spectral reflectance level appeared in comparison to [START_REF] Salisbury | Thermal infrared (2.5 to 13.5 µm) directional and hemispherical of leaves[END_REF]- [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], [START_REF] Da | Spectral reflectance and emissivity features of broad leaf plants: Prospects for remote sensing in the thermal infrared (8.0-14.0 µm)[END_REF], and [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF] data. For example, sorghum spectra had similar levels as those of Indian grass [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], Zea Mays [START_REF] Da | Propriedades espectrais das plantas no infravermelho termal (2,5-14 µm)[END_REF], and Rye grass [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF] in the 8-15-µm range, but lower levels at shorter wavelengths. In the case of sunflower, the reflectance levels were always different from those found in the literature. These differences might be explained by variations in the chemical composition of leaves, particularly in the surface layers, and/or by variations in the structure of the leaves. Such changes might be related to differences in plant cultivar, as well as differences in leaf age and growth conditions.

To our knowledge, our study is the first one to report the impact of leaf drying on reflectance in the 3-15-µm range. A previous study at canopy level by Olioso et al. [START_REF] Olioso | Evidences of low land surface thermal infrared emissivity in presence of dry vegetation[END_REF] suggested the impact of such drying but no direct evidence was given at leaf level. This letter showed that the drying impact may be different according to species, spectral domain, and LWC level. The first evidence in our results was that the investigated domain could be split into two parts: one above and the other below 5.5 µm. All samples had the same spectral behavior in the 3-5.5-µm range, whereas it changed according to species in the 8-15-µm range. For all leaf samples, the reflectance at a given wavelength in 3-5.5 µm had an exponential behavior according to water content decrease. A quite similar behavior was noticed for most of the samples in the 8-15 µm, except for the adaxial side cherry tree leaf which was almost unaffected by drying. In this last range, the observed changes in reflectance were limited and occurred mostly for LWC lower than 30%. At higher LWC, reflectances were almost unaffected, while in the 3-5.5-µm range, reflectance changes were continuous.

Only measurements of leaf directional-hemispherical reflectances were presented in this letter, and no attempt was done for deriving leaf emissivity. If the validity of the Kirchhoff law is not challenged, the derivation of leaf spectral emissivities would require the measurement of leaf directionalhemispherical transmittance in order to derive the absorptance of directional radiation. Indeed, the works in [START_REF] Gerber | Modeling infrared optical properties of leaves to improve water content estimation[END_REF] and [START_REF] Olioso | Adaptation of the leaf optical property model PROSPECT to thermal infrared[END_REF] highlighted that drying leaves might not be opaque in the thermal infrared domain (in particular, in the spectral domains characterized by high leaf reflectances).

Our study was the first step for considering the impact of vegetation water content on spectral properties in the midand thermal-infrared domains. It only considered the leaf level and could not be used per se for the direct interpretation of land surface spectral properties and satellite signal. Many other factors would have to be accounted for, such as sensor performances, geometry of acquisition, vegetation architecture, soil optical properties, atmospheric properties, and the relative contributions of reflected radiations (solar and atmospheric) and of emitted radiations according to wavelength. The assessment of all these factors would require the integration of various specific models, including atmospheric radiative transfer models, sensor models, and vegetation canopy radiative transfer models. These last ones have to be used for analyzing the impact of geometry of acquisition, vegetation architecture (Leaf Area Index, leaf angle distribution, etc.), and soil spectral properties on land surface spectral emissivities or spectral radiances [START_REF] Salisbury | Emissivity of terrestrial materials in the 3-5 µm atmospheric window[END_REF], [START_REF] Tanguy | Intercomparaison, évaluation et amélioration de l'étalonnage d'instruments de mesure de réflectance spectrale en infrarouge[END_REF], [START_REF] Olioso | Adaptation of the leaf optical property model PROSPECT to thermal infrared[END_REF], [START_REF] Lesaignoux | Soil moisture impact on lab measured reflectance of bare soils in the optical domain 0.4-15 µm[END_REF]. It is important to notice that coherent modeling of directional radiances requires solving leaf and soil energy balances in order to account for the distribution of surface temperatures inside the canopy and at the soil surface, which also require the knowledge of leaf optical properties over the whole spectral range. Another point to consider in the future concerns the differences between drying conditions in the laboratory and in vivo. This drying process difference might affect the relationship between LWC and spectral properties.

V. C ONCLUSION

We have presented the laboratory measurements of spectral directional-hemispherical reflectances acquired on leaf samples in the 3-15-µm domain. Three species have been investigated. Our measurements were in good agreement with existing data published in literature for fresh and fully dried leaves (even if the reflectance levels differ). The originality of our work consists not only in measuring spectral signatures for fresh and dry samples but also for intermediate LWCs. Our measurement analysis has led to the following significant points.

1) The investigated spectral domain might be divided into two parts: 3-5.5 and 5.5-15 µm. In the first domain (mid-wave infrared), the reflectance was very sensitive to LWC, whereas in the LWIR domain, this dependence remains very weak for wavelengths higher than 10 µm.

2) The considered species and leaf side had a different behavior according to the water content amount. In future laboratory experiments, these measurements will be completed by increasing the set of leaf samples in order to validate our first results. Moreover, it will be interesting to measure both reflectance and transmittance in order to describe the complete optical property behavior, including emissivity, according to water content amount. The experiment measurement geometry and observation conditions were not representative of those related to remote sensing observations. Further analyses are in progress for including the impact of other factors such as atmospheric conditions, canopy architecture [START_REF] Olioso | Test of the SAIL-Thermique radiative transfer model for simulating thermal infrared emissivity and emissivity spectra of plant canopies[END_REF], or soil water content [START_REF] Lesaignoux | Soil moisture impact on lab measured reflectance of bare soils in the optical domain 0.4-15 µm[END_REF].

Fig. 2 .

 2 Fig. 2. Spectral reflectances of leaves for several water contents in the thermal domain. (a) Adaxial and (b) abaxial sides of cherry tree, (c) adaxial side of sunflower, and (d) sorghum leaf samples. Arrow labels locate a specific behavior, detailed in the document.

  (b) → 1] and at 9.67 µm [Fig. 2(a) → 2] due to the absorption of triterpenes existing in resin, absorption peaks at 5.75 µm [Fig. 2(a) → 3] due to the C-O double bond, and between 8.7-and 9.7-µm [Fig. 2(a) → 1] reflectance peaks due to leaf compounds other than cellulose. In fact, the very subtle reflectance maxima of cellulose in this range are normally accompanied by high transparency between 10.5 and 12 µm [Fig. 2(b) → 2] and by a volume scattering feature at 11.15 µm [Fig. 2(b) → 3] [5], [7].

  ) → 1] due to water [Fig. 1] and 9.5 µm [Fig. 2(d) → 2], the maximum of reflectance at 9.62 µm [Fig. 2(d) → 3] ap r i o r irelated to COH bending, and the presence of hemicellulose.

Fig. 3 .

 3 Fig. 3. Leaf reflectances at (a) 4.1 and (b) 11 µm as a function of LWC (wavelengths chosen in atmospheric window only).

ACKNOWLEDGMENT

The authors would like to thank S. De Gang, H. Mourin (ERM Belgium), and B. Tanguy (ONERA) for the help during the experimental instrumentation manipulation.