This protocol presents a method to perform rheology characterization of mucus that resides on gill rakers (GRs) of the silver carp. Viscoelastic characteristics of GR-mucus, obtained by measuring viscosity, storage and loss moduli, are evaluated for the apparent yield stress to understand the filter feeding mechanism in GRs.
The silver carp, Hypophthalmichthys molitrix, is an invasive planktivorous filter feeder fish that infested the natural waterways of the upper Mississippi River basin due to its highly efficient filter feeding mechanism. The characteristic organs called gill rakers (GRs), found in many such filter feeders, facilitate the efficient filtration of food particles such as phytoplankton that are of a few microns in size.
The motivation to investigate the rheology of the GR mucus stems from our desire to understand its role in aiding the filter feeding process in the silver carp. The mucus-rich fluid, in a ‘thick and sticky’ state may facilitate the adhesion of food particulates. The permeation and transport through the GR membrane are facilitated by the action of external shear forces that induce varying shear strain rates. Therefore, mucus rheology can provide a vital clue to the tremendous outcompeting nature of the silver carp within the pool of filter feeding fish. Based on this it was posited that GR mucus may provide an adhesive function to food particles and act as a transport vehicle to assist in the filter feeding process.
The main objective of the protocol is to determine the yield stress of the mucus, attributed to the minimum shear stress required to initiate flow at which irreversible plastic deformation is first observed across a structured viscoelastic material. Accordingly, rheological properties of the GR mucus, i.e., viscosity, storage, and loss moduli, were investigated for its non-Newtonian, shear-thinning nature using a rotational rheometer.
A protocol presented here is employed to analyze the rheological properties of mucus extracted from the gill rakers of a silver carp, fished at Hart Creek location of the Missouri River. The protocol aims to develop an effective strategy for rheological testing and material characterization of mucus assumed to be a structured viscoelastic material.
The silver carp, Hypophthalmichthys molitrix, is a planktivorous filter feeder and an invasive species that has infiltrated several natural waterways in the United States. This species was initially introduced in the upper Mississippi River basin to control algal blooms1,2,3. The silver carp is an extremely efficient feeder. Typically, its consumable food particle sizes range from 4 to 20 μm to larger zooplankton that are around 80 μm3,4,5. This species has outcompeted other native fish and can potentially cause enormous damage to native waterways by limiting available resources1,2,6. Thus, filter feeding fish such as the silver carp and the bighead carp pose a major threat to the Great Lakes1,2,6,7,8.
Filter feeding fish possess special organs called the gill rakers (GRs) with a thin layer of mucus residing on their surface. These organs improve the efficiency of filtration and aggregation of small particles from the incoming fluid. The goal of the protocol presented herein is to characterize the non-Newtonian, shear thinning material property and yield stress of the GR mucus acquired from the inner surface of the gill rakers in the silver carp. The value of yield stress of the GR-mucus, ascertained using a rotational rheometer, is of interest in this study. The measured yield stress also called the “apparent yield stress” depends on the testing methods such as steady shear rate- or dynamic oscillatory strain-type9,10. The shear-thinning, ‘yield-stress fluid,’ undergoes a transition from solid-like to liquid-like behavior at a critical applied stress9,11. The apparent yield stress is the minimum shear stress required to initiate flow or that at which irreversible plastic deformation is first observed when the mucus transitions from a gel-like material to a fluid-like material. This behavior can be observed in structured viscoelastic materials. The transition from gel-like to fluid-like behavior of the GR mucus entails two functions i.e., an adhesive role to gather food particulates and a transport vehicle role to assist in the particulate delivery and filtration process. The extended function of the mucus includes creating diffusion barriers in disease resistance and respiration, providing controlled release of nutritional factors, toxic components and excretion, creating metabolic pathways for feeding and nesting, helping in predator protection, and producing boundary layer modifications that improve the locomotion and propulsive efficiency12,13,14.
Unlike simple fluids, complex fluids like the mucus possess properties that vary with flow conditions and require additional measurement parameters to define their bulk scale physical behavior. To monitor the viscosity and yield stress of GR mucus, rheological measurements are performed using a rotational rheometer. The rotational rheometer applies a steady or oscillatory shear stress or strain by means of a rotating disk in contact with the fluid sample and measures its response. The rationale behind using this instrument and technique is that the rheometer can provide a set of measurements to describe the material properties of the GR mucus of the silver carp, which cannot be defined by viscosity alone.
The mucus is a viscoelastic material and its mechanical response to an imposed deformation is between that of a pure solid (governed by Hooke’s law of elasticity) and that of a pure liquid (governed by Newton’s law of viscosity)15,16. The complex macromolecular network contained within the mucus can stretch and reorient in response to external forces or deformation. A rotational rheometer is comprised of a cone geometry and a Peltier plate as shown in Figure 1 and Figure 2 (see Table 1 for instrumentation specifications). The objective of this study was to develop a protocol to determine the rheological properties of the GR mucus. An advantage of the rotational rheometer over a viscometer is its ability to make dynamic measurements using small sample volumes. The GR mucus sample volume in this study was approximately 1.4 mL. The viscometer, on the other hand, is limited to constant shear rates and requires large sample volumes.
The rheological properties of the mucus are expected to vary greatly within the silver carp anatomy. For example, the properties of the mucus residing on the GR surfaces may be different from the epibranchial organ. To account for the potential variability of mucus properties in different regions of the fish, the acquired GR mucus sample was diluted, and solutions of three concentrations were created and tested using the rotational rheometer.. The data and results regarding mucus rheology reported after executing the protocol demonstrated the efficacy of the measurement technique. The illustrative data presented in this paper are not meant to be generalized across the entire silver carp population. The protocol presented herein can be extended to investigate mucus rheology across larger sample sets to test other hypotheses.
The purpose of this study is to demonstrate the variation of rheological properties of GR mucus rheology with three different mucus concentrations (400 mg/mL, 200 mg/mL and 100 mg/mL). The 400 mg/mL concentration represents the raw mucus sample harvested from the fish GRs. Deionized water (DI) was used to dilute the raw mucus sample into 200 mg/mL and 100 mg/mL concentrations. Diluting the mucus samples allowed for the evaluation of the degree of shear thinning and apparent yield stress as a function of concentration and the determination of the concentration at which the GR mucus transitions to non-Newtonian behavior. A shaker was used to break down any large clumps of mucus in the samples to mitigate errors in the rheological data due to inhomogeneity.
In most vertebrates, including fish, the predominant mucus-forming macromolecules are glycoproteins (mucins) that tend to swell in water by entanglements or chemical cross-linking and create a gel-like material12,13,17,18,19,20. The high-molecular-weight, gel-forming macromolecules and high-water content reflects the slipperiness in the mucus13. A high degree of inter-macromolecular interactions leads to gel-formation whereas lower levels of inter-macromolecular interactions or broken bonds result in high-viscosity fluids21.
The processes of food particulate filtration in filter feeding fish are aided by GR mucus-related properties such as cohesion and viscosity that determine its potential for adhesion and tack22. The strength of mucus-based adhesion depends on specific intermolecular, electrostatic or hydrophobic interactions23. Sanderson et al.24 conducted a suspension-feeding study in blackfish wherein they found the evidence for mucus-based adhesion. They stated that the adhesion of suspended food particulates with a mucosal surface is followed by the transport of aggregated clumps of particles bound together with mucus by directed water-flow acting on it24. The mucus exposed to shear strain rates generated from water-flow facilitates the delivery of food particulates to digestive organs. Endoscopic techniques were used to observe filtered particles24.
Literature on the range of shear rates and practical limits in the rheological testing of GR mucus is scarce. Therefore, guidance was sought from rheological studies on gastric, nasal, cervical and lung mucus, salmon skin mucus, hagfish slime, and bone-joint surface lubricant wherein the rheological characterization and non-Newtonian attributes were studied11,12,25,26,27,28,29,30,31. More recently, the effect of fish skin mucus on locomotion and propulsive efficiency has been studied using constant shear rate viscometry. Skin mucus rheology studies (without any dilution or homogenization) pertaining to seabream, sea bass and meagre demonstrated non-Newtonian behavior at typically low shear rates14. In another related study, the raw skin mucus samples from dorsal and ventral sides of the Senegalese sole were found to exhibit non-Newtonian behavior, indicating a higher viscosity of the ventral mucus at all shear rates considered32. Other rheological protocols pertaining to the hydrogel scaffold development and for highly concentrated suspensions using a constant shear rate viscometer have also been reported in the literature33,34.
In this study, the GR mucus properties were investigated using a strain rate controlled, rotational rheometer that has been widely used in rheology experiments on complex biological fluids25. For Newtonian fluids, the apparent viscosity remains constant, is shear-rate-independent and the shear stresses vary linearly with shear strain rates (Figure 3A, B). For non-Newtonian fluids (such as shear-thinning fluids) viscosity is shear-rate-dependent or deformation-history-dependent (Figure 3A, B). The loss modulus (G”) represents the extent to which the material resists the tendency to flow and is representative of fluid viscosity (Figure 4). The storage modulus (G’) represents the tendency of the material to recover its original shape following stress-induced deformation and is equivalent to elasticity (Figure 4). The phase angle (δ) or loss tangent value, is calculated from the inverse tangent of G”/G’. It represents the balance between energy loss and storage and is also a common parameter for characterizing viscoelastic materials (δ = 0° for a Hookean solid; δ = 90° for a viscous liquid; δ < 45° for a viscoelastic solid and δ > 45° for a viscoelastic liquid) (Figure 4)25. The apparent yield stress (σy) in structured fluids represents a change of state that can be observed in rheological data from steady state sweep and dynamic stress-strain sweeps10. If the external applied stress is less than the apparent yield stress, the material will deform elastically. When the stress exceeds the apparent yield stress (marked as “average stress” in Figure 3B), the material will transition from elastic to plastic deformation and begin to flow in its liquid state35. Measuring the storage modulus (G’) and loss modulus (G”) in the mucus-sample under oscillatory stress (or strain) conditions quantifies the change in the material state from gel-like to viscoelastic liquid-like behavior.
The types of rheometer tests performed to monitor data pertaining to storage modulus (G’), loss modulus (G”) and apparent viscosity (η) are described here. The dynamic oscillation tests (strain sweeps and frequency sweeps) monitored G’ and G” under controlled oscillation of cone geometry. The dynamic strain sweep tests determined the linear viscoelastic region (LVR) of the mucus by monitoring the intrinsic material response (Figure 4). Strain sweeps were used to determine the yielding behavior at constant oscillation frequency and temperature. The dynamic frequency sweep tests monitored the material response to increasing frequency (rate of deformation) at a constant amplitude (strain or stress) and temperature. Strain was maintained in the linear viscoelastic region (LVR) for the dynamic frequency sweep tests. The steady-state shear rate tests monitored the apparent viscosity (η) under steady rotation of the cone geometry. The GR mucus was subjected to incremental stress steps and apparent viscosity (η, Pa.s) was monitored for varying shear rate (ý, 1/s).
The protocol presented in this paper treats the GR mucus as a complex structured material of unknown viscoelasticity with a certain linear viscoelastic response range. The fish mucus was extracted from the GRs of the silver carp during a fishing expedition at the Hart creek location in the Missouri River by Professor L. Patricia Hernandez (Department of Biological Sciences, The George Washington University) 1,2,36. An array of GRs inside the mouth of a Silver carp is shown in Figure 5A and a schematic drawing is presented in Figure 5B. An excised GR is shown in Figure 5C. The extraction of mucus from GRs of the silver carp is presented as an example in the schematic drawings, Figure 5D, E. All the rheometer tests were performed under a constant, controlled temperature of 22 ± 0.002 °C, the temperature recorded at the fishing site1,2,36. Each mucus sample was tested three times with the rheometer, and the averaged results are presented along with the statistical error bars.
1. Preparation of the mucus solutions of various concentrations
NOTE: Three concentrations of the mucus solution (400 mg/mL, 200 mg/mL and 100 mg/mL with approximate volumes, 1 mL, 1 mL, and 2 mL, respectively) are prepared for this experiment.
2. Measurements and data acquisition using a rheometer
NOTE: The software used in this protocol for instrument control and data acquisition with rheometer are noted in the Table of Materials. This software will be referred as ‘rheometer instrument control software’.
3. Repeat the protocol for other concentrations of mucus solutions of 200 mg/mL and 400 mg/mL.
4. Graphical representation and data analysis
NOTE: The code provided in the supplemental code file performs data averaging and generates repeatability-errors, overlays the data from all experiments. The standard-deviation calculation features are not available in the rheometer instrument control software. The code is written in a programming language for data analysis, post-processing and graphical representation (see Table of Materials for details).
In this section, we present the results of the experiments on GR mucus using a rotational rheometer with a cone geometry (40 mm diameter, 1° 0’ 11’’) and a Peltier plate. The experiments helped in characterizing the non-Newtonian, shear-thinning behavior of the GR mucus and the apparent yield stress depicting the mucus transition from a gel-like material to a fluid-like material. The representative results entail quantitative descriptions of low-torque limits and secondary flow effects of the rotational rheometer instrumentation. The instrumentation limits and steady-state and dynamic strain rate measurements helped in accurately ascertaining the viscoelastic behavioral trends and apparent yield stress of GR mucus. The apparent yield stress measurements provided a means to observe the minimum stress required for irreversible plastic deformation of the GR mucus and to initiate the flow. The flow initiation tendency of the GR mucus could be attributed to the food particulate adhesion and transport functions. The adhesion and transport functions of GR mucus were macroscopic material attributes that were informed by the rheological measurements in the protocol experiments. Therefore, macro-rheology characterization of the GR mucus was performed with this protocol.
The mucus used in the experiment was acquired from several gill rakers of up to three silver carp and did not have any visible traces of blood1,2. The acquired sample was diluted into two additional samples, as described in the protocol. All measurements were made at the controlled temperature of 22 ± 0.002 °C36. This temperature was maintained on the Peltier plate of the rheometer. The cone geometry was chosen for its versatility in measuring a wide range of viscosities in biological materials such as the GR mucus. The minimum torque under steady state shear conditions (10 x 10-9 Nm), minimum torque under oscillatory conditions (2 x 10-9 Nm) along with the cone angle (1° 0’ 11”) and the summary of rheometer specifications required for assessing low-torque and secondary flow regimes are presented in Table 1. We report replicates of final three samples for the characterization and comparison of non-Newtonian and shear thinning behavior.
Broad inferences after successful execution of the protocol
The successful execution of the protocol and analysis resulted in the characterization of (macro) rheological properties entailing non-Newtonian, shear thinning behavior of mucus extracted from the gill rakers of the silver carp, Hypophthalmichthys molitrix. Particularly, the yielding phenomenon was resolved and an apparent yield stress of the mucus (400 mg/mL concentration, closest to the actual extracted mucus consistency) was ascertained (σy = 0.2736 Pa). The protocol was well-suited for measurements involving very small sample volume (approximately, 1.4 mL) of mucus. Due to scarcity of literature pertaining to GR mucus characterization, these data will aid in analytical modeling and extended rheometric studies.
Results of dynamic sweep experiments
The results of the dynamic frequency and amplitude sweep experiments are presented in this section. These are the results of the procedures created in steps 2.8.2 – 2.8.5. The low-torque limits for oscillatory frequency sweep and amplitude sweep of the GR mucus with 400 mg/mL concentration is marked in Figure 6A, B.
The frequency sweep data (Figure 6A) were acquired for an angular frequency range, 6.28 ≤ ω ≤ 125.66 rad s-1 at a constant oscillatory strain amplitude of 0.01. The angular frequency value, 6.28 rad/s (1 Hz) was chosen as an approximate frequency of motion of palatal folds in the interstices of the gill rakers and marked in Figure 6A. The choice of the strain amplitude value was derived from the linear viscoelastic region ascertained in protocol step 2.7. In Figure 6A, two potential lower limits of the low-torque regime were calculated using strain amplitudes of 0.01 and 0.001 (γ0), and minimum torques of 2×10-9 Nm and 10 x 10-9 Nm (Tmin, See Table 1.), respectively. The data presented in Figure 6A start at the approximated frequency of palatal fold motion (1 Hz or 6.28 rad/s) and increase to higher angular frequencies that are beyond the scope for physical interpretation in this study. Hence, these data were not analyzed further, as they require a more detailed parametric investigation of strain amplitude and palatal fold motion frequency.
The amplitude sweep data (Figure 6B) were acquired at a fixed angular frequency (ω) of 6.28 rad/s (1 Hz). It should be noted that the amplitude sweep data were not affected by the low-torque regime of the instrumentation (Figure 6B). Hence, these data were further analyzed for all three mucus concentrations (100 mg/mL, 200 mg/mL and 400 mg/mL) to determine the extent of viscoelasticity and yielding behavior.
The graphical representation shown in Figure 4 was used as a guideline for extended analysis of the amplitude sweep experiments. The results of three mucus solutions with concentrations 100 mg/mL, 200 mg/mL and 400 mg/mL are discussed below.
The result for 100 mg/mL mucus concentration (Figure 7A), shows that at low oscillation stresses (0.01 ≤ ≤ 0.1 Pa) the storage and loss moduli (G’ and G”) overlapped significantly. At oscillation stresses greater than 0.1 Pa, the storage modulus declines, indicating lower elasticity. The loss modulus, representative of viscosity, remains constant in the full range of oscillation stress (0.01 ≤ ≤ 0.5 Pa). This phenomenon can be attributed to a Newtonian fluid-like behavior and is in agreement with the constant apparent viscosity of the 100 mg/mL mucus concentration (Figure 7A and Figure 8A,B). The corresponding phase angle (δ) data show that at moderate and high oscillation stresses (0.05 ≤ ≤ 0.3 Pa), the values vary between 55° and 70° (Figure 7D). It can therefore be inferred that the 100 mg/mL mucus solution demonstrates fluid-like behavior, with negligible apparent yield stress.
As observed in Figure 7B, 200 mg/mL concentration at low oscillation stresses (0.02 ≤ ≤ 0.04 Pa), the storage modulus (G’) decreases but still remains greater than the loss modulus (G”). Within the oscillation stress range (0.04 ≤ ≤ 0.07 Pa), there was a “crossover” region wherein the G’ and G” values remain approximately equal. This region is marked in Figure 7B with dashed lines and corresponding oscillation stress values were noted (0.04193 ≤ ≤ 0.06467 Pa). Beyond this region, G” attained a higher value than G’ suggesting a transition to a fluid-like behavior. However, G” (representing viscosity) remained constant within the full range of oscillation stress (0.01 ≤ ≤ 0.5 Pa). The phase angle data presented in Figure 7E show a higher degree of variance, especially in the oscillation stress range (0.04193 ≤ ≤ 0.06467 Pa). From Figures 7B,E one can infer that there was transitional behavior of the fluid from linear viscoelastic to non-linear viscoelastic region. Further, the 200 mg/mL mucus concentration represented non-Newtonian characteristics and propensity to yield with in the oscillation stress range, (0.04193 ≤ ≤ 0.06467 Pa). The non-Newtonian, shear-thinning fluid-like behavior agreed with the apparent viscosity data presented in Figure 8A and corresponding stress variations in Figure 8B.
The 400 mg/mL mucus concentration data are presented in Figure 7C,F. The G’ and G” trends in Figure 7C clearly demonstrate a yielding phenomenon with a crossover point between G” and G’. The apparent yield stress (σy) value was recorded as 0.2736 Pa indicating a clear change in the state of the mucus from gel-like to a non-Newtonian fluid-like state. The phase angle data presented in Figure 7F shows a sharp increase at the apparent yield stress (σy = 0.2736 Pa) from approximately 20° to 65°. Such a sharp increase in the phase angle can occur when the material undergoes yielding and starts to flow like a fluid. The non-Newtonian fluid-like behavior can be further supported with the results of steady state shear tests reported in Figure 8A,8B. The apparent yield stress as reported in steady state shear tests was 0.2272 Pa (Figure 8B).
Results of steady state shear rate experiments
The results of the steady state shear rate experiments are presented for three mucus solutions with concentrations 100 mg/mL, 200 mg/mL and 400 mg/mL in this section using graphical representations as a guideline (Figures 3A,B). These results correspond to the procedure initialized in protocol steps 2.8.6-2.8.8.
In Figure 8A for the 100 mg/mL mucus concentration, the apparent viscosity data with high variance at low-shear rates (1 ≤ ≤ 4 s-1) are shown along with the slope, -1.4. The location of low-torque regime is also marked. The high variance of 100 mg/mL data within that range (1 ≤ ≤ 4 s-1), is assumed to be an effect of the (shaded) low-torque regime. In Figure 8B, the corresponding stress variation with shear rate data indicated a small range of shear rates where the sample attained a ‘stress-plateau’ (or the flat region). This region is neglected for yield stress estimation as the corresponding viscosity data are subject to low-torque effects. Within the high-shear rate range (2500 ≤ ≤ 10000 s-1) the apparent viscosity data were affected by the secondary flow regime. The 100 mg/mL mucus solution, therefore, behaves as a Newtonian fluid that is independent of the shear rate outside the low-torque and secondary flow regimes, and with a constant apparent viscosity of 0.00088 Pa s (±1.656 x10-5 Pa s).
As observed in Figure 8A, the 200 mg/mL concentration of mucus remained unaffected by low-torque limits and demonstrated a shear-thinning effect in the shear rate range, 1 ≤ ≤ 15 s-1. The zero-shear strain rate viscosity (ηo) was noted as 0.032 Pa s (±0.024 Pa s) and the infinite-shear strain rate viscosity (η∞) at shear rate (), 1995 s-1, was noted as 0.00085 Pa s (±2.495 x 10-5 Pa s). The shear thinning effect of the fluid was demonstrated with a slope of -1.8 within the shear rate range, 1 ≤ ≤ 4 s-1. The corresponding stress variation in Figure 8B, demonstrates a ‘stress-plateau’ that represents a yielding phenomenon with the average yield stress of 0.1446 Pa (±0.0037 Pa).
The 400 mg/mL concentration of mucus is the least diluted and consequently, the closest in material consistency to the actual extracted GR mucus. In Figure 8A, note that the shear-thinning characteristic is well-defined for the 400 mg/mL mucus concentration when compared to the 200 mg/mL mucus concentration. The zero-shear (ηo) and infinite-shear strain rate viscosities (η∞) at shear rate (), 1995 s-1, were 0.137 Pa s (±0.032 Pa s) and 0.00099 Pa s (±9.323 x 10-5 Pa s), respectively. In addition, the slope of the shear thinning region was established as -0.91 within the shear rate range, 1 ≤ ≤ 32 s-1. The corresponding ‘stress-plateau’ from the stress variations with shear rate observed in Figure 8B, represents an apparent yield stress of 0.2272 Pa (±0.0948 Pa).
Component | Model/Part No./Version | Parameter Description | Specifications |
Rheometer | DHR-2 | Frequency Range | 1 x 10−7 – 100 Hz |
Maximum Angular Velocity | 300 rad/s | ||
Minimum Torque under Steady shear | 10 nN.m | ||
Minimum Torque under Oscillation | 2 nN.m | ||
Maximum Torque | 200 μN *m | ||
Torque Resolution | 0.1 nN.m | ||
Shear Rate Range | 5.73 x 10−6 to 1.72 x 104 [1/s] | ||
Maximum Normal Force | 50 N | ||
Normal Force Resolution | 0.5 mN | ||
Geometry | 513404.905 | Dimensions | 40 mm diameter |
1° 0’ 11” Cone angle | |||
Peltier Plate | 533210.901 | Temperature Range | -40°C to 200°C ± 0.1°C |
Table 1: Specifications of Rheometer
Figure 1: CAD rendering of rheometer components. (A) 40 mm 1° cone geometry, (B) Peltier Plate attachment. The cone geometry should be attached to the shaft of the rheometer, and the Peltier Plate should be connected to the base of the rotational rheometer. Please click here to view a larger version of this figure.
Figure 2: Placement of fluid onto the Peltier plate. The fluid sample should be placed in the center of the Peltier plate to ensure an even spread of fluid throughout the plate when the geometry is lowered. Please click here to view a larger version of this figure.
Figure 3: Graphical representation of steady rheological properties. Variation of (A) Apparent Viscosity (η) and (B) Shear stress () with Shear Strain rate. Please click here to view a larger version of this figure.
Figure 4: Graphical representation of dynamic rheological properties. Variation of Storage (G’) and Loss (G”) Modulii and Phase Angle (δ) with Oscillation Stress. Please click here to view a larger version of this figure.
Figure 5: Representative images and schematic drawings of gill rakers (GR). (A) View of gill raker array and palatal folds (B) Schematic drawing of the gill raker array and palatal folds (C) Excised gill raker (D) Schematic drawing of the gill raker with salient features (E) Location of the mucus extraction in the gill raker. Images 5A and 5C were taken during a dissection performed by Professor L. Patricia Hernandez of the Department of Biological Sciences at The George Washington University. Please click here to view a larger version of this figure.
Figure 6: Low-torque effects. Storage and loss modulus variation for 400 mg/mL, mucus concentration with (A) Frequency sweeps at strain amplitude = 0.01 and (B) Amplitude sweeps at oscillation frequency = 1 Hz (or 6.28 rad/s). Please click here to view a larger version of this figure.
Figure 7: Amplitude sweeps for three concentrations of Silver Carp mucus. Tests were completed at f = 1 Hz or ω = 6.28 rad/s (A) Storage and Loss modulus variation for mucus concentration, 100 mg/mL (B) Storage and Loss modulus variation for mucus concentration, 200 mg/mL (C) Storage and Loss modulus variation for GR mucus concentration, 400 mg/mL (D) Phase angle variation for GR mucus concentration, 100 mg/mL (E) Phase angle variation for mucus concentration, 200 mg/mL (F) Phase angle variation for mucus concentration, 400 mg/mL. Please click here to view a larger version of this figure.
Figure 8: Variation of apparent viscosity (η) and stress (σ) with shear rate () for all three concentrations of fish mucus. (A) Apparent viscosity variation with shear rate for mucus concentrations, 400 mg/mL, 200 mg/mL and 100 mg/mL along with regimes of low-torque effects and secondary flow effects (B) Stress variation with shear rate for GR mucus concentrations, 400 mg/mL, 200 mg/mL and 100 mg/mL, marking the ‘stress plateaus’ (or flat region) with dashed-lines. Dashed lines represent the average apparent yield stress values for the three GR mucus concentrations. Please click here to view a larger version of this figure.
Supplementary Figure 1: Three concentrations of silver carp GR mucus. From left to right: 400 mg/mL, 200 mg/mL, 100 mg/mL. The initial concentration of 400 mg/mL was chosen with the criteria that once diluted the sample would contain a reasonable amount of fish mucus while also providing a large enough volume to run several tests. The two subsequent concentrations were diluted by 50% DI water by volume. Please click here to download this file.
Supplementary Figure 2: Launch of rheometer instrument control software. This software must be launched only after the machine is turned on. Otherwise, the instrument may not be calibrated correctly. Please click here to download this file.
Supplementary Figure 3: Calibration of instrument. Inertial calibration is the only calibration needed for the instrument. There are other calibrations performed after the geometry is installed. Please click here to download this file.
Supplementary Figure 4: Smart swap toggle. This option is for enabling or disabling smart swap. Smart swap is a rheometer instrument control software feature that automatically detects a geometry once it is installed on the rheometer shaft. Please click here to download this file.
Supplementary Figure 5: Measurement Gap Options The “Gap” tab options were accessed to set the conditions for zero-gap mode and traverse velocity of the measurement head. An axial contact force between the geometry and the Peltier plate was set to 1 N to ensure the zero-gap reference, i.e., the contact between cone geometry and the surface of the Peltier plate. The measurement head was then made to accurately traverse to the measurement gap of 24 μm between the 40 mm 1° cone geometry and the Peltier plate. Please click here to download this file.
Supplementary Figure 6: Calibration of geometry attachment. Upon the installation of the 40 mm 1° cone geometry and it’s detection by the rheometer instrument control software, the geometry was calibrated in the same way as the instrument calibration to ensure accurate measurements during an experiment. Please click here to download this file.
Supplementary Figure 7: Zero Gap icon. Zero-gap initialization is performed using this icon. Once zero-gap initialization is complete, the rheometer can accurately reference the position of the shaft, with the geometry attached to it, as it is raised or lowered. Please click here to download this file.
Supplementary Figure 8: CAD rendering of cone geometry and the Peltier plate after zero-gap reference is established. The geometry is set to establish the zero-gap when an axial contact force of 1 N is generated as it makes contact with the Peltier plate. Please click here to download this file.
Supplementary Figure 9: Step 2.5.1 of the protocol. The figure represents the way sample naming and file and data output are set. Please click here to download this file.
Supplementary Figure 10: Step 2.5.2 of the protocol. The figure represents how geometry settings, such as sample volume, geometry gap, and trim gap can be set. For some geometries, namely the cone used in this experiment, these settings are unchangeable and defined based on the geometry. Please click here to download this file.
Supplementary Figure 11: Step 2.5.3 of the protocol. The figure represents how test setting and conditioning are set in this procedural step. Please click here to download this file.
Supplementary Figure 12: Trim gap icon. The trim gap was set so that the geometry could be lowered enough to trim excess fluid that leaks out of the area between the Peltier plate and cone geometry. The gap is dependent upon the geometry in use. For the 40mm, 1° cone geometry used in the protocol, the trim gap was 28 µm. Please click here to download this file.
Supplementary Figure 13: Motor velocity icon. The motor settings were used to adjust rotational velocity of the shaft and minimize the geometry inertia. Please click here to download this file.
Supplementary Figure 14: Geometry gap icon. The geometry gap lowers geometry to a specific distance above the Peltier Plate as specified by the cone-plate geometry. For the 40 mm, 1° cone geometry used in the protocol, the geometry gap is 24 µm. Please click here to download this file.
Supplementary Figure 15: Start icon. The start button initiates the entire sequence of procedures that were set up earlier. Please click here to download this file.
Supplementary Figure 16: Changing the Graph Variables. The figure represents the variables that can be defined for data presentation when running the procedures. Particularly, oscillation strain and oscillation stress are of importance during the dynamic sweep experiments in the protocol. Please click here to download this file.
Supplementary Figure 17: Exporting of rheometer instrument control software files to a spreadsheet software. After the files are exported as spreadsheets, the data analysis was made feasible using other programming software. Please click here to download this file.
Supplemental Code File: Post-processing of data files using data analysis program. Please click here to download this file.
One of the main objectives of developing this protocol is to establish that it is well-suited for rheological characterization of GR mucus when very small sample volumes are available. We acknowledge that more samples from a school of silver carp are needed to fully characterize the rheological properties of the GR mucus and the data presented herein are not a generalization across the entire silver carp population. Our technique is justified because of its efficacy with rheological characterization of small sample volumes and with extended investigations involving larger ensembles of mucus samples.
The critical steps within the protocol are the preparation of mucus solutions of various concentrations, measurements and data acquisition using a rotational rheometer, and graphical representation and data analysis for physical insights.
Physical insights into GR mucus data are drawn from graphical representations shown in Figure 3 and Figure 4, that are annotated with attributes of the expected material behavior. Zero-shear strain rate viscosity (ηo) values can be observed at low-shear strain rates where mobility of the material molecules dominates (Figure 3A and Figure 8A). Infinite-shear strain viscosity (η∞) values in non-Newtonian fluids are orders of magnitudes lower than the zero-shear strain rate viscosity. These data can be noticed at high shear rates where there is little or no dependence on intermolecular interactions (Figure 3A and Figure 8A). For non-Newtonian fluids, apparent viscosities progressively decrease as the shear rates increase and attain a constant low value (Figure 3A and Figure 8A). Yielding behavior in the GR mucus under steady state measurements can be represented with slope as shown in Figure 3A and presented in Equation 1., where ηa represents the apparent viscosity, σy is the (constant) yield stress and is the shear strain rate.
Figure 3A and Figure 8B are presented on a log-log scale and therefore, Equation 1 attains the following form:
where k – represents the apparent yield stress. On a log-log scale, the apparent viscosity decreases with a slope of ‘-1’ indicates the material yield as shown in Figure 3A10. The 200 mg/mL and 400 mg/mL mucus concentrations possessed slopes of -1.8 and -0.91, respectively, and demonstrate the yielding behavior (Figure 8A). Under dynamic oscillation measurements, the viscoelastic characteristics are independent of the strain amplitude in the Linear Viscoelastic Region (LVR) (Figure 4). The yielding behavior in the GR mucus under dynamic oscillation measurements can be observed as the viscoelastic material (GR mucus) enters the non-linear viscoelastic region (NLVR) as the storage modulus (G’) decreases (Figure 4). In the NLVR regime the viscoelastic material will demonstrate solid-gel-like behavior if the storage modulus is greater than the loss modulus (G’ > G”). When the loss modulus exceeds the storage modulus (G’ < G”), a “crossover” between G’ and G” data occur. As shown in Figure 7B,C, the 200 mg/mL and 400 mg/mL GR mucus concentrations demonstrated fluid-like behavior marked by the “crossover” between G’ and G” data. The apparent yield stress under steady state measurements is represented as the average value of stress until an inflection point is reached (Figure 3B). Thereafter, the stress begins to increase sharply with an increase in the shear strain rate as shown in Figure 3B and Figure 8B. The GR mucus data (200 mg/mL and 400 mg/mL concentrations) showed shear-thinning fluid behavior until the material begins to yield (Figure 8A,B). The apparent yield stress was observed clearly in 200 mg/mL and 400 mg/mL mucus concentrations due to their non-Newtonian characteristics (Figure 8B). The apparent yield stress under dynamic oscillation measurements are shown in Figure 4 and Figure 7B,C as the “crossover” region between G’ and G” data, followed by G” values exceeding G’. The 400 mg/mL GR mucus data showed shear-thinning, non-Newtonian behavior. The onset point of material yield was observed with an apparent yield stress of approximately 0.2736 Pa (Figure 7C). The hydrogel-to-fluid like transition with phase angle (δ = tan-1 (G”/G’)) changes are presented in Figures 4 and 7D-F. The extrema in the phase angle is associated with a Hookean solid at 0˚ and viscous fluid at 90˚ as shown in Figure 4. The phase angle values around 45˚ were attributed to transition of gel-like behavior of the material to a fluid-like behavior. The 400 mg/mL mucus concentration clearly showed a change in the material characteristic from hydrogel to fluid like behavior through the process of yielding with an apparent yield stress of ~ 0.2736 Pa (Figure 7F).
Understanding the measurement limitations and avoiding data unsuitable for physical interpretation is a challenge with complex and soft biological fluids, especially in studies involving small sample volumes11. The data generated under low-torque and secondary flow effects are unsuitable for physical interpretation and are dependent on the geometry used in the rheometer (such as cone and plate in this study). These regimes were identified to avoid any misrepresentation of experimental data suffering from instrument resolution and measurement artifacts due to momentum diffusion. Low-torque limits (Figure 6A and Figure 8A) are functions of geometry and minimum torque generated by the instrument (Table 1). Under steady shear measurement conditions, the criterion for rejecting data affected by the low-torque limit for a cone-plate geometry of radius (R) with minimum torque (Tmin = 10 x 10-9 Nm, Table 1) has been discussed by Ewoldt et al. and is presented below11:
where is the shear strain rate. Unlike the 100 mg/mL GR mucus concentrations, the 200 mg/ml and 400 mg/mL GR mucus concentrations were unaffected by low-torque effects clearly demonstrate non-Newtonian, shear thinning behavior with high zero-shear strain rate viscosities at low shear strain rates. The criterion for minimum measurable viscoelastic moduli under dynamic oscillation measurements has been discussed by Ewoldt et al. and presented below (Equation 4)11. In Equation 4, for a cone-plate geometry of radius (R) the minimum torque under oscillatory shear (Tmin = 2 x 10-9 Nm, Table 1).
where Gmin is the storage modulus (G’) or loss modulus (G”) and is the shear strain rate. The regimes of instrumentation limitation governed by low-torque effects are marked in Figure 6A and 6B. The secondary flow regime under steady state measurements is governed by an inward momentum diffusion of the fluid by means of an eddy residing within the rotational cone and plate geometry11. The secondary flow pattern increases torque incorrectly making the fluid appear to be shear-thickening (Figure 8A). The secondary flow limit, proposed by Ewoldt et al. in Figure 8A was drawn using the following relation11:
where L = βR, β is the cone angle, R is the cone radius, ρ = 1000 kg m-3, Recrit = 4 and is the shear rate. This regime helped in estimating the infinite-shear strain viscosity (η∞) values in GR mucus samples.
A modification of the protocol can be made by using a flat-plate geometry instead of the cone-plate geometry as shown in protocol presented herein. The flat-plate tests should be performed with a parametric variation of the measurement gap in the rotational rheometer to reveal the dependence of apparent yield stress on the measurement gap and geometry. The suggested improvements of the protocol presented in this paper are described below. A parametric variation of the strain amplitude in the linear viscoelastic regime (LVR) and oscillation frequency should be performed. ‘Tack and peel’ rheology tests should be performed to develop a full understanding of the adhesivity of the GR mucus. Rheology characteristics of GR mucus should be performed on larger sample volume ensembles along with studies to measure any traces of blood cells to account its effect on the overall GR rheological properties.
The limitations of the protocol are described below. The intricacies of the GR mucus extraction-procedures and the presence of blood cells or tissue fragments in the mucus samples may influence the rheology of the mucus. However, it should be noted the mucus used in the protocol did not have any visible traces of blood. The GR mucus sample is a heterogeneous material and can possess different rheological properties due to the variance in location of and the conditions post-extraction. This limitation was addressed by sufficiently homogenizing the GR mucus using a shaker to breakdown any large clumps of mucus and tissue presence. Another important limitation is the very small GR mucus sample volumes (approximately 1.4 mL), harvested for analyses that constrain a generalization of GR mucus properties.
The significance of this protocol is that it allows for an accurate rheological characterization of non-Newtonian, biological fluids such as the mucus. The protocol presented herein paves the way for investigating other similar biological fluids associated with human, animal and plant secretions. In addition, synthetic fluids or polymer-based solutions that are analogs of biological fluids can be testing using this protocol to understand material properties under varying stresses, oscillation frequencies, and temperature. The protocol is well-suited for rheological characterization of biological fluids when very small sample volumes are made available.
The extended outcome of the protocol is that the apparent viscosity and apparent yield stress of GR mucus will facilitate the creation of analytical models to interpret results from fundamental hydrodynamic investigation of filter feeding and advance technologies requiring and involving crossflow and membrane filtration.
The macro-rheological study posits that the mucus in contact with food particulates is initially, in a gel-like state that serves as an adhesive. Upon initiation of flow and shear forces the mucus attains an apparent yield stress and undergoes plastic deformation. The protocol execution using a rotational rheometer helped in characterization of the transition of mucus from gel-like to fluid-like behavior. This transition was experimentally observed, and the apparent yield stress was recorded at 0.2736 Pa in rotational rheometer experiments. When the external stresses on the mucus are less than the apparent yield stress, the mucus will demonstrate gel-like behavior to facilitate adhesion of food particulates. When the external stresses exceed the apparent yield stress, mucus will demonstrate shear-thinning behavior that will facilitate transport of agglomerated food particulates to the digestive organs in the silver carp.
The authors have nothing to disclose.
The authors acknowledge support and funding from the GW Center for Biomimetics and Bioinspired Engineering. We thank Professor L. Patricia Hernandez of the Department of Biological Sciences at The George Washington University for inspiring the investigation and ongoing collaboration, providing biological expertise on the physiology of the silver carp and providing the mucus samples. We thank the students, Mr. David Palumbo, Ms. Carly Cohen, Mr. Isaac Finberg, Mr. Dominick Petrosino, Mr. Alexis Renderos, Ms. Priscilla Varghese, Mr. Carter Tegen and Mr. Raghav Pajjur for help in the laboratory and Mr. Thomas Evans and Mr. James Thomas of TA Instruments, New Castle, DE for support with training and maintenance of the rheometer. Images for Figures 5A,C were taken during a dissection performed by Professor L. Patricia Hernandez of the Department of Biological Sciences at The George Washington University.
Materials | |||
Kim Wipes | VWR | 470224-038 | To clean Sample from plate |
Gloves | VWR | 89428-750 | To prevent contamination of sample |
Pipette | VWR | 89079-974 | To transport sample from vial to rheometer |
Pipette Tips | Thermo Scientific | 72830-042 | To transport sample from vial to rheometer |
Shaker | VWR | 89032-094 | To homogenously mix sample of mucus |
Vials | VWR | 66008-710 | Contains measured sample volumes |
Weigh Scale | Ohaus | Scout –SPX Balances | To weigh mass of mucus samples |
Chemical Reagents | |||
De-Ionized Water (H20) | – | – | Liquid |
Sterile 70% Isopropanol (C3H8O) | VWR | 89108-162 | Liquid |
GR Mucus | |||
100 mg/mL concentration, 2mL | – | – | Viscoelastic Material |
400 mg/mL concentration, 1mL | – | – | Viscoelastic Material |
200 mg/mL concentration, 1mL | – | – | Viscoelastic Material |
Software | |||
MATLAB | Mathworks | R2017a | Data analysis, post-processing and graphical representation |
Trios | TA Instruments | v4.5.042498 | Rheometer instrument control and analysis software |