The Causal Effects of Lipid Profiles on Sleep Apnea
Hongyi Tang
Qing Zhou
Fu Zheng1 
Tong Wu
Yi-Da Tang
Jiuhui Jiang- 1Department of Orthodontics, Peking University School and Hospital of Stomatology, Beijing, China
- 2Department of Cardiology, State Key Laboratory of Cardiovascular Disease, National Center for Cardiovascular Diseases, Fuwai Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 3Department of Cardiology, Graduate School of Peking Union Medical College, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 4Department of Cardiology and Institute of Vascular Medicine, Peking University Third Hospital, Beijing, China
- 5Key Laboratory of Molecular Cardiovascular Science, Ministry of Education, Beijing, China
Introduction: Observational studies have suggested that lipid profiles were associated with risk of sleep apnea (SA). However, the specific lipid types and whether this relationship has a causal effect are uncertain. This study conducted two-sample Mendelian randomization (MR) and multivariable Mendelian randomization (MVMR) to investigate the potential causal relationship between lipid profiles and risk of SA.
Materials and Methods: We used the largest genome-wide association study (GWAS) on European participants on the UK Biobank. After a rigorous single nucleotide polymorphism screening process to remove confounding effects, we performed MR and MVMR to explore the causal relationship between lipid profiles and SA risk.
Results: Both MR and MVMR showed causal effects of increased triglyceride on SA risk [MR: per 10 units, odds ratio (OR): 1.0156; 95% CI: 1.0057–1.0257; P value = 0.002; MVMR: per 10 units, OR: 1.0229; 95% CI: 1.0051–1.0411; P value = 0.011]. The sensitivity analysis including Cochran’s Q test, MR-Egger intercept, and MR pleiotropy residual sum and outlier (MR-PRESSO) test indicated that our findings were robust. The causal effects of triglyceride on SA did not change after adjusting for potential confounders (obesity, age, sex, and airway obstruction).
Conclusion: Genetically increased triglyceride levels have independent causal effects on risk of sleep apnea without the confounding effects of obesity, suggesting that lowering triglyceride concentrations may help to reduce the risk of sleep apnea.
Introduction
Sleep apnea (SA) is defined as the absence of inspiratory airflow for at least 10 s during sleep (1) and is associated with cardiovascular disease (2), stroke (3), diabetes (4), increased traffic accidents (5), lost workdays (6), and even death (7). Thus, prevention of SA is necessary, as it might cause huge social and personal burden.
Lipid profiles are an important component of human metabolism. Recent investigations have focused on the potential link between lipid profiles and SA (8–11). Obstructive sleep apnea (OSA) is the most prevalent type of SA (2). A meta-regression analysis reported that patients with OSA appeared to have elevated levels of total cholesterol, triglyceride and low-density lipoprotein (LDL), and decreased level of high-density lipoprotein (HDL) (8). A twin study demonstrated that the co-occurrence of OSA and hypertriglyceridemia had a genetic influence (9). Another cross-sectional study reported that the association between dyslipidemia and SA was limited to severe patients (10). Despite previous studies suggesting there was an association between lipid profiles and SA, the results were mainly based on observational studies and meta-analyses, which were prone to systematic bias and could not explore the causal relationship.
The causal relationship between lipid profiles and SA needs to be further explored. There are challenges in carrying out this exploration. The causal relationship between lipid profiles and SA might be confounded by other factors. Obesity is associated with lipid profiles (12) and is also considered to be a major factor contributing to SA (11, 13). This makes it difficult to explain whether the association is a causal relationship between lipid profiles and SA or the relationship comes from their common link with obesity. Additionally, common lipid profile test items include triglyceride, LDL, HDL, apolipoprotein A-1 (ApoA-1), and apolipoprotein B (ApoB). Their mutual influences may also interfere with the exploration of specific components significantly related to SA and lipid profiles. Mendelian randomization (MR) has become a reliable method to estimate for causal relationship (14) and is less prone to conventional confounding issues because of random assortment of alleles (15). Recently, multivariable Mendelian randomization (MVMR) has emerged as a method that allows for simultaneous assessment of separate but correlated exposures (16).
In this study, we extracted genetic instruments using summary statistics from a large-scale genome-wide association study (GWAS) on lipid profiles (triglyceride, LDL, HDL, ApoA-1, and ApoB) and SA. We conducted two-sample MR to investigate the causal relationship between lipid profiles and risk of SA. We carried out MVMR to explore more strong evidence for an independent causal effect of each component of lipid profiles on SA.
Materials and Methods
Selection of Data Sources
In this study, we performed two-sample MR as well as MVMR analysis of lipid profiles (triglyceride, LDL, HDL, ApoA-1, and ApoB) and SA using summary statistics from a large-scale genome-wide association study (GWAS). All original GWASs received ethical permission from corresponding ethics committees, and related participants signed informed consent (17). However, ethics permission was not required for the present study because it was derived from summary statistical data.
Genetic datasets of lipid profiles (triglyceride, HDL, LDL, ApoA-1, and ApoB) were retrieved from the largest GWAS (sample size: triglyceride: N = 441,016, HDL: N = 403,943, LDL: N = 440,546, ApoA-1: N = 393,193, and ApoB: N = 439,214), whose participants were of white European ancestry, on the UK Biobank (UKB) (17, 18). We selected all single nucleotide polymorphisms (SNPs) that reached genome-wide significance (P < 5 × 10–8), removed relative SNPs with stringent linkage disequilibrium (LD) (19) (default of the LD_clump function: R2 = 0.001, clumping window = 10,000), and obtaining the primary instrumental variables: triglyceride, HDL, LDL, ApoA-1, and ApoB. The primary instrumental variables, triglyceride, HDL, LDL, ApoA-1, and ApoB, consisted of 313, 362, 177, 299, and 198 SNPs, respectively. Because a large number of studies have shown that obesity is one of the main factors causing SA (11, 13, 20, 21) and that the occurrence of obesity is closely related to lipid metabolism (12, 22), we pruned all SNPs related to obesity from the primary instrumental variables to remove the confounding effects of obesity. We searched the traits of all the preliminary screened SNPs from PhenoScanner, which is a database of human genotype-phenotype associations.1 The following traits were selected as potential confounders because they are associated with obesity: body mass index, whole body fat mass, weight, and body fat percentage. SNPs that were not found on PhenoScanner were also removed (no search results). After the removal of SNPs associated with the potential confounders, 255, 305, 158, 237, and 169 SNPs remained in the instrumental variables triglyceride, HDL, LDL, ApoA-1, and ApoB, respectively. The screened SNPs associated with the potential confounders can be viewed online in Supplementary Table 1.
We used the data of SA from the MR-data database (23).2 We searched using “sleep apnea” as the keyword, selecting data that used participants of European descents in the database. We chose the largest sample size datasets of SA in the MRC Integrative Epidemiology Unit (MRC-IEU) consortium3 from the UKB (18) (dataset ID: ukb-b-7853, N = 463,010, diagnosis: sleep apnea). Apnea-hypopnea index (AHI) is defined as the frequency of obstructive or mixed apnea or hypopnea per hour and is the disease-defining threshold for laboratory diagnosis of sleep apnea obtained by laboratory polysomnography test (24, 25). The diagnosis of SA is determined by AHI > 5/h for adults (26). Pediatric SA is diagnosed as AHI ≥ 1/h or obstructive hypoventilation, manifested by Pa CO2 > 50 mm Hg for 25% of sleep time, coupled with snoring, paradoxical thoracoabdominal movement, or flattening of nasal airway pressure waveform (27).
We extracted the SNPs of the instrumental variables of lipid profiles (triglyceride, HDL, LDL, ApoA-1, and ApoB) from the dataset of sleep apnea and excluded palindromic SNPs, acquiring 106, 120, 63, 94, and 59 SNPs, respectively. The whole process of extracting SNPs is shown in Figure 1. The final extracted SNPs can be viewed online in Supplementary Table 2.
Figure 1. Whole process of extracting single nucleotide polymorphisms (SNPs).
Mendelian Randomization Analysis
We conducted five independent two-sample MR analyses to evaluate the potential causal relationship between lipid profiles (triglyceride, HDL, LDL, ApoA-1, and ApoB) and risk of SA. MR analysis is based on the following three assumptions (28, 29): (1) the genetic variables must be closely related to the exposure, (2) the variables must affect the outcome only through their effects on the exposure, and (3) the variables must be independent of any confounding factors of the association between the exposure and the outcome (Figure 2).
Figure 2. Schematic representation of two-sample Mendelian randomization.
We performed an inverse variance-weighted (IVW) meta-analysis of each Wald ratio for each SNP (15). Applying fixed effects IVW requires the assumption that there is no heterogeneity or horizontal pleiotropy in the SNPs (15). We also used maximum likelihood method to estimate the causal effect. The method maximized the likelihood given the SNP-exposure and the SNP-outcome effects directly and assumed a linear regression between the exposure and the outcome (15). Leave-one-out, scatter, forest, and funnel plots were produced to illustrate the results.
We conducted multivariate Mendelian randomization (MVMR) to eliminate the interaction among the selected exposures (triglyceride, HDL, LDL, ApoA-1, and ApoB) that could be confounding factors (16) (Figure 3). For the MVMR analysis, we excluded any SNPs with LD (R2 > 0.001, clumping window = 10,000) and SNPs being palindromic with intermediate allele frequencies, retaining 433 SNPs. In addition, it should be considered that SA is also closely related to age, gender, and airway obstruction. Therefore, after finding the specific lipid type with significant causal effects on SA from the lipid profiles, we conducted a multivariable Mendelian randomization analysis to analyze its effects on SA adjusted for age, gender, and airway obstruction. As no available GWAS on age could be found, we used telomere length as substitute for age, because a wide range of studies have manifested that telomere length is largest at birth and decreases with age, and that it is a biomarker of aging (30). Therefore, the traits of GWAS representing age, gender, and airway were telomere length, genetic sex, and chronic obstructive airway disease, respectively. The GWASs mentioned above were from the MR-data database (see text footnote 2) (23). The dataset IDs of the three GWASs are as follows: ieu-b-4879 (telomere length), ukb-d-is_female (sex), and ukb-b-13447 (obstructive airway disease). All the analyses were performed using R version 4.0.5 (RStudio).
Figure 3. Schematic representation of multivariable Mendelian randomization.
Sensitivity Analysis
To further assess the robustness of the findings, we performed a Cochran’s Q test to assess the heterogeneity between the genetic variables (31). Horizontal pleiotropy occurs when genetic instruments are associated with more than one independent biological pathway, which can lead to a violation of the fundamental assumption of the MR study (i.e., the variable is not related to the outcome other than via the confounding factors) (32). We conducted an MR–Egger intercept test to evaluate the horizontal pleiotropy (33). An MR pleiotropy residual sum and outlier (MR-PRESSO) test was performed to provide outlier removal, verify the results, and assess the horizontal pleiotropy (34). We applied F statistic to evaluate instrument strength (35). F statistic was calculated with the formula (35) F = (N − k − 1)/N × R2/(1 − R2) (N = sample size, k = the number of selected SNPs, R2 represents the phenotype variance induced by the SNPs.) When R2 is not available, we used the formula R2 = 2 × MAF × (1 − MAF) × (beta/SD)2 (beta = the effect value of the genetic variant in the exposure (35), MAF = the effect allele frequency (30), SD (standard deviation) = SE × √N, SE = the standard error of the genetic variant in the exposure, N = sample size).
Results
Mendelian Randomization Estimates
The evaluation for the association between lipid profiles and SA risk is shown in Table 1. The two-sample MR analysis found that increased triglyceride levels showed a causal association with risk of SA based on inverse-variance weighting. The OR (odds ratio) indicated that a 10-unit increase in triglyceride was causally associated with a 1.56% increase in SA risk (Table 1, triglyceride: N = 106 SNPs, OR: 1.0156, 95% CI: 1.0057–1.0257, P = 0.002). No causal relationship was found between the other lipids (HDL, LDL, ApoA-1, and ApoB) and SA risk (Table 1, HDL: N = 120 SNPs, OR: 1.0031, 95% CI: 0.9931–1.0132, P = 0.544; LDL: N = 63 SNPs, OR: 1.0043, 95% CI: 0.992–1.0168, P = 0.494; ApoA-1: N = 94 SNPs, OR: 1.0051, 95% CI: 0.9957–1.0146, P = 0.289; ApoB: N = 59 SNPs, OR: 1.0031, 95% CI: 0.9908–1.0156, P = 0.623). The maximum likelihood method showed similar results (Table 1, triglyceride: OR = 1.0157, 95% CI: 1.0059–1.0256, P = 0.001; HDL: OR = 1.0031, 95% CI: 0.9939–1.0124, P = 0.508; LDL: OR = 1.0044, 95% CI: 0.9929–1.016, P = 0.456; ApoA-1: OR = 1.0051, 95% CI: 0.9961–1.0141, P = 0.267; ApoB: OR = 1.0031, 95% CI: 0.9907–1.0156, P = 0.623). In the MVMR analysis, increased triglyceride still showed a causal effect on SA risk and the OR of the MVMR indicated that a 10-unit increase in triglyceride was causally associated with a 2.29% increase in SA risk (Table 2, triglyceride, OR: 1.0229, 95% CI: 1.0051–1.0411, P = 0.011; HDL, OR: 1.0128, 95% CI: 0.9708–1.0565, P = 0.556; LDL, OR: 1.004, 95% CI: 0.9463–1.0651, P = 0.895; ApoA-1, OR: 0.9902, 95% CI: 0.9498–1.0323, P = 0.642; ApoB, OR: 0.9866, 95% CI: 0.932–1.0444, P = 0.641). In the MVMR analysis related to age, sex, and airway obstruction, the causal effects of triglyceride on SA adjusted for age, gender, and upper airway obstruction were still significant (triglyceride adjusted for telomere length, OR: 1.0142, 95% CI: 1.005–1.0234; triglyceride adjusted for genetic sex, OR: 1.0169, 95% CI: 1.0082–1.0256; triglyceride adjusted for obstructive airway disease, OR: 1.0171, 95% CI: 1.0086–1.0256). The scatter plots illustrated the causal relationship between triglyceride and SA risk (Figure 4). The leave-one-out plots of triglyceride illustrated that even if any SNP in the instrumental variables was removed, the rest of the data could still achieve a significant causal effect (Figure 5). The forest plots showed causal effects of triglycerides on risk of SA, and all the funnel plots showed no asymmetry (Supplementary Figures 1, 2). In addition, the F statistics of all the instrumental variables were greater than 10.
Table 1. Mendelian randomization estimates lipid profiles and sleep apnea.
Table 2. Multivariable Mendelian randomization estimates lipid profiles and sleep apnea.
Figure 4. Scanner plot. (A) Triglyceride, (B) high-density lipoprotein, (C) low-density lipoprotein, (D) apolipoprotein A-1, and (E) apolipoprotein B.
Figure 5. Leave-one-out plot. (A) Triglyceride, (B) high-density lipoprotein, (C) low-density lipoprotein, (D) apolipoprotein A-1, and (E) apolipoprotein B.
Sensitivity Analysis
The Q variable showed no heterogeneity for the analysis of triglyceride, HDL, LDL, ApoA-1, and ApoB and SA (Table 3, MR-Egger Q, triglyceride: P = 0.427, HDL: P = 0.068, LDL: P = 0.159, ApoA-1: P = 0.224, and ApoB: P = 0.659). The MR-Egger intercept did not find significant horizontal pleiotropy for all the analyses mentioned above (Table 3, Egger intercept, triglyceride: P = 0.086, HDL: P = 0.458, LDL: P = 0.705, ApoA-1: P = 0.611, and ApoB: P = 0.949). The MR-PRESSO showed no outliers in all the instrumental variables. The estimates of the MR-PRESSO found causal effects of increased triglycerides on SA risk, in agreement with the results of the IVW and maximum likelihood methods (Table 4, triglyceride: P = 0.019, HDL: P = 0.494, LDL: P = 0.702, ApoA-1: P = 0.292, and ApoB: P = 0.597).
Table 3. The estimates of Q test and MR-Egger intercept.
Table 4. Mendelian randomization pleiotropy residual sum and outlier test (MR-PRESSO) estimates between lipid profiles and SA.
Discussion
In this study, we applied both two-sample Mendelian randomization (MR) and multivariable Mendelian randomization (MVMR) to demonstrate a causal relationship between genetically increased triglyceride and risk of sleep apnea (SA). The estimates of heterogeneity and horizontal pleiotropy showed that our results were robust. The finding of the MR-PRESSO was in agreement with the inverse variance weighting (IVW) and maximum likelihood, further strengthening the robustness of our results. The MR-PRESSO estimates found no outliers, which indicated the stability of our instrumental variables. The F statistics of all the instrumental variables were greater than 10, showing a strong association between the SNPs and the exposure with sufficient statistical power (35).
Many previous studies have demonstrated the association between risk of SA and lipid profiles including LDL, HDL and triglyceride (8–10, 36, 37). Continuing from previous relative studies, our findings showed causal effects of triglycerides on risk of SA on the basis of MR analysis. To our knowledge, this study is the first to investigate the relationship between lipid profiles and SA risk from the perspective of causality.
Our study revealed that the specific lipid causing SA was triglyceride, which is consistent with a variety of previous studies. For example, a study evaluating the association between triglyceride-glucose index (TyG index) and OSA found that higherTyG index was independently associated with increased OSA risk (OR = 3.348, 95% CI = 1.081–10.372, P < 0.05) (38). TyG was calculated using the following equation: ln [fasting triglyceride (mg/dl) × fasting glucose (mg/dl)/2]. The study also suggested a predictive role of TyG in OSA onset. Another study evaluated the visceral adipose tissue (VAT) and subcutaneous adipose tissue (SAT) of 41 participants with OSA (22 men and 19 women) and 39 controls (20 men and 19 women). In the men, apnea was associated with VAT, whereas in the women it was associated with subcutaneous, visceral, and total fats (39). These findings suggested the need for sex-specific therapeutic strategies such as reduction of visceral fat through exercise or pharmacological treatment in men and weight loss in women. In another study, Lebkuchen et al. used metabolomic and lipidomic strategies to select potential biomarkers for OSA. From 22 lipids initially selected, glycerophosphoethanolamines, sphingomyelin, and lyso-phosphocholines proved to be best associated with OSA and were considered to be potential early biomarkers in OSA screening (40). In terms of mechanisms linking lipids with OSA, a study assessed the AHI, inflammatory factors, and precise body fat of 392 adolescents to explore the mediation effects of systemic inflammation in the linkage between visceral adiposity and incident OSA. The study found that 42% of the association between visceral fat and OSA in the adolescents was mediated by interleukin-6 (IL-6) (P = 0.03), and that 82% of the association was mediated by C-reactive protein (CRP) (P = 0.01), suggesting that inflammation is a mediator in the causal relationship between body fat and SA (41). The studies mentioned above conducted various laboratory tests and manifested the link between lipids and SA, which makes our conclusion more convincing and robust, and pointed out directions for future research.
The underlying mechanisms linking triglyceride and SA can be classified into two parts: anatomical pathways and pathophysiological pathways (Figure 6). Triglyceride is any of a class of compounds that consist of an ester of glycerol with three fatty acids. Most natural fats and oils are triglycerides. When blood triglyceride level is elevated, there is an excessive triglyceride deposit on various parts of the human body. The accumulation of triglyceride from blood is the key point in the linkage between blood triglyceride level and SA. Three kinds of triglyceride accumulation constitute the anatomical pathways. (1) Oral and pharyngeal triglyceride accumulations narrow the airway directly (42–45). The fat tissue under the mandible, in the tongue, soft palate, and uvula largely compresses the anatomical space for respiration in the oral cavity and pharynx. Patients are able to compensate for the upper airway narrowing by increasing the activity of upper airway muscles during wakefulness. However, this protective effect is lost during sleep because of relaxation of muscles (42). (2) Abdominal fat accumulation including subcutaneous fat and visceral fat that surround visceral organs in the abdominal cavity markedly increases abdominal pressure, leading to reduction in lung volume. Reduction of lung volume may decrease longitudinal tracheal traction forces and pharyngeal wall tension, which is predisposed to narrowing of the airway (46). (3) Neck fat accumulation increases the collapsibility of the upper airway, decreases the efficiency of dilator muscle contraction (47), and leads to sarcopenia, denervation, and skeletal muscle dysfunction. These changes are also linked with upper airway narrowing and SA (48).
Figure 6. Underlying mechanisms linking triglyceride with sleep apnea.
Visceral and vascular triglyceride accumulations are crux of pathophysiological pathways. Higher levels of visceral fat have been observed in both obese and non-obese men with OSA, even when matched with controls for age and BMI (39, 49). When blood triglyceride level is elevated, adipocytes in visceral fat tissues grow quite large before dividing (50), and the size of lipid droplets in adipocytes increases a lot. The change makes visceral adipose tissues have more glucocorticoid and adrenergic and androgen receptors, and become more susceptible to glucocorticoid, catecholamine, and androgen signaling (51). In addition, adipocytes with abnormal large lipid droplets uptake glucose at a high rate and develop to insulin resistance (52). Insulin resistance represents an important mechanism of cardiometabolic dysfunction, and SA has long been considered as a common epiphenomenon of cardiometabolic dysfunction. In addition, along with weight gain, adipocytes themselves begin secreting tumor necrosis factor α (TNF-α). TNF-α stimulates surrounding endothelial cells to produce monocyte chemoattract protein-1 (MCP-1), which promotes macrophage recruitment. Macrophages secret many kinds of cytokines including IL-6 and IL-1. As a result, plasma CRP and other kinds of acute phase reactive proteins are elevated, triggering systemic inflammation (53, 54). Inflammatory response plays a causal role in metabolic dysregulation and subsequent development of SA. Excessive blood triglyceride can also deposit in blood vessels directly, especially where endothelia are destroyed. This kind of deposition is a main trigger of atherosclerotic plaque formation, leading to atherosclerosis and arterial stiffness. Atherosclerosis exacerbates organ ischemia and hypoxia, which is one of the main causes of various cardiovascular diseases. Again, cardiometabolic dysregulation is triggered in this pathway, and then SA happens (55, 56).
We have clarified how triglycerides play a causal role in the pathogenesis of SA and why LDL, HDL, ApoA-1, and ApoB have no causal effects on SA needs to be explained. The reasons can be listed as follows: 1. Triglycerides are hydrophobic, non-polar neutral molecules. This kind of structure makes triglycerides water-insoluble and easy to get together and deposit in extracellular fluid. However, HDL, LDL, ApoA-1, and ApoB have hydrophilic polar groups. HDL and LDL have a central hydrophobic core of non-polar lipids, primarily cholesterol esters and triglycerides. The hydrophobic core is surrounded by a hydrophilic membrane consisting of phospholipids, free cholesterol, and apolipoproteins directed outward (57). ApoA-1 and ApoB also have hydrophilic amino acid sequences (58). The hydrophilic parts make them soluble in the salt-water-based internal environment. Therefore, even if the levels are elevated, it is hard for them to deposit and accumulate in the oral cavity, pharynx, neck, and abdomen to narrow the airway, and neither can they form fat tissues to cause metabolic dysregulation like triglycerides. In brief, their hydrophilic structures suggest that they are not able to have an obvious direct causal effect on SA like triglycerides. 2. The function of LDL and HDL is transporting cholesterol instead of triglycerides. ApoA-1 and ApoB are major apolipoproteins of HDL and LDL, respectively (59). Fluctuation in levels of LDL and HDL has a direct influence on the level of cholesterol but hardly affects triglycerides. In terms of lipid component, triglycerides only account for 12% of total lipid content of both LDL and HDL, while cholesterol esters account for 59 and 40% of total lipid content of LDL and HDL, respectively (60). Cholesterol, the major core lipid of LDL and HDL, is an essential structural component of the cell membrane, where it is required to establish proper membrane permeability and fluidity (61). In addition, cholesterol is an important component for the manufacture of steroid hormones, bile acids, and vitamin D (62). Cholesterol also has hydrophilic group (-OH) and is not as easy as to accumulate as triglycerides. Therefore, the physiological function of HDL, LDL, ApoA-1, and ApoB suggests that they are unlikely to cause SA indirectly through triglyceride pathways.
From the perspective of diet, the content of triglycerides in the body is directly related to the intake of exogenous glucose and fats, which suggests that reducing the risk of SA can be achieved by decreasing the intake of exogenous triglycerides. Relative studies reported that the optimal diet consisting of not more than 50–60% carbohydrate sources, comprising mostly complex carbohydrates such as whole grains and fruits and vegetables, was an effective management to lower triglyceride concentrations (63). Besides optimal diet, regular aerobic exercise, avoiding alcohol abuse, and drug treatment including fibrates, omega-3 fatty acids, niacin, and statins were also important factors in the prevention and treatment of hypertriglyceridemia (63).
It should be noted that obesity is a major factor leading to the occurrence of SA, which has been confirmed by a large number of studies (11, 64, 65). Furthermore, a study on an MR analysis reported causal effects of obesity on risk of SA (66). Many studies reported a strong association between lipids levels and obesity (67, 68). Two reviews pointed out that obesity was an important factor increasing the prevalence of hypertriglyceridemia (67, 68). Therefore, we considered the impact of obesity in the causal effects of triglycerides on the risk of SA. Before the MR analysis, we searched all the traits of SNPs in the primary instrumental variables from PhenoScanner and excluded all SNPs related to obesity. The removed SNPs are shown online in Supplementary Table 1. Through this procedure, we demonstrated that the causal impact of triglycerides on risk of SA was significant without the confounding effects of obesity.
The interaction between lipid profiles should also be considered. MVMR analysis removes interfering factors effectively. By MVMR, we also found significant causal effects of genetically increased triglycerides on risk of SA, but no significant causal relationship was found with HDL, LDL, ApoA-1, and ApoB. The results provided stronger support for the causal relationship we found. In addition, we also considered that SA might be influenced by age, sex, and airway obstruction. The related multivariable Mendelian randomization analysis indicated that the causal effects of triglycerides on SA were not affected by the confounding factors, namely, age, sex, and airway obstruction.
One limitation of our study was that although the sensitivity analysis did not detect horizontal pleiotropy, confounding and pleiotropic factors may still exist (31). For example, we only removed SNPs related to the trait obesity, because we only found a causal relationship between obesity and SA risk in previous studies based on MR analysis (66). However, SA might be associated with many systemic diseases (2–4), and the incorporated SNPs in the lipid profile instruments might still impact SA risk through other ways. Another limitation was that our study was mainly based on a GWAS, whose participants were European, on the UK-Biobank, and different types of studies based on other populations need to be further analyzed.
Conclusion
The present study demonstrates the independent causal effects of genetically increased triglycerides on risk of sleep apnea without the confounding effects of obesity by two-sample Mendelian randomization and multivariable Mendelian randomization, and suggests that lowering triglyceride concentration contributes to the reduction in risk of sleep apnea.
Data Availability Statement
The original contributions presented in this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
HT contributed in software application, resources, and wrote the manuscript. QZ contributed in methodology and data curation. FZ and TW contributed in preliminary investigation. JJ contributed to editing the manuscript and supervision. Y-DT contributed to editing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by National Key R&D Program of China (2020YFC2004705), Research Unit of Medical Science Research Management/Basic and Clinical Research of Metabolic Cardiovascular Diseases from Chinese Academy of Medical Sciences (2021RU003), National Natural Science Foundation of China (81825003 and 91957123), and Science and Technology Project of Xicheng District Finance (XCSTS-SD2021-01).
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.
Acknowledgments
We would like to thank the UK Biobank for the GWAS summary statistics.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.910690/full#supplementary-material
Supplementary Figure 1 | Forest plot. (A) Triglyceride, (B) high-density lipoprotein, (C) low-density lipoprotein, (D) apolipoprotein A-1, and (E) apolipoprotein B.
Supplementary Figure 2 | Funnel plot. (A) Triglyceride, (B) high-density lipoprotein, (C) low-density lipoprotein, (D) apolipoprotein A-1, and (E) apolipoprotein B.
Supplementary Table 1 | Removal of single nucleotide polymorphisms (SNPs) associated with potential confounders.
Supplementary Table 2 | Final extracted SNPs.
Footnotes
- ^ http://www.phenoscanner.medschl.cam.ac.uk/
- ^ https://gwas.mrcieu.ac.uk/
- ^ http://www.bristol.ac.uk/integrative-epidemiology/
References
1. Javaheri S, Barbe F, Campos-Rodriguez F, Dempsey JA, Khayat R, Javaheri S, et al. Sleep apnea: types, mechanisms, and clinical cardiovascular consequences. J Am Coll Cardiol. (2017) 69:841–58. doi: 10.1016/j.jacc.2016.11.069
2. Nieto FJ, Young TB, Lind BK, Shahar E, Samet JM, Redline S, et al. Association of sleep-disordered breathing, sleep apnea, and hypertension in a large community-based study. Sleep Heart Health Study JAMA. (2000) 283:1829–36. doi: 10.1001/jama.283.14.1829
3. Lau DH, Nattel S, Kalman JM, Sanders P. Modifiable risk factors and atrial fibrillation. Circulation. (2017) 136:583–96. doi: 10.1161/CIRCULATIONAHA.116.023163
4. Strausz S, Havulinna AS, Tuomi T, Bachour A, Groop L, Mäkitie A, et al. Obstructive sleep apnoea and the risk for coronary heart disease and type 2 diabetes: a longitudinal population-based study in Finland. BMJ Open. (2018) 8:e022752. doi: 10.1136/bmjopen-2018-022752
5. Terán-Santos J, Jiménez-Gómez A, Cordero-Guevara J. The association between sleep apnea and the risk of traffic accidents. Cooperative Group Burgos-Santander. N Engl J Med. (1999) 340:847–51. doi: 10.1056/NEJM199903183401104
6. Sjösten N, Vahtera J, Salo P, Oksanen T, Saaresranta T, Virtanen M, et al. Increased risk of lost workdays prior to the diagnosis of sleep apnea. Chest. (2009) 136:130–6. doi: 10.1378/chest.08-2201
7. The National Commission on Sleep Disorders Research. Wake Up America: A National Sleep Alert. Washington, DC: US Government Printing Office (2002).
8. Nadeem R, Singh M, Nida M, Waheed I, Khan A, Ahmed S, et al. Effect of obstructive sleep apnea hypopnea syndrome on lipid profile: a meta-regression analysis. J Clin Sleep Med. (2014) 10:475–89. doi: 10.5664/jcsm.3690
9. Meszaros M, Tarnoki AD, Tarnoki DL, Kovacs DT, Forgo B, Lee J, et al. Obstructive sleep apnea and hypertriglyceridaemia share common genetic background: results of a twin study. J Sleep Res. (2020) 29:e12979. doi: 10.1111/jsr.12979
10. Martínez-Cerón E, Casitas R, Galera R, Sánchez-Sánchez B, Zamarrón E, Garcia-Sanchez A, et al. Contribution of sleep characteristics to the association between obstructive sleep apnea and dyslipidemia. Sleep Med. (2021) 84:63–72. doi: 10.1016/j.sleep.2021.05.012
11. Wu WT, Tsai SS, Shih TS, Lin MH, Chou TC, Ting H, et al. The association between obstructive sleep apnea and metabolic markers and lipid profiles. PLoS One. (2015) 10:e0130279. doi: 10.1371/journal.pone.0130279
12. Iyer A, Fairlie DP, Prins JB, Hammock BD, Brown L. Inflammatory lipid mediators in adipocyte function and obesity. Nat Rev Endocrinol. (2010) 6:71–82. doi: 10.1038/nrendo.2009.264
13. Kasasbeh A, Kasasbeh E, Krishnaswamy G. Potential mechanisms connecting asthma, esophageal reflux, and obesity/sleep apnea complex–a hypothetical review. Sleep Med Rev. (2007) 11:47–58. doi: 10.1016/j.smrv.2006.05.001
14. Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32:1–22. doi: 10.1093/ije/dyg070
15. Georgiopoulos G, Ntritsos G, Stamatelopoulos K, Tsioufis C, Aimo A, Masi S, et al. The relationship between blood pressure and risk of atrial fibrillation: a mendelian randomization study. Eur J PrevCardiol. (2021) 2021:zwab005. doi: 10.1093/eurjpc/zwab005
16. Burgess S, Thompson SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. (2015) 181:251–60. doi: 10.1093/aje/kwu283
17. Richardson TG, Sanderson E, Palmer TM, Ala-Korpela M, Ference BA, Davey Smith G, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable Mendelian randomisation analysis. PLoS Med. (2020) 17:e1003062. doi: 10.1371/journal.pmed.1003062
18. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. (2015) 12:e1001779. doi: 10.1371/journal.pmed.1001779
19. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. doi: 10.1086/519795
20. Strobel RJ, Rosen RC. Obesity and weight loss in obstructive sleep apnea: a critical review. Sleep. (1996) 19:104–15. doi: 10.1093/sleep/19.2.104
21. Meurling IJ, Shea DO, Garvey JF. Obesity and sleep: a growing concern. Curr Opin Pulm Med. (2019) 25:602–8. doi: 10.1097/MCP.0000000000000627
22. Pan DA, Hulbert AJ, Storlien LH. Dietary fats, membrane phospholipids and obesity. J Nutr. (1994) 124:1555–65. doi: 10.1093/jn/124.9.1555
23. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. (2018) 7:e34408. doi: 10.7554/eLife.34408
24. Cade BE, Hassan SM, Dashti HS, Kiernan M, Pavlova MK, Redline S, et al. Sleep apnea phenotyping and relationship to disease in a large clinical biobank. JAMIA Open. (2022) 5:ooab117. doi: 10.1093/jamiaopen/ooab117
25. Raynes-Greenow CH, Hadfield RM, Cistulli PA, Bowen J, Allen H, Roberts CL. Sleep apnea in early childhood associated with preterm birth but not small for gestational age: a population-based record linkage study. Sleep. (2012) 35:1475–80. doi: 10.5665/sleep.2192
26. Veasey SC, Rosen IM. Obstructive sleep apnea in adults. N Engl J Med. (2019) 380:1442–9. doi: 10.1056/NEJMcp1816152
27. Sateia MJ. International classification of sleep disorders-third edition: highlights and modifications. Chest. (2014) 146:1387–94. doi: 10.1378/chest.14-0970
28. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. (2008) 27:1133–63. doi: 10.1002/sim.3034
29. Egger M, Smith GD, Phillips AN. Meta-analysis: principles and procedures. BMJ. (1997) 315:1533–7. doi: 10.1136/bmj.315.7121.1533
30. Rizvi S, Raza ST, Mahdi F. Telomere length variations in aging and age-related diseases. Curr Aging Sci. (2014) 7:161–7. doi: 10.2174/1874609808666150122153151
31. Bowden J, Del Greco MF, Minelli C, Zhao Q, Lawlor DA, Sheehan NA, et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. (2019) 48:728–42. doi: 10.1093/ije/dyy258
32. Gormley M, Dudding T, Sanderson E, Martin RM, Thomas S, Tyrrell J, et al. A multivariable Mendelian randomization analysis investigating smoking and alcohol consumption in oral and oropharyngeal cancer. Nat Commun. (2020) 11:6071. doi: 10.1038/s41467-020-19822-6
33. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. (2017) 36:1783–802. doi: 10.1002/sim.7221
34. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
35. You R, Chen L, Xu L, Zhang D, Li H, Shi X, et al. High Level of Uromodulin Increases the Risk of Hypertension: a Mendelian Randomization Study. Front Cardiovasc Med. (2021) 8:736001. doi: 10.3389/fcvm.2021.736001
36. Kono M, Tatsumi K, Saibara T, Nakamura A, Tanabe N, Takiguchi Y, et al. Obstructive sleep apnea syndrome is associated with some components of metabolic syndrome. Chest. (2007) 131:1387–92. doi: 10.1378/chest.06-1807
37. McArdle N, Hillman D, Beilin L, Watts G. Metabolic risk factors for vascular disease in obstructive sleep apnea: a matched controlled study. Am J Respir Crit Care Med. (2007) 175:190–5. doi: 10.1164/rccm.200602-270OC
38. Kang HH, Kim SW, Lee SH. Association between triglyceride glucose index and obstructive sleep apnea risk in Korean adults: a cross-sectional cohort study. Lipids Health Dis. (2020) 19:182. doi: 10.1186/s12944-020-01358-9
39. Kritikou I, Basta M, Tappouni R, Pejovic S, Fernandez-Mendoza J, Nazir R, et al. Sleep apnoea and visceral adiposity in middle-aged male and female subjects. Eur Respir J. (2013) 41:601–9. doi: 10.1183/09031936.00183411
40. Lebkuchen A, Carvalho VM, Venturini G, Salgueiro JS, Freitas LS, Dellavance A, et al. Metabolomic and lipidomic profile in men with obstructive sleep apnoea: implications for diagnosis and biomarkers of cardiovascular risk. Sci Rep. (2018) 8:11270. doi: 10.1038/s41598-018-29727-6
41. Gaines J, Vgontzas AN, Fernandez-Mendoza J, Calhoun SL, He F, Liao D, et al. Inflammation mediates the association between visceral adiposity and obstructive sleep apnea in adolescents. Am J Physiol Endocrinol Metab. (2016) 311:E851–8. doi: 10.1152/ajpendo.00249.2016
42. Schwab RJ, Pasirstein M, Pierson R, Mackley A, Hachadoorian R, Arens R, et al. Identification of upper airway anatomic risk factors for obstructive sleep apnea with volumetric magnetic resonance imaging. Am J Respir Crit Care Med. (2003) 168:522–30. doi: 10.1164/rccm.200208-866OC
43. Horner RL, Shea SA, McIvor J, Guz A. Pharyngeal size and shape during wakefulness and sleep in patients with obstructive sleep apnoea. Q J Med. (1989) 72:719–35.
44. Shelton KE, Woodson H, Gay S, Suratt PM. Pharyngeal fat in obstructive sleep apnea. Am Rev Respir Dis. (1993) 148:462–6. doi: 10.1164/ajrccm/148.2.462
45. Stauffer JL, Buick MK, Bixler EO, Sharkey FE, Abt AB, Manders EK, et al. Morphology of the uvula in obstructive sleep apnea. Am Rev Respir Dis. (1989) 140:724–8. doi: 10.1164/ajrccm/140.3.724
46. Isono S. Obesity and obstructive sleep apnoea: mechanisms for increased collapsibility of the passive pharyngeal airway. Respirology. (2012) 17:32–42. doi: 10.1111/j.1440-1843.2011.02093.x
47. Deegan PC, McNicholas WT. Pathophysiology of obstructive sleep apnoea. Eur Respir J. (1995) 8:1161–78. doi: 10.1183/09031936.95.08071161
48. Poelkens F, Lammers G, Pardoel EM, Tack CJ, Hopman MT. Upregulation of skeletal muscle inflammatory genes links inflammation with insulin resistance in women with the metabolic syndrome. Exp Physiol. (2013) 98:1485–94. doi: 10.1113/expphysiol.2013.072710
49. Vgontzas AN, Papanicolaou DA, Bixler EO, Hopper K, Lotsikas A, Lin HM, et al. Sleep apnea and daytime sleepiness and fatigue: relation to visceral obesity, insulin resistance, and hypercytokinemia. J Clin Endocrinol Metab. (2000) 85:1151–8. doi: 10.1210/jcem.85.3.6484
50. Ibrahim MM. Subcutaneous and visceral adipose tissue: structural and functional differences. Obes Rev. (2010) 11:11–8. doi: 10.1111/j.1467-789X.2009.00623.x
51. Rebuffé-Scrive M, Lundholm K, Björntorp P. Glucocorticoid hormone binding to human adipose tissue. Eur J Clin Invest. (1985) 15:267–71. doi: 10.1111/j.1365-2362.1985.tb00182.x
52. Mårin P, Andersson B, Ottosson M, Olbe L, Chowdhury B, Kvist H, et al. The morphology and metabolism of intraabdominal adipose tissue in men. Metabolism. (1992) 41:1242–8. doi: 10.1016/0026-0495(92)90016-4
53. Xu H, Barnes GT, Yang Q, Tan G, Yang D, Chou CJ, et al. Chronic inflammation in fat plays a crucial role in the development of obesity-related insulin resistance. J Clin Invest. (2003) 112:1821–30. doi: 10.1172/JCI19451
54. Wellen KE, Hotamisligil GS. Obesity-induced inflammatory changes in adipose tissue. J Clin Invest. (2003) 112:1785–8. doi: 10.1172/JCI20514
55. Drager LF, Togeiro SM, Polotsky VY, Lorenzi-Filho G. Obstructive sleep apnea: a cardiometabolic risk in obesity and the metabolic syndrome. J Am Coll Cardiol. (2013) 62:569–76. doi: 10.1016/j.jacc.2013.05.045
56. Gaines J, Vgontzas AN, Fernandez-Mendoza J, Bixler EO. Obstructive sleep apnea and the metabolic syndrome: the road to clinically-meaningful phenotyping, improved prognosis, and personalized treatment. Sleep Med Rev. (2018) 42:211–9. doi: 10.1016/j.smrv.2018.08.009
57. Ginsberg HN. Lipoprotein physiology. Endocrinol Metab Clin North Am. (1998) 27:503–19. doi: 10.1016/s0889-8529(05)70023-2
58. Mahley RW, Innerarity TL, Rall SC Jr., Weisgraber KH. Plasma lipoproteins: apolipoprotein structure and function. J Lipid Res. (1984) 25:1277–94.
59. Illingworth DR. Lipoprotein metabolism. Am J Kidney Dis. (1993) 22:90–7. doi: 10.1016/s0272-6386(12)70173-7
60. Feingold KR. Introduction to lipids and lipoproteins. In: Feingold KR, Anawalt B, Boyce A, Chrousos G, de Herder WW, Dhatariya K, et al. editors. Endotext [Internet]. South Dartmouth, MA: MDText.com, Inc (2021). 2000 p.
61. Luo J, Yang H, Song BL. Mechanisms and regulation of cholesterol homeostasis. Nat Rev Mol Cell Biol. (2020) 21:225–45. doi: 10.1038/s41580-019-0190-7
62. Schade DS, Shey L, Eaton RP. Cholesterol review: a metabolically important molecule. EndocrPract. (2020) 26:1514–23. doi: 10.4158/EP-2020-0347
64. Shechter A. Obstructive sleep apnea and energy balance regulation: a systematic review. Sleep Med Rev. (2017) 34:59–69. doi: 10.1016/j.smrv.2016.07.001
65. Ischander MM, Lloyd RD Jr. Severe paediatric obesity and sleep: a mutual interactive relationship! J Sleep Res. (2021) 30:e13162. doi: 10.1111/jsr.13162
66. Strausz S, Ruotsalainen S, Ollila HM, Karjalainen J, Kiiskinen T, Reeve M, et al. Genetic analysis of obstructive sleep apnoea discovers a strong association with cardiometabolic health. Eur Respir J. (2021) 57:2003091. doi: 10.1183/13993003.03091-2020
67. Laufs U, Parhofer KG, Ginsberg HN, Hegele RA. Clinical review on triglycerides. Eur Heart J. (2020) 41:99c–109c. doi: 10.1093/eurheartj/ehz785
Keywords: genetics, sleep apnea, lipid profiles, multivariable Mendelian randomization, triglyceride
Citation: Tang H, Zhou Q, Zheng F, Wu T, Tang Y-D and Jiang J (2022) The Causal Effects of Lipid Profiles on Sleep Apnea. Front. Nutr. 9:910690. doi: 10.3389/fnut.2022.910690
Received: 01 April 2022; Accepted: 26 May 2022;
Published: 21 June 2022.
Edited by:
Justyna Godos, University of Catania, ItalyReviewed by:
Jian Yin Zou, Shanghai Jiao Tong University, ChinaYanzhong Li, Shandong University, China
Copyright © 2022 Tang, Zhou, Zheng, Wu, Tang 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: Yi-Da Tang, tangyida@bjmu.edu.cn; Jiuhui Jiang, drjiangw@163.com
†These authors have contributed equally to this work and share first authorship