Research Paper Volume 15, Issue 16 pp 8026—8043

NDUFB11 and NDUFS3 play a role in atherosclerosis and chronic stress

Yin Yang1, *, , Jing-Hui Li1, *, , Bo-Chen Yao1, , Qing-Liang Chen1, , Nan Jiang1, , Lian-Qun Wang1, , Zhi-Gang Guo1, ,

  • 1 Clinical School of Thoracic, Tianjin Medical University, Tianjin Chest Hospital, Jinnan, Tianjin 300222, P.R. China
* Equal contribution

Received: May 11, 2023       Accepted: July 18, 2023       Published: August 28, 2023      

https://doi.org/10.18632/aging.204947
How to Cite

Copyright: © 2023 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Objective: Atherosclerosis is characterized by the formation of fibrofatty plaques in the intima of arteries, resulting in thickening of the vessel wall and narrowing of the lumen. Chronic stress refers to individuals in a state of long-term chronic stress. However, the relationship between NDUFB11 and NDUFS3 and atherosclerosis and chronic stress is unclear.

Method: The atherosclerosis with chronic stress group data file was used. DEGs were screened and WGCNA was performed. Construction and analysis of PPI Network. Functional enrichment analysis, GSEA, gene expression heatmap, immune infiltration analysis and mRNA analysis were performed. CTD was used to find diseases most related to core genes. WB was performed. TargetScan was used to screen miRNAs of DEGs.

Results: 1708 DEGs were identified. According to GO analysis, they were mainly enriched in catabolic processes, organic acid metabolism processes, carboxylic acid metabolism processes. KEGG analysis showed that they were mainly enriched in metabolic pathways, fatty acid metabolism, pentose phosphate pathway, glycolysis / gluconeogenesis, fructose and mannose metabolism. Gene expression heatmap showed that the core genes (NDUFB11, NDUFS3) were lowly expressed in samples of those with atherosclerosis accompanied by chronic stress and highly expressed in the normal samples. NDUFB11 and NDUFS3 were associated with necrosis, hyperplasia, inflammation, renal disease, weight loss, memory impairment, and cognitive impairment. WB showed that the expression level of NDUFS3 in atherosclerosis and chronic stress was lower than that in control group.

Conclusions: NDUFB11 and NDUFS3 are underexpressed in atherosclerosis and chronic stress; the lower NDUFB11 and NDUFS3 levels, the worse the prognosis.

Introduction

Atherosclerosis is a chronic disease caused by the deposition of lipids, cholesterol, and other substances on the inner wall of blood vessels, as well as endothelial cell damage due to inflammatory reactions [1]. Atherosclerosis usually develops over age of 40, and men are at higher risk of developing atherosclerosis [2]. Atherosclerosis is mainly characterized by lipid deposition and inflammatory reactions on the intimal layer of blood vessels, which progressively lead to the thickening and loss of elasticity of the intimal layer and eventually the formation of arterial plaques. Atherosclerosis is a chronic progressive disease. The clinical manifestations of atherosclerosis may be asymptomatic in the early stages, but as the disease progresses, cardiovascular, cerebrovascular, peripheral vascular, renal lesions may develop [3, 4]. Chronic stress refers to a state in which an individual is stimulated by persistent psychological or physical stress over a period of time [5]. Studies have shown that young people under 20 years of age, those with high intensity, long hours and high professional responsibilities, and those with low income and educational levels are more susceptible to chronic stress [6]. Chronic stress is a chronically accumulated stress that may affect multiple aspects of physical and mental health and may affect different individuals [7]. The pathological features of chronic stress are increased or decreased secretion of neuroendocrine hormones causing changes in the cardiovascular, immune, and metabolic systems, decreased immune function or dysregulated immune responses, increased oxidative stress damaging cells and tissues, affecting energy metabolism, and increased inflammatory responses [8]. The causes of atherosclerosis and chronic stress are not clear, and genetic factors, chromosomal abnormalities, and gene fusions may contribute to this disease. Therefore, it is particularly important to deeply investigate molecular mechanisms of atherosclerosis and chronic stress.

Bioinformatics is an interdisciplinary field that involves computer science, mathematics, biology, and statistics [9]. The development of bioinformatics technology has greatly assisted biological research, accelerating the interpretation and understanding of biomolecules such as genomes, proteins, and metabolomes. Bioinformatics technology includes sequence analysis, structure analysis, functional prediction, systems biology, genomics, and proteomics. The advantages of bioinformatics technology are mainly reflected in its efficiency, accuracy, visualization, and reproducibility [10, 11].

However, the relationship between NDUFB11 and NDUFS3 and atherosclerosis and chronic stress is still unclear. Therefore, we wanted to use bioinformatics to mine core genes between atherosclerosis and chronic stress and normal tissues, perform enrichment analysis, pathway analysis. Validation of significant role of NDUFB11 and NDUFS3 in atherosclerosis and chronic stress using public datasets. And the basal cell experiment was applied to verify it.

Results

Differential gene analysis

In the present study, following the set cut-off value, 1708 DEGs were identified based on identification of differentially expressed genes in the atherosclerosis accompanying chronic stress matrix (Figure 1).

Analysis of differentially expressed genes. 1708 DEGs were identified.

Figure 1. Analysis of differentially expressed genes. 1708 DEGs were identified.

Functional enrichment analysis

DEGs

We analyzed DEGs by GO and KEGG. According to GO analysis, they were mainly enriched in catabolic processes, organic acid metabolism processes, carboxylic acid metabolism processes. Through KEGG analysis, they were mainly enriched in metabolic pathways, fatty acid metabolism, pentose phosphate pathway, glycolysis / gluconeogenesis, fructose and mannose metabolism (Figure 2A2D).

Functional enrichment analysis. (A–D) DEGs (E–H) GSEA.

Figure 2. Functional enrichment analysis. (AD) DEGs (EH) GSEA.

GSEA

A genome-wide GSEA was performed to find possible enrichment terms among the non-differentially expressed genes and to validate the differentially expressed genes. The intersection of the enriched terms with the GO KEGG enriched terms of the differentially expressed genes is shown in, and was mainly enriched in the pentose phosphate pathway, fatty acid metabolism, glycolysis / gluconeogenesis (Figure 2E2H).

Metascape enrichment analysis

Purine nucleotide metabolic process, monocarboxylate metabolic process, amino acid metabolism were found in GO enrichment items in the enrichment items of metascape (Figure 3A), meanwhile we also exported the enrichment network colored by enrichment term and p-value (Figures 3B, 3C, 4) to visualize the association and confidence representing each enrichment item.

Metascape enrichment analysis. (A) Purine nucleotide metabolic process, monocarboxylate metabolic process, amino acid metabolism were found in GO enrichment items in the enrichment items of metascape (B) output the enrichment network colored by enrichment terms (C) output the enrichment network colored by p-value.

Figure 3. Metascape enrichment analysis. (A) Purine nucleotide metabolic process, monocarboxylate metabolic process, amino acid metabolism were found in GO enrichment items in the enrichment items of metascape (B) output the enrichment network colored by enrichment terms (C) output the enrichment network colored by p-value.

Metascape enrichment analysis. Visualize the association and confidence representing each enrichment item. (A) Protein-protein interact by the Metascape. (B) The significant modules.

Figure 4. Metascape enrichment analysis. Visualize the association and confidence representing each enrichment item. (A) Protein-protein interact by the Metascape. (B) The significant modules.

WGCNA

The network topology is analyzed and the soft threshold power of WGCNA is set to 9, which was the lowest power for a scale-free topology fit index of 0.9 (Figure 5A, 5B). A hierarchical clustering tree of all genes was constructed, and 3 significant modules were generated (Figure 5C). The interaction between modules was analyzed (Figure 5D). And generated a module to phenotype correlation Heatmap (Figure 5E). GS to MM correlation scatter plot of the associated hub genes (Figure 5F).

WGCNA. (A) β=6, 0.70 (B) β=6, 163.48 (C) A hierarchical clustering tree of all genes was constructed, and 3 significant modules were generated (D) Interactions between these modules were then analyzed (E) Generated a module to phenotype correlation Heatmap (F) GS to MM correlation scatter plot of the associated hub genes.

Figure 5. WGCNA. (A) β=6, 0.70 (B) β=6, 163.48 (C) A hierarchical clustering tree of all genes was constructed, and 3 significant modules were generated (D) Interactions between these modules were then analyzed (E) Generated a module to phenotype correlation Heatmap (F) GS to MM correlation scatter plot of the associated hub genes.

Construction and analysis of protein-protein interaction (PPI) networks

The PPI network of DEGs was constructed from the STRING online database and analyzed by Cytoscape software, we selected the enrichment items in the gokegg analysis related to adrenal gland and input the list of differential genes enriched on these enrichment items into STRING, constructed PPI network (Figure 6A), hub genes were identified using two different algorithms (Figure 6B, 6C), which were merged with the Venn diagram (Figure 6D), Fourteen core genes (UQCRB, COX5B, NDUFB11, NDUFA8, NDUFS6, COX6B1, NDUFA4, UQCRQ, 2NDUFA2, NDUFB9, NDUFV2, NDUFA10, COX4I1, NDUFS3) were obtained.

Construction and analysis of protein-protein interaction (PPI) networks. (A) The PPI network (B) DMNC algorithm was adopted to identify the core genes (C) MCC algorithm was adopted to identify the core genes (D) Venn diagram.

Figure 6. Construction and analysis of protein-protein interaction (PPI) networks. (A) The PPI network (B) DMNC algorithm was adopted to identify the core genes (C) MCC algorithm was adopted to identify the core genes (D) Venn diagram.

Gene expression heatmap

A heat map of the expression of core genes in the samples was plotted (Figure 7), and we found that the core genes (NDUFB11, NDUFS3) were lowly expressed in the atherosclerosis accompanying chronic stress samples and highly expressed in the normal samples.

Gene expression heatmap. Visualized the expression Heatmap of the core genes in the samples.

Figure 7. Gene expression heatmap. Visualized the expression Heatmap of the core genes in the samples.

Immune infiltration analysis

We used the cibersort package to analyze the atherosclerosis accompanying chronic stress matrix file, and at 95% confidence, obtained the proportion results of immune cells from the full gene expression matrix (Figure 8A) and the immune cell expression Heatmap in the dataset (Figure 8B), and also performed the correlation analysis on infiltrated immune cells, resulting in a plot of co expression patterns among immune cell components (Figure 8C).

Immune infiltration analysis. (A) Obtained the proportion results of immune cells from the full gene expression matrix (B) the immune cell expression Heatmap in the dataset (C) performed the correlation analysis on infiltrated immune cells, resulting in a plot of co expression patterns among immune cell components.

Figure 8. Immune infiltration analysis. (A) Obtained the proportion results of immune cells from the full gene expression matrix (B) the immune cell expression Heatmap in the dataset (C) performed the correlation analysis on infiltrated immune cells, resulting in a plot of co expression patterns among immune cell components.

The messenger RNA (mRNA) analysis

We obtained autophagy related mRNA differentially expressed genes by extracting autophagy related genes in the matrix of atherosclerosis accompanying chronic stress, intersecting with the differential genes obtained from the screening.

We then performed KEGG analysis on these differentially expressed genes to obtain relevant enrichment items KEGG analysis of DEGs was performed to obtain related enrichment entries (Figure 9), and entered these differentially expressed genes into the STRING website to create a PPI network (Figure 10A), using two algorithms (MCC, MNC) to identify hub genes (Figure 10B), followed by a Venn diagram to obtain a Union (Figure 10C), Core genes related to autophagy mRNAs were obtained (HSPA5, CANX, LAMP1, mTOR, RAB7A, LAMP2, RAB5A, HSPA8, MAPK3).

The mRNA analysis. Performed KEGG analysis on these differentially expressed genes to obtain relevant enrichment items.

Figure 9. The mRNA analysis. Performed KEGG analysis on these differentially expressed genes to obtain relevant enrichment items.

The mRNA analysis (A) PPI network (B, C) using two algorithms (MCC, MNC) to identify hub genes (D) Venn diagram to obtain a Union.

Figure 10. The mRNA analysis (A) PPI network (B, C) using two algorithms (MCC, MNC) to identify hub genes (D) Venn diagram to obtain a Union.

CTD analysis

Core genes was entered into CTD to find diseases related to core genes. The core genes NDUFB11 and NDUFS3 were found to be associated with necrosis, hyperplasia, inflammation, renal disease, weight loss, memory impairment, and cognitive impairment (Figure 11).

CTD analysis. The core genes NDUFB11 and NDUFS3 were found to be associated with necrosis, hyperplasia, inflammation, renal disease, weight loss, memory impairment, and cognitive impairment.

Figure 11. CTD analysis. The core genes NDUFB11 and NDUFS3 were found to be associated with necrosis, hyperplasia, inflammation, renal disease, weight loss, memory impairment, and cognitive impairment.

WB

The expression level of NDUFS3, E-cadherin, BCL2, P21 and SOX9 in atherosclerosis and chronic stress were lower than that in control group. They were higher in AS+CS/NDUFS3-OE group than in AS+CS group. The expression was lower in the AS+CS/NDUFS3-OE group than in the AS+CS group.

The expression level of Vimentin, N-cadherin, Snail1, ZEB1, PI3K, p-AKT, p-mTOR and BAX in atherosclerosis and chronic stress were higher than that in control group. They were lower in AS+CS/NDUFS3-OE group than in AS+CS group. The expression was higher in the AS+CS/NDUFS3-OE group than in the AS+CS group (Figure 12).

The expression level of NDUFS3 in the atherosclerosis and chronic stress groups was lower than that in the control group.

Figure 12. The expression level of NDUFS3 in the atherosclerosis and chronic stress groups was lower than that in the control group.

The miRNAs prediction and functional annotation associated with core genes

The hub genes were entered into TargetScan to search for relevant miRNA (Table 1). The related miRNAs of NDUFS3 gene are mmu-mir-125b-5p, mmu-mir-6394, mmu-mir-6367.

Table 1. A summary of miRNAs that regulate hub genes.

GeneMIRNA
1Ndufs3mmu-miR-125b-5pmmu-miR-6394mmu-miR-6367
2Ndufb11none

Discussion

Atherosclerosis is a chronic and progressive disease in which plaques that form become occluded within the arterial wall reduce blood flow, increase blood pressure, and trigger complications such as vascular rupture, thrombosis, and embolization, leading to diseases such as glomerulosclerosis, renal insufficiency, and renal failure [12]. Chronic stress may lead to mental health problems such as mood disorders, depression, anxiety and insomnia, resulting in decreased function of the immune system and increased risk of infections, inflammatory and autoimmune diseases, metabolic disorders, peptic ulcer inflammatory bowel disease and other digestive diseases, infertility, sexual dysfunction and other reproductive diseases [13, 14]. Several molecular mechanisms are involved in the development and progression of atherosclerosis, including oxidative stress, inflammatory response, platelet activation, abnormal cholesterol metabolism, and apoptosis [1517]. The molecular mechanisms of chronic stress involve disturbances in several biological systems [1820]. The main results of this study are that NDUFB11 and NDUFS3 are underexpressed in atherosclerosis and chronic stress, and the lower the NDUFB11 and NDUFS3, the worse the prognosis.

NDUFB11 is one of the subunits of complex I of mitochondrial respiratory chain and is a key factor in the electron transfer process [21]. Mitochondrial respiratory chain complex I, contains multiple subunits and cofactors and is required for maintaining the energy metabolism of mitochondria. NDUFB11 is mainly involved in complex I assembly and function stability, but also in regulation of respiratory function, apoptosis, and oxidative stress response in mitochondria [22]. Recent studies have shown that mutations in NDUFB11 can be associated with the development of many diseases. For example, NDUFB11 deficiency was found to be a rare mitochondrial disorder known as NDUFB11 deficient mitochondriopathy. Patients with this disease usually present with a variety of symptoms, including muscle weakness, neurological disorders, and metabolic disorders, among others. NDUFB11 expression levels may also be involved in the development of other diseases, such as tumors, neurodegenerative and cardiovascular diseases, among others [23]. NDUFB11 mutations all have an established role in the mitochondrial respiratory chain [24]. The NDUFB11 gene is associated with neurogenetic disorders [25].

NDUFS3 is a subunit of mitochondrial respiratory chain complex I [26, 27]. Mitochondrial respiratory chain complex I, a complex in the mitochondria, is involved in the process of oxidative phosphorylation in cells. NDUFS3 is one of the key subunits of this complex responsible for important steps in the electron transport process, and its absence or aberrant function is associated with disorders associated with mitochondrial respiratory chain disorders, such as muscular, neurodegenerative, and metabolic diseases, among others. The study of NDUFS3 may provide insight into mitochondrial function and oxidative phosphorylation processes and shed new light on the treatment of related diseases [27, 28]. Therefore, we speculate that NDUFB11 and NDUFS3 may play important roles in intracellular energy metabolism in atherosclerosis and chronic stress.

Although this paper has carried out rigorous bioinformatics analysis, there are still some shortcomings. Animal experiments with overexpression or knockdown of the gene were not performed in this study to further verify the function.

In summary, NDUFB11 and NDUFS3 are underexpressed in atherosclerosis and chronic stress, and the lower NDUFB11 and NDUFS3, the worse the prognosis. NDUFB11 and NDUFS3 may serve as molecular targets for the precision treatment of atherosclerosis and chronic stress, thus providing a foundation of direction for mechanistic studies of atherosclerosis and chronic stress.

Materials and Methods

Atherosclerosis accompanied by chronic stress dataset

In this study, the atherosclerosis with chronic stress group data file including 3 atherosclerosis with chronic stress samples and 3 normal samples.

Screening of differentially expressed genes (DEGs)

Probe aggregation and background correction of merge matrix of atherosclerotic concomitant chronic stress using R package “limma”. P value were adjusted using Benjamini-Hochberg method. The fold change (FC) is calculated using false discovery rate (FDR). The cutoff value of DEG is p less than 0.05 and FC is greater than 1.5. And make a visual representation of the volcano.

Weighted gene co expression network analysis (WGCNA)

First of all, use de-batch and post-merge matrix of atherosclerosis accompanying chronic stress to calculate median absolute deviation (MAD) of each gene. Outlier genes and samples were removed by good sample gene method of WGCNA in R package, and scale-free co-expression network was constructed. We calculated characteristic gene differences of modules, and selected tangent line for module tree view, incorporated part of modules.

Construction and analysis of protein-protein interaction (PPI) networks

Search Tool for the Retrieval of Interacting Genes (STRING) is a search system for known and predicted PPI. STRING database also contains the predicted results using bioinformatics methods. The differential genes were input into STRING to construct PPI network and predict core genes. PPI network was visualized, core genes are predicted by Cytoscape software. First of all, we import PPI network into the Cytoscape, and then find module with the best correlation through MCODE, and MCC and MNC were used to calculate the best correlated genes. Finally, the list of core genes was obtained after visualization.

Functional enrichment analysis

Gene Ontology (GO) analysis is a computational method to evaluate gene functions and biological pathways, and it is a key step to endow sequence information with practical biological significance. Kyoto Encyclopedia of Gene and Genome (KEGG) is an online database dedicated to collecting information on genomes, molecular interaction networks, enzyme catalytic pathways, and biochemical products. The genomic information and gene function were linked, and gene function was systematically analyzed. The list of differential genes screened by Wayne map was input into KEGG rest API obtained latest KEGG Pathway gene annotation. Gene set enrichment results were obtained using R package cluster Profiler.

Metascape (http://metascape.org/) can realize cognition of gene or protein function, and can be visually exported. We used Metascape database to analyze functional enrichment of the above differential gene list and derive it.

Gene set enrichment analysis (GSEA)

The samples were divided into two groups according to the atherosclerosis accompanying chronic stress and normal groups, and the data were GSEA is based on level-specific gene probes that evaluate data from microarrays and is a way to uncover genomic expression data through fundamental knowledge. The samples were divided into atherosclerosis accompanying chronic stress and normal groups. 5 is minimum gene set and 5000 is maximum gene set, 1000 resampling times.

Gene expression heatmap

We use R-packet heatmap to map expression of core genes found in PPI network, and to visualize difference of core gene expression between atherosclerosis accompanying chronic stress and normal samples.

Immune infiltration analysis

The CIBERSORT (http://CIBERSORT.stanford.edu/) is a very common method for calculating immune cell infiltration. We applied the integrated bioinformatics method, used the CIBERSORT software package to analyze the de-batch merging matrix of atherosclerotic concomitant chronic stress, and immune cell abundance was estimated by deconvoluting the expression matrix of immune cell subtypes by linear support vector regression principle. At the same time, the samples with sufficient confidence were selected by using confidence P < 0.05as the truncation criterion.

The messenger RNA (mRNA) analysis

We intersected with the screening derived differential genes by extracting autophagy related genes in the matrix of atherosclerosis with chronic stress for identification of differentially expressed genes related to autophagy.

Comparative toxicogenomics database (CTD)

The CTD (https://ctdbase.org/) is a powerful public database that links toxicological information on chemicals, genes, phenotypes, diseases, and exposures to advance understanding of human health. CTD incorporates literature based, manually curated interactions to create a knowledge base that orchestrates chemical exposures and their biological effects across heterogeneous data across species. CTD identifies integrated chemical diseases, chemical genes, and gene-disease interactions that are used to predict disease-related gene/protein relationships as well as new associations and generate extended networks. We input core gene into CTD, find disease most related to core gene, and drew an expression difference radar plot for each gene with Excel.

WB

Western blotting, also known as immunoblotting, is a method to detect the expression of a certain protein in complex samples according to the specific binding of antigens and antibodies, and can qualitatively and semi-quantitatively analyze proteins. Total protein was extracted and the protein content was determined. After SDS-PAGE electrophoresis and membrane transfer, the protein samples were blocked with 5% skim milk for 1 h at room temperature, shaken with Tris Buffered Saline Tween at high speed on a shaker, washed for 5 min, and repeated three times. The primary antibody was added and incubated overnight at 4° C, followed by TBST shaking 3 times (5 min each time) and TBST shaking 3 times (5 min each time). The results were analyzed after chemiluminescence development.

miRNA

TargetScan (https://www.targetscan.org) can predict and analyze miRNA and target genes. Screening of miRNAs regulating central DEGs was performed using TargetScan in this study.

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

DEGs: Differentially expressed genes; WGCNA: weighted gene co-expression network analysis; PPI: protein-protein interaction; GSEA: gene set enrichment analysis; CTD: comparative toxicogenomics database; GO: gene ontology; KEGG: Kyoto Encyclopedia of Gene and Genome; FC: fold change; FDR: false discovery rate; MAD: median absolute deviation; STRING: Search Tool for the Retrieval of Interacting Genes; mRNA: messenger RNA.

Author Contributions

Yin Yang and Jing-Hui Li participated in the work concept and design of the study to performed experiments mentioned, Bo-Chen Yao and Qing-Liang Chen collected the data, Nan Jiang and Lian-Qun Wang performed the statistical analysis of the data, Zhidang Guo drafted the manuscript, Yin Yang, Jing-Hui Li, and Zhi-Gang Guo revised the main content of the manuscript. All authors read and agree on the manuscript.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Ethical Statement

The Ethics Committee of Tianjin Medical University approved the study.

Funding

Early warning method and application of adverse events in the perioperative period of off-pump coronary artery bypass grafting based on artificial intelligence data analysis (Project number: 20JCZDJC00810).

References

  • 1. Fan J, Watanabe T. Atherosclerosis: Known and unknown. Pathol Int. 2022; 72:151–60. https://doi.org/10.1111/pin.13202 [PubMed]
  • 2. Man JJ, Beckman JA, Jaffe IZ. Sex as a Biological Variable in Atherosclerosis. Circ Res. 2020; 126:1297–319. https://doi.org/10.1161/CIRCRESAHA.120.315930 [PubMed]
  • 3. Vuong JT, Stein-Merlob AF, Nayeri A, Sallam T, Neilan TG, Yang EH. Immune Checkpoint Therapies and Atherosclerosis: Mechanisms and Clinical Implications: JACC State-of-the-Art Review. J Am Coll Cardiol. 2022; 79:577–93. https://doi.org/10.1016/j.jacc.2021.11.048 [PubMed]
  • 4. Foks AC, Bot I. Preface: Pathology and Pharmacology of Atherosclerosis. Eur J Pharmacol. 2017; 816:1–2. https://doi.org/10.1016/j.ejphar.2017.10.052 [PubMed]
  • 5. Yao BC, Meng LB, Hao ML, Zhang YM, Gong T, Guo ZG. Chronic stress: a critical risk factor for atherosclerosis. J Int Med Res. 2019; 47:1429–40. https://doi.org/10.1177/0300060519826820 [PubMed]
  • 6. Novais A, Monteiro S, Roque S, Correia-Neves M, Sousa N. How age, sex and genotype shape the stress response. Neurobiol Stress. 2016; 6:44–56. https://doi.org/10.1016/j.ynstr.2016.11.004 [PubMed]
  • 7. Golbidi S, Frisbee JC, Laher I. Chronic stress impacts the cardiovascular system: animal models and clinical outcomes. Am J Physiol Heart Circ Physiol. 2015; 308:H1476–98. https://doi.org/10.1152/ajpheart.00859.2014 [PubMed]
  • 8. Radley JJ, Herman JP. Preclinical Models of Chronic Stress: Adaptation or Pathology? Biol Psychiatry. 2023; 94:194–202. https://doi.org/10.1016/j.biopsych.2022.11.004 [PubMed]
  • 9. Fu Y, Ling Z, Arabnia H, Deng Y. Current trend and development in bioinformatics research. BMC Bioinformatics. 2020 (Suppl 9); 21:538. https://doi.org/10.1186/s12859-020-03874-y [PubMed]
  • 10. Azad RK, Shulaev V. Metabolomics technology and bioinformatics for precision medicine. Brief Bioinform. 2019; 20:1957–71. https://doi.org/10.1093/bib/bbx170 [PubMed]
  • 11. Chen C, Hou J, Tanner JJ, Cheng J. Bioinformatics Methods for Mass Spectrometry-Based Proteomics Data Analysis. Int J Mol Sci. 2020; 21:2873. https://doi.org/10.3390/ijms21082873 [PubMed]
  • 12. Dong Z, Zhao Y, Chen Y, Liu Z, Song H, Li H, Shi D, Zhou C, Zhou J, Liu R. Evaluating Atherosclerosis of the Abdominal Aorta in Rabbits Using 2-D Strain Imaging. Ultrasound Med Biol. 2022; 48:2199–206. https://doi.org/10.1016/j.ultrasmedbio.2022.05.029 [PubMed]
  • 13. Cherpak CE. Mindful Eating: A Review Of How The Stress-Digestion-Mindfulness Triad May Modulate And Improve Gastrointestinal And Digestive Function. Integr Med (Encinitas). 2019; 18:48–53. [PubMed]
  • 14. Osborne MT, Shin LM, Mehta NN, Pitman RK, Fayad ZA, Tawakol A. Disentangling the Links Between Psychosocial Stress and Cardiovascular Disease. Circ Cardiovasc Imaging. 2020; 13:e010931. https://doi.org/10.1161/CIRCIMAGING.120.010931 [PubMed]
  • 15. Simion V, Zhou H, Haemmig S, Pierce JB, Mendes S, Tesmenitsky Y, Pérez-Cremades D, Lee JF, Chen AF, Ronda N, Papotti B, Marto JA, Feinberg MW. A macrophage-specific lncRNA regulates apoptosis and atherosclerosis by tethering HuR in the nucleus. Nat Commun. 2020; 11:6135. https://doi.org/10.1038/s41467-020-19664-2 [PubMed]
  • 16. Bakogiannis C, Sachse M, Stamatelopoulos K, Stellos K. Platelet-derived chemokines in inflammation and atherosclerosis. Cytokine. 2019; 122:154157. https://doi.org/10.1016/j.cyto.2017.09.013 [PubMed]
  • 17. Zhu Y, Xian X, Wang Z, Bi Y, Chen Q, Han X, Tang D, Chen R. Research Progress on the Relationship between Atherosclerosis and Inflammation. Biomolecules. 2018; 8:80. https://doi.org/10.3390/biom8030080 [PubMed]
  • 18. Zefferino R, Di Gioia S, Conese M. Molecular links between endocrine, nervous and immune system during chronic stress. Brain Behav. 2021; 11:e01960. https://doi.org/10.1002/brb3.1960 [PubMed]
  • 19. Herman JP, McKlveen JM, Ghosal S, Kopp B, Wulsin A, Makinson R, Scheimann J, Myers B. Regulation of the Hypothalamic-Pituitary-Adrenocortical Stress Response. Compr Physiol. 2016; 6:603–21. https://doi.org/10.1002/cphy.c150015 [PubMed]
  • 20. Noushad S, Ahmed S, Ansari B, Mustafa UH, Saleem Y, Hazrat H. Physiological biomarkers of chronic stress: A systematic review. Int J Health Sci (Qassim). 2021; 15:46–59. [PubMed]
  • 21. Amate-García G, Ballesta-Martínez MJ, Serrano-Lorenzo P, Garrido-Moraga R, González-Quintana A, Blázquez A, Rubio JC, García-Consuegra I, Arenas J, Ugalde C, Morán M, Guillén-Navarro E, Martín MA. A Novel Mutation Associated with Neonatal Lethal Cardiomyopathy Leads to an Alternative Transcript Expression in the X-Linked Complex I NDUFB11 Gene. Int J Mol Sci. 2023; 24:1743. https://doi.org/10.3390/ijms24021743 [PubMed]
  • 22. Reinson K, Kovacs-Nagy R, Õiglane-Shlik E, Pajusalu S, Nõukas M, Wintjes LT, van den Brandt FC, Brink M, Acker T, Ahting U, Hahn A, Schänzer A, Haack TB, et al. Diverse phenotype in patients with complex I deficiency due to mutations in NDUFB11. Eur J Med Genet. 2019; 62:103572. https://doi.org/10.1016/j.ejmg.2018.11.006 [PubMed]
  • 23. Petruzzella V, Tessa A, Torraco A, Fattori F, Dotti MT, Bruno C, Cardaioli E, Papa S, Federico A, Santorelli FM. The NDUFB11 gene is not a modifier in Leber hereditary optic neuropathy. Biochem Biophys Res Commun. 2007; 355:181–7. https://doi.org/10.1016/j.bbrc.2007.01.140 [PubMed]
  • 24. Indrieri A, Franco B. Linear Skin Defects with Multiple Congenital Anomalies (LSDMCA): An Unconventional Mitochondrial Disorder. Genes (Basel). 2021; 12:263. https://doi.org/10.3390/genes12020263 [PubMed]
  • 25. Gurok U, Bork K, Nuber U, Spörle R, Nöhring S, Horstkorte R. Expression of Ndufb11 encoding the neuronal protein 15.6 during neurite outgrowth and development. Gene Expr Patterns. 2007; 7:370–4. https://doi.org/10.1016/j.modgep.2006.07.004 [PubMed]
  • 26. Lou X, Shi H, Wen S, Li Y, Wei X, Xie J, Ma L, Yang Y, Fang H, Lyu J. Correction: A novel NDUFS3 mutation in a Chinese patient with severe Leigh syndrome. J Hum Genet. 2018; 63:1283–4. https://doi.org/10.1038/s10038-018-0517-9 [PubMed]
  • 27. Samanta D, Veerapandiyan A, Burrow TA, Gokden M. Mitochondrial Ultrastructural Defects in NDUFS3-Related Disorder. J Pediatr Neurosci. 2021; 16:299–302. https://doi.org/10.4103/jpn.JPN_182_20 [PubMed]
  • 28. Pereira CV, Peralta S, Arguello T, Bacman SR, Diaz F, Moraes CT. Myopathy reversion in mice after restauration of mitochondrial complex I. EMBO Mol Med. 2020; 12:e10674. https://doi.org/10.15252/emmm.201910674 [PubMed]