Elasticity and stress relaxation of a very small vocal fold

Across mammals many vocal sounds are produced by airflow induced vocal fold oscillation. We tested the hypothesis that stress-strain and stress-relaxation behavior of rat vocal folds can be used to predict the fundamental frequency range of the species' vocal repertoire. In a first approximation vocal fold oscillation has been modeled by the string model but it is not known whether this concept equally applies to large and small species. The shorter the vocal fold, the more the ideal string law may underestimate normal mode frequencies. To accommodate the very small size of the tissue specimen, a custom-built miniaturized tensile test apparatus was developed. Tissue properties of 6 male rat vocal folds were measured. Rat vocal folds demonstrated the typical linear stress-strain behavior in the low strain region and an exponential stress response at strains larger than about 40%. Approximating the rat's vocal fold oscillation with the string model suggests that fundamental frequencies up to about 6 kHz can be produced, which agrees with frequencies reported for audible rat vocalization. Individual differences and time-dependent changes in the tissue properties parallel findings in other species, and are interpreted as universal features of the laryngeal sound source.     

Optical measurements of vocal fold tensile properties: implications for phonatory mechanics

In voice research, in vitro tensile stretch experiments of vocal fold tissues are commonly employed to determine the tissue biomechanical properties. In the standard stretch-release protocol, tissue deformation is computed from displacements applied to sutures inserted through the thyroid and arytenoid cartilages, with the cartilages assumed to be rigid. Here, a non-contact optical method was employed to determine the actual tissue deformation of vocal fold lamina propria specimens from three excised human larynges in uniaxial tensile tests. Specimen deformation was found to consist not only of deformation of the tissue itself, but also deformation of the cartilages, as well as suture alignment and tightening. Stress-stretch curves of a representative load cycle were characterized by an incompressible Ogden model. The initial longitudinal elastic modulus was found to be considerably higher if determined based on optical displacement measurements than typical values reported in the literature. The present findings could change the understanding of the mechanics underlying vocal fold vibration. Given the high longitudinal elastic modulus the lamina propria appeared to demonstrate a substantial level of anisotropy. Consequently, transverse shear could play a significant role in vocal fold vibration, and fundamental frequencies of phonation should be predicted by beam theories accounting for such effects.     

Adapted to roar: functional morphology of tiger and lion vocal folds

Vocal production requires active control of the respiratory system, larynx and vocal tract. Vocal sounds in mammals are produced by flow-induced vocal fold oscillation, which requires vocal fold tissue that can sustain the mechanical stress during phonation. Our understanding of the relationship between morphology and vocal function of vocal folds is very limited. Here we tested the hypothesis that vocal fold morphology and viscoelastic properties allow a prediction of fundamental frequency range of sounds that can be produced, and minimal lung pressure necessary to initiate phonation. We tested the hypothesis in lions and tigers who are well-known for producing low frequency and very loud roaring sounds that expose vocal folds to large stresses. In histological sections, we found that the Panthera vocal fold lamina propria consists of a lateral region with adipocytes embedded in a network of collagen and elastin fibers and hyaluronan. There is also a medial region that contains only fibrous proteins and hyaluronan but no fat cells. Young's moduli range between 10 and 2000 kPa for strains up to 60%. Shear moduli ranged between 0.1 and 2 kPa and differed between layers. Biomechanical and morphological data were used to make predictions of fundamental frequency and subglottal pressure ranges. Such predictions agreed well with measurements from natural phonation and phonation of excised larynges, respectively. We assume that fat shapes Panthera vocal folds into an advantageous geometry for phonation and it protects vocal folds. Its primary function is probably not to increase vocal fold mass as suggested previously. The large square-shaped Panthera vocal fold eases phonation onset and thereby extends the dynamic range of the voice.     

Mechanical stress during phonation in a self-oscillating finite-element vocal fold model

The stress information during phonation in the vocal folds is helpful in understanding the etiologies of vocal trauma and its related vocal diseases, such as nodules. In this paper, a self-oscillating finite-element model, which combines aerodynamic properties, tissue mechanics, airflow-tissue interactions, and vocal fold collisions, was used to simulate the vocal fold vibration during phonation. The spatial and temporal characteristics of mechanical stress in the vocal folds were predicted by this model. Temporally, it was found that mechanical stress periodically undulates with vibration of the vocal folds and that vocal fold impact causes a jump in the normal stress value. Spatially, the normal stress is significantly higher on the vocal fold surface than inside of the vocal folds. At the midpoint of the medial surface, the peak-to-peak amplitude of the normal stress reaches its maximum value. Using different lung pressures (0-1.5kPa) to drive the self-oscillating model, we found that lower lung pressure can effectively decrease the mechanical stress in the vocal folds. This study supports the fatigue damage hypothesis of vocal trauma. With this hypothesis and the numerical simulation in this study, the clinical observations of vocal fold trauma risk can be explained. This implies the mechanical stress predicted by this self-oscillating model could be valuable for predicting, preventing, and treating vocal fold injury.     

Phase relationship between dynamics of the subglottic pressure and oscillatory movement of the vocal folds. I. Sustained phonation

The phase relationship between subglottic pressure and vocal fold length has been studied during sustained phonation in five subjects with normal larynx. Pressure was measured by tracheal puncture and vocal fold length was deduced from simultaneous measurement of translaryngeal impedance in the horizontal plane and transglottal light flux in the vertical plane. The pressure sine wave shows a phase lead of slightly less than 90 degrees relative to the length sine wave. Thus during sustained phonation the vocal apparatus behaves like a harmonic oscillator; the frequency of oscillation is determined by the mechanical parameters of the vibrating system; the source of periodic energy supply is the subglottal pressure wave.     

Modeling viscous dissipation during vocal fold contact: the influence of tissue viscosity and thickness with implications for hydration

The mechanics of vocal fold contact during phonation is known to play a crucial role in both normal and pathological speech production, though the underlying physics is not well understood. Herein, a viscoelastic model of the stresses during vocal fold contact is developed. This model assumes the cover to be a poroelastic structure wherein interstitial fluid translocates in response to mechanical squeezing. The maximum interstitial fluid pressure is found to generally increase with decreasing viscous dissipation and/or decreasing tissue elasticity. A global minimum in the total contact stress, comprising interstitial fluid pressure and elastic stress in the tissue, is observed over the studied dimensionless parameter range. Interestingly, physiologically reasonable estimates for the governing parameters fall within this global minimum region. The model is validated against prior experimental and computational work, wherein the predicted contact stress magnitude and impact duration agree well with published results. Lastly, observations of the potential relationship between vocal fold hydration and increased risk of tissue damage are discussed based upon model predictions of stress as functions of cover layer thickness and viscosity.     

Changes in expression of extracellular matrix genes, fibrogenic factors, and actin cytoskeletal organization in retinol treated and untreated vocal fold stellate cells.

The regulation of extracellular matrix (ECM) constituency is critical in maintaining vocal cord biomechanical viscoelasticity required for phonation. Recently our laboratory successfully isolated and cultured a novel cell called a vocal fold stellate cell (VFSC), thought to play a central role in laryngeal ECM metabolism, aging, scarring and cancer. Our laboratory has shown that these cells undergo transdifferentiation that is partially reversed by exposure to all-trans retinol (ATROH). Here we make the first report on the expression of various ECM components, MMPs, TIMPs, pro-fibrogenic cytokines, and other ECM modulators in transdifferentiated and deactivated VFSCs. We show that VFSCs maintain an ECM expression pattern similar to laryngeal cancer and scars but distinct from tracheal fibroblasts. Exposure to ATROH differentially affects the VFSC expression of ECM components, matrix-regulating enzymes, and fibrogenic factors suggesting that the inhibitory effects of this synthetic cofactor should be studied further in laryngeal fibrosis and scarring. We also show that increased exposure to retinol induces sequential reorganization of the actin cytoskeleton in activated VFSCs. Our findings demonstrate that VFSCs are capable of regulating vocal fold ECM constituency important throughout normal laryngeal development. Furthermore, our results implicate VFSC activation in ECM misregulation which is a hallmark of several laryngeal pathologies.     

Non-invasive in vivo measurement of the shear modulus of human vocal fold tissue

Voice is the essential part of singing and speech communication. Voice disorders significantly affect the quality of life. The viscoelastic mechanical properties of the vocal fold mucosa determine the characteristics of the vocal folds oscillations, and thereby voice quality. In the present study, a non-invasive method was developed to determine the shear modulus of human vocal fold tissue in vivo via measurements of the mucosal wave propagation speed during phonation. Images of four human subjects' vocal folds were captured using high speed digital imaging (HSDI) and magnetic resonance imaging (MRI) for different phonation pitches, specifically fundamental frequencies between 110 and 440 Hz. The MRI images were used to obtain the morphometric dimensions of each subject's vocal folds in order to determine the pixel size in the high-speed images. The mucosal wave propagation speed was determined for each subject and at each pitch value using an automated image processing algorithm. The transverse shear modulus of the vocal fold mucosa was then calculated from a surface (Rayleigh) wave propagation dispersion equation using the measured wave speeds. It was found that the mucosal wave propagation speed and therefore the shear modulus of the vocal fold tissue were generally greater at higher pitches. The results were in good agreement with those from other studies obtained via in vitro measurements, thereby supporting the validity of the proposed measurement method. This method offers the potential for in vivo clinical assessments of vocal folds viscoelasticity from HSDI.     

Hyaluronic acid-based microgels and microgel networks for vocal fold regeneration

Vocal fold scarring disrupts the viscoelastic properties of the lamina propria that are critical for normal phonation. There is a clinical need for the development of advanced biomaterials that approximate the mechanical properties of the lamina propria for in vivo vocal fold regeneration. We have developed hyaluronic acid (HA)-based microgels and cross-linked microgel networks with tunable degradation and mechanical properties. HA microgels were prepared by cross-linking HA derivatives carrying hydrazide (HAADH) and aldehyde (HAALD) functionalities within the inverse emulsion droplets. Alternatively, poly(ethylene glycol) dialdehyde (PEGDiALD) was employed in place of HAALD. Microgels based on HAADH/HAALD are more resistant to enzymatic degradation than those generated from HAADH/PEGDiALD. In vitro cytotoxicity studies using vocal fold fibroblasts indicate that microgels synthesized from HAADH/HAALD are essentially nontoxic, whereas microgels derived from HAADH/PEGDiALD exhibit certain adverse effects on the cultured cells at high concentration (> or =2 mg/mL). These microgels exhibit residual functional groups that can be used as reactive handles for covalent conjugation of therapeutic molecules. The presence of residual functional groups also allows for subsequent cross-linking of the microgels with other reactive polymers, giving rise to doubly cross-linked networks (DXNs) with tunable viscoelasticity. Mechanical measurements using a torsional wave apparatus indicate that HA-based DXNs exhibit elastic moduli that are similar to those of vocal fold lamina propria at frequencies close to the range of human phonation. These HA-based microgel systems are promising candidates for the treatment of vocal fold scarring, not just as biocompatible filler materials, but as smart entities that can repair focal defects, smooth the vocal fold margin, and potentially soften and dissolve scar tissue.     

Mechanotransduction of vocal fold fibroblasts and mesenchymal stromal cells in the context of the vocal fold mechanome

The design of cell-based therapies for vocal fold tissue engineering requires an understanding of how cells adapt to the dynamic mechanical forces found in the larynx. Our objective was to compare mechanotransductive processes in therapeutic cell candidates (mesenchymal stromal cells from adipose tissue and bone marrow, AT-MSC and BM-MSC) to native cells (vocal fold fibroblasts-VFF) in the context of vibratory strain. A bioreactor was used to expose VFF, AT-MSC, and BM-MSC to axial tensile strain and vibration at human physiological levels. Microarray, an empirical Bayes statistical approach, and geneset enrichment analysis were used to identify significant mechanotransductive pathways associated with the three cell types and three mechanical conditions. Two databases (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes) were used for enrichment analyses. VFF shared more mechanotransductive pathways with BM-MSC than with AT-MSC. Gene expression that appeared to distinguish the vibratory strain condition from polystyrene condition for these two cells types related to integrin activation, focal adhesions, and lamellipodia activity, suggesting that vibratory strain may be associated with cytoarchitectural rearrangement, cell reorientation, and extracellular matrix remodeling. In response to vibration and tensile stress, BM-MSC better mimicked VFF mechanotransduction than AT-MSC, providing support for the consideration of BM-MSC as a cell therapy for vocal fold tissue engineering. Future research is needed to better understand the sorts of physical adaptations that are afforded to vocal fold tissue as a result of focal adhesions, integrins, and lamellipodia, and how these adaptations could be exploited for tissue engineering.     

Mechanism of initiation of oscillatory motion in human glottis

Phonatory onsets of 17 normal subjects under usual speech conditions were investigated by measuring the cross-sectional area of the glottic aperture, using a photoelectric device. During a normal soft or breathy onset of phonation, the vocal fold vibration begins with an open glottis; the pattern of the optical signal is a sine wave of increasing amplitude with one to eight cycles before steady state oscillation is achieved. The first deflection of the base-line is either towards adduction or towards abduction. The classically assumed--since van den Berg et al. (1957)--retro-aspiration phenomenon, consisting in a narrowing of the glottic chink due to the accelerated air flow, according to Bernoulli's law, is incompatible with these observations. An aerodynamic study, with accurate calibration of the photoelectric transducer in one trained subject, as well as flux and subglottic pressure measurements, suggests that the gas flux (air as well as a helium-oxygen mixture) reaches the condition of turbulence at the level of the glottic nozzle, just before vocal folds are set into oscillatory motion. The setting in motion of the free edge of the vocal folds in normal soft or breathy onset of phonation can be explained by a sudden modification of flow conditions within the expiratory gas : the flow is laminar in the trachea and suddenly becomes turbulent at the level of the glottic nozzle. On the other hand, approximately normal atmospheric pressure values are attained due to a Bernoulli-effect at that level, allowing the vocal folds to oscillate sinusoidally about their vibration axes, corresponding to their virtual resting position, like a forced oscillator.     

Modeling mechanical stresses as a factor in the etiology of benign vocal fold lesions

Vocal fold tissue lesions such as nodules and polyps are thought to develop in response to mechanical stress that occurs during vocal fold collision. Two computational models of vocal fold collision during voice production are used to investigate this hypothesis. A one-dimensional lumped mass model, whose parameters are derived from vocal fold tissue dimensions and material properties, predicts stress perpendicular to the direction of impact (normal stress). A previously published three-dimensional finite element model that incorporates the same dimensions and properties predicts the entire stress tensor. The hypothesis is supported by predictions from the finite element model that three components of normal stress and one component of shear stress are increased during collision in the typical location of lesions (i.e. the center of the superior medial edge of the vocal fold in the middle of the vibrating and contact region). The lumped mass model predicts that mechanical stress is negatively correlated with mucosal thickness (increased by voice warm-up and hydration), is positively correlated with driving force (proportional to voice intensity), and is affected by voice production method. These relationships are consistent with clinical observations of vocal fold lesion risk factors and have implications for improving prevention and treatment of benign vocal fold lesions.     

Material characterization of the pig kidney in relation with the biomechanical analysis of renal trauma

The objective of this study was an investigation of the material properties of the fresh pig kidney and parametric characterization of its elastic and inelastic material behavior. The material investigation included density measurements, uniaxial as well as three-dimensional compression tests, tensile tests. and shear tests on the samples extracted from the fresh pig kidney. For comparison, density measurements on a number of soft synthetic materials were also performed. Compression tests on the radial and the tangential specimens from the cortex tissue were performed at various loading rates. Three-axial compression tests were performed on the cortex tissues placed in a compression chamber. Shear tests were performed by punching a cylinder into a slice of the cortex. Tensile tests were carried out on the outer capsule. For characterization of the material behavior, a non-linear theoretical simulation based on a two parameter Blatz model was used. For characterization of the time-dependent behavior of the pig kidney cortex, a four-parameter linear viscoelastic model was employed. From the present experimental and theoretical studies, a number of conclusions were derived: (1) The general behavior of the pig kidney cortex samples under compression showed the general non-linear features typical of the soft tissues; the stress strain diagram was composed of a very flat part at very low stress level to about 30% relative deformation which was followed by a steeply rising stiffening leading to the radial rupture of samples marked by a maximum nominal rupture strain of about 50%. (2) The uniaxial compression tests on the radial and the tangential samples from the cortex tissue showed an increase of the rupture stress with the increase in the loading rate, but a decrease in the related rupture strain. (3) The long-term uniaxial compression tests on the cortex specimens under sustained constant load showed an instantaneous deformation followed by a creep response which eventually approached an asymptote. (4) Simulation of the non-linear material behavior of the cortex tissue under uniaxial compression by the Blatz model gave two pairs of material parameters for the cortex in the radial and the tangential directions. Furthermore, fitting of the assumed four-parameter linear viscoelastic model with the experimental data resulted in the viscoelastic material parameters.     

Effect of the ventricular folds in a synthetic larynx model

Within the human larynx, the ventricular folds serve primarily as a protecting valve during swallowing. They are located directly above the sound-generating vocal folds. During normal phonation, the ventricular folds are passive structures that are not excited to periodical oscillations. However, the impact of the ventricular folds on the phonation process has not yet been finally clarified.An experimental synthetic human larynx model was used to investigate the effect of the ventricular folds on the phonation process. The model includes self-oscillating vocal fold models and allows the comparison of the pressure distribution at multiple locations in the larynx for configurations with and without ventricular folds.The results indicate that the ventricular folds increase the efficiency of the phonation process by reducing the phonation threshold level of the pressure below the vocal folds. Two effects caused by the ventricular folds could be identified as reasons: (1) a decrease in the mean pressure level in the region between vocal and ventricular folds (ventricles) and (2) an increase in the glottal flow resistance.The reason for the first effect is a reduction of the pressure level in the ventricles due to the jet entrainment and the low static pressure in the glottal jet. The second effect results from an increase in the glottal flow resistance that enhances the aerodynamic energy transfer into the vocal folds. This effect reduces the onset threshold of the pressure difference across the glottis.     

Cross-sample validation provides enhanced proteome coverage in rat vocal fold mucosa

The vocal fold mucosa is a biomechanically unique tissue comprised of a densely cellular epithelium, superficial to an extracellular matrix (ECM)-rich lamina propria. Such ECM-rich tissues are challenging to analyze using proteomic assays, primarily due to extensive crosslinking and glycosylation of the majority of high M(r) ECM proteins. In this study, we implemented an LC-MS/MS-based strategy to characterize the rat vocal fold mucosa proteome. Our sample preparation protocol successfully solubilized both proteins and certain high M(r) glycoconjugates and resulted in the identification of hundreds of mucosal proteins. A straightforward approach to the treatment of protein identifications attributed to single peptide hits allowed the retention of potentially important low abundance identifications (validated by a cross-sample match and de novo interpretation of relevant spectra) while still eliminating potentially spurious identifications (global single peptide hits with no cross-sample match). The resulting vocal fold mucosa proteome was characterized by a wide range of cellular and extracellular proteins spanning 12 functional categories.     

A numerical analysis of phonation using a two-dimensional flexible channel model of the vocal folds.

A two-dimensional flexible channel model of the vocal folds coupled with an unsteady one-dimensional flow model is presented for an analysis of the mechanism of phonation. The vocal fold is approximated by springs and dampers distributed in the main flow direction that are enveloped with an elastic cover. In order to approximate three-dimensional collision of the vocal folds using the two-dimensional model, threshold values for the glottal width are introduced. The numerical results show that the collision plays an important role in speech sound, especially for higher resonant frequency components, because it causes the source sound to include high-frequency components.     

Synthetic,Multi-Layer,Self-Oscillating Vocal Fold Model Fabrication

Sound for the human voice is produced via flow-induced vocal fold vibration. The vocal folds consist of several layers of tissue, each with differing material properties ¹. Normal voice production relies on healthy tissue and vocal folds, and occurs as a result of complex coupling between aerodynamic, structural dynamic, and acoustic physical phenomena. Voice disorders affect up to 7.5 million annually in the United States alone ² and often result in significant financial, social, and other quality-of-life difficulties. Understanding the physics of voice production has the potential to significantly benefit voice care, including clinical prevention, diagnosis, and treatment of voice disorders.Existing methods for studying voice production include in vivo experimentation using human and animal subjects, in vitro experimentation using excised larynges and synthetic models, and computational modeling. Owing to hazardous and difficult instrument access, in vivo experiments are severely limited in scope. Excised larynx experiments have the benefit of anatomical and some physiological realism, but parametric studies involving geometric and material property variables are limited. Further, they are typically only able to be vibrated for relatively short periods of time (typically on the order of minutes).Overcoming some of the limitations of excised larynx experiments, synthetic vocal fold models are emerging as a complementary tool for studying voice production. Synthetic models can be fabricated with systematic changes to geometry and material properties, allowing for the study of healthy and unhealthy human phonatory aerodynamics, structural dynamics, and acoustics. For example, they have been used to study left-right vocal fold asymmetry ^3,4, clinical instrument development ⁵, laryngeal aerodynamics ^6-9, vocal fold contact pressure ¹⁰, and subglottal acoustics ¹¹ (a more comprehensive list can be found in Kniesburges et al. ¹²)Existing synthetic vocal fold models, however, have either been homogenous (one-layer models) or have been fabricated using two materials of differing stiffness (two-layer models). This approach does not allow for representation of the actual multi-layer structure of the human vocal folds ¹ that plays a central role in governing vocal fold flow-induced vibratory response. Consequently, one- and two-layer synthetic vocal fold models have exhibited disadvantages ^3,6,8 such as higher onset pressures than what are typical for human phonation (onset pressure is the minimum lung pressure required to initiate vibration), unnaturally large inferior-superior motion, and lack of a "mucosal wave" (a vertically-traveling wave that is characteristic of healthy human vocal fold vibration).In this paper, fabrication of a model with multiple layers of differing material properties is described. The model layers simulate the multi-layer structure of the human vocal folds, including epithelium, superficial lamina propria (SLP), intermediate and deep lamina propria (i.e., ligament; a fiber is included for anterior-posterior stiffness), and muscle (i.e., body) layers ¹. Results are included that show that the model exhibits improved vibratory characteristics over prior one- and two-layer synthetic models, including onset pressure closer to human onset pressure, reduced inferior-superior motion, and evidence of a mucosal wave.