Next Article in Journal
Longitudinal Data to Assess Relationships among Plasma Folate, Vitamin B12, Non-esterified Fatty Acid, and β-Hydroxybutyrate Concentrations of Holstein Cows during the Transition Period
Next Article in Special Issue
Metabolomic Profiling of Bipolar Disorder by 1H-NMR in Serbian Patients
Previous Article in Journal
Metabolic Inheritance and the Competition for Calories between Mother and Fetus
Previous Article in Special Issue
Association of Metabolomic Biomarkers with Sleeve Gastrectomy Weight Loss Outcomes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 13
Issue 4
10.3390/metabo13040546
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Profiling in Mouse Model of Menopause-Associated Asthma

by
William P. Pederson
1,2,
Laurie M. Ellerman
3,
Yan Jin
4,
Haiwei Gu
4 and
Julie G. Ledford
3,5,*
1
Physiological Sciences GIDP, University of Arizona, Tucson, AZ 85724, USA
2
Department of Physiology, University of Arizona, Tucson, AZ 85724, USA
3
Asthma and Airway Disease Research Center, Tucson, AZ 85724, USA
4
Center for Translational Science, Florida International University, Port St. Lucie, FL 34987, USA
5
Department of Cellular and Molecular Medicine, University of Arizona, Tucson, AZ 85724, USA
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(4), 546; https://doi.org/10.3390/metabo13040546
Submission received: 24 March 2023 / Revised: 6 April 2023 / Accepted: 7 April 2023 / Published: 11 April 2023
(This article belongs to the Special Issue Advanced Metabolomics and Lipidomics Approaches in Studying Human Diseases)
Browse Figures
Versions Notes

Abstract

:
Menopause-associated asthma impacts a subset of women, tends to be more severe, and is less responsive to current treatments. We recently developed a model of menopause-associated asthma using 4-Vinylcyclohexene Diepoxide (VCD) and house dust mites (HDM). The goal of this study was to uncover potential biomarkers and drivers of menopause-onset asthma by assessing serum and bronchoalveolar lavage fluid (BALF) samples from mice with and without menopause and HDM challenge by large-scale targeted metabolomics. Female mice were treated with VCD/HDM to model menopause-associated asthma, and serum and BALF samples were processed for large-scale targeted metabolomic assessment. Liquid chromatography–tandem mass spectrometry (LC-MS/MS) was used to examine metabolites of potential biological significance. We identified over 50 individual metabolites, impacting 46 metabolic pathways, in the serum and BALF that were significantly different across the four study groups. In particular, glutamate, GABA, phosphocreatine, and pyroglutamic acid, which are involved in glutamate/glutamine, glutathione, and arginine and proline metabolisms, were significantly impacted in the menopausal HDM-challenged mice. Additionally, several metabolites had significant correlations with total airway resistance including glutamic acid, histamine, uridine, cytosine, cytidine, and acetamide. Using metabolic profiling, we identified metabolites and metabolic pathways that may aid in discriminating potential biomarkers for and drivers of menopause-associated asthma.

Graphical Abstract

1. Introduction

Asthma is a common respiratory disease, characterized by episodic, and reversible, airway obstruction, affecting nearly 8% of the US population [1]. Symptoms include wheezing, coughing, and shortness of breath, all of which contribute to an inability to breathe [2,3,4]. However, asthma is a heterogeneous disease, with prevalence and severity varying wildly across age, sex, and ethnicity, emphasizing environmental and genetic factors that contribute to disease pathogenesis [5]. This heterogeneity has led researchers to identify and characterize specific asthma endotypes based on symptoms and cellular and molecular phenotypes [5,6]. The most common of these endotypes is allergic asthma, which involves eosinophilic inflammation in the airways in response to allergenic stimuli [6]. Allergic asthma, along with several other common asthma endotypes, has been well studied, and there are a variety of treatments and therapies available. However, one subtype of asthma that remains less well defined is menopause-associated asthma [7,8].
Women experiencing menopause-associated asthma often encounter more severe symptoms, which often require hospitalization and advanced treatment beyond the standard of care [9]. As such, late-onset asthma in menopausal women often goes uncontrolled, since standard asthma treatments are relatively ineffective [7]. Recently, our lab developed a mouse model of menopause-associated asthma in order to study the mechanisms driving disease onset associated with menopause and identify differences from non-menopause onset asthma [8]. Utilizing the 4-Vinylcyclohexene Diepoxide (VCD) mouse model of menopause, in combination with the house dust mite (HDM) mouse model of asthma, we were able to demonstrate significantly increased airway hyper-responsiveness and non-Th2 airway inflammation in menopausal, HDM-challenged mice [8]. To the best of our knowledge, this is the first animal model to mimic clinical phenotypes observed in women with menopause-associated asthma.
Despite the success of the VCD/HDM model in replicating the severe airway function and non-Th2 inflammation seen in women, the underlying molecular mechanisms remain unknown. It is crucial to have a comprehensive understanding of cellular and molecular mechanisms of asthma in order to identify biomarkers unique to specific endotypes and to discover novel targets for new therapies. In recent years, metabolomics, or quantitative analysis of small molecules in biological samples, has become a powerful tool for identifying mechanisms of complex diseases, including asthma [10,11,12]. Metabolomic studies performed on human samples, including blood, lung tissue, sputum, urine, and exhaled breath condensates, have reported a wide range of metabolites and metabolic pathways that distinguish asthma subtypes based on age, asthma severity, obesity, and other factors [12,13]. Some of the most commonly reported metabolic pathways affected by asthma include lipid metabolism, amino acid metabolism, and oxidative stress [11]. While there is an abundance of metabolomic data describing potential mechanisms for asthma, these studies also highlight the severe gap in our understanding of the complexity of the signaling pathways at play. Additional metabolomic studies and comprehensive review of reported metabolic pathways is necessary to develop a stronger understanding of asthma endotypes and their respective cellular and molecular mechanisms.
In this study, we performed large-scale targeted metabolomics on samples from mice from two compartments, bronchoalveolar lavage fluid (BALF) and serum, across four experimental conditions: control, HDM challenged, menopausal non-HDM challenged, and menopausal HDM challenged. These findings contribute to the current metabolomic profiles being developed for asthma endotypes in humans and shed light on potential biomarkers that could be used to discriminate menopause-associated asthma from non-menopause asthma, in addition to providing novel therapeutic targets for a group of patients with limited treatment options.

2. Materials and Methods

Mouse Models: All experiments utilizing animal models were performed in accordance with University of Arizona Institutional Animal Care and Use Committee-approved animal protocols. Wild-Type (WT) female mice on a C57BL/6J background were ordered from Jackson Laboratories. Mice were age-matched (7 weeks) at the time of delivery, and experiments began when age reached 8 weeks. As shown in Figure 1A, intraperitoneal injections of the VCD compound (160 mg/kg, Sigma, St. Louis, MO, USA, #94956) in sesame oil (Sigma, #S3547), or sesame oil alone as vehicle control, were performed daily for twenty consecutive days [8]. Two weeks following the final IP injection, cyclicity was performed daily via a vaginal lavage, the cytology of which was monitored daily to determine the onset of menopause (ovarian failure), defined as 10 consecutive days of diestrus. Sensitization and challenge involved intranasal administration of HDM (Greer Laboratories, Lenoir, NC, USA, #NC1468783), or saline vehicle control, once a week for a total of three weeks. Mice that received HDM were given 50 μL of solution at a concentration of 2 mg/mL [8].
Sample Collection: Twenty-four hours after the third HDM challenge, mice were weighed, and lung function tests were performed on the flexiVent system (SciReq, Montreal, QC, Canada) [8]. Bronchoalveolar Lavage Fluid (BALF) was obtained by gently flushing the lungs with 1.5 mL of lavage fluid (PBS, 0.1 mM EDTA) three times. Mice were dissected, and blood was collected via cardiac puncture. Blood samples were centrifuged to separate the serum, and serum was transferred to a new tube. Cells suspended in solution in the BALF were separated via centrifugation, and the cell-free BALF was aliquoted into new tubes. BALF and serum samples were stored at −80 °C from the 4 experimental groups as detailed in Figure 1B until further analysis.
Metabolomics Analysis: Before metabolomic analysis was performed, protein precipitation and metabolite extraction were performed on all samples. To begin, 50 μL of either BALF and serum sample were transferred to a 2 mL Eppendorf tube after being thawed overnight at 4 °C. Then, 500 μL MeOH and 50 μL internal standard solution (containing 1810.5 μM 13C3-lactate and 142 μM 13C5-glutamic acid) was added to each sample, followed by vortexing for 10 s, storage at −20 °C for 30 min, and finally centrifugation at 14,000 RPM for 10 min at 4 °C. The supernatants, typically around 450 μL, were then transferred to a new Eppendorf tube, dried using a CentriVap Concentrator (Labconco, Fort Scott, KS, USA), and reconstituted in 150 μL of 40% PBS/60% ACN. A quality control (QC) sample was created by pooling all plasma samples together. We used a targeted LC-MS/MS method for metabolomic assessment that was modeled after methods used regularly in a large number of studies [14,15,16,17,18,19,20]. All LC-MS/MS experiments were performed utilizing the Agilent 1290 UPLC-6495 QQQ-MS (Santa Clara, CA, USA) system. Each sample was analyzed using negative ionization mode (using 10 μL of total sample) and positive ionization mode (using 4 μL of total sample). Hydrophilic interaction chromatography (HILIC) mode on a Waters XBridge BEH Amide column (150 × 2.1 mm, 2.5 µm particle size, Waters Corporation, Milford, MA, USA) was used for both chromatographic separations. The auto-sampler temperature was kept at 4 °C, the column compartment was set at 40 °C, and the flow rate was 0.3 mL/min. Solvents A (10 mM ammonium acetate, 10 mM ammonium hydroxide in 95% H2O/5% ACN) and B (10 mM ammonium acetate, 10 mM ammonium hydroxide in 95% ACN/5% H2O) made up the mobile phase. The percentage of Solvent B decreased to 40% at t = 11 min after the initial 1 min isocratic elution of 90% B. The composition of Solvent B was then maintained at 40% for 4 min (t = 15 min), before gradually increasing to 90% to prepare for the next injection. An electrospray ionization (ESI) source was equipped to the mass spectrometer and targeted data acquisition was performed in multiple-reaction-monitoring (MRM) mode. Agilent Masshunter Workstation software (version 10.1, Santa Clara, CA, USA) controller the entire LC-MS system. Agilent MassHunter Quantitative Data Analysis (Santa Clara, CA, USA) was used to integrate the extracted MRM peaks [14,15,16,17,18,19,20].
Statistical Analysis: Statistical analysis of previously published airway function data [8] was account for variations of recorded data on separate experimental days, and statistical significance was assessed using unpaired t-tests, according to previous publications [8,21]. Analysis of metabolites and pathway analyses were performed using MetaboAnalyst 5.0 software. Metabolite data were normalized using log transformation, and statistical analysis was performed using an ANOVA, for comparison of all four groups together, or an unpaired t-test, for comparison of only two groups. For analysis using ANOVA, significance was achieved with p ≤ 0.05, using Tukey’s Honest Significant Difference. For analysis using unpaired t-test, significance was achieved with p ≤ 0.05, and a fold change threshold of 1.5. Pathway analysis was performed comparing each experimental group to one another. Partial least squares-discriminant analysis (PLS-DA) was performed comparing all four experimental groups.

3. Results

3.1. Impact of Menopause on Metabolomic Profile

As previously published, the HDM and menopause HDM experimental groups had a robust eosinophilic response compared to their saline controls but were not different from one another; both having >70% eosinophils and averaging ~6 × 105 Eos/mL BALF [8]. The menopause HDM group had approximately twice the number and percentage of neutrophils in BALF as compared to the HDM only group [8]; however, neutrophils were much lower compared to the eosinophils, as was expected of the HDM model of allergic airways disease in mice.
By comparing the menopausal group to the controls, we observed how menopause affects metabolomic changes in the BALF and serum. There were no statistically significant upregulated metabolites in the BALF of menopausal mice compared to control, but 4-aminobutyric acid (GABA) showed an increasing trend (Table 1).
Menopausal mice had a significant upregulation of several metabolites in the serum, including pyroglutamic acid (Figure 2D), glucosamine, 3-hydroxyisovaleric acid, myristic acid, nicotinuric acid, 4-pyrodoxic acid, oxoglutaric acid, 2-methylglutaric acid and indole 3-acetic acid compared to non-menopausal mice (Table 2).
Conversely, menopause resulted in a downregulation of several metabolites in the BALF, including TMAO, 2-hydroxybenzoic acid (salicylic acid), threonine, homoserine, hippuric acid, hydroxyproline, 5-aminolevulinic acid (Table 1). Interestingly, menopause had an even greater impact on the metabolomic profile in the serum. Menopause resulted in a significant downregulation of many metabolites in the serum, including 2-hydroxybenzoic acid (salicylic acid), oxidized glutathione, mucic acid, 5-aminolevulinic acid, acetylcholine, ribose, methyl guanidine, xylose, hydroxyproline, carnosine, uridine, and GA3P (Table 2).
Pathway analysis was performed, which groups individual metabolites into known metabolic pathways. In BALF, menopause did not significantly impact any metabolic pathways in the KEGG database (Figure 3A, Table 3). However, in serum, menopause significantly impacted 16 different known metabolic pathways, the most significantly affected being glutathione metabolism, glycine, serine, and threonine metabolism, pentose and glucuronate interconversions, pentose phosphate pathway, alanine, aspartate, and glutamate metabolism, d-glutamine and d-glutamate metabolism, and tryptophan metabolism (Figure 3E, Table 4). PLS-DA analysis of the BALF metabolites in all groups showed some overlap between menopausal and control groups, although discrimination of these two groups can be observed (Figure 4A). PLS-DA analysis of the serum metabolites shows very distinct discrimination between menopausal mice and control mice (Figure 4B).

3.2. Impact of HDM on Metabolomic Profile

Comparing the HDM-challenged group to the non-HDM-challenged controls, we assessed how HDM, our asthma model, affected metabolomic changes in the BALF and serum. Asthma resulted in a significant downregulation of several metabolites in the BALF, including fructose, ribose, and TMAO (Table 1), while metabolites 9-octadecynoic acid (Figure 2A), acetamide (Figure 2B), and glycocyamine were significantly upregulated in the asthma group (Table 1). When assessing the serum samples, asthma led to significant downregulation of metabolites 3-methyladipic acid, methyl-D-mannopyranoside, histamine, and carnosine (Table 2). Using a pathway analysis, we discovered that asthma significantly impacted arginine and proline metabolism, amino sugar and nucleotide sugar metabolism, and glycolysis gluconeogeneis pathways in the BALF (Figure 3B, Table 3). Additionally, asthma significantly affected the glycerophospholipid metabolism, selenocompound metabolism, histidine metabolism, and the glycine, serine, and threonine metabolism pathways in the serum (Figure 3F, Table 4). Asthma and control groups overlap entirely in the PLS-DA analysis of serum metabolites (Figure 4A). However, these groups do begin to separate when PLS-DA is performed on BALF (Figure 4B).

3.3. Impact of HDM in Menopausal Mice on Metabolomic Profile

By comparing the menopausal, HDM-challenged group to the menopausal control group, we assessed the impact of HDM challenge in menopausal mice. In menopausal mice, HDM challenge resulted in a significant upregulation of many metabolites in the BALF, including cytosine (Figure 2C), glycocyamine, 9-octadecynoic acid (Figure 2A), d-ribose 5-phosphate, cytidine, and phenylpyruvic acid (Table 1). In the serum, HDM challenge in menopausal mice only resulted in a significant upregulation of two metabolites: uridine and acetylcholine (Table 2). Using pathway analysis, HDM challenge in menopausal mice significantly impacted fatty acid elongation and degradation metabolism, purine metabolism, pyrimidine metabolism, pentose phosphate metabolism, and biosynthesis of unsaturated fatty acids in the BALF (Figure 3C, Table 3). In the serum, we saw significant impacts on glycerophospholipid metabolism, arachidonic metabolism, pyrimidine metabolism, and selenocompound metabolism in the menopausal HDM-challenged mice compared to the menopausal non-HDM-challenged mice (Figure 3G, Table 4). Menopausal HDM-challenged and menopause groups overlap entirely in the PLS-DA analysis of serum and BALF metabolites (Figure 4A,B).

3.4. Impact of Menopause in HDM-Challenged Mice on Metabolomic Profile

To distinguish menopause-associated asthma from non-menopausal asthma, we compared the metabolomic profiles of the menopausal HDM-challenged mice to the HDM-challenged control mice. In the BALF, we found significant upregulation of cytosine (Figure 2C), cytidine, glutamic acid, and d-ribose 5 phosphate in menopausal HDM-challenged mice compared to non-menopausal HDM-challenged mice (Table 1). In the serum, we observed significant upregulation of a single metabolite, phosphocreatine (Table 2). Pathway analysis of the BALF between these groups revealed significant impacts on pyrimidine metabolism, pentose phosphate metabolism, nitrogen metabolism, d-glutamine and d-glutamate metabolism, porphyrin and chlorophyll metabolism, arginine and proline metabolism, butanoate metabolism, and purine metabolism (Figure 3D, Table 3). Interestingly, a single pathway was significantly impacted in the serum, arginine and proline metabolism (Figure 3H, Table 4). PLS-DS analysis in the menopausal HDM-challenged and non-menopausal HDM-challenged group shows distinct separation of populations in the serum and BALF (Figure 4A,B).

3.5. Associations between Metabolites and Airway Hyper-Responsiveness

We next wanted to assess whether any of the discovered metabolites were associated with airway function in the same mice across all experimental groups. The most significantly different metabolites were paired with the maximum total airway resistance measured at the highest dose of methacholine, which was obtained during airway function studies published previously [8]. In BALF, we found that cytosine, cytidine, acetamide, and glutamic acid were all significantly and positively correlated with total airway resistance at the 100mg/mL dose of methacholine (Figure 5B,C,E,F), indicating that when these metabolites are elevated, airway resistance is also high, and lung function is worse. In serum, we found that histamine had a significant negative correlation with total airway resistance at the 100mg/mL dose of methacholine (Figure 5A), suggesting that when histamine levels are low in the serum, airway resistance is high, and lung function is worse. Additionally, serum uridine levels had a significant positive correlation with total airway resistance (Figure 5D), indicating that when uridine is high in the serum, airway resistance is also high and lung function is worse.

4. Discussion

Asthma is a heterogeneous disease that consists of complex inflammatory pathways that contribute to prevalence and severity of disease. Distinct mechanistic pathways have been described to differentiate asthma endotypes and allow for the development of more precise treatments [22]. However, menopause-associated asthma has yet to be investigated in depth, and has the potential to be a distinct asthma endotype, with a unique inflammatory signature, and may require specific treatments to maintain asthma control. The lack of knowledge surrounding menopause-associated asthma demonstrates an urgent need to study and understand the molecular mechanisms driving menopause-associated asthma. In this study, we have utilized large-scale targeted metabolomics to describe the metabolomic profiles of mice undergoing the VCD/HDM model of menopause-associated asthma, with the goal of discovering novel metabolites that may serve as potential biomarkers for the identification and treatment of menopause-associated asthma.
Metabolomics offers a powerful tool to uncover novel mechanisms of complex diseases and involves the quantitative analysis of small molecules in biological samples and their relation to pathophysiological changes in organisms [10,11,12]. Metabolomic profiling can also lead to the discovery of novel biomarkers that can help identify diseases [23]. Metabolomics has been used extensively to study asthma and airway inflammation, describing a wide variety of changes in metabolomic profiles, some of which include lipid metabolism, amino acid metabolism, energy metabolism disorders, and oxidation and reduction imbalance [11,12,13]. With this study, we aimed to describe the metabolomic profile in BALF and serum of menopause-associated asthma in the VCD/HDM model to distinguish it from non-menopausal asthma.
The metabolic pathways that are significantly impacted in the BALF of the menopausal HDM-challenged group compared to the non-menopausal HDM-challenged group share a specific metabolite: glutamic acid is a player in nitrogen, d-glutamine and d-Glutamate, porphyrin and chlorophyll, arginine and proline, and butanoate metabolic pathways. In this case, glutamic acid is the main metabolite driving the dysregulation of these metabolic pathways. glutamic acid, also known as l-glutamate has been identified in our study as significantly upregulated in the BALF of menopausal HDM-challenged mice compared to non-menopausal HDM-challenged mice. In addition, glutamate shows a positive correlation with airway resistance, indicating that higher levels of glutamate in the BALF correlate with higher levels of airway resistance. Glutamate is one of the most abundant excitatory neurotransmitters in the body, with many synapses using a variety of different glutamate receptors [24]. Additionally, glutamate is the metabolic precursor of y-aminobutyric acid, also known as GABA [25]. Previously, GABA has been shown to be secreted by pulmonary neuroendocrine cells to induce goblet cell hyperplasia in primate models [26]. GABA has also been shown to be produced by pulmonary neuroendocrine cells in murine lungs [27]. While not quite achieving statistical significance, GABA was also detected in the metabolomic assessment of BALF when comparing menopausal mice to non-menopausal mice, demonstrating a trending upregulation of the metabolite (p = 0.0796, FC 5.42). This suggests that menopause alone may impact neuroendocrine signaling in the lung, and in turn, affect airway epithelial cell differentiation. With increased goblet cell hyperplasia, we would also expect to see mucus hypersecretion, a key contributor to the airway obstruction seen in HDM challenged airways. However, more studies will be needed to validate these findings.
The arginine and proline metabolism pathway was the only one identified as significantly impacted in the serum analysis of the menopausal HDM-challenged mice compared to the non-menopausal HDM-challenged mice. While glutamic acid is still involved in this pathway, it was not significantly different between these two groups when individual metabolites were assessed in serum. Phosphocreatine, which was upregulated, was the only metabolite significantly different between the menopausal HDM-challenged and non-menopausal HDM-challenged groups. Interestingly, studies have shown that glutamate uptake in synaptic vesicles is not only dependent on ATP but can also be regulated by phosphocreatine [28]. While this phenomenon was demonstrated in isolated bovine brain synapses, it is possible that phosphocreatine may influence glutamate uptake in pulmonary neuroendocrine synapses as well, which could influence airway epithelial cell differentiation similar to GABA, as previously described [26]. Phosphocreatine was also upregulated in the serum of non-human primates sensitized and challenged with Ascaris suum extract and was returned to control levels with experimental treatment [29]. Having been detected in other animal models of allergic asthma, serum phosphocreatine holds potential for being a biomarker for severe asthma. However, other researchers suggest that elevated serum phosphocreatine during airway exacerbations may simply be the product of the increased work of the respiratory muscles [30].
Our previous work suggests that non-Th2 inflammatory pathways, specifically Th17, may be contributing to the severe asthmatic phenotype seen in menopausal HDM-challenged mice [8]. Non-Th2 driven asthma, including Th1 and Th17, is often more severe and difficult to treat [31,32,33]. Other studies have demonstrated that Th1 and Th17 T cell differentiation is dependent on the conversion of glutamine to glutamate via the glutaminase enzyme [34,35]. The authors concluded that glutamine metabolism has a unique role in promoting Th17 cell differentiation but restricting Th1 and CTL effector cell differentiation. With the significant changes in “glutathione”, “alanine, aspartate, and glutamate”, and “d-glutamine and d-glutamate” metabolism observed in the serum of menopausal mice compared to control, it is possible that increased glutamic acid may be affecting the T-helper cell differentiation, which then influences the inflammatory response with allergen challenge.
Pyroglutamic acid, an intermediate in glutathione metabolism that can be converted to glutamate by 5-oxoprolinase [36], is another metabolite upregulated in the serum of menopausal mice. A recent study detected significantly increased pyroglutamic acid, as well as glutamine, in the plasma of ovalbumin sensitized mice, which was then reduced with experimental treatment [37]. In addition, both metabolites were significantly associated with Penh values (measure of lung function) at various concentrations of methacholine challenge, eosinophil numbers, and Th17 cell percentages. Interestingly, pyroglutamic acid was not significantly upregulated in the non-menopausal HDM-challenged group compared to control in our study. Instead, we observed significant increases in the menopausal group compared to control and non-menopausal HDM challenged, as well as the menopausal HDM challenged compared to control, suggesting that menopause is the driving factor in the dysregulation of pyroglutamic acid in the VCD/HDM model.
A previous study identified significantly impaired glutathione metabolism in the lung during an endoplasmic reticulum (ER) stress mouse model [38]. This same study also noted TMAO significantly downregulated in the ER stress model compared to control [38]. TMAO has also been identified in plasma and urine of asthma patients and is being considered as a possible biomarker for asthma exacerbation [39,40]. This is in line with our findings, in which we discovered significantly downregulated TMAO in the BALF of menopausal and HDM-challenged mice compared to control, and downregulated oxidized glutathione in serum of menopausal mice compared to control. Glutathione metabolism is believed to be involved in the regulation of redox homeostasis and relieving ER stress in the lung [38] and may a significant player in the pathogenesis of menopause-associated asthma.
Histamine was significantly lower in the serum of HDM-challenged mice compared to control (Table 2) and had a significant negative correlation with airways hyper-responsiveness (Figure 5A). This contrasts with previous studies that demonstrated high histamine levels in the serum of sensitized animals as asthmatic patients [41]. This discrepancy may be due to the short half-life of histamine [42], or the timepoint at which collection occurred did not align with peak histamine release. Uridine, which is involved in nucleotide and pyrimidine metabolism, was significantly lower in the serum of menopausal mice compared to control, but significantly higher in menopausal HDM-challenged mice compared to menopausal mice (Table 2). Additionally, uridine had a significant correlation with airway hyper-responsiveness (Figure 5D). Previous literature has described uridine as being anti-inflammatory in models of lung inflammation [43].
We also identified several novel metabolites and metabolic pathways that have not been previously linked to asthma or lung disease, which warrant further investigation. Cytosine, which is one of the four base nucleotides that make up DNA, and cytidine, which is formed when cytosine is attached to a ribose ring in RNA, were both significantly upregulated in the BALF of menopausal HDM-challenged mice compared to all other groups, with the exception of menopausal HDM challenged vs. non-menopausal HDM challenged with cytosine (Figure 2C, Table 1). Both of these metabolites had significant positive correlations with airway hyper-responsiveness (Figure 5B,C). While not well studied, cytosine was elevated in patients with COVID-19 compared to healthy controls [44]. Additionally, researchers have shown that elevated cytosine levels are involved in the tryptophan/nicotinamide metabolic pathway, which has been linked to inflammatory signaling [45,46]. Coincidentally, nicotinuric acid was upregulated in serum of menopausal mice compared to control in our study. Cytidine was found to be upregulated in the sputum on patients with bronchial asthma [47,48].
9-octadecynoid acid, also known as stearolic acid, and acetamide were significantly elevated in the BALF of both HDM-challenged groups compared to controls, with the menopausal HDM-challenged mice having the highest concentrations (Figure 2A,B, Table 1). Acetamide had a significant positive correlation with airway hyper-responsiveness (Figure 5E). 9-Octadecynoic acid has only been described in a handful of metabolomics studies [49,50,51], none of which were related to the lung. Also very understudied, acetamide has been identified as a metabolomic biomarker for pregnancy complications [52,53,54,55]. With the significant differences in our studies and their strong associations with lung function, these two relatively understudied metabolites in regard to lung function and asthma serve as the most novel potential biomarkers for menopause-associated asthma and warrant additional study.
The major strength of this study lies in the VCD/HDM model of menopause in mice. This model results in robust and consistent changes in airway function in the menopausal HDM-challenged mice compared to non-menopausal mice, measures that typically yield high variability when using the Flexivent system. Additionally, we were able to assess metabolomic profiles from two compartments: the serum, which allowed us to understand the global effects of menopause and/or asthma on the animal, and BALF, which gave us a more precise view of the lungs and airways. The primary limitation of this study is the sample size. While n = 5 per group is sufficient to observe metabolomic changes, an increased sample size would aid in narrowing down the most significant metabolites and pathways that may be involved in disease severity and progression. Moreover, our samples were only collected at a single timepoint after the onset of menopause, so we are unable to infer longitudinal changes in metabolites across the menopausal transition. Additionally, a better understanding of the underlying changes in the central nervous system during menopause, related to hormonal fluctuation and the impact on lung function, remains a gap in knowledge that warrants further exploration in future studies.
In conclusion, our metabolomic analysis on samples from mice allowed us to better describe the metabolic landscape in the serum and BALF in the menopausal HDM-challenged mice in comparison to the other experimental groups. We discovered a large variety of significantly impacted metabolites and metabolic pathways in both compartments impacted by both menopause and asthma status. Additionally, PLS-DA analysis of the metabolites of the experimental groups was able to discriminate menopausal HDM-challenged and non-menopausal HDM-challenged groups. With this new information, we proposed possible mechanisms for disease, and potential novel biomarkers, which may be used to aid in the development of new therapies for a group of patients that currently have very few options. While additional studies are needed to validate results in humans and delve further into mechanisms of action driving asthma onset and reduced lung function in menopausal asthmatics, our study provides the foundation for detailing metabolomic changes in the lung and in circulation that may impact the pathological phenotypes manifested in menopause-associated asthma.

Author Contributions

Conceptualization, W.P.P., H.G. and J.G.L.; methodology, W.P.P., L.M.E., Y.J. and H.G.; software, W.P.P., Y.J. and H.G.; validation, Y.J. and H.G.; formal analysis, W.P.P.; investigation, W.P.P., L.M.E. and J.G.L.; resources, H.G. and J.G.L.; data curation, Y.J. and H.G.; writing—original draft preparation, W.P.P.; writing—review and editing, L.M.E., Y.J., H.G. and J.G.L.; visualization, W.P.P.; supervision, J.G.L.; project administration, H.G. and J.G.L.; funding acquisition, J.G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH, grant numbers HL142769, 5T32AG058503-02.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Animal Care and Use Committee of University of Arizona (protocol # 15-575, 2 May 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. Ongoing studies with this particular data set are still in progress.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. National Current Asthma Prevalence. Centers for Disease Control and Prevention (CDC). 2019. Available online: https://www.cdc.gov/asthma/most_recent_national_asthma_data.htm (accessed on 1 January 2023).
  2. Castillo, J.R.; SPPeters Busse, W.W. Asthma Exacerbations: Pathogenesis, Prevention, and Treatment. J. Allergy Clin. Immunol. Pract. 2017, 5, 918–927. [Google Scholar] [CrossRef]
  3. Maslan, J.; Mims, J.W. What is asthma? Pathophysiology, demographics, and health care costs. Otolaryngol. Clin. N. Am. 2014, 47, 13–22. [Google Scholar] [CrossRef] [PubMed]
  4. Mims, J.W. Asthma: Definitions and pathophysiology. Int. Forum Allergy Rhinol. 2015, 5 (Suppl. 1), S2–S6. [Google Scholar] [CrossRef]
  5. Wenzel, S.E. Asthma phenotypes: The evolution from clinical to molecular approaches. Nat. Med. 2012, 18, 716–725. [Google Scholar] [CrossRef] [PubMed]
  6. Desai, M.; Oppenheimer, J. Elucidating asthma phenotypes and endotypes: Progress towards personalized medicine. Ann. Allergy Asthma Immunol. 2016, 116, 394–401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Balzano, G.; Fuschillo, S.; De Angelis, E.; Gaudiosi, C.; Mancini, A.; Caputi, M. Persistent airway inflammation and high exacerbation rate in asthma that starts at menopause. Monaldi Arch. Chest Dis. 2007, 67, 135–141. [Google Scholar] [CrossRef]
  8. Pederson, W.P.; Ellerman, L.M.; Sandoval, E.C.; Boitano, S.; Frye, J.B.; Doyle, K.P.; Brooks, H.L.; Polverino, F.; Ledford, J.G. Development of a Novel Mouse Model of Menopause-associated Asthma. Am. J. Respir. Cell Mol. Biol. 2022, 67, 605–609. [Google Scholar] [CrossRef]
  9. Balzano, G.; Fuschillo, S.; Melillo, G.; Bonini, S. Asthma and sex hormones. Allergy 2001, 56, 13–20. [Google Scholar] [CrossRef]
  10. Luxon, B.A. Metabolomics in asthma. Adv. Exp. Med. Biol. 2014, 795, 207–220. [Google Scholar]
  11. Wang, C.; Jiang, S.; Zhang, S.; Ouyang, Z.; Wang, G.; Wang, F. Research Progress of Metabolomics in Asthma. Metabolites 2021, 11, 567. [Google Scholar] [CrossRef] [PubMed]
  12. Xu, S.; Panettieri, R.A., Jr.; Jude, J. Metabolomics in asthma: A platform for discovery. Mol. Asp. Med. 2022, 85, 100990. [Google Scholar] [CrossRef] [PubMed]
  13. Pite, H.; Morais-Almeida, M.; Rocha, S.M. Metabolomics in asthma: Where do we stand? Curr. Opin. Pulm. Med. 2018, 24, 94–103. [Google Scholar] [CrossRef] [PubMed]
  14. Carroll, P.A.; Diolaiti, D.; McFerrin, L.; Gu, H.; Djukovic, D.; Du, J.; Cheng, P.F.; Anderson, S.; Ulrich, M.; Hurley, J.B.; et al. Deregulated Myc requires MondoA/Mlx for metabolic reprogramming and tumorigenesis. Cancer Cell 2015, 27, 271–285. [Google Scholar] [CrossRef] [Green Version]
  15. Eghlimi, R.; Shi, X.; Hrovat, J.; Xi, B.; Gu, H. Triple Negative Breast Cancer Detection Using LC-MS/MS Lipidomic Profiling. J. Proteome Res. 2020, 19, 2367–2378. [Google Scholar] [CrossRef]
  16. Gu, H.; Carroll, P.A.; Du, J.; Zhu, J.; Neto, F.C.; Eisenman, R.N.; Raftery, D. Quantitative Method to Investigate the Balance between Metabolism and Proteome Biomass: Starting from Glycine. Angew. Chem. Int. Ed. Engl. 2016, 55, 15646–15650. [Google Scholar] [CrossRef] [PubMed]
  17. Gu, H.; Zhang, P.; Zhu, J.; Raftery, D. Globally Optimized Targeted Mass Spectrometry: Reliable Metabolomics Analysis with Broad Coverage. Anal. Chem. 2015, 87, 12355–12362. [Google Scholar] [CrossRef] [Green Version]
  18. Jasbi, P.; Mitchell, N.M.; Shi, X.; Grys, T.E.; Wei, Y.; Liu, L.; Lake, D.F.; Gu, H. Coccidioidomycosis Detection Using Targeted Plasma and Urine Metabolic Profiling. J. Proteome Res. 2019, 18, 2791–2802. [Google Scholar] [CrossRef]
  19. Jasbi, P.; Shi, X.; Chu, P.; Elliott, N.; Hudson, H.; Jones, D.; Serrano, G.; Chow, B.; Beach, T.G.; Liu, L.; et al. Metabolic Profiling of Neocortical Tissue Discriminates Alzheimer’s Disease from Mild Cognitive Impairment, High Pathology Controls, and Normal Controls. J. Proteome Res. 2021, 20, 4303–4317. [Google Scholar] [CrossRef]
  20. Shi, X.; Wang, S.; Jasbi, P.; Turner, C.; Hrovat, J.; Wei, Y.; Liu, J.; Gu, H. Database-Assisted Globally Optimized Targeted Mass Spectrometry (dGOT-MS): Broad and Reliable Metabolomics Analysis with Enhanced Identification. Anal. Chem. 2019, 91, 13737–13745. [Google Scholar] [CrossRef]
  21. Addison, K.J.; Morse, J.; Robichaud, A.; Daines, M.O.; Ledford, J.G. A Novel in vivo System to Test Bronchodilators. J. Infect. Pulm. Dis. 2017, 3, 10-16966. [Google Scholar]
  22. Kuruvilla, M.E.; Lee, F.E.; Lee, G.B. Understanding Asthma Phenotypes, Endotypes, and Mechanisms of Disease. Clin. Rev. Allergy Immunol. 2019, 56, 219–233. [Google Scholar] [CrossRef] [PubMed]
  23. Villasenor, A.; Eguiluz-Gracia, I.; Moreira, A.; Wheelock, C.E.; Escribese, M.M. Metabolomics in the Identification of Biomarkers of Asthma. Metabolites 2021, 11, 346. [Google Scholar] [CrossRef] [PubMed]
  24. Pal, M.M. Glutamate: The Master Neurotransmitter and Its Implications in Chronic Stress and Mood Disorders. Front. Hum. Neurosci. 2021, 15, 722323. [Google Scholar] [CrossRef] [PubMed]
  25. Petroff, O.A. GABA and glutamate in the human brain. Neuroscientist 2002, 8, 562–573. [Google Scholar] [CrossRef]
  26. Barrios, J.; Kho, A.T.; Aven, L.; Mitchel, J.A.; Park, J.A.; Randell, S.H.; Miller, L.A.; Tantisira, K.G.; Ai, X. Pulmonary Neuroendocrine Cells Secrete gamma-Aminobutyric Acid to Induce Goblet Cell Hyperplasia in Primate Models. Am. J. Respir. Cell Mol. Biol. 2019, 60, 687–694. [Google Scholar] [CrossRef]
  27. Reader, J.R.; Tepper, J.S.; Schelegle, E.S.; Aldrich, M.C.; Putney, L.F.; Pfeiffer, J.W.; Hyde, D.M. Pathogenesis of mucous cell metaplasia in a murine asthma model. Am. J. Pathol. 2003, 162, 2069–2078. [Google Scholar] [CrossRef] [Green Version]
  28. Xu, C.J.; Klunk, W.E.; Kanfer, J.N.; Xiong, Q.; Miller, G.; Pettegrew, J.W. Phosphocreatine-dependent glutamate uptake by synaptic vesicles. A comparison with atp-dependent glutamate uptake. J. Biol. Chem. 1996, 271, 13435–13440. [Google Scholar] [CrossRef] [Green Version]
  29. Camateros, P.; Kanagaratham, C.; Najdekr, L.; Holub, D.; Vrbkova, J.; Cote, L.; Fournier, J.; Gourdon, J.; Creery, D.; Olivenstein, R.; et al. Toll-Like Receptor 7/8 Ligand, S28463, Suppresses Ascaris suum-induced Allergic Asthma in Nonhuman Primates. Am. J. Respir. Cell Mol. Biol. 2018, 58, 55–65. [Google Scholar] [CrossRef]
  30. Burki, N.K.; Diamond, L. Serum creatine phosphokinase activity in asthma. Am. Rev. Respir. Dis. 1977, 116, 327–331. [Google Scholar]
  31. Al Heialy, S.; Ramakrishnan, R.K.; Hamid, Q. Recent advances in the immunopathogenesis of severe asthma. J. Allergy Clin. Immunol. 2022, 149, 455–465. [Google Scholar] [CrossRef]
  32. Boonpiyathad, T.; Sozener, Z.C.; Satitsuksanoa, P.; Akdis, C.A. Immunologic mechanisms in asthma. Semin. Immunol. 2019, 46, 101333. [Google Scholar] [CrossRef]
  33. Newcomb, D.C.; Peebles, R.S., Jr. Th17-mediated inflammation in asthma. Curr. Opin. Immunol. 2013, 25, 755–760. [Google Scholar] [CrossRef] [Green Version]
  34. Johnson, M.O.; Wolf, M.M.; Madden, M.Z.; Andrejeva, G.; Sugiura, A.; Contreras, D.C.; Maseda, D.; Liberti, M.V.; Paz, K.; Kishton, R.J.; et al. Distinct Regulation of Th17 and Th1 Cell Differentiation by Glutaminase-Dependent Metabolism. Cell 2018, 175, 1780–1795.e19. [Google Scholar] [CrossRef] [Green Version]
  35. Miao, Y.; Zheng, Y.; Geng, Y.; Yang, L.; Cao, N.; Dai, Y.; Wei, Z. The role of GLS1-mediated glutaminolysis/2-HG/H3K4me3 and GSH/ROS signals in Th17 responses counteracted by PPARgamma agonists. Theranostics 2021, 11, 4531–4548. [Google Scholar] [CrossRef]
  36. Schilling, S.; Wasternack, C.; Demuth, H.U. Glutaminyl cyclases from animals and plants: A case of functionally convergent protein evolution. Biol. Chem. 2008, 389, 983–991. [Google Scholar] [CrossRef]
  37. Zhou, Y.; Zhao, H.; Wang, T.; Zhao, X.; Wang, J.; Wang, Q. Anti-Inflammatory and Anti-asthmatic Effects of TMDCT Decoction in Eosinophilic Asthma Through Treg/Th17 Balance. Front. Pharmacol. 2022, 13, 819728. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, Z.; Ma, P.; Wang, Y.; Hou, B.; Zhou, C.; Tian, H.; Li, B.; Shui, G.; Yang, X.; Qiang, G.; et al. Untargeted metabolomics and transcriptomics identified glutathione metabolism disturbance and PCS and TMAO as potential biomarkers for ER stress in lung. Sci. Rep. 2021, 11, 14680. [Google Scholar] [CrossRef]
  39. Chiu, C.Y.; Cheng, M.L.; Chiang, M.H.; Wang, C.J.; Tsai, M.H.; Lin, G. Metabolomic Analysis Reveals Distinct Profiles in the Plasma and Urine Associated with IgE Reactions in Childhood Asthma. J. Clin. Med. 2020, 9, 887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Loureiro, C.C.; Duarte, I.F.; Gomes, J.; Carrola, J.; Barros, A.S.; Gil, A.M.; Bousquet, J.; Bom, A.T.; Rocha, S.M. Urinary metabolomic changes as a predictive biomarker of asthma exacerbation. J. Allergy Clin. Immunol. 2014, 133, 261–263. [Google Scholar] [CrossRef] [PubMed]
  41. Chen, S.; Chen, G.; Shu, S.; Xu, Y.; Ma, X. Metabolomics analysis of baicalin on ovalbumin-sensitized allergic rhinitis rats. R. Soc. Open Sci. 2019, 6, 181081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Hashimoto, T.; Ohata, H.; Honda, K. Lysophosphatidic acid (LPA) induces plasma exudation and histamine release in mice via LPA receptors. J. Pharmacol. Sci. 2006, 100, 82–87. [Google Scholar] [CrossRef] [Green Version]
  43. Trevethick, M.A. Is uridine a treatment for asthma? Clin. Exp. Allergy 2010, 40, 1436–1438. [Google Scholar] [CrossRef]
  44. Roberts, I.; Wright Muelas, M.; Taylor, J.M.; Davison, A.S.; Xu, Y.; Grixti, J.M.; Gotts, N.; Sorokin, A.; Goodacre, R.; Kell, D.B. Untargeted metabolomics of COVID-19 patient serum reveals potential prognostic markers of both severity and outcome. Metabolomics 2021, 18, 6. [Google Scholar] [CrossRef] [PubMed]
  45. Blasco, H.; Bessy, C.; Plantier, L.; Lefevre, A.; Piver, E.; Bernard, L.; Marlet, J.; Stefic, K.; Benz-de Bretagne, I.; Cannet, P.; et al. The specific metabolome profiling of patients infected by SARS-COV-2 supports the key role of tryptophan-nicotinamide pathway and cytosine metabolism. Sci. Rep. 2020, 10, 16824. [Google Scholar] [CrossRef] [PubMed]
  46. Tounta, V.; Liu, Y.; Cheyne, A.; Larrouy-Maumus, G. Metabolomics in infectious diseases and drug discovery. Mol. Omics 2021, 17, 376–393. [Google Scholar] [CrossRef] [PubMed]
  47. Takhar, P.; Corrigan, C.J.; Smurthwaite, L.; O’Connor, B.J.; Durham, S.R.; Lee, T.H.; Gould, H.J. Class switch recombination to IgE in the bronchial mucosa of atopic and nonatopic patients with asthma. J. Allergy Clin. Immunol. 2007, 119, 213–218. [Google Scholar] [CrossRef]
  48. Tian, M.; Chen, M.; Bao, Y.L.; Xu, C.D.; Qin, Q.Z.; Zhang, W.X.; He, Y.T.; Shao, Q. Sputum metabolomic profiling of bronchial asthma based on quadruple time-of-flight mass spectrometry. Int. J. Clin. Exp. Pathol. 2017, 10, 10363–10373. [Google Scholar]
  49. Berry, D.E.; Chan, J.A.; MacKenzie, L.; Hecht, S.M. 9-Octadecynoic acid: A novel DNA binding agent. Chem. Res. Toxicol. 1991, 4, 195–198. [Google Scholar] [CrossRef]
  50. Meyer, F.; Bloch, K. Metabolism of Stearolic Acid in Yeast. J. Biol. Chem. 1963, 238, 2654–2659. [Google Scholar] [CrossRef]
  51. Thomasson, H.J. Stearolic acid, an essential fatty acid? Nature 1954, 173, 452. [Google Scholar] [CrossRef]
  52. Bahado-Singh, R.O.; Akolekar, R.; Mandal, R.; Dong, E.; Xia, J.; Kruger, M.; Wishart, D.S.; Nicolaides, K. Metabolomics and first-trimester prediction of early-onset preeclampsia. J. Matern. Fetal Neonatal Med. 2012, 25, 1840–1847. [Google Scholar] [CrossRef] [PubMed]
  53. Bahado-Singh, R.O.; Akolekar, R.; Mandal, R.; Dong, E.; Xia, J.; Kruger, M.; Wishart, D.S.; Nicolaides, K. First-trimester metabolomic detection of late-onset preeclampsia. Am. J. Obstet. Gynecol. 2013, 208, 58.e1–58.e7. [Google Scholar] [CrossRef] [PubMed]
  54. Bahado-Singh, R.O.; Akolekar, R.; Mandal, R.; Dong, E.; Xia, J.; Kruger, M.; Wishart, D.S.; Nicolaides, K. Metabolomic analysis for first-trimester Down syndrome prediction. Am. J. Obstet. Gynecol. 2013, 208, 371.e1–371.e8. [Google Scholar] [CrossRef] [PubMed]
  55. Bahado-Singh, R.O.; Ertl, R.; Mandal, R.; Bjorndahl, T.C.; Syngelaki, A.; Han, B.; Dong, E.; Liu, P.B.; Alpay-Savasan, Z.; Wishart, D.S.; et al. Metabolomic prediction of fetal congenital heart defect in the first trimester. Am. J. Obstet. Gynecol. 2014, 211, 240.e1–240.e14. [Google Scholar] [CrossRef] [PubMed]
Metabolites 13 00546 g001 550
Figure 1. Experimental setup. (A) Timeline of the VCD/HDM model denoting times of VCD dosing, cyclicity checking, HDM dosing, Flexivent and tissue harvest. (B) Table denoting the four experimental groups, as well as n values per group.
Figure 1. Experimental setup. (A) Timeline of the VCD/HDM model denoting times of VCD dosing, cyclicity checking, HDM dosing, Flexivent and tissue harvest. (B) Table denoting the four experimental groups, as well as n values per group.
Metabolites 13 00546 g001
Metabolites 13 00546 g002 550
Figure 2. Significantly different metabolites across the four experimental groups by one-way ANOVA with multiple comparisons. (A) 9-Octadecynoic acid in BALF; (B) acetamide in BALF; (C) cytosine in BALF; (D) pyroglutamic acid in serum. * p ≤ 0.05, ** p ≤ 0.005, *** p ≤ 0.0005, **** p ≤ 0.00005.
Figure 2. Significantly different metabolites across the four experimental groups by one-way ANOVA with multiple comparisons. (A) 9-Octadecynoic acid in BALF; (B) acetamide in BALF; (C) cytosine in BALF; (D) pyroglutamic acid in serum. * p ≤ 0.05, ** p ≤ 0.005, *** p ≤ 0.0005, **** p ≤ 0.00005.
Metabolites 13 00546 g002
Metabolites 13 00546 g003 550
Figure 3. Pathway analysis plots for group comparison in BALF (A–D) and serum (E,F). (A,E) Menopause vs. Control; (B,F) HDM vs. Control, (C,G) Menopause HDM vs. Menopause; (D,H) Menopause HDM vs. HDM. The log p-values of the significantly different metabolites are plotted within a given metabolic pathway by the impact score that those metabolites have on the pathway. The size of the circle correlates to the number of significantly different metabolites detected for that pathway. Numbers next to circles correspond to numbers next to pathways listed in Table 3 and Table 4.
Figure 3. Pathway analysis plots for group comparison in BALF (A–D) and serum (E,F). (A,E) Menopause vs. Control; (B,F) HDM vs. Control, (C,G) Menopause HDM vs. Menopause; (D,H) Menopause HDM vs. HDM. The log p-values of the significantly different metabolites are plotted within a given metabolic pathway by the impact score that those metabolites have on the pathway. The size of the circle correlates to the number of significantly different metabolites detected for that pathway. Numbers next to circles correspond to numbers next to pathways listed in Table 3 and Table 4.
Metabolites 13 00546 g003
Metabolites 13 00546 g004 550
Figure 4. PLS-DA analysis of control, HDM-challenged, menopausal, and menopausal HDM-challenged mice. Score plots of BALF (A) and serum (B) using all significant metabolites.
Figure 4. PLS-DA analysis of control, HDM-challenged, menopausal, and menopausal HDM-challenged mice. Score plots of BALF (A) and serum (B) using all significant metabolites.
Metabolites 13 00546 g004
Metabolites 13 00546 g005 550
Figure 5. Correlation plots of metabolites with total airway resistance (Rrs). Airway resistance from the 100mg/mL methacholine dose was assessed using the Flexivent system and association analysis performed for the most significantly different metabolites detected. (A) Serum histamine concentration with a negative correlation to Rrs, * p = 0.044. (B) BALF cytosine concentration with a positive correlation to Rrs, *** p = 0.0005. (C) BALF cytidine with a positive correlation to Rrs, * p = 0.013. (D) Serum uridine with a positive correlation with Rrs, * p = 0.039. (E) BALF acetamide with a positive correlation to Rrs, ** p = 0.0076. (F) BALF glutamic acid with a positive correlation to Rrs, * p = 0.0163. All statistical significance was achieved using Pearson’s correlation coefficient.
Figure 5. Correlation plots of metabolites with total airway resistance (Rrs). Airway resistance from the 100mg/mL methacholine dose was assessed using the Flexivent system and association analysis performed for the most significantly different metabolites detected. (A) Serum histamine concentration with a negative correlation to Rrs, * p = 0.044. (B) BALF cytosine concentration with a positive correlation to Rrs, *** p = 0.0005. (C) BALF cytidine with a positive correlation to Rrs, * p = 0.013. (D) Serum uridine with a positive correlation with Rrs, * p = 0.039. (E) BALF acetamide with a positive correlation to Rrs, ** p = 0.0076. (F) BALF glutamic acid with a positive correlation to Rrs, * p = 0.0163. All statistical significance was achieved using Pearson’s correlation coefficient.
Metabolites 13 00546 g005
Table
Table 1. Significantly different individual metabolites in the BALF.
Table 1. Significantly different individual metabolites in the BALF.
ComparisonChangeMetabolitep ValueFold Change
Menopause
vs.
Control		DownregulatedTMAO0.0070.11
2-HydroxybenzoicAcid 0.0120.30
Threonine0.0140.59
Homoserine0.0160.60
Hippuric Acid0.0380.46
Hydroxyproline0.0460.63
5-Aminolevulinic Acid0.0490.67
Upregulated4-Aminobutyric Acid0.0805.42
HDM
vs.
Control		DownregulatedFructose0.0410.60
Ribose0.0450.61
TMAO0.0480.25
Upregulated9-Octadecynoic Acid0.0013.31
Acetamide0.0051.55
Glycocyamine0.0122.15
Menopause HDM
vs.
Menopause		UpregulatedCytosine0.00033.96
GlycoCyamine0.0051.92
9-Octadecynoic Acid0.0112.24
D-Ribose 5-phosphate0.0132.98
Cytidine0.0232.21
Phenylpyruvic Acid0.0401.86
Menopause HDM
vs.
HDM		UpregulatedCytosine0.000052.17
Cytidine0.0152.05
Glutamic Acid0.0332.44
D-Ribose 5 Phosphate0.0471.96
Statistical significance acquired using t-test with a p-value cutoff of 0.05 and a fold change cutoff of 1.5. Italicized metabolites did not reach statistical significance of p ≤ 0.05 but were of particular interest.
Table
Table 2. Significantly different individual metabolites in the serum.
Table 2. Significantly different individual metabolites in the serum.
ComparisonChangeMetabolitep ValueFold Change
Menopause
vs.
Control		Downregulated2-Hydroxybenzoic Acid 0.0070.12
Oxidized Glutathione0.0070.61
Mucic Acid0.0080.38
5-aminolevulinic acid0.0080.67
Acetylcholine0.0080.67
Ribose0.0110.22
Methyl Guanidine0.0110.34
Xylose0.0110.19
Hydroxyproline0.0120.65
Carnosine0.0320.59
Uridine0.0310.47
GA3P0.0350.33
UpregulatedPyroglutamic Acid0.00032.00
Glucosamine0.0061.63
3-hydroxyisovaleric Acid0.0071.63
Myristic Acid0.0151.66
Nicotinuric Acid0.0201.83
4-Pyridoxic Acid0.0241.73
Oxoglutaric Acid0.0321.82
2-Methylglutaric Acid0.0371.77
Indole 3-acetic Acid0.0422.33
HDM
vs.
Control		Downregulated3-Methyladipic Acid0.0050.66
Methyl-D-Mannopyranoside0.0230.61
Histamine0.0260.40
Carnosine0.0450.64
Menopause HDM
vs.
Menopause		UpregulatedUridine0.0282.27
Acetylcholine0.0341.66
Menopause HDM
vs.
HDM		UpregulatedPhosphocreatine0.0141.86
Statistical significance acquired using t-test with a p-value cutoff of 0.05 and a fold change cutoff of 1.5.
Table
Table 3. Significantly impacted metabolic pathways in the BALF.
Table 3. Significantly impacted metabolic pathways in the BALF.
ComparisonPathway NameMatch Statusp ValueImpact
Menopause
vs.
Control		1. Arginine & Proline Metabolism9/380.0940.30
2. Amino Sugar & Nucleotide Sugar Metabolism1/370.1000.0
3. Phenylalanine Metabolism3/120.1100.62
HDM
vs.
Control		1. Arginine & Proline Metabolism9/380.0070.30
2. Amino Sugar & Nucleotide Sugar Metabolism1/370.0410.0
3. Glycolysis/Gluconeogenesis1/260.0590.0
Menopause HDM
vs.
Menopause		1. Fatty Acid Elongation1/390.0030.0
2. Fatty Acid Degradation1/390.0030.0
3. Purine Metabolism1/660.0130.01
4. Pyrimidine Metabolism1/390.0230.01
5. Pentose Phosphate Pathway3/220.0240.17
6. Biosynthesis of Unsaturated Fatty Acids2/360.0410.0
Menopause HDM
vs.
HDM		1. Pyrimidine Metabolism1/390.0150.01
2. Pentose Phosphate Pathway3/220.0250.17
3. Nitrogen Metabolism1/60.0330.0
4. D-Glutamine & D-Glutamate metabolism2/60.0350.5
5. Porphyrin and Chlorophyll Metabolism3/300.0350.03
6. Arginine and Proline Metabolism9/380.0450.30
7. Butanoate Metabolism5/150.0450.03
8. Purine Metabolism1/660.0470.01
Numbers in front of pathway names correspond to numbers next to circles in Figure 3A–D. Statistical significance acquired using t-test with a p-value cutoff of 0.05. Italicized metabolic pathways did not reach statistical significance of p ≤ 0.05 but were of particular interest.
Table
Table 4. Significantly impacted metabolic pathways in the serum.
Table 4. Significantly impacted metabolic pathways in the serum.
ComparisonPathway NameMatch Statusp ValueImpact
Menopause
vs.
Control		1. Glutathione Metabolism7/280.00050.43
2. Fatty Acid Elongation1/390.0040.0
3. Fatty Acid Degradation1/390.0040.0
4. Fatty Acid Biosynthesis5/470.0060.02
5. Glycine, serine, & Threonine Metabolism11/340.0070.54
6. Glycerophospholipid Metabolism2/360.0090.03
7. Pentose & Glucuronate Interconversions5/180.0110.38
8. Porphyrin and Chlorophyll Metabolism3/300.0140.03
9. Pentose Phosphate Pathway3/220.0140.17
10. Alanine, aspartate, & Glutamate Metabolism9/280.0190.59
11. Vitamin B6 Metabolism1/90.0240.0
12. Butanoate Metabolism5/150.0330.0
13. D-Glutamine & D-Glutamate metabolism3/60.0360.5
14. Beta-Alanine Metabolism3/210.0380.06
15. Pyrimidine Metabolism5/390.0380.09
16. Tryptophan Metabolism4/410.0500.25
HDM
vs.
Control		1. Glycerophospholipid Metabolism2/360.0230.03
2. Selenocompound Metabolism1/200.0240.0
3. Histidine Metabolism9/160.0260.71
4. Glycine, Serine, & Threonine Metabolism11/340.0410.54
Menopause HDM
vs.
Menopause		1. Glycerophospholipid metabolism2/360.0330.03
2. Arachidonic Acid Metabolism1/360.0350.0
3. Pyrimidine Metabolism5/390.0420.09
4. Selenocompound Metabolism1/200.0450.0
Menopause HDM
vs.
HDM		1. Arginine and Proline Metabolism8/380.0150.38
Numbers in front of pathway names correspond to numbers next to circles in Figure 3E–H. Statistical significance acquired using t-test with a p-value cutoff of 0.05.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pederson, W.P.; Ellerman, L.M.; Jin, Y.; Gu, H.; Ledford, J.G. Metabolomic Profiling in Mouse Model of Menopause-Associated Asthma. Metabolites 2023, 13, 546. https://doi.org/10.3390/metabo13040546

AMA Style

Pederson WP, Ellerman LM, Jin Y, Gu H, Ledford JG. Metabolomic Profiling in Mouse Model of Menopause-Associated Asthma. Metabolites. 2023; 13(4):546. https://doi.org/10.3390/metabo13040546

Chicago/Turabian Style

Pederson, William P., Laurie M. Ellerman, Yan Jin, Haiwei Gu, and Julie G. Ledford. 2023. "Metabolomic Profiling in Mouse Model of Menopause-Associated Asthma" Metabolites 13, no. 4: 546. https://doi.org/10.3390/metabo13040546

APA Style

Pederson, W. P., Ellerman, L. M., Jin, Y., Gu, H., & Ledford, J. G. (2023). Metabolomic Profiling in Mouse Model of Menopause-Associated Asthma. Metabolites, 13(4), 546. https://doi.org/10.3390/metabo13040546

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop