- Research
- Open access
- Published:
Proteomic analysis reveals that aging rabbit vocal folds are more vulnerable to changes caused by systemic dehydration
BMC Genomics volume 23, Article number: 762 (2022)
Abstract
Background
Older adults are more prone to develop systemic dehydration. Systemic dehydration has implications for vocal fold biology by affecting gene and protein expression. The objective of this study was to quantify vocal fold protein changes between two age groups and hydration status, and to investigate the interaction of age and hydration status on protein expression, which has not been investigated in the context of vocal folds before. Comparative proteomics was used to analyze the vocal fold proteome of 6.5-month-old and > 3-year-old rabbits subjected to water ad libitum or water volume restriction protocol.
Results
Young and older adult rabbits (n = 22) were either euhydrated (water ad libitum) or dehydrated by water volume restriction. Dehydration was confirmed by body weight loss of − 5.4% and − 4.6% in young and older groups, respectively, and a 1.7-fold increase of kidney renin gene expression in the young rabbits. LC-MS/MS identified 2286 proteins in the rabbit vocal folds of young and older adult rabbits combined. Of these, 177, 169, and 81 proteins were significantly (p ≤ 0.05) affected by age, hydration status, or the interaction of both factors, respectively. Analysis of the interaction effect revealed 32 proteins with opposite change patterns after dehydration between older and young rabbit vocal folds, while 31 proteins were differentially regulated only in the older adult rabbits and ten only in the young rabbits in response to systemic dehydration. The magnitude of changes for either up or downregulated proteins was higher in the older rabbits. These proteins are predominantly related to structural components of the extracellular matrix and muscle layer, suggesting a disturbance in the viscoelastic properties of aging vocal fold tissue, especially when subjected to systemic dehydration.
Conclusions
Water restriction is a laboratory protocol to assess systemic dehydration-related changes in the vocal fold tissue that is translatable to human subjects. Our findings showed a higher number of proteins differentially regulated with a greater magnitude of change in the vocal folds of older adult rabbits in the presence of systemic dehydration compared to younger rabbits. The association of these proteins with vocal fold structure and biomechanical properties suggests that older human subjects may be more vulnerable to the effects of systemic dehydration on vocal function. The clinical implications of these protein changes warrant more investigation, but age should be taken into consideration when evaluating vocal treatment recommendations that interfere with body fluid balance.
Background
Vocal folds are specialized viscoelastic tissue located in the larynx responsible for vibrating the airstream to produce sound. The larynx and vocal folds age as people age, resulting in physiological speech and voice changes [1,2,3]. The process of physiological voice aging is referred to as presbyphonia [4]. Several studies have described the changes associated with this process [2, 3, 5,6,7]. Generally, age-related alterations are attributed to the anatomy and physiology of the structural support of vocal folds, such as muscle and vocal fold atrophy, stiffness of the collagen layer, thinning of elastin fibers, and decreased mucous production [3, 4, 8, 9]. These alterations may result in poor vocal projection, reduced loudness, hoarseness, change in pitch, decreased phonation time, tremor, breathiness, and decreased accuracy of motor control [3, 4, 6, 8, 10]. Understanding the effects of aging on the systems that support voice and speech is critical to discriminating normal physiology from pathology and allows for studies focused on the effects of concomitant aggravating factors on voice health, such as dehydration, a common condition in older individuals [11, 12].
Maintaining water homeostasis is essential to preserve life, as virtually all biological and chemical reactions in the body require an aqueous medium [13,14,15]. Additionally, water serves as a carrier of nutrients required to maintain viable cells, eliminates toxic compounds from the body, and provides structural support for tissues [14, 15]. Thus, water and electrolyte balance is tightly regulated [14]. Total body water decreases with age [14, 16, 17]; in humans, the water body changes range from 75% in infants to 55% in the elderly [13] and 45% in older women, making older individuals more prone to developing dehydration [14]. There are multiple factors contributing to the development of dehydration as we age [13, 14], including the increased proportion of body fat over lean body mass [14, 17], decreased renal function impairing efficient water conservation [14, 17], decreased sense of thirst [18], and the presence of comorbidities [14, 19]. Although several studies have focused on fluid homeostasis, dehydration, and contributing factors in the aging process, the impact of systemic dehydration on older subjects in the context of vocal fold biology has not yet been investigated.
Vocal fold hydration is a typical recommendation to maintain vocal health; however, many aspects of the underlying mechanisms of its benefits remain unclear. Studies focused on systemic (de)hydration – manipulating water consumption and cellular absorption – in humans have shown that dehydration negatively affects vocal fold physiology and function. However, due to ethical concerns, the use of human subjects is limited to mechanistic studies of laryngeal imaging, acoustic assessments and aerodynamic measures of voice, or postmortem studies [20,21,22,23,24]. In addition, some human studies combined systemic and surface dehydration as challenges, making it difficult to determine the specific biological effects of systemic (de)hydration in the vocal folds [21, 22, 25]. Nevertheless, these studies suggest that dehydration negatively impacts voice production by altering the viscoelasticity properties of the vocal folds. Conversely, animal models allow for ex vivo and in vivo research, using different dehydration protocols. Studies using excised canine larynges have shown that osmotically induced dehydration increases vocal fold stiffness and viscosity [26] and that changes in viscosity affect the vocal fold vibration [27]. Moreover, a decrease in hyaluronan content in the lamina propria was observed in vocal folds from water-restricted rats [28]; hyaluronan is a critical component for maintaining hydration and optimal tissue viscosity [24]. These results support the premise that hydration status alters vocal fold viscoelasticity, as observed in human studies.
In addition to clinical and experimental studies, research focused on the underlying molecular mechanisms of hydration has shown that systemic dehydration affected gene and protein expression in the vocal folds of different animal models [28,29,30]. Transcriptome studies of the vocal folds from furosemide-induced systemically dehydrated rabbits revealed hundreds of differentially expressed genes compared with euhydrated rabbits [29]. In that study, epithelial and cell adhesion-related genes were downregulated, suggesting systemic dehydration may perturb the vocal fold epithelium. Proteomic studies using the same rabbit model revealed dozens of differentially regulated proteins and suggested novel molecular evidence for increased susceptibility to oxidative stress [30] as observed in rats after acute vocal fold injury [31, 32]. Interestingly, treatment with antioxidants prevented age-associated changes in rat vocal folds [9]. Collectively, these findings represent strong evidence that systemic dehydration adversely affects vocal fold biology and that systemic hydration is an effective practice to improve voice quality; however, the interaction between hydration status and aging is an unexplored topic in voice research.
To start filling this gap, our study aimed to identify protein signatures and related mechanisms associated with the intersection of age and hydration status in vocal fold biology. We used comparative label-free quantitative proteomics to analyze the alterations of protein expression patterns in the vocal folds of 6-month-old and > 3-year-old rabbits in response to water ad libitum or water restriction-dehydration protocols. Proteomics-based technology was utilized given its high sensitivity and high-throughput capacity [33]. We used the rabbit as in vivo model given the similarities to human vocal fold anatomy and biology [34, 35]. Furthermore, the rabbit is a well-validated model in studies focused on the vocal folds [29, 30, 36,37,38,39,40,41,42,43,44,45]. Water volume restriction was chosen as the systemic dehydration method to mimic a more physiologically-realistic state that can be better translated to human behavior [46,47,48]. Our comparative proteomics analysis showed a higher number of proteins, related to tissue structure and viscoelastic properties, differentially regulated and with a greater magnitude of change in the vocal folds of older adult rabbits due to systemic dehydration compared to younger rabbits, suggesting a higher vulnerability of aging vocal folds to systemic dehydration.
Results
Water intake
The water intake protocol was divided into two stages: baseline (5 days, day 1 to 5) and water restriction-dehydration period (5 days, day 6 to 10) (Fig. 1). The water intake of euhydrated-control rabbits (young and older adults), which received water ad libitum during the entire experiment, did not change throughout the study, although the baseline water intake was 30% lower for older rabbits than younger rabbits (Fig. 2a). For animals in the dehydrated group, water volume intake was reduced by 50% during the water restriction period compared to the volume consumed in the baseline period, as our water volume restriction protocol intended (Fig. 2b).
Systemic dehydration markers
The body weights of all rabbits were monitored daily during the entire experiment as a marker of systemic dehydration. Systemic dehydration was verified at the end of the water restriction period by body weight loss of − 5.4% in the young-dehydrated rabbits and − 4.6% in the older-dehydrated group, while young and older-control rabbits gained 2.2 and 1.5% body weight, respectively (Fig. 2c–d). The percentage changes in packed cell volume (PCV) and total solids (TS), from baseline to the end of the dehydration period, were compared between control and dehydrated rabbits in both age groups as blood markers of dehydration. PCV and TS did not change significantly between control and dehydrated rabbits in either of the age groups (Fig. 2e–f). Finally, at the end of the experiment, renal renin gene expression was compared between controls and dehydrated rabbits in the young and older adult groups. Renin expression is expected to increase due to systemic dehydration, which was verified for young rabbits (1.7-fold increase) but not for the older adult rabbits compared to their respective controls (Fig. 2g–h).
Comparative proteomics of vocal folds of young and old rabbits under euhydration (control condition) or water restriction-dehydration
Label-free quantitative (LFQ) proteomics was used to determine the proteome changes related to age, hydration status, and interaction between age and hydration in rabbit vocal folds. Proteins identified at 1% FDR, with at least one unique peptide and MS/MS counts ≥2, were considered true identifications. These data were uploaded to the Perseus platform, and after removal of “contaminants”, “reverse”, and “only identified by site”, a total of 3667 proteins were identified. Then, the data were Log2 transformed and filtered to keep only proteins with a minimum of three valid LFQ (non-zero) values in at least one of the four groups: young-control (YC), young-dehydrated (YD), older-control (OC), or older-dehydrated (OD) rabbits. After filtering and imputation of values based on normal distribution, 2286 proteins were identified and used in further analyses. Two-way ANOVA identified 177, 169, and 81 proteins significantly (p ≤ 0.05) affected by age, hydration status, or interaction of both factors, respectively (Additional file 1: Table S1). Principal component analysis (PCA) was conducted in parallel with ANOVA to further characterize the collective contribution of proteins to the factors investigated. We considered the correlation of proteins affected by age and hydration status to principal component 1 (PC1) as an alternative to adjusted statistical significance as a criterion in defining proteins of interest (Additional file 2: Tables S2). These subsets of proteins with p ≤ 0.05 were then used for heat maps to show protein expression patterns across groups, and pathway and process enrichment analysis to define functional clusters. Analysis of the set of proteins that showed a significant interaction effect was performed to identify the pattern of changes in young and older rabbits’ vocal folds due to systemic dehydration. Results from these analyses are described below.
Principal component analysis (PCA) based on the vocal fold proteome was performed for each main effect of age and hydration status, and interaction effect. Different subsets of proteins representing the top 500 (for age and hydration status), 250 (for interaction), and all significant proteins combined (p ≤ 0.05) arranged by ascending p-values for the respective factor were used to identify the separation between groups and the overall variance collectively provided by these proteins. Young (YC, YD) and older (OC, OD) age groups had a linear separation with the top 500 proteins (Fig. 3a), explaining 32.8% of the overall variance (PC1 and PC2 combined). The separation between hydration status groups (control: YC, OC vs. dehydrated: YD, OD) was observed for the top 500 proteins providing an overall explained variance of 31.5% (Fig. 3b). The distinction of groups considering the interaction effect was observed with the top 250 proteins explaining an overall variance of 31.8% (Fig. 3c). Finally, all significant proteins (p ≤ 0.05) identified for age, hydration status, and interaction effects provided an overall explained variance of 33%, with clear differentiation between the four groups and a more pronounced distinction of the older-dehydrated from the other groups (Fig. 3d). The focus of this study are the proteins with p ≤ 0.05 affected by interaction and those affected by age and hydration status that showed moderate to high correlation (≥ 0.6) with the respective first principal component indicating their contribution to the separation of biological groups (Additional file 2: Tables S2).
Hierarchical clustering of the significant proteins (p ≤ 0.05) into heat maps shows the pattern of up or downregulation of each protein subset by age, hydration status, and interaction (Fig. 4a–c). It is noteworthy that proteins affected by interaction show opposite expression patterns between young and older adult rabbits in response to dehydration (YD x OD) (Fig. 4c). A Venn diagram shows the overlap between the significant proteins affected by each individual factor and by the interaction (Fig. 4d). The majority of proteins are unique to each subset: age (142 proteins), hydration status (134 proteins), and interaction (49 proteins). Twenty-three common proteins were significantly affected by both main factors of age and hydration status. Twenty common proteins were significantly affected by age and interaction, and 20 were commonly found differentially regulated by hydration status and interaction. Only eight proteins were commonly found across the three subsets. The sets of unique proteins affected by only one factor or by the interaction highlights changes in protein expression patterns that are age- or hydration-dependent, with 49 proteins differentially affected in young and older adult rabbit vocal folds depending on the hydration level.
Enrichment analysis for biological pathways and processes of significant proteins affected by age, hydration status, or interaction
Enrichment analysis was performed using Metascape [49], and enriched terms were based on Gene Ontology Biological Processes (BP). The top 20 enriched terms for age, hydration status, and interaction are shown in Fig. 5a, c, and e, respectively. The Molecular Complex Detection (MCODE) was applied to identify functional clusters based on protein-protein interaction networks (Fig. 5b,d,f). Age-related top enrichment BP terms are associated with “muscle contraction”, “cell redox homeostasis”, “protein degradation”, and “extracellular matrix organization”, among others (Fig. 5a). The protein-protein interaction (PPI) based on the MCODE algorithm identified five functional clusters shown in Fig. 5b. These functional clusters are formed by networks of proteins related to “extracellular matrix organization”, “muscle contraction”, “cell proliferation”, “cellular oxidant detoxification”, and “carboxylic acid catabolism”, which highlights the effect of age on the vocal fold tissue structure and oxidative stress. Protein members of these functional categories include: The top BP enriched terms related to the proteins affected by hydration status include “protein degradation”, “muscle structure development and system process”, “translation”, “regulation of DNA biosynthesis”, and interestingly, “aging” (Fig. 5c). Five PPI MCODE clusters were identified with the largest network (25 proteins) associated with “protein degradation”, followed by “endocytosis”, “translation”, “muscle structure,” and “extracellular matrix organization” (Fig. 5d). These functional clusters suggest an effect of hydration/dehydration on the vocal fold structure and fundamental cellular processes. Among the top 20 enriched terms for the vocal fold proteins with an interaction effect, “striated muscle contraction”, “positive regulation of sodium ion transmembrane transport”, “extracellular matrix organization”, “tissue homeostasis”, and “aging” are listed (Fig. 5e). The MCODE clusters are associated with “muscle contraction” and “translation” process in the vocal fold tissue (Fig. 5f).
The protein members of the combined enriched GO BP terms forming each functional cluster related to age, hydration, and interaction are listed in Table 1.
Proteins differentially regulated in the rabbit vocal folds due to age, hydration status, and interaction
The age-related protein changes, i.e., proteins differentially regulated by the main effect of age, are associated with downregulation of muscle and extracellular matrix structural proteins in the vocal folds of older adult rabbits compared to young, such as Myosin light chains 2, 3, 4 and 6B, Myosin heavy chains 6, 7 and 10, Troponin I, Fibulin 5, Collagens type XIV and XVIII alpha 1 chain, Serpin H1, and Versican. Desmin, a muscle contraction-related protein, was upregulated in older vocal folds compared to young. Changes related to hydration status resulted in upregulation of proteins involved in muscle function such as Myosin light chain 2 and Tropomyosin alpha-1 chain, and proteins of the extracellular matrix structure including Fibulins 1 and 5, Collagen type XVIII alpha 1 chain, Fibrillin 1, Microfibril associated proteins 4 and 5, in the presence of dehydration. The transcription factor NIF3-like protein 1, with no function previously described in vocal folds, was also upregulated by dehydration. Upregulation of the intermediate filament Lamin A/C, involved in muscle structure development and related to tissue stiffness, was observed as an effect of dehydration, although the magnitude of increase was below the mean difference cut-off used in our study. Among downregulated proteins affected by dehydration, Sarcolemmal membrane-associated protein involved with muscle contraction was identified. The changes in the protein expression patterns related to an interaction effect were analyzed, and the number of proteins upregulated, downregulated, or unchanged (below the cut-off of 1.3-fold change; i.e., |Log2 mean difference| = 0.38) in young and older adult rabbits due to systemic dehydration are shown in Fig. 6a. The fold change cut-off of 1.3 was based on previous studies and considerable changes for large structural proteins [50,51,52]. Amongst the 81 proteins differentially regulated by an interaction effect, 32 showed an opposite pattern of change in the young and older rabbit vocal folds in response to systemic dehydration, with some proteins considered unchanged in one of the age groups (Additional file 3: Table S3). Seventeen proteins were downregulated in older adults and upregulated in young rabbits, while 15 were upregulated in older adults and downregulated in young rabbits. These opposite change patterns were seen for Myosin light chains 2 and 3, and Troponin I, which were increased by 18.2, 9.6, and 4.1-fold, respectively, in older, and decreased by − 1.5, − 1.8, and − 2-fold, respectively, in young rabbit vocal folds due to dehydration (Fig. 6b). Forty-one proteins presented a change in the expression pattern due to dehydration in only one of the age groups. Thirty-one proteins changed solely in older adult rabbits, while ten changed in the young rabbits alone (Additional file 3: Table S3). Fibulin 5, Collagen type XVIII alpha 1 chain, and Laminin subunit gamma 1 were upregulated by 5, 4.3, and 1.5-fold, respectively, in older adult rabbits due to dehydration, while they remained unchanged in young rabbits (fold change ≤1.3) (Fig. 6c). Proteins involved in translation as the Eukaryotic translation initiation factor 4B and Ribosomal protein S15 were downregulated by − 3.9 and − 2.5-fold, respectively, in older adult rabbits and unchanged in young rabbits, while the Small nuclear ribonucleoprotein U1 was downregulated by − 2.3-fold in older adult and increased by 1.3-fold in young rabbits (Fig. 6d). Interestingly, the overall magnitude of change of differentially regulated proteins was higher in the older adult rabbits compared to young, suggesting that older vocal folds are more susceptible to protein changes caused by systemic dehydration. In addition, a higher number of proteins were considered unchanged (below 1.3-fold change) in the young rabbits compared to 18 in the older rabbits, indicating that more proteins are affected by systemic dehydration in aging vocal folds (Additional file 3: Table S3).
RT-qPCR of selected genes based on differentially regulated proteins
The genes tested by RT-qPCR were chosen based on a subset of proteins identified as differentially regulated by either hydration status or interaction effect. These genes were also selected based on their functional enrichment terms related to extracellular matrix and muscle functions. The transcription factor, NIF3L1, involved in neuron differentiation that was upregulated by dehydration, was also tested. Eleven genes were tested; five based on proteins differentially regulated by hydration status (FBN1, MFAP5, SLMAP, LMNA, NIF3L1) (Fig. 7a-e), and six selected from the interaction (FBLN5, LAMC1, COL18A1, MFAP4, MYL2, DMD) (Fig. 8a-f) results. Interestingly, the expression of FBN1, LMNA, and NIF3L1 genes was significantly affected by hydration status (p ≤ 0.05) like the respective proteins, but while the proteins were upregulated, the genes were downregulated by systemic dehydration (Fig. 7a, d, and e, respectively). MFAP5 and SLMAP gene expression (Fig. 7b and c, respectively) was not significantly affected by hydration status. FBLN5 gene was significantly affected by age (p ≤ 0.05) and was decreased in older rabbit vocal folds (Fig. 8a), matching the protein expression pattern. MYL2 gene was not significantly affected by interaction (Fig. 8e) despite having a similar pattern of expression observed for the protein (Fig. 6b). None of the genes tested were significantly affected by interaction, unlike the correspondent proteins.
Discussion
This study aimed to investigate the proteomic changes in the vocal folds of young and older adult rabbits subjected to the same protocol of water restriction-dehydration: five days of 50% water volume restriction. A water restriction protocol was chosen over pharmacologically-induced dehydration or complete water deprivation to provide a more physiologically-relevant translation to human behavior. A significant proportion of the population of adults and older adults consume about half of the water intake recommended by various governmental institutes worldwide [46,47,48, 53, 54]. Additionally, we sought to achieve only mild to moderate dehydration to provide further physiological validity to the model. After water restriction, young and older adult rabbits showed body weight losses of − 5.4% and − 4.6%, respectively, consistent with mild to moderate systemic dehydration [55, 56].
LC-MS/MS analysis quantified 2286 proteins in the vocal folds of young and older adult rabbits combined after applying the filtering criteria for this study. These proteins include structural components of the extracellular matrix such as collagens, fibrillin 1 and 2, fibulins, elastins, laminins, proteoglycans, decorin and fibronectin; and muscle proteins such as myosins and actins. These proteins are also present in human vocal folds [57], corroborating our protein extraction and fractionation methods and the LC-MS/MS approach used to analyze the proteome of rabbit vocal folds.
Subsets of proteins were identified as significantly affected by age (young vs. older adult rabbits), hydration status (water ad lib-control vs. water restriction-dehydrated), or interaction of these factors. Interestingly, the majority of proteins affected by the individual factors or the interaction are unique, with only a few common proteins between subsets. These results indicate that the vocal folds have independent mechanisms of protein regulation related to the aging process, response to hydration status, and a combination of age and level of hydration (interaction). Despite the identification of unique proteins regulated by the individual effects considered in this study, the enrichment pathway analysis showed a common association of these proteins to extracellular matrix organization and muscle structure/contraction, suggesting that age and hydration, as single or interacting factors, interfere with the physiological biomechanical and viscoelastic properties of the vocal folds responsible for structural support and voice production [57,58,59,60,61,62].
Our study revealed several changes related to muscle proteins, especially in older adult rabbits. When comparing the effect of age on vocal fold protein regulation, the downregulation of essential muscle proteins, such as Myosin heavy and light chains and Troponin I, was observed in older adult rabbits. The Myosin heavy chains 6, 7 (β-subunit) and 10, and the Myosin light chains 2, 3, and 4 were downregulated in older rabbit vocal folds, an effect of age alone. Myosin heavy chain β-subunit is expressed in muscle, rich in slow-twitch type 1 fibers; the Myosin light chains downregulated in older rabbit vocal folds are responsible for modulatory and regulatory functions of muscle conformational changes and contraction [63]. Likewise, Troponin I regulates contraction by modulating myofibril’s calcium sensitivity, causing relaxation. The downregulation of these proteins corroborates the presence of muscle decline expected in aged individuals [8]. Studies in humans have shown the decreased synthesis of muscle proteins in middle-aged and older adults compared to their young counterparts [64]. Specifically, Myosin heavy chain decreased with age. Myosin is crucial for muscle contraction by hydrolyzing ATP to generate muscle contraction [63], and its decrease may also contribute to the loss of muscle mass and muscle strength. An opposite interaction effect was observed for Myosin light chains 2, 3, and 6B, which were upregulated in older adult rabbits and downregulated in young rabbits in response to dehydration, suggesting that dehydration also plays a role in the expression of these proteins. There is a direct connection between the effect of hydration/dehydration in muscle cells; water mediates interactions of myosin and actin filaments, thus, interfering with muscle contraction [65,66,67]. Interestingly, the intermediate muscle filament Desmin was upregulated in older in comparison to young rabbit vocal folds, indicating an effect of aging. Desmin has an important role in cellular homeostasis in response to chronic stress or muscle atrophy [68, 69], and increased Desmin expression has been observed in the muscle of aging rats [70]. However, the impact of these changes have yet to be investigated in the context of the vocal folds and voice function. Nevertheless, these data provide evidence of age alone and interactions between age and hydration status in the laryngeal or vocal fold muscle.
Various extracellular matrix-related proteins were differentially regulated herein. Extracellular matrix proteins provide structural support, mechanical properties, and strength of tissues, including vocal folds, playing a pivotal role in phonation [62, 71, 72]. Collagens XIV, XVIII, and Fibulin 5 were downregulated in older rabbit vocal folds compared to young tissue. To our knowledge, these specific collagen types have not been investigated in depth in vocal fold tissue; however, studies suggest that the changes in the collagen fiber density and arrangements within the lamina propria may affect phonation [73, 74]. Collagen type IV is exclusive to extracellular matrix basal membranes [75] and is present in the human vocal fold basal membrane providing support to epithelial and endothelial cells [76]. Collagen type IV was upregulated in older rabbit vocal folds compared to young, an effect of aging observed in our study. The relationship between Collagen type IV and aging is not well established. Increased accumulation of Collagen type IV is reported in the basal lamina of cerebral microvessels in humans [77] but decreased in the skin of older adults [78]. Conversely, several extracellular matrix proteins were upregulated, including Collagen type XVIII and Fibulin 5, in the presence of dehydration when observing the effect of hydration status alone. These protein changes may be related to the remodeling of the extracellular matrix [79] in response to dehydration. Moreover, the accumulation of collagens and the decrease of elastins may result in extracellular matrix stiffness in aging larynx and other organs [59, 79]. Finally, Lamin A was upregulated by dehydration, by a smaller magnitude, especially when observing the mean difference within the young groups. Previous data has identified that Lamin proteins A and C are important for imparting the nucleus with its stiffness, and their expression has been reported to scale with tissue stiffness [80]. Thus, upregulation of this protein due to dehydration may be related to tissue stiffness in the vocal fold of rabbits.
We tested the relative expression of a small subset of genes based on the proteins differentially regulated in the vocal folds. Four out of the 11 (36%) tested genes followed the direction of the expression pattern of the respective proteins; however, their relative expression was not statistically significant as their proteins. The correlation between gene and protein expression can be as little as 40% depending on the system being investigated [81]. Multiple factors contribute to the inconsistency between mRNA and protein abundance for a given gene, such as mRNA and protein stability, lower rate of mRNA transcription compared to protein translation, and post-translational modifications [81, 82]. Ideally, comparing the complete proteome to the transcriptome is a more desirable approach when the study design and biological replicates allow it.
The number of proteins identified in this study represents an increment of 602 vocal fold proteins compared to our previous study using 6-month-old rabbits (young) subjected to systemic dehydration with furosemide, which identified 1684 proteins [30]. This difference can be due to the different number of biological replicates used in each study, eight rabbits in the furosemide study versus 22 rabbits in the present study, and the filtering criteria applied in each study. Regardless of the differences between studies, 1307 proteins were detected in both, including main extracellular matrix and muscle components and cellular metabolism related-proteins. While this study identified a higher number of vocal fold proteins compared to our previous study [30], the number of differentially regulated proteins due to age, dehydration, or interaction was still relatively low. This might suggest that many vocal fold proteins are regulated at the post-translational level, and measurement of total proteome may only confer a partial view of the overall changes in vocal fold biology. Therefore, future studies that focus on identifying site-specific changes in protein phosphorylation and other modifications may provide a deeper and complete understanding of the regulation of vocal fold proteome under different conditions such as age or hydration status. Our study shows that proteomic analysis can provide critical information and contribute to our understanding of vocal fold biology under different physiological conditions.
The life expectancy of White New Zealand rabbits may reach ten years with proper care, spaying and neutering, but rabbits in research settings have a lower life expectancy due to stress [83]. The young and older rabbits used herein correspond to sexually mature young adults and retired breeders representing older adults, respectively. Studies including elderly rabbits (> 6 years of age) may provide additional age-related protein expression changes. However, the effects of aging were observed in this study in the older rabbits group. Interestingly, as seen in humans [14, 17], the older rabbits in this study had a lower water consumption than young rabbits when provided with water ad libitum. Among the indicators of dehydration, packed cell volume (PCV) and total solids (TS) did not change significantly after water restriction; however, the average percentage change of TS in the dehydrated groups is higher than those for the control groups, but the high variability within groups likely prevents a significant statistical difference. Given the lack of sensitivity of these indicators of dehydration, our group has developed a method to quantify renin gene expression in the kidneys [28]. Renin is released by juxtaglomerular cells in the kidney in response to decreased blood pressure, and its transcription is expected to be increased in the case of dehydration due to decreased blood volume [84]. In this study, there was a significant increase in renin expression in young rabbits subjected to water restriction compared with same-age rabbits receiving water ad libitum; conversely, there was no significant difference in renin expression in older adult rabbits between treatments. Interestingly, the impact of age on the kidneys, including a decrease in renin activity, is well documented in aging humans [13, 16, 85]; these factors contribute to the increased susceptibility to dehydration as we age. Thus, it is possible that renin expression was affected by the age in the older rabbits in this study. Nevertheless, the rabbits in both groups were systemically dehydrated based on imposed water restriction and decreased body weight, which is commonly used as a consistent indicator of dehydration in humans [86,87,88] and animal models [28, 30, 89,90,91,92]. Clinical translation to human subjects is preliminary and awaits further investigation. However, given the similarity in rabbit laryngeal anatomy to humans [34, 35], it is appropriate to posit that systemic dehydration in older speakers, as induced by reduced water intake, may impact the extracellular matrix which scaffolds the vocal folds. This structural impact could negatively prime these tissues for trauma from mechanical vibration. The similarities between older adult rabbits and aging humans regarding decreased water intake and renin expression serve as additional support for using the rabbit as an animal model of dehydration.
Conclusions
Our study successfully used the rabbit as an in vivo model to investigate the impact of aging and hydration status on the proteome of vocal fold tissue. The proteomic data generated in this study provide a robust resource for the laryngological research community. Our findings suggest changes in the relative abundance of extracellular matrix components and muscle-related proteins in young and older rabbit vocal folds with a higher impact of systemic dehydration on the older tissue. The extracellular matrix is directly related to the viscoelastic properties and vibration of the vocal folds, while the muscle fibers scaffold the tissue and are also involved in phonation. The clinical impact of the magnitude of protein changes observed in the rabbit vocal folds, the translation to human vocal folds, and implications for vocal output await further investigation.
Methods
Rabbits and water restriction protocol
Experiments were approved by the Purdue Animal Care and Use Committee (protocol #1606001428) and conducted following the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals [93]. Twenty-two male New Zealand White rabbits (Envigo, Indianapolis, IN, USA) were used in the study. The rabbits were acquired as two age groups: the young group comprised 6.5-month-old rabbits (n = 11), and the older group comprised rabbits between 3.2 and 3.7 years of age (average = 3.3 years old) (n = 11). The 6.5-month-old rabbits are sexually mature young adult rabbits, and the ~ 3.5-year-old rabbits are retired breeders and represent older adults [94]. Male rabbits were used to avoid the effects of female sex hormones on body fluid regulation that occur in hormonally-cycling humans and rodents [95,96,97]. Animals were acclimatized for one week prior to the experiment. After acclimation, rabbits from both age groups were randomly assigned to either a water ad libitum-control (n = 5 young, and n = 5 older adults) or water restriction-dehydrated (n = 6 young, and n = 6 older adults) group. The experiment consisted of ten days divided into two periods: the baseline period, when all rabbits received water ad libitum and had individual water intake measured daily for 5 days, and the water restriction period. Rabbits assigned to the water restriction-dehydrated group received half volume of their baseline-average water intake, and control rabbits had free access to water for 5 days (Fig. 1). Body weights and water intake were recorded daily for all rabbits throughout the experiment. Whole blood was collected from all rabbits via ear vein puncture into heparinized microhematocrit capillary tubes to assess the packed cell volume (PCV) and total solids (TS) on day 5 of baseline and day 10 of water restriction. All rabbits had unrestricted access to food and hay during the entire study. At the end of the water restriction period, rabbits were euthanized using 1.0 mL IV Beuthanasia-D Special (Schering Plough Animal Health Corp., Union, NJ, USA).
Systemic dehydration markers
Percentage changes in body weight, packed cell volume (PCV), and total solids (TS) were assessed as markers of systemic dehydration: [final value after water restriction dehydration (day 10) – baseline value (day 5)]/baseline value*100. Additionally, the relative expression of renin gene in the kidney was compared between water ad lib-control and water restriction-dehydrated rabbits in the young and older adult groups.
Tissue collection
Following euthanasia, kidney cortex tissue and larynges from all rabbits were collected. The kidney tissue specimens were immediately flash frozen in liquid nitrogen, and stored at − 80 °C until processing. Immediately after collecting the larynx, the vocal fold tissue was microdissected under a dissection microscope and flash frozen as described before [30]. The full thickness of the vocal fold tissue, identified as the soft tissue at the level of the arytenoid cartilages at the transverse level of the glottis, was collected bilaterally (left and right vocal folds). Each bilateral section included both mucosa and thyroarytenoid muscle. The left vocal fold was processed for proteomic analysis and the right vocal fold for RT-qPCR.
Vocal fold protein extraction and preparation for mass spectrometry
Protein extraction and mass spectrometry analysis were performed at the Purdue Proteomics Facility as previously described [30, 98]. Briefly, one vocal fold specimen per rabbit (n = 22) was thawed on ice and washed twice with 25 mM ammonium bicarbonate (ABC) before centrifugation at 6500 rpm for 90 sec. Next, tissues were lysed with 100 mM ABC and 1.4 mm ceramic beads and homogenized using a Precellys tissue homogenizer (Bertin Technologies, Rockville, MD, USA). Homogenates were transferred to new tubes and protein concentration was measured using the bicinchoninic acid assay (BCA) [99].
To increase proteome coverage, protein fractionation was performed for each sample, and the soluble and insoluble fractions were analyzed independently in the mass spectrometer. This procedure was not applied to define the subcellular localization of proteins. The fractionation protocol is described elsewhere [30]. Briefly, fractionation was performed using 50 μg of proteins per sample across all 22 samples. Samples were centrifuged in an Optima MAX-XP ultracentrifuge (Beckman Coulter, Indianapolis, IN, USA) at 55 K rpm for 40 min, and the supernatant containing the soluble fraction was recovered and mixed with 4 volumes of cold acetone and incubated overnight at − 20 °C for protein precipitation. The pellets containing the insoluble fractions were kept at − 20 °C. Following precipitation, soluble fractions were centrifuged at 13,500 rpm for 15 min at 4 °C and pellets were collected and dried in a vaccum centrifuge. Pellets from both soluble and insoluble fractions were incubated with 10 mM dithiothreitol (DTT) and 8 M urea in 25 mM ABC for 1 h at 37 °C for reduction, and 2% iodoethanol, 0.5% triethylphosphine in acetonitrile (ACN) for 1 h at 37 °C for alkylation. Next, samples were dried in a vacuum centrifuge.
Proteins in each sample (22 soluble and 22 insoluble) were digested with trypsin and Lys-C mix (Promega, Madison, WI, USA) at a 1:25 enzyme to substrate ratio, at high pressure using a NEP2320 Barocycler (Pressure Biosciences, South Easton, MA, USA) at 50 °C with 60 cycles of 50 seconds at 20 kpsi, and 10 seconds at atmospheric pressure. MicroSpin C18 silica columns (The Nest Group, Inc., Southborough, MA, USA) were used to desalt the digested peptides and elute with 80% ACN and 0.1% formic acid (FA). Purified peptides were dried in a vacuum centrifuge and stored at − 80 °C until LC-MS/MS. Both soluble and insoluble peptides were resuspended in 3% ACN, 0.1% FA in water at a final concentration of 1 μg/μL, and 1 μL was loaded to the HPLC system.
LC-MS/MS data analysis
LC-MS/MS data were collected at the Purdue Proteomics Facility. Peptides were analyzed in a Dionex UltiMate™ 3000 RSLCnano System (Thermo Fisher Scientific, Waltham, MA, USA) coupled online to an ETD-enabled Orbitrap Fusion Lumos Mass Spectrometer (Thermo Fisher Scientific) as previously described [30, 98]. Briefly, reverse phase peptide separation was accomplished using a trap column (300 μm ID × 5 mm) packed with 5 μm 100 Å PepMap C18 medium coupled to a 50-cm long × 75 μm inner diameter analytical column packed with 2 μm 100 Å PepMap C18 silica (Thermo Fisher Scientific). The column temperature was maintained at 50 °C. Samples were loaded to the trap column in a loading buffer (2% ACN, 0.1% FA) at a flow rate of 5 μL/min for 5 min and eluted from the analytical column at a flow rate of 150 nL/min using a 280-min LC gradient as follows: linear gradient of 8 to 30% of solvent B in 225 min, 30–60% of B in next 35 min, before reverting back to 2% of B, and hold at 2% of B for next 15 min for column equilibration. The column was further washed and equilibrated by using three 30-min LC gradients before injecting the following sample. All data were acquired in the Orbitrap Fusion Lumos Mass Spectrometer (Thermo Fisher Scientific). Orbitrap mass analyzer and data were collected using a higher-energy collisional dissociation (HCD) fragmentation scheme. For MS scans, the scan range was from 375 to 1500 m/z at a resolution of 120,000, the automatic gain control (AGC) target was set at 4 × 105, maximum injection time 50 ms, dynamic exclusion 60s, and intensity threshold 5.0 × 103. MS data were acquired in data-dependent mode with a cycle time of 3 s/scan. MS/MS data were collected at a resolution of 7500. The mass spectrometry data underlying this study were deposited in the Center for Computational Mass Spectrometry database under MassIVE ID: MSV000089683 (https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp).
LC-MS/MS data were analyzed using MaxQuant software (version 1.6.3.3) [100] against rabbit (Oryctolagus cuniculus) database downloaded from the UniProt (www.uniprot.org) containing 41,458 protein sequences. Soluble and insoluble fractions acquired data were merged for each correspondent sample before database searches. The following parameters were edited for the searches: precursor mass tolerance of 10 ppm; enzyme specificity of trypsin/Lys-C enzyme allowing up to 2 missed cleavages; oxidation of methionine (M) as a variable modification and iodoethanol (C) as a fixed modification. False discovery rate (FDR) of peptide spectral match (PSM) and protein identification was set to 1%. The unique plus razor peptides (non-redundant, non-unique peptides assigned to the protein group with most other peptides) were used for peptide quantitation. Only proteins detected with at least one unique peptide and MS/MS ≥ 2 (spectral counts) were considered valid identifications. Label-free quantitation (LFQ) intensity values were used for relative protein abundance comparisons.
Statistical analysis for comparative proteomics
Perseus software platform (v1.6.12.0) [101] was used to identify vocal fold proteins differentially regulated between the four groups: young-control (YC), young-dehydrated (YD), older-control (OC), and older-dehydrated (OD). First, MaxQuant data was uploaded to Perseus and “contaminants”, “reverse”, and “only identified by site” proteins were removed. Next, LFQ values were Log2 transformed, and data was filtered to have a minimum of three LFQ valid values (non-zero LFQ values in at least three replicates) in at least one group. After filtering, imputation for missing values was performed sample-wise based on normal distribution. Statistical analysis was completed with Two-way ANOVA to compare the LFQ means of each protein between groups for each given factor and interaction. Significant results were considered for the main effect of age (young vs. older adults; regardless of hydration status), the main effect of hydration status (control vs. dehydration; regardless of age), and the interaction of both factors. Proteins with a p-value ≤0.05 were considered significant between age groups, hydration status or interaction. Proteins with a |Log2(LFQ mean difference)| ≥ 0.38 were considered differentially regulated [30, 50].
Principal component analysis
Principal component analysis (PCA) was completed for each set of Two-way ANOVA results based on age, hydration status, and interaction. The three sets share the same list of proteins, but with different p-values depending on the main effect or interaction analyzed (Additional file 1: Table S1). Common proteins between age and interaction, and between hydration status and interaction were removed before the PCA since an interaction effect precludes interpretation of the main effects directly. Proteins in each list were organized by ascending p-value and PCA plots were created based on the top 500 proteins (for age and hydration status), top 250 proteins (for interaction), and all significant proteins (p ≤ 0.05) for each main effect and interaction combined. Correlations between proteins with p ≤ 0.05 for age and hydration status with principal component 1 (PC1) were calculated. PCA and correlation analysis were completed with R (v4.1.2) using RStudioTM version 1.4.1717 (RStudio, PBC, Boston, MA, http://www.rstudio.com).
Pathway and process enrichment analysis of significant proteins
Three enrichment analyses were performed for the protein subsets identified in the vocal folds as significantly affected by age, hydration status, or interaction of both. Uniprot IDs from these significant subsets of proteins were converted into gene symbols using the UniProt Retrieve/ID mapping tool (https://www.uniprot.org/uploadlists/), which were used as the input for the enrichment analysis with Metascape [49] (https://www.metascape.org). For each given subset of genes, pathway and process enrichment analysis was conducted selecting GO Biological Processes as the ontology source. Enrichment analysis default parameters were applied: minimum overlap of 3 gene counts, p-value cut-off of 0.01, and enrichment factor > 1.5 (ratio between the observed counts and the counts expected by chance). All genes in the genome were used as the enrichment background. Enriched ontology terms are grouped into clusters based on membership similarities with the most statistically significant term within a cluster representing that cluster. Metascape calculates p-values based on the cumulative hypergeometric distribution [102] and q-values using the Benjamini-Hochberg procedure to account for multiple testing [103]. Kappa scores [104] are used as the similarity metric for hierarchical clustering of the enriched terms, and clusters are determined based on sub-trees with a similarity > 0.3.
Protein-protein interaction enrichment analysis based on significant proteins
Protein-protein interaction (PPI) enrichment analysis was carried out for each subset of genes affected by age, hydration status, or interaction by applying the databases STRING [105], BioGrid [106], OmniPath [107], InWeb_IM [108]. Protein networks were created and the Molecular Complex Detection (MCODE) algorithm [109] was applied to identify densely connected network components and associated functional descriptions. Pathway and process enrichment analysis was applied to each MCODE component independently, and the three best-scoring terms by p-value were retained as the functional description of the corresponding components, which are shown along with networks as a color-coded list.
RNA extraction and RT-qPCR
Total RNA was extracted from renal cortex (n = 22) and vocal fold tissue (n = 22; one vocal fold per rabbit) using RNeasy Plus Mini Kit and RNeasy Fibrous Tissue Mini Kit (QIAGEN, Hilden, Germany) following manufacturer’s protocol, respectively. Concentration and purity of RNA was assessed by spectrophotometry (NanoDrop™, Thermo Fisher Scientific). The cDNA synthesis was performed with 1.0 μg and 250 ng of total RNA extracted from renal cortex and vocal fold tissue, respectively, using SuperScript™ IV VILO™ Master Mix (Invitrogen™ by Thermo Fisher Scientific). Target vocal fold genes were selected based on proteins differentially regulated identified in the proteome associated with gene ontology enrichment terms (based on the enrichment analysis described above) related to extracellular matrix and muscle functions. The primers for qPCR were designed using Primer-BLAST tool [110], and at least one primer of each pair span an exon-exon junction to prevent residual genomic DNA amplification (Table 2). Each sample was run in triplicates for each gene in 96-well PCR plates with Power SYBR® Green PCR Master Mix (Applied Biosystems™ by Thermo Fisher Scientific), and a final concentration of 100 nM of each primer per reaction. Reactions were carried out in a QuantStudio™ 3 Real-Time PCR System (Applied Biosystems™) with the thermal cycling conditions of 95 °C for 10 min; 40 cycles of 95 °C for 15 sec, and 60 °C for 1 min; and melt curve stage of 95 °C for 15 sec, 60 °C for 1 min, 95 °C for 1 sec. The relative expression level of kidney renin was analyzed with the 2−ΔΔCT method [111] using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as normalizer, and the relative expression of the vocal fold target genes was calculated with 2−ΔCT method [112] using hypoxanthine phosphoribosyl-transferase 1 (HPRT1) as normalizer.
Statistical analysis
The data for baseline water intake, % changes in body weight, PCV and TS, and kidney renin gene expression were tested for normal distribution with Shapiro-Wilk test (α = 0.05). For comparisons between two groups, unpaired t-test was applied when data passed the normality test, and Mann-Whitney test, when data did not pass normality. Two-way ANOVA was applied for comparisons with two factors, such as water intake between the baseline and water restriction period (Two-way repeated-measures ANOVA followed by Sidak’s test), and relative gene expression of vocal fold genes between rabbits (ordinary Two-way ANOVA). GraphPad Prism 9.3.1 software was used for the statistical analyses and graphics. Significance was defined at α = 0.05.
Availability of data and materials
The proteomic data generated and analyzed during the current study are available in the MassIVE repository under ID MSV000089683:
https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=fff675b02b354a669eb130c4d6923962
References
Vaca M, Mora E, Cobeta I. The aging voice: influence of respiratory and laryngeal changes. Otolaryngol Neck Surg. 2015;153(3):409–13.
Santos AO, Godoy J, Silverio K, Brasolotto A. Vocal changes of men and women from different age decades: an analysis from 30 years of age. J Voice [Internet]. 2021; [cited 2022 May 30];0(0). Available from: https://www.jvoice.org/article/S0892-1997(21)00185-5/fulltext.
Torre P, Barlow JA. Age-related changes in acoustic characteristics of adult speech. J Commun Disord. 2009;42(5):324–33.
Kendall K. Presbyphonia: a review. Curr Opin Otolaryngol Head Neck Surg. 2007;15(3):137–40.
Gorham-Rowan MM, Laures-Gore J. Acoustic-perceptual correlates of voice quality in elderly men and women. J Commun Disord. 2006;39(3):171–84.
Pessin ABB, Tavares ELM, Gramuglia ACJ, De Carvalho LR, Martins RHG. Voice and ageing: clinical, endoscopic and acoustic investigation. Clin Otolaryngol. 2017;42(2):330–5.
Taylor S, Dromey C, Nissen SL, Tanner K, Eggett D, Corbin-Lewis K. Age-related changes in speech and voice: spectral and cepstral measures. J Speech Lang Hear Res. 2020;63(3):647–60.
Martins RHG, Benito Pessin AB, Nassib DJ, Branco A, Rodrigues SA, Matheus SMM. Aging voice and the laryngeal muscle atrophy: vocal muscle atrophy in the aging voice. Laryngoscope. 2015;125(11):2518–21.
Mizuta M, Hirano S, Hiwatashi N, Kobayashi T, Tateya I, Kanemaru S, et al. Effect of AST on age-associated changes of vocal folds in a rat model. Laryngoscope. 2014;124(10):E411–7.
Rapoport SK, Menier J, Grant N. Voice changes in the elderly. Otolaryngol Clin N Am. 2018;51(4):759–68.
Xiao H, Barber J, Campbell ES. Economic burden of dehydration among hospitalized elderly patients. Am J Health-Syst Pharm AJHP Off J Am Soc Health-Syst Pharm. 2004;61(23):2534–40.
Bryant H. Dehydration in older people: assessment and management: Hannah Bryant relates the importance of assessing and managing dehydration in older people. Emerg Nurse. 2007;15(4):22–7.
Casado Á, Ramos P, Rodríguez J, Moreno N, Gil P. Types and characteristics of drinking water for hydration in the elderly. Crit Rev Food Sci Nutr. 2015;55(12):1633–41.
Sheehy CM, Perry PA, Cromwell SL. Dehydration: biological considerations, age-related changes, and risk factors in older adults. Biol Res Nurs. 1999;1(1):30–7.
Armstrong-Esther CA, Browne KD, Armstrong-Esther DC, Sander L. The institutionalized elderly: dry to the bone! Int J Nurs Stud. 1996;33(6):619–28.
Luckey AE, Parsa CJ. Fluid and electrolytes in the aged. Arch Surg Chic Ill. 2003;138(10):1055–60.
Benelam B, Wyness L. Hydration and health: a review. Nutr Bull. 2010;35(1):3–25.
Phillips PA, Rolls BJ, Ledingham JG, Forsling ML, Morton JJ, Crowe MJ, et al. Reduced thirst after water deprivation in healthy elderly men. N Engl J Med. 1984;311(12):753–9.
Masot O, Miranda J, Santamaría AL, Paraiso Pueyo E, Pascual A, Botigué T. Fluid intake recommendation considering the physiological adaptations of adults over 65 years: A critical review. Nutrients. 2020;12(11):E3383.
Sivasankar M, Erickson E, Schneider S, Hawes A. Phonatory effects of airway dehydration: preliminary evidence for impaired compensation to oral breathing in individuals with a history of vocal fatigue. J Speech Lang Hear Res JSLHR. 2008;51(6):1494–506.
Sivasankar M, Leydon C. The role of hydration in vocal fold physiology. Curr Opin Otolaryngol Head Neck Surg. 2010;18(3):171–5.
Solomon NP, DiMattia MS. Effects of a vocally fatiguing task and systemic hydration on phonation threshold pressure. J Voice Off J Voice Found. 2000;14(3):341–62.
Alves M, Krüger E, Pillay B, van Lierde K, van der Linde J. The effect of hydration on voice quality in adults: A systematic review. J Voice Off J Voice Found. 2019;33(1):125.e13–28.
Chan RW, Gray SD, Titze IR. The importance of hyaluronic acid in vocal fold biomechanics. Otolaryngol Neck Surg. 2001;124(6):607–14.
Verdolini-Marston K, Titze IR, Druker DG. Changes in phonation threshold pressure with induced conditions of hydration. J Voice. 1990;4(2):142–51.
Chan RW, Tayama N. Biomechanical effects of hydration in vocal fold tissues. Otolaryngol--Head Neck Surg Off J Am Acad Otolaryngol-Head Neck Surg. 2002;126(5):528–37.
Nakagawa H, Fukuda H, Kawaida M, Shiotani A, Kanzaki J. Lubrication mechanism of the larynx during phonation: an experiment in excised canine larynges. Folia Phoniatr Logop Off Organ Int Assoc Logop Phoniatr IALP. 1998;50(4):183–94.
Duan C. Do Nascimento NC, calve S, Cox A, Sivasankar MP. Restricted water intake adversely affects rat vocal fold biology. Laryngoscope. 2021;131(4):839–45.
Cannes do Nascimento N, dos Santos AP, Sivasankar MP, Cox A. Unraveling the molecular pathobiology of vocal fold systemic dehydration using an in vivo rabbit model, PloS One. Jetté M, 2020;15(7):e0236348.
do Nascimento NC, Dos Santos AP, Mohallem R, Aryal UK, Xie J, Cox A, et al. Furosemide-induced systemic dehydration alters the proteome of rabbit vocal folds. J Proteome. 2022;10(252):104431.
Mizuta M, Hirano S, Hiwatashi N, Tateya I, Kanemaru SI, Nakamura T, et al. Effect of astaxanthin on vocal fold wound healing. Laryngoscope. 2014;124(1):1–7.
Mizuta M, Hirano S, Ohno S, Tateya I, Kanemaru SI, Nakamura T, et al. Expression of reactive oxygen species during wound healing of vocal folds in a rat model. Ann Otol Rhinol Laryngol. 2012;121(12):804–10.
Aslam B, Basit M, Nisar MA, Khurshid M, Rasool MH. Proteomics: technologies and their applications. J Chromatogr Sci. 2017;55(2):182–96.
Ajlan AM, Al-Khatib T, Al-Sheikah M, Jastaniah S, Salih A, Althubaiti A, et al. Helical computed tomography scanning of the larynx and upper trachea in rabbits. Acta Vet Scand. 2015;57(1):1–6.
Okhovat S, Milner TD, Clement WA, Wynne DM, Kunanandam T. Validation of animal models for simulation training in pediatric Laryngotracheal reconstruction. Ann Otol Rhinol Laryngol. 2020;129(1):46–54.
Branski RC, Rosen CA, Verdolini K, Hebda PA. Biochemical markers associated with acute vocal fold wound healing: A rabbit model. J Voice. 2005;19(2):283–9.
Döllinger M, Kniesburges S, Berry DA, Birk V, Wendler O, Dürr S, et al. Investigation of phonatory characteristics using ex vivo rabbit larynges. J Acoust Soc Am. 2018;144(1):142–52.
Kimball EE, Sayce L, Powell M, Gartling GJ, Brandley J, Rousseau B. Different vibratory conditions elicit different structural and biological vocal fold changes in an in-vivo rabbit model of phonation. J Voice. 2019;35(2):216-25.
Kojima T, Valenzuela CV, Novaleski CK, Van Deusen M, Mitchell JR, Garrett CG, et al. Effects of phonation time and magnitude dose on vocal fold epithelial genes, barrier integrity, and function. Laryngoscope. 2014;124(12):2770–8.
Kolosova K, Gao Q, Tuznik M, Bouhabel S, Kost KM, Wang H, et al. Characterizing vocal fold injury recovery in a rabbit model with three-dimensional virtual histology. Laryngoscope. 2021;131(7):1578–87.
Maytag AL, Robitaille MJ, Rieves AL, Madsen J, Smith BL, Jiang JJ. Use of the rabbit larynx in an excised larynx setup. J Voice Off J Voice Found. 2013;27(1):24–8.
Novaleski CK, Kimball EE, Mizuta M, Rousseau B. Acute exposure to vibration is an apoptosis-inducing stimulus in the vocal fold epithelium. Tissue Cell. 2016;48(5):407–16.
Rousseau B, Kojima T, Novaleski CK, Kimball EE, Valenzuela CV, Mizuta M, et al. Recovery of vocal fold epithelium after acute Phonotrauma. Cells Tissues Organs. 2017;204(2):93–104.
Rousseau B, Suehiro A, Echemendia N, Sivasankar M. Raised intensity phonation compromises vocal fold epithelial barrier integrity. Laryngoscope. 2011;121(2):346–51.
Rousseau B, Ge PJ, Ohno T, French LC, Thibeault SL. Extracellular matrix gene expression after vocal fold injury in a rabbit model. Ann Otol Rhinol Laryngol. 2008;117(8):598–603.
Le Bellego L, Jean C, Jiménez L, Magnani C, Tang W, Boutrolle I. Understanding fluid consumption patterns to improve healthy hydration. Nutr Today. 2010;45(6):S22.
Liska D, Mah E, Brisbois T, Barrios PL, Baker LB, Spriet LL. Narrative review of hydration and selected health outcomes in the general population. Nutrients. 2019;11(1):70.
Chang T, Ravi N, Plegue MA, Sonneville KR, Davis MM. Inadequate hydration, BMI, and obesity among US adults: NHANES 2009–2012. Ann Fam Med. 2016;14(4):320–4.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun [Internet]. 2019;10(1):1–10. https://doi.org/10.1038/s41467-019-09234-6.
Mohallem R, Aryal UK. Regulators of TNFα mediated insulin resistance elucidated by quantitative proteomics. Sci Rep. 2020;10(1):1–15.
Kim SQ, Mohallem R, Franco J, Buhman KK, Kim KH, Aryal UK. Multi-omics approach reveals dysregulation of protein phosphorylation correlated with lipid metabolism in mouse non-alcoholic fatty liver. Cells. 2022;11(7):1172.
Mohallem R, Aryal UK. Quantitative proteomics and Phosphoproteomics reveal TNF-α-mediated protein functions in hepatocytes. Molecules. 2021;26(18):5472.
Pross N, Demazières A, Girard N, Barnouin R, Metzger D, Klein A, et al. Effects of changes in water intake on mood of high and low drinkers. PLoS One. 2014;9(4):e94754.
Vergne S. Methodological aspects of fluid intake records and surveys. Nutr Today. 2012;47(4):S7.
Benton D, Young HA. Do small differences in hydration status affect mood and mental performance? Nutr Rev. 2015;73(suppl_2):83–96.
Popkin BM, D’Anci KE, Rosenberg IH. Water, hydration, and health. Nutr Rev. 2010;68(8):439–58.
Grossmann T, Darnhofer B, Birner-Gruenberger R, Kirsch A, Gugatschka M. Descriptive proteomics of paired human vocal fold and buccal mucosa tissue. PROTEOMICS – Clin Appl. 2022;16(2):2100050.
Ohno T, Hirano S. Treatment of aging vocal folds: novel approaches. Curr Opin Otolaryngol Head Neck Surg. 2014;22(6):472–6.
Branco A, Todorovic Fabro A, Gonçalves TM, Garcia Martins RH. Alterations in extracellular matrix composition in the aging larynx. Otolaryngol Neck Surg. 2015;152(2):302–7.
Miri AK, Li NYK, Avazmohammadi R, Thibeault SL, Mongrain R, Mongeau L. Study of extracellular matrix in vocal fold biomechanics using a two-phase model. Biomech Model Mechanobiol. 2015;14(1):49–57.
Kim JM, Heo HS, Shin SC, Kwon HK, Lee JC, Sung ES, et al. Increased calcium channel in the lamina propria of aging rat. Aging. 2019;11(20):8810–24.
Tang SS, Mohad V, Gowda M, Thibeault SL. Insights into the role of collagen in vocal fold health and disease. J Voice. 2017;31(5):520–7.
Cooper GM. Actin, Myosin, and Cell Movement. Cell Mol Approach 2nd 2000 [cited 2022 Jun 12]; Available from: https://www.ncbi.nlm.nih.gov/books/NBK9961/
Balagopal P, Rooyackers OE, Adey DB, Ades PA, Nair KS. Effects of aging on in vivo synthesis of skeletal muscle myosin heavy-chain and sarcoplasmic protein in humans. Am J Physiol-Endocrinol Metab. 1997;273(4):E790–800.
Oplatka A. The role of water in the mechanism of muscular contraction. FEBS Lett. 1994;355(1):1–3.
Lorenzo I, Serra-Prat M, Yébenes JC. The role of water homeostasis in muscle function and frailty: A review. Nutrients. 2019;11(8):E1857.
Yoo H, Nagornyak E, Das R, Wexler AD, Pollack GH. Contraction-induced changes in hydrogen bonding of muscle hydration water. J Phys Chem Lett. 2014;5(6):947–52.
Agnetti G, Herrmann H, Cohen S. New roles for desmin in the maintenance of muscle homeostasis. FEBS J. 2022;289(10):2755–70.
Bornemann A, Schmalbruch H. Desmin and vimentin in regenerating muscles. Muscle Nerve. 1992;15(1):14–20.
Russ DW, Grandy JS. Increased desmin expression in hindlimb muscles of aging rats. J Cachexia Sarcopenia Muscle. 2011;2(3):175.
Bühler RB, Sennes LU, Tsuji DH, Mauad T, Ferraz da Silva L, Saldiva PN. Collagen type I, collagen type III, and Versican in vocal fold Lamina Propria. Arch Otolaryngol Neck Surg. 2011;137(6):604–8.
The R-BS, Family C. Cold Spring Harb Perspect Biol. 2011;3(1):a004978.
Hammond TH, Gray SD, Butler JE. Age- and gender-related collagen distribution in human vocal folds. Ann Otol Rhinol Laryngol. 2000;109(10 Pt 1):913–20.
Madruga de Melo EC, Lemos M, Aragão Ximenes Filho J, Sennes LU, Nascimento Saldiva PH, Tsuji DH. Distribution of collagen in the lamina propria of the human vocal fold. Laryngoscope. 2003;113(12):2187–91.
Khoshnoodi J, Pedchenko V, Hudson BG. Mammalian collagen IV. Microsc Res Tech. 2008;71(5):357–70.
Tateya T, Tateya I, Bless DM. Collagen subtypes in human vocal folds. Ann Otol Rhinol Laryngol. 2006;115(6):469–76.
Uspenskaia O, Liebetrau M, Herms J, Danek A, Hamann GF. Aging is associated with increased collagen type IV accumulation in the basal lamina of human cerebral microvessels. BMC Neurosci. 2004;24(5):37.
Vázquez F, Palacios S, Alemañ N, Guerrero F. Changes of the basement membrane and type IV collagen in human skin during aging. Maturitas. 1996;25(3):209–15.
Bonnans C, Chou J, Werb Z. Remodelling the extracellular matrix in development and disease. Nat Rev Mol Cell Biol. 2014;15(12):786–801.
Swift J, Ivanovska IL, Buxboim A, Harada T, Dingal PCDP, Pinter J, et al. Nuclear Lamin-A scales with tissue stiffness and enhances matrix-directed differentiation. Science. 2013;341(6149):1240104.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13(4):227–32.
Liu Y, Beyer A, Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell. 2016;165(3):535–50.
Sengupta P, Dutta S. Mapping the age of laboratory rabbit strains to human. Int J Prev Med. 2020;11(12):44–51.
Persson PB. Renin: origin, secretion and synthesis. J Physiol. 2003;552(Pt 3):667–71.
Weidmann P, De Myttenaere-Bursztein S, Maxwell MH, de Lima J. Effect of aging on plasma renin and aldosterone in normal man. Kidney Int. 1975;8(5):325–33.
Gorelick MH, Shaw KN, Murphy KO. Validity and reliability of clinical signs in the diagnosis of dehydration in children. Pediatrics. 1997;99(5):e6.
Victora CG, Fuchs SC, Kirkwood BR, Lombardi C, Barros FC. Low body weight: A simple Indicator of the risk of dehydration among children with Diarrhoea. J Diarrhoeal Dis Res. 1997;15(1):7–11.
Cheuvront SN, Kenefick RW. Dehydration: physiology, assessment, and performance effects. Compr Physiol. 2014;4(1):257–85.
Islam S, Abély M, Alam NH, Dossou F, Chowdhury AKA, Desjeux JF. Water and electrolyte salvage in an animal model of dehydration and malnutrition. J Pediatr Gastroenterol Nutr. 2004;38(1):27–33.
Schwarz C, Hentschke H, Butovas S, Haiss F, Stüttgen MC, Gerdjikov TV, et al. The head-fixed behaving rat—procedures and pitfalls. Somatosens Mot Res. 2010;27(4):131–48.
Guo ZV, Hires SA, Li N, O’Connor DH, Komiyama T, Ophir E, et al. Procedures for Behavioral Experiments in Head-Fixed Mice. PLoS One, Simon SA, 2014;9(2):e88678.
Urai AE, Aguillon-Rodriguez V, Laranjeira IC, Cazettes F, Mainen ZF, Churchland AK. Citric Acid Water as an Alternative to Water Restriction for High-Yield Mouse Behavior. eNeuro. 2021;8(1) ENEURO.0230–20.2020.
Guide for the care and use of laboratory animals. Washington, D.C.: National Academies Press; 2011. Available from: http://www.nap.edu/catalog/12910
Dutta S, Sengupta P. Rabbits and men: relating their ages. J Basic Clin Physiol Pharmacol. 2018;29(5):427–35.
Santollo J, Daniels D. Control of fluid intake by estrogens in the female rat: role of the hypothalamus. Front Syst Neurosci. 2015;9:22–3.
Giersch GEW, Charkoudian N, Stearns RL, Casa DJ. Fluid balance and hydration considerations for women: review and future directions. Sports Med Auckl NZ. 2020;50(2):253–61.
Giersch GEW, Colburn AT, Morrissey MC, Butler CR, Pruchnicki ML, Kavouras SA, et al. Effects of sex and menstrual cycle on volume-regulatory responses to 24-h fluid restriction. Am J Physiol - Regul Integr Comp Physiol. 2020;319(5):R560–5.
Bailey TW, Pires dos Santos A, Cannes do Nascimento N, Xie J, Preeti Sivasankar M, Cox A. Recurring exposure to low humidity induces transcriptional and protein level changes in the vocal folds of rabbits. Sci Rep. 2021;11:24180.
Walker JM. The bicinchoninic acid (BCA) assay for protein quantitation. Methods Mol Biol. 1994;32:5-8.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72.
Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods. 2016;13(9):731–40.
Zar JH. Biostatistical analysis. 4th ed. Upper Saddle River, N.J: Prentice Hall; 1999. p. 1.
Hochberg Y, Benjamini Y. More powerful procedures for multiple significance testing. Stat Med. 1990;9(7):811–8.
Cohen J. A coefficient of agreement for nominal scales. Educ Psychol Meas. 1960;20(1):37–46.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Stark C. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34(90001):D535–9.
Türei D, Korcsmáros T, Saez-Rodriguez J. OmniPath: guidelines and gateway for literature-curated signaling pathway resources. Nat Methods. 2016;13(12):966–7.
Li T, Wernersson R, Hansen RB, Horn H, Mercer J, Slodkowicz G, et al. A scored human protein–protein interaction network to catalyze genomic interpretation. Nat Methods. 2017;14(1):61–4.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;27:1–27.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13(1):1–11.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25(4):402–8.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.
Acknowledgements
The authors thank the Purdue Proteomics Facility, where the proteomic experiments were conducted.
Funding
Funding was provided from R01DC015545 (National Institutes of Health/National Institute on Deafness and other Communication Disorders). The NIH/NIDCD did not participate in any component of the study design, execution, or analysis.
Author information
Authors and Affiliations
Contributions
NCN was responsible for the conception and design of the study, data acquisition, data analysis and interpretation, and manuscript writing. TWB was involved in the design of the study, data analysis, and reviewing and editing the manuscript. APS was involved in the design of the study, provided resources, reviewed and edited the manuscript. CD was involved in the design of the study and reviewing the manuscript. RM was involved in the data acquisition and analysis, and reviewing the manuscript. JF was involved in the data acquisition and analysis, and reviewing the manuscript. UKA was involved in the data acquisition and analysis, and reviewing the manuscript. JX was involved in supervising the statistical analysis, and reviewing the manuscript. AC was responsible for conception and design of the study, supervision of the study, funding, and reviewing and editing the manuscript. MPS was responsible for conception and design of the study, supervision of the study, funding, reviewing and editing the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Experiments were approved by the Purdue Institutional Animal Care and Use Committee (IACUC) under the Protocol # 1606001428 and were conducted following the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals.
The study is reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
do Nascimento, N.C., Bailey, T.W., Santos, A.P. et al. Proteomic analysis reveals that aging rabbit vocal folds are more vulnerable to changes caused by systemic dehydration. BMC Genomics 23, 762 (2022). https://doi.org/10.1186/s12864-022-08975-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-022-08975-x