Abstract
Background: Along with acceleration of population aging, the increasing prevalence of sarcopenia has posed a heavy burden on families as well as society. In this context, it is of great significance to diagnose and intervene sarcopenia as early as possible. Recent evidence has indicated the role of cuproptosis in the development of sarcopenia. In this study, we aimed to seek the key cuproptosis-related genes that can be used for identification and intervention of sarcopenia.
Methods: The GSE111016 dataset was retrieved from GEO. The 31 cuproptosis-related genes (CRGs) were obtained from previous published studies. The differentially expressed genes (DEGs) and Weighed gene co-expression network analysis (WGCNA) were subsequently analyzed. The core hub genes were acquired by the intersection of DEGs, WGCNA and CRGs. Through logistic regression analysis, we established a diagnostic model of sarcopenia based on the selected biomarkers and was validated in muscle samples from GSE111006 and GSE167186. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis were performed on these genes. Furthermore, the gene set enrichment analysis (GSEA), and immune cell infiltration were also conducted on the identified core genes. Finally, we screened the potential drugs targeting the potential biomarkers of sarcopenia.
Results: A total of 902 DEGs and WGCNA containing 1,281 significant genes were preliminarily selected. Intersection of DEGs, WGCNA and CRGs yielded four core genes (PDHA1, DLAT, PDHB, and NDUFC1) as potential biomarkers for the prediction of sarcopenia. The predictive model was established and validated with high AUC values. KEGG pathway and Gene Ontology biological analysis indicated these core genes may play a crucial role in energy metabolism in mitochondria, oxidation process, and aging-related degenerative diseases. In addition, the immune cells may be involved in the development of sarcopenia through mitochondrial metabolism. Finally, metformin was identified as a promising strategy of sarcopenia treatment via targeting NDUFC1.
Conclusion: The four cuproptosis-related genes PDHA1, DLAT, PDHB and NDUFC1 may be the diagnostic biomarkers for sarcopenia, and metformin holds great potential to be developed as a therapy for sarcopenia. These outcomes provide new insights for better understanding of sarcopenia and innovative therapeutic approaches.
Introduction
Sarcopenia is an age-related disease characterized by progressive loss of skeletal muscle mass and strength with decline in physical function (). The aging individuals with sarcopenia are at higher risk of adverse outcomes, including frailty, impaired quality of life, falls, fracture, disability, and mortality (; ; ). The prevalence of sarcopenia increases with advancing age, which was approximately 10% for older people aged 60–70 years, and 30% for those over 80 years old (; ). Sarcopenia is connected to skeletal muscle fibers loss and muscle atrophy. The muscle fibers are generally divided into slow-twitch type I muscle fiber and fast-twitch type II muscle fiber based on expression of different myosin isoforms and their contractile properties (). On the muscle fiber level, aging does not affect the muscle fiber types equivalently, and the loss of muscle fibers is not homogeneous (Tosetti et al., 2020). Skeletal muscle aging predominantly causes fast-twitch type II muscle fiber atrophy and fiber necrosis with decline in muscle strength, whereas type I fibers are less affected due to transition of muscle fibers from type II to type I with age (). The reduction in the numbers of both two types of muscle fibers and specific type II muscle fiber atrophy, leading to the ageing-related phenotype known as sarcopenia. In addition, reduced number of type II satellite cells (Verdijk et al., 2014), together with progressive infiltration of fat and non-contractile tissue into muscle are also involved in the development of sarcopenia (Visser et al., 2005). Despite the rapidly increasing translation research in the therapeutic field of sarcopenia, exercise remains the primary treatment of sarcopenia (), and no specific drug has been approved for the treatment of sarcopenia (). Hence, it is of great significance to fully understand the potential mechanism and to develop new effective strategy for treating sarcopenia.
Copper is one of the essential trace elements in the human body playing several vital roles in living organisms. As an enzyme cofactor, copper is considerably involved in a broad range of biological processes, such as mitochondrial respiration, antioxidant defense, and iron metabolism (). Copper also exerts its function as a dynamic signaling element with effects on the processes of lipolysis, cytotoxicity, cell proliferation, oxidative stress. In this regard, copper homeostasis is crucial for cell survival and metabolism. Copper is a redox active metal, the imbalance in its concentration and metabolism, is connected with the development of many diseases, including cancer, neurodegenerative diseases, metabolic syndrome, anemia and Wilson disease (). Accumulating evidence suggests that copper toxicity disrupts mitochondrial metabolic enzymes, which induces a specific regulatory cell death mechanism defined as cuproptosis (). Distinguishable from oxidative stress-related cell death, such as apoptosis, necroptosis and ferroptosis, cuproptosis is ignited by mitochondrial stress, especially the aggregation of mitochondrial enzymes lipoylation. The overload of copper can directly bind to lipoylated components of the tricarboxylic acid cycle in mitochondrial respiration, leading to the accumulating protein toxic stress and cell death (Tsvetkov et al., 2022). Current data has revealed the role of cuproptosis ignited by aging and genetic mutation in the progression of pathological sequelae, including neurodegenerative diseases, cancer, and inflammation (; ; Zhang et al., 2022; Zhao et al., 2022). The development of sarcopenia, characterized by age-associated with loss of skeletal muscle mass and strength, has been shown to be caused by dysregulated mitochondrial function, increased ROS generation and denervation (; ). However, at present, there are few studies on the cuproptosis and their association with sarcopenia. Therefore, in this study, we aim to explore the cuproptosis related biomarkers that can be used as prognostic indicators for sarcopenia, to further understand the potential mechanism of cuproptosis in the development of sarcopenia, which may improve the therapeutic strategies for prevention, diagnosis and treatment of sarcopenia in the future.
Methods
Data sources
The expression profiling generated by the high-throughput sequencing dataset GSE111016 was downloaded from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). The GSE111016 expression profiling contained 40 muscle biopsy samples, including 20 subjects with sarcopenia and 20 healthy controls. The cuproptosis-related genes were extracted from previous published studies, which contained 31 cuproptosis-related genes (Supplementary Table S1) (; ; Tsvetkov et al., 2022).
Data processing
The raw count of GSE111016 is normalized by “edgeR”, and DEGs between sarcopenia and healthy control were identified via using the “limma” package [adjusted p value < 0.05 and log2 (fold change (FC)]≥ 0.3785. The DEGs with log FC < 0 were considered as down-regulated, while log FC > 0 were regarded as up-regulated. The volcano plots of DEGs and cluster analysis heat map were visualized via using “ggplot2” and “heatmap” packages in R software.
Weighed gene co-expression network analysis
Weighed gene co-expression network analysis (WGCNA) is a bioinformatic method for the identification of the trait-related genes with importance (). The co-expression network for sarcopenia was constructed with “WGCNA” package in R. We removed the abnormal expressed genes from RNA-seq data of the training set by using goodSampleGenes method. The spearman’s correlation matrices and average linkage method were performed for all the pair-wise genes, and the similarity matrix was constructed by calculating the correlation coefficient between gene pairs. A suitable soft threshold power was set to transform the similarity matrix into an adjacency matrix to ensure the construction of a scale-free network. Subsequently, the adjacency was transformed into a topological overlap matrix (TOM) to measure the mean network connectivity of each gene, defined as the sum of its adjacency with all other genes for network generation, and then, the corresponding dissimilarity (1-TOM) was calculated to form clusters. According to the TOM-based dissimilarity measure with a minimum size of 30 for the genes dendrogram, average linkage hierarchical clustering was performed to classify co-expressed genes into different gene modules. To further analyze the module, we calculated the dissimilarity of module eigen (ME) genes, chose a cut line for module dendrogram and merged some module. We combined the modules with distance of less than 0.25, and finally obtained 14 co-expression modules shown in different colors. The module with the highest absolute value of the correlation coefficient was identified as the hub gene candidates for further analysis. The correlation coefficient representing the association between a gene module and sarcopenia phenotype, is considered significant with p value < 0.05.
Hub genes identification and enrichment analysis
To further identify the sarcopenia related hub genes, we used the “VennDiagram” package to intersect the genes obtained by WGCNA and DEGs, as well as the cuproptosis-related genes. To explore the expression differences of the hub genes, the expression of each hub gene in two groups were validated and shown in box plots. To further understand the biological mechanism of the hub genes that are associated with sarcopenia, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis were performed and visualized by “clusterProfiler” and “enrichplot” packages in R, the adjust p value < 0.05 was considered statistically significant.
Verification of hub genes
In this study, according to the expression profile data of hub genes, a receiver operating characteristic curve (ROC) was established via “pROC” in R. The predictive sensitivity and specificity of the diagnostic model were evaluated by the area under the ROC curve (AUC). To further evaluate the diagnostic value of the hub genes in the training GSE111016 dataset, two independent datasets (GSE111006 and GSE167186) were retrieved from GEO to verify the robustness of the established model.
Gene-set enrichment analysis
To investigate the biological role of the four hub genes in sarcopenia, the patients were categorized into high and low expression groups based on the median expression levels of these hub genes. The background gene set was obtained from MSigDB (). The outcomes of the analysis were presented using the “enrichplot” package in the R software for visualization.
Immune cell infiltration and correlation analysis
To explore the immune microenvironment in sarcopenia patients, the immune infiltrations of immune cells were analyzed using the CIBERSORT algorithm (). A total of 17 immune cells was analyzed, including T cells, B cells, NK cells, monocytes, macrophages, neutrophils, and dendritic cells. Box plots were drawn to show the expression levels of 17 immune infiltrating cells. Spearman correlations were calculated to explore the association of different type of immune cells and were visualized via the ggplot2 package.
Target drugs analysis
The Drug Gene Interaction database (DGIdb, https://dgidb.org/) was utilized to screen the potential drugs targeting cuproptosis-related genes in sarcopenia (). The 3D chemical structures of corresponding drugs were visualized and downloaded from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) (Wang et al., 2012). 3D structures of target proteins were retrieved from the RCSB Protein Data Bank (https://www.rcsb.org/) to explore the docking mode of drugs (). The molecular docking with receptor protein and ligands acquired from the above steps was performed with ChimeraX 1.5 software (). To further verify the binding properties of ligand to the target protein, AutoDock vina 1.1.2 and PyMOL 4.3.0 tools were used for perform docking simulation and visualization. The docking binding energy < 0 kcal/mol indicates that the receptor and ligand can bind spontaneously, and binding energy < −5.0 kcal/mol is considered to have an excellent binding activity.
Results
Identification of sarcopenia related RNAs
The gene expression profiles in vastus lateralis muscle of sarcopenia men and their age and race-matched corresponding healthy controls were retrieved from the GSE111016 dataset, which was further analyzed to investigate the biological functions of cuproptosis in the progression of sarcopenia. The age of healthy control was 69.95 ± 4.15 while the age of sarcopenia group was 72.75 ± 4.28, however, there was no statistical difference between the two groups. The detailed information of subjects in the GSE111016 was presented in Supplementary Table S2. A detailed flowchart of the study was showed in Figure 1. A total of 31,623 RNA-seq were obtained from GSE111016 database, and 902 DEGs were preliminary screened according to p value < 0.05 and FC ≥ 1.3, of which 288 were upregulated and 614 were downregulated. The DEGs are shown in the volcano plot (Figure 2A), the top 30 DEGs were plotted with a heatmap (Figure 2B).
FIGURE 1
FIGURE 2
Construction and analysis of gene co-expression network
We applied WGCNA analysis on the data of GSE111016 dataset to build the co-expression network. In WGCNA analysis, the soft thresholding power (β = 4) was select to ensure a relatively balanced scale independence and mean connectivity (Figures 3A,B). No outlier samples were excluded by sample cluster analysis (Figure 3C). In total, 14 gene modules distinct were generated through hierarchical clustering tree, in which each tree branch represented a module, and each leaf constitutes a gene in the dendrogram (Figure 3E). The correlations between modules and sarcopenia were computed. Among the modules, the turquoise module was selected as the most significant module correlated with sarcopenia (r = 0.33 and p = 0.03) (Figure 3D). Furthermore, the 1,281 genes were obtained in the turquoise module, which had the most remarkable correlation with both module and sarcopenia phenotype, the correlation was 0.49 and p value was less than 0.001, which indicated a strong relationship between these genes and clinical traits (Figure 3F).
FIGURE 3
Identification and functional enrichment analysis of hub genes
A total of 31 cuproptosis-related genes (CRGs) were summarized from previous published studies. To identify the hub genes associated with sarcopenia and cuproptosis, a venn diagram was utilized to screen the genes between the DEGs, WGCNA and CRGs. As shown in Figure 4A, a total of four overlapping genes were extracted as the hub genes involved in sarcopenia and cuproptosis, which were PDHA1, DLAT, PDHB and NDUFC1. Further, the functional enrichment analysis was conducted to explore the functions and pathways of the four hub genes via clusterProfiler R package. Pathway enriched by KEGG were mainly related to RNAs with citrate cycle, glucagon and metabolic pathway, HIF-signaling pathway, retrograde endocannabinoid signaling, and several ageing related degenerative diseases (Figure 4B). GO enrichment analysis revealed that the four overlapping genes were enriched in multiple biological processes (BP) of metabolism, including ribonucleotide metabolic process, glucose metabolic process, generation of precursor metabolites and energy, oxidation reduction process and oxoacid metabolic process. With regards to cellular components (CC), the hub genes were associated with mitochondrion, mitochondrial matrix, mitochondrial pyruvate dehydrogenase complex and respiratory chain activity. Concerning molecular functions (MF), they were involved in oxidoreductase activity, catalytic activity and NADH dehydrogenase activity (Figure 4C).
FIGURE 4
Construction and verification of multigene diagnostic model
The expression levels of the four hub genes were demonstrated via box plots. As is showed in Figure 5A, the expression levels of PDHA1 (p = 2.4e-3), DLAT (p = 2.2e-3), PDHB (p = 3.3e-3) and NDUFC1 (p = 1.4e-3) in the sarcopenia group were significantly decreased than these in the healthy controls. Then, the ROC curve analysis of multigene prediction model of sarcopenia was constructed via logistic regression. To evaluate the sensitivity and specificity for the prediction model, the area under the curve (AUC) values of the four hub genes were calculated with the value of 0.83, which indicated the high predictive value of the hub genes for sarcopenia (Figure 5B). To further confirm the utility of this model, the diagnostic value of the four hub genes was validated in two independent datasets subsequently, with AUC values of 0.87 in GSE111006 dataset and 0.64 in GSE167186 dataset (Figures 5C,D). The sensitivity and specificity indicated the four genes had good diagnostic efficiency in sarcopenia.
FIGURE 5
GSEA analysis of hub genes
The samples were divided into high/low expression groups of PHDA1, DLAT, PDHB and NDUFC1 for GSEA analysis, respectively. GSEA function and pathway analysis indicated that the four hub genes were significantly associated with several neurodegenerative diseases, which included Huntington diseases, Parkinson diseases, Alzheimer diseases. In addition, pathways related to several biological metabolic processes, such as oxidative phosphorylation, peroxisome, and pyrimidine metabolism, were enriched in the high expression levels of hub genes (Figure 6).
FIGURE 6
Immune cell infiltration and correlation analysis
To investigate the immune microenvironment of sarcopenia, the relative abundance of several immune cell subtypes was calculated according to the transcriptome of all samples. The result of CIBERORT indicated that M2 macrophages, CD4+ memory T cells, and resting mast cells were the main infiltrative immune cells (Figure 7A). The infiltration of immune cells was compared between sarcopenia and healthy control groups (Figure 7B). In comparison with healthy control group, the sarcopenia group had higher fraction of M0 macrophages (p = 0.04), however, the proportions of other immune cell subsets showed no significant difference between the two groups. Subsequently, we evaluated the correlation between the expression of hub genes and infiltration levels of different cell subsets (Figure 7C). The results illustrated that PDHA1, DLAT, PDHB and NDUFC1 were negatively correlated with M0 macrophages and resting NK cells, were positively correlated with CD8+ T cells. Among them, NDUFC1 had the strongest association with M0 macrophages (r = −0.41, p = 0.0085). Next, the correlation between the immune cell populations was analyzed (Figure 7D). There was a strong negative correlation between M2 macrophages and monocytes (r = −0.65). Besides, the negative association between M2 macrophages and resting CD4+T cells was obviously observed (r = −0.59). The resting CD4+T cells were positively correlated with both activated NK cells and monocytes (r = 0.41).
FIGURE 7
Prediction of potential therapeutic drugs
To explore the potential anti-sarcopenia medicines, drugs that may target on the cuproptosis-related hub genes were screened based on the DGIdb database. According to the findings, metformin hydrochloride was found to be the therapeutic medicine of NDUFC1, while no potential drugs was found for the other three genes (PDHA1, DLAT and PDHB). The 3D chemical structure of metformin hydrochloride was obtained from PubChem database for drug docking analysis. Simultaneously, the structure of target protein NDUFC1were retrieved from PCSB PDB database (PDB ID: 5XTC). Then, the molecular docking with receptor and ligand was carried out through ChimeraX 1.5 software (Figure 8). The docking binding energy was −43.968 kcal/mol, which displayed an excellent binding affinity of drug to the target protein. In order to investigate the binding mechanism of the small molecule metformin to the NDUFC1 protein, Autodock and PyMOL tools were used to perform molecular docking simulations. The binding energy between protein NDUFC1 and metformin is −5.3 kcal/mol, which revealed metformin exhibit a strong affinity for NDUFC1. In the two-dimensional plot, hydrogen bonds are depicted as green dashed lines (Figure 9). The results showed that metformin formed two hydrogen bonds with the 106th SER amino acid residue of the 5j3j protein, with bond lengths of 3.10 Å and 2.80 Å, as well as two hydrogen bonds with the 286th THR amino acid residue of the same protein, with bond lengths of 2.79 Å and 3.10 Å. These interactions, which were revealed through the analysis of the crystal structure, were critical for the stable binding of metformin to NDUFC1. These findings provide novel insight into the molecular interactions that enable metformin to modulate the activity of target protein.
FIGURE 8
FIGURE 9
Discussion
Copper is an essential catalytic cofactor in the human body, which plays a crucial role in several biological processes ranging from oxygen metabolism to cell growth (). In general, intracellular copper is maintained at relatively low concentration, the overload of copper can cause cytotoxicity and ultimately lead to cuproptosis, which is a recently identified form of regulated cell death (RCD) (; Tsvetkov et al., 2022; Wang et al., 2022b). Distinct from the traditional oxidative stress-related cell death such as apoptosis, necroptosis and ferroptosis, cuproptosis is triggered by mitochondrial proteotoxic stress when excessive copper tightly binds to lipoylated proteins of tricarboxylic acid (TCA) cycle, followed by the aggregation of lipoylated mitochondrial enzymes in TCA and loss of Fe-S cluster proteins (; ; Tsvetkov et al., 2022). The current findings have indicated that the copper imbalance contributed to skeletal muscle mitochondrial defects, which accelerate the development of muscle atrophy and weakness, and is a major initiator of sarcopenia (). However, the connections between cuproptosis and sarcopenia remains unknown. Therefore, we aimed to elucidate the specific mechanism of cuproptosis in this aging muscle phenotypes. In addition, gene signatures were used to identify the predictive model and potential therapeutic strategy with drugs to combat sarcopenia.
In this study, we comprehensively analyzed the expression profiles of muscle tissue between sarcopenia patients and healthy individuals. Totally, 902 DEGs were identified to be associated with the development of sarcopenia. WGCNA was employed to explore the highly correlated mRNA with clinical traits. To further obtain the cuproptosis-related hub genes, we took intersections of DEGs, key WGCNA modules, and cuproptosis related genes. A total of four possible hub genes of sarcopenia were screened, namely, PDHA1, DLAT, PDHB, and NDUFC1. The GO enrichment analysis revealed that these hub genes were predominantly involved in the energy and nucleotide metabolism related biological processes. Similarly, GO cellular component analysis showed the high enrichment of these genes for mitochondrial structures and mitochondrial respiration. The GO molecular function terms also highlighted their important roles in oxidoreductase activity, catalytic activity as well as NADH dehydrogenase activity, which is consistent with the perspective of earlier studies that dysfunctional mitochondria play a key role in initiating and exacerbating sarcopenia (). It is well acknowledged that mitochondrial dysfunction is one the hallmarks of aging with the presence of decline in mitochondrial respiration and increase in oxidative stress (). In our study, the GSEA analysis confirmed that the four hub genes were also involved in the age-related neurodegenerative diseases and cellular mitochondrial functions in skeletal muscles. This finding is in line with the previous evidence that the aging myoblasts isolated from old individuals showed reduced mitochondrial respiration with increased glycolysis, and the alternations in energy metabolism in muscle tissues during aging may be the leading causes of unhealthy aging and aging-related diseases ().
Previous studies have suggested that these genes were functionally enriched in processed related to metabolism and cell activities such as cellular respiration and oxidation-reduction process (Yang et al., 2021). KEGG analysis showed that they were involved in biological processes related to mitochondrial functions in skeletal muscles, including central carbon metabolism, citrate cycle, glycolysis, and oxidative phosphorylation (). In this study, we also discovered that the four hub genes were enriched in the age-related diseases, such as Alzheimer’s disease and Parkinson’s diseases. Previous researches have suggested that cuproptosis plays a pathogenic role in the progression of neurodegenerative diseases in the context of aging (; ). Cognitive impairment and declining physical performance resulting from neurodegenerative diseases have been linked to an increased prevalence and severity of sarcopenia (Tay et al., 2018; Yazar et al., 2018; Yun et al., 2021). It is widely recognized that mitochondria play a crucial role in regulating energy metabolism, cell differentiation and death. Mitochondria function in muscle cells and motor neurons deteriorate with aging (; ), and crosstalk between these two compartments may promote the death of both ().
The four hub genes (PDHA1, DLAT, PDHB, and NDUFC1) were included to construct a diagnostic model for sarcopenia. GSE111006 dataset was subsequently utilized to validate the predictive model. The AUC values in the ROC analysis were 0.83 in the training cohort and 0.87 in the external validation cohort, demonstrating the reliable and efficient prediction performance of model. According to the results of ROC analysis, PDHA1, DLAT, PDHB, and NDUFC1 were selected as the potential promising diagnostic biomarkers of sarcopenia. These four hub genes, are involved in metabolic pathways known as cellular respiration, playing crucial roles in cellular energy metabolism. PDHA1 and PDHB are major subunits of the pyruvate dehydrogenase complex, which can convert pyruvate, the product of glycolysis, into acetyl-CoA during oxidative phosphorylation. Pyruvate can enter the citric acid cycle to generate ATP, the primary energy source for the cells (). PDHA1 and PDHB were downregulated in sarcopenia patients, which indicated the dysfunction of mitochondrial respiration in sarcopenia patients. Accumulating evidence has implicated that mitochondrial dysfunction may underlie a growing number of neurodegenerative diseases and even the process of ageing (; ; ). In this study, we also found that expression levels of DLAT was downregulated in the sarcopenia group. DLAT is an enzyme component of pyruvate dehydrogenase complex in mitochondria that acts as an acetyltransferase catalyzing the transfer of an acetyl group from acetyl-CoA to lipoamide, which ensures the energy generated from pyruvate to be properly utilized (). The mutations in DLAT genes lead to energy deficiency in the cell, has been demonstrated to propagate various human diseases, such as movement and metabolic disorders (). NDUFC1 is an essential component of NADH dehydrogenase (Complex I) in the electron transport chain required for the structural and functional stability of the complex. NDUFC1, a key gene related to cuproptosis, is also downregulated among the sarcopenia subjects. NDUFC1 plays a crucial role in the transfer of electrons from NADH to ubiquinone, contributing to maintain the energy balance in the cells. According to the previous evidence, the knockdown of NDUFC1 can inhibit complex I function and enhance ROS generation (), which may promote the accumulation of mtDNA mutations. The mtDNA damage to sarcopenia and corresponding impaired function may be implicated in the significantly reduced complex I activity in aging skeletal muscles (; ). The present study inferred the significance of mitochondrial dysfunction in the progression of sarcopenia, and the four cuproptosis-related genes, PDHA1, DLAT, PDHB, and NDUFC1, may serve as the potential targets for further therapeutic intervention. Currently, the precise mechanism of the aberrant expression of the four hub genes in promoting sarcopenia progression remains to be addressed.
The intramuscular immune cells may be involved in the progression of sarcopenia (Villalta et al., 2009; Wang et al., 2022a). In this study, we further analyzed the immune cell infiltration and the correlation of hub genes with immune factors. For most immune cell infiltration, there were no significant differences except for the quantity of M0 macrophages. The correlation between PDHA1, DLAT, PDHB, NDUFC1 and immune cells suggested that these hub genes were associated with CD8+ T cells, M0 macrophages, and resting NK cells. Evidence has indicated that the immune senescence plays a role in the onset of sarcopenia (). The immune system undergoes senescence with drastic changes, characterized by thymic involution, the accumulation of senescent T cells, and impaired function of innate immune cells, such as macrophages and NK cells. Age-dependent thymic atrophy significantly reduces the output of thymocytes from thymus, ultimately leading to impairment of T cell-mediated immunity and an increase in the proportion of CD8+ T cells, a major cause of sarcopenia (). Macrophages that infiltrate skeletal muscle are responsible for the regeneration and inflammatory responses of skeletal muscle tissues. M0 macrophages are non-differentiated immune cells that can polarize into both pro-inflammatory M1 macrophages and anti-inflammatory M2 macrophages (). In aged muscle, macrophages immunophenotyping reveals divergent muscle polarization. The dysregulated macrophage response during recovery from disuse corresponded to impaired muscle growth in old adults (). Consequently, the impaired activity of immunity in old age may contribute to both skeletal muscle atrophy and inflammatory responses. Additional research is required to fully comprehend the intricate interplay between immune cells and skeletal muscle cells.
Molecular docking, currently performed in structure-based virtual screening campaigns via rationalizing ligands activity towards a targeted protein, is a powerful tool in drug discovery (). In this study, the molecular docking results suggested that metformin could be a potential therapeutic drug for sarcopenia by targeting NDUFC1. Metformin, a well-established anti-diabetes drug for decades, has recently shown the anti-aging effects extend longevity in different organism (). Metformin has been identified as an inhibitor of mitochondrial oxidative phosphorylation through inhibiting complex I, which can reduce the production of ROS and DNA damage in senescent cells (). Metformin partially inhibiting mitochondrial respiratory complex I activity has been demonstrated to prevent mitochondrial dysfunction (). A growing body of epidemiological evidence has indicated that T2DM patients receiving metformin had less decline in walking speed than those not taking (). Similarly, the T2DM patients using metformin lost less muscle mass than controls (), suggesting metformin may play a beneficial role in sarcopenia (). In this study, the drug analysis showed that metformin is a potential therapeutic agent on sarcopenia via binding to NDUFC1, which may have significant implications for possible therapeutic interventions on sarcopenia, however, the precise mechanism by which metformin act remains further research.
As with all studies, this study has several limitations. First, the transcriptome datasets of sarcopenia were available in only a limited number. In the present study, we performed the analyses with one training dataset and one validation set, which is the main limitation of the study. Second, despite we performed detailed analysis with application of bioinformatics technology, such as differential gene expression analysis, WGCNA, and biological pathway enrichment analysis, the key genes that may drive the progression of sarcopenia still be missed. Third, the causal mechanism of cuproptosis-related hub genes and immune cell infiltration that underlie these associations and lead to sarcopenia, still needed a large number of experiments to demonstrate.
Conclusion
Our findings indicate the cuproptosis-related genes, PDHA1, DLAT, PDHB, and NDUFC1 are the diagnostic biomarker for sarcopenia, which provide evidence concerning the role of cuproptosis in sarcopenia. We also discover that metformin, the novel anti-aging drug, may be a potential therapeutic agent for the treatment of sarcopenia.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
HJ initiated and designed the study. YZ, XC, and SG carried out the bioinformatics and statistical analyses. YZ and QL drew figures and drafted the manuscript. YZ, YL, and HY edited the manuscript. All authors contributed to the study and approved the manuscript.
Funding
This work was supported by grants from Shanghai Pudong New District Health and Family Planning Commissions (No. PW2022B-13, No. PW2019E-4).
Acknowledgments
The datasets in this study were retrieved from GEO database. We would like to thank all those who contributed and participated in establishment and management of the database.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1136763/full#supplementary-material
References
1
Alvarez-CalderonF.GregoryM. A.Pham-DanisC.DeRyckereD.StevensB. M.ZaberezhnyyV.et al (2015). Tyrosine kinase inhibition in leukemia induces an altered metabolic state sensitive to mitochondrial perturbations. Clin. Cancer Res.21 (6), 1360–1372. 10.1158/1078-0432.Ccr-14-2146
2
AlwayS. E.MohamedJ. S.MyersM. J. (2017). Mitochondria initiate and regulate sarcopenia. Exerc Sport Sci. Rev.45 (2), 58–69. 10.1249/jes.0000000000000101
3
BandmannO.WeissK. H.KalerS. G. (2015). Wilson's disease and other neurological copper disorders. Lancet Neurol.14 (1), 103–113. 10.1016/s1474-4422(14)70190-5
4
BermanH. M.WestbrookJ.FengZ.GillilandG.BhatT. N.WeissigH.et al (2000). The protein Data Bank. Nucleic Acids Res.28 (1), 235–242. 10.1093/nar/28.1.235
5
BurtéF.CarelliV.ChinneryP. F.Yu-Wai-ManP. (2015). Disturbed mitochondrial dynamics and neurodegenerative disorders. Nat. Rev. Neurol.11 (1), 11–24. 10.1038/nrneurol.2014.228
6
CalvaniR.JosephA. M.AdhihettyP. J.MiccheliA.BossolaM.LeeuwenburghC.et al (2013). Mitochondrial pathways in sarcopenia of aging and disuse muscle atrophy. Biol. Chem.394 (3), 393–414. 10.1515/hsz-2012-0247
7
ChaeS.KimS. J.Do KooY.LeeJ. H.KimH.AhnB. Y.et al (2018). A mitochondrial proteome profile indicative of type 2 diabetes mellitus in skeletal muscles. Exp. Mol. Med.50 (9), 129–214. 10.1038/s12276-018-0154-6
8
ChenJ.JiangY.ShiH.PengY.FanX.LiC. (2020). The molecular mechanisms of copper metabolism and its roles in human diseases. Pflugers Arch.472 (10), 1415–1429. 10.1007/s00424-020-02412-2
9
ChenL.MinJ.WangF. (2022). Copper homeostasis and cuproptosis in health and disease. Signal Transduct. Target Ther.7 (1), 378. 10.1038/s41392-022-01229-y
10
ChoM. R.LeeS.SongS. K. (2022). A review of sarcopenia pathophysiology, diagnosis, treatment and future direction. J. Korean Med. Sci.37 (18), e146. 10.3346/jkms.2022.37.e146
11
CobineP. A.BradyD. C. (2022). Cuproptosis: Cellular and molecular mechanisms underlying copper-induced cell death. Mol. Cell82 (10), 1786–1787. 10.1016/j.molcel.2022.05.001
12
CooperC.DereW.EvansW.KanisJ. A.RizzoliR.SayerA. A.et al (2012). Frailty and sarcopenia: Definitions and outcome parameters. Osteoporos. Int.23 (7), 1839–1848. 10.1007/s00198-012-1913-1
13
Cruz-JentoftA. J.SayerA. A. (2019). Sarcopenia. Lancet393 (10191), 2636–2646. 10.1016/s0140-6736(19)31138-9
14
De BuyserS. L.PetrovicM.TaesY. E.ToyeK. R.KaufmanJ. M.LapauwB.et al (2016). Validation of the FNIH sarcopenia criteria and SOF frailty index as predictors of long-term mortality in ambulatory older men. Age Ageing45 (5), 602–608. 10.1093/ageing/afw071
15
DentE.MorleyJ. E.Cruz-JentoftA. J.AraiH.KritchevskyS. B.GuralnikJ.et al (2018). International clinical practice guidelines for sarcopenia (ICFSR): Screening, diagnosis and management. J. Nutr. Health Aging22 (10), 1148–1161. 10.1007/s12603-018-1139-9
16
FreshourS. L.KiwalaS.CottoK. C.CoffmanA. C.McMichaelJ. F.SongJ. J.et al (2021). Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res.49 (D1), D1144–d1151. 10.1093/nar/gkaa1084
17
GengY.Hernández VillanuevaA.OunA.Buist-HomanM.BlokzijlH.FaberK. N.et al (2020). Protective effect of metformin against palmitate-induced hepatic cell death. Biochim. Biophys. Acta Mol. Basis Dis.1866 (3), 165621. 10.1016/j.bbadis.2019.165621
18
GottliebR. A.BernsteinD. (2016). Mitochondrial remodeling: Rearranging, recycling, and reprogramming. Cell Calcium60 (2), 88–101. 10.1016/j.ceca.2016.04.006
19
GromadzkaG.TarnackaB.FlagaA.AdamczykA. (2020). Copper dyshomeostasis in neurodegenerative diseases-therapeutic implications. Int. J. Mol. Sci.21 (23), 9259. 10.3390/ijms21239259
20
HanF.LiuJ.ChuH.CaoD.WuJ.FuH.et al (2022). Knockdown of NDUFC1 inhibits cell proliferation, migration, and invasion of hepatocellular carcinoma. Front. Oncol.12, 860084. 10.3389/fonc.2022.860084
21
HitachiK.HondaM.TsuchidaK. (2022). The functional role of long non-coding RNA in myogenesis and skeletal muscle atrophy. Cells11 (15), 2291. 10.3390/cells11152291
22
HoodD. A.MemmeJ. M.OliveiraA. N.TrioloM. (2019). Maintenance of skeletal muscle mitochondria in health, exercise, and aging. Annu. Rev. Physiol.81, 19–41. 10.1146/annurev-physiol-020518-114310
23
HuangX.ZhouS.TóthJ.HajduA. (2022). Cuproptosis-related gene index: A predictor for pancreatic cancer prognosis, immunotherapy efficacy, and chemosensitivity. Front. Immunol.13, 978865. 10.3389/fimmu.2022.978865
24
IslamM. T. (2017). Oxidative stress and mitochondrial dysfunction-linked neurodegenerative disorders. Neurol. Res.39 (1), 73–82. 10.1080/01616412.2016.1251711
25
JohnsonM. L.RobinsonM. M.NairK. S. (2013). Skeletal muscle aging and the mitochondrion. Trends Endocrinol. Metab.24 (5), 247–256. 10.1016/j.tem.2012.12.003
26
KulkarniA. S.GubbiS.BarzilaiN. (2020). Benefits of metformin in attenuating the hallmarks of aging. Cell Metab.32 (1), 15–30. 10.1016/j.cmet.2020.04.001
27
LaiY.LinC.LinX.WuL.ZhaoY.LinF. (2022). Identification and immunological characterization of cuproptosis-related molecular clusters in Alzheimer's disease. Front. Aging Neurosci.14, 932676. 10.3389/fnagi.2022.932676
28
LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma.9, 559. 10.1186/1471-2105-9-559
29
LeeC. G.BoykoE. J.Barrett-ConnorE.MiljkovicI.HoffmanA. R.Everson-RoseS. A.et al (2011). Insulin sensitizers may attenuate lean mass loss in older men with diabetes. Diabetes Care34 (11), 2381–2386. 10.2337/dc11-1032
30
LeeC. G.SchwartzA. V.YaffeK.HillierT. A.LeBlancE. S.CawthonP. M.et al (2013). Changes in physical performance in older women according to presence and treatment of diabetes mellitus. J. Am. Geriatr. Soc.61 (11), 1872–1878. 10.1111/jgs.12502
31
LiS. R.BuL. L.CaiL. (2022). Cuproptosis: Lipoylated TCA cycle proteins-mediated novel cell death pathway. Signal Transduct. Target Ther.7 (1), 158. 10.1038/s41392-022-01014-x
32
LiberzonA.BirgerC.ThorvaldsdóttirH.GhandiM.MesirovJ. P.TamayoP. (2015). The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst.1 (6), 417–425. 10.1016/j.cels.2015.12.004
33
LinS.HuangH.LingM.ZhangC.YangF.FanY. (2022). Development and validation of a novel diagnostic model for musculoskeletal aging (sarcopenia) based on cuproptosis-related genes associated with immunity. Am. J. Transl. Res.14 (12), 8523–8538.
34
LiuH. (2022). Pan-cancer profiles of the cuproptosis gene set. Am. J. Cancer Res.12 (8), 4074–4081.
35
López-OtínC.BlascoM. A.PartridgeL.SerranoM.KroemerG. (2023). Hallmarks of aging: An expanding universe. Cell186 (2), 243–278. 10.1016/j.cell.2022.11.001
36
MarzettiE.LeesH. A.WohlgemuthS. E.LeeuwenburghC. (2009). Sarcopenia of aging: Underlying cellular mechanisms and protection by calorie restriction. Biofactors35 (1), 28–35. 10.1002/biof.5
37
MassiminoE.IzzoA.RiccardiG.Della PepaG. (2021). The impact of glucose-lowering drugs on sarcopenia in type 2 diabetes: Current evidence and underlying mechanisms. Cells10 (8), 1958. 10.3390/cells10081958
38
MéneretA.GarcinB.FrismandS.LannuzelA.MarianiL. L.RozeE. (2021). Treatable hyperkinetic movement disorders not to Be missed. Front. Neurol.12, 659805. 10.3389/fneur.2021.659805
39
MijnarendsD. M.KosterA.ScholsJ. M.MeijersJ. M.HalfensR. J.GudnasonV.et al (2016). Physical activity and incidence of sarcopenia: The population-based AGES-reykjavik study. Age Ageing45 (5), 614–620. 10.1093/ageing/afw090
40
MurrayP. J. (2017). Macrophage polarization. Annu. Rev. Physiol.79, 541–566. 10.1146/annurev-physiol-022516-034339
41
NelkeC.DziewasR.MinnerupJ.MeuthS. G.RuckT. (2019). Skeletal muscle as potential central link between sarcopenia and immune senescence. EBioMedicine49, 381–388. 10.1016/j.ebiom.2019.10.034
42
NewmanA. M.LiuC. L.GreenM. R.GentlesA. J.FengW.XuY.et al (2015). Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods12 (5), 453–457. 10.1038/nmeth.3337
43
NgoiN. Y.LiewA. Q.ChongS. J. F.DavidsM. S.ClementM. V.PervaizS. (2021). The redox-senescence axis and its therapeutic targeting. Redox Biol.45, 102032. 10.1016/j.redox.2021.102032
44
OhH. J.JinH.NahS. Y.LeeB. Y. (2021). Gintonin-enriched fraction improves sarcopenia by maintaining immune homeostasis in 20- to 24-month-old C57BL/6J mice. J. Ginseng Res.45 (6), 744–753. 10.1016/j.jgr.2021.07.006
45
PestronkA.KeelingR.ChoksiR. (2017). Sarcopenia, age, atrophy, and myopathy: Mitochondrial oxidative enzyme activities. Muscle Nerve56 (1), 122–128. 10.1002/mus.25442
46
PettersenE. F.GoddardT. D.HuangC. C.MengE. C.CouchG. S.CrollT. I.et al (2021). UCSF ChimeraX: Structure visualization for researchers, educators, and developers. Protein Sci.30 (1), 70–82. 10.1002/pro.3943
47
PinziL.RastelliG. (2019). Molecular docking: Shifting paradigms in drug discovery. Int. J. Mol. Sci.20 (18), 4331. 10.3390/ijms20184331
48
ReidyP. T.McKenzieA. I.MahmassaniZ. S.PetrocelliJ. J.NelsonD. B.LindsayC. C.et al (2019). Aging impairs mouse skeletal muscle macrophage polarization and muscle-specific abundance during recovery from disuse. Am. J. Physiol. Endocrinol. Metab.317 (1), E85–e98. 10.1152/ajpendo.00422.2018
49
RuizL. M.LibedinskyA.ElorzaA. A. (2021). Role of copper on mitochondrial function and metabolism. Front. Mol. Biosci.8, 711227. 10.3389/fmolb.2021.711227
50
RygielK. A.GradyJ. P.TurnbullD. M. (2014). Respiratory chain deficiency in aged spinal motor neurons. Neurobiol. Aging35 (10), 2230–2238. 10.1016/j.neurobiolaging.2014.02.027
51
SchaapL. A.van SchoorN. M.LipsP.VisserM. (2018). Associations of sarcopenia definitions, and their components, with the incidence of recurrent falling and fractures: The longitudinal aging study amsterdam. J. Gerontol. A Biol. Sci. Med. Sci.73 (9), 1199–1204. 10.1093/gerona/glx245
52
SchiaffinoS.ReggianiC. (2011). Fiber types in mammalian skeletal muscles. Physiol. Rev.91 (4), 1447–1531. 10.1152/physrev.00031.2010
53
SebastiánD.Acín-PérezR.MorinoK. (2016). Mitochondrial health in aging and age-related metabolic disease. Oxid. Med. Cell Longev.2016, 5831538. 10.1155/2016/5831538
54
SuY.ClaflinD. R.HuangM.DavisC. S.MacphersonP. C. D.RichardsonA.et al (2021). Deletion of neuronal CuZnSOD accelerates age-associated muscle mitochondria and calcium handling dysfunction that is independent of denervation and precedes sarcopenia. Int. J. Mol. Sci.22 (19), 10735. 10.3390/ijms221910735
55
TangD.ChenX.KroemerG. (2022). Cuproptosis: A copper-triggered modality of mitochondrial cell death. Cell Res.32 (5), 417–418. 10.1038/s41422-022-00653-7
56
TatsutaT.LangerT. (2008). Quality control of mitochondria: Protection against neurodegeneration and ageing. Embo J.27 (2), 306–314. 10.1038/sj.emboj.7601972
57
TayL.LeungB. P.WeeS.TayK. S.AliN.ChanM.et al (2018). Association of nutrition and immune-endocrine dysfunction with muscle mass and performance in cognitively impaired older adults. Arch. Gerontol. Geriatr.75, 20–27. 10.1016/j.archger.2017.11.008
58
TosettiB.BrodesserS.BrunnA.DeckertM.BlüherM.DoehnerW.et al (2020). A tissue-specific screen of ceramide expression in aged mice identifies ceramide synthase-1 and ceramide synthase-5 as potential regulators of fiber size and strength in skeletal muscle. Aging Cell19 (1), e13049. 10.1111/acel.13049
59
TsvetkovP.CoyS.PetrovaB.DreishpoonM.VermaA.AbdusamadM.et al (2022). Copper induces cell death by targeting lipoylated TCA cycle proteins. Science375 (6586), 1254–1261. 10.1126/science.abf0529
60
VerdijkL. B.SnijdersT.DrostM.DelhaasT.KadiF.van LoonL. J. (2014). Satellite cells in human skeletal muscle; from birth to old age. Age (Dordr)36 (2), 545–547. 10.1007/s11357-013-9583-2
61
VillaltaS. A.NguyenH. X.DengB.GotohT.TidballJ. G. (2009). Shifts in macrophage phenotypes and macrophage competition for arginine metabolism affect the severity of muscle pathology in muscular dystrophy. Hum. Mol. Genet.18 (3), 482–496. 10.1093/hmg/ddn376
62
VisserM.GoodpasterB. H.KritchevskyS. B.NewmanA. B.NevittM.RubinS. M.et al (2005). Muscle mass, muscle strength, and muscle fat infiltration as predictors of incident mobility limitations in well-functioning older persons. J. Gerontol. A Biol. Sci. Med. Sci.60 (3), 324–333. 10.1093/gerona/60.3.324
63
WangY.XiaoJ.SuzekT. O.ZhangJ.WangJ.ZhouZ.et al (2012). PubChem's BioAssay database. Nucleic Acids Res.40, D400–D412. 10.1093/nar/gkr1132
64
WangY.WelcS. S.Wehling-HenricksM.KongY.ThomasC.Montecino-RodriguezE.et al (2022a). Myeloid cell-specific mutation of Spi1 selectively reduces M2-biased macrophage numbers in skeletal muscle, reduces age-related muscle fibrosis and prevents sarcopenia. Aging Cell21 (10), e13690. 10.1111/acel.13690
65
WangY.ZhangL.ZhouF. (2022b). Cuproptosis: A new form of programmed cell death. Cell Mol. Immunol.19 (8), 867–868. 10.1038/s41423-022-00866-1
66
YangY.YanJ.FanX.ChenJ.WangZ.LiuX.et al (2021). The genome variation and developmental transcriptome maps reveal genetic differentiation of skeletal muscle in pigs. PLoS Genet.17 (11), e1009910. 10.1371/journal.pgen.1009910
67
YazarT.YazarH. O.ZayimoğluE.ÇankayaS. (2018). Incidence of sarcopenia and dynapenia according to stage in patients with idiopathic Parkinson's disease. Neurol. Sci.39 (8), 1415–1421. 10.1007/s10072-018-3439-6
68
YunJ. H.KimD. H.ChangM. C. (2021). A simple bedside exercise method to enhance lower limb muscle strength in moderate alzheimer's disease patients with sarcopenia. Healthc. (Basel)9 (6), 680. 10.3390/healthcare9060680
69
ZhangZ.ZengX.WuY.LiuY.ZhangX.SongZ. (2022). Cuproptosis-related risk score predicts prognosis and characterizes the tumor microenvironment in hepatocellular carcinoma. Front. Immunol.13, 925618. 10.3389/fimmu.2022.925618
70
ZhaoJ.GuoS.SchrodiS. J.HeD. (2022). Cuproptosis and cuproptosis-related genes in rheumatoid arthritis: Implication, prospects, and perspectives. Front. Immunol.13, 930278. 10.3389/fimmu.2022.930278
Summary
Keywords
sarcopenia, cuproptosis, mitochondrial metabolism, metformin, bioinformatics
Citation
Zhu Y, Chen X, Geng S, Li Q, Li Y, Yuan H and Jiang H (2023) Identification of the cuproptosis-related hub genes and therapeutic agents for sarcopenia. Front. Genet. 14:1136763. doi: 10.3389/fgene.2023.1136763
Received
03 January 2023
Accepted
28 February 2023
Published
17 March 2023
Volume
14 - 2023
Edited by
Jiang-An Yin, University of Zurich, Switzerland
Reviewed by
Amit Kumar, Albert Einstein College of Medicine, United States
Zhiyong Shao, Fudan University, China
Updates
Copyright
© 2023 Zhu, Chen, Geng, Li, Li, Yuan and Jiang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hua Jiang, huajiang2013@tongji.edu.cn
This article was submitted to Genetics of Aging, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.