Statins for the Treatment of Pulmonary Hypertension in Patients with Chronic Obstructive Pulmonary Disease

Background: Patients with chronic obstructive pulmonary disease (COPD) are at risk for pulmonary hypertension (PH). The aim of our study was to investigate the benefit of statins for PH in patients with COPD.

Methods: The study enrolled 23 million individuals from Taiwan’s population database from January 1, 2002, to December 31, 2017. COPD patients who met the inclusion criteria were enrolled, and patients with lung cancer, less than one year of observation, specific drug therapy for PH and lung transplantation were excluded.

Results: A total of 643,131 COPD patients were included in the study, and only 12,308 patients developed PH during follow-up. Based on the inclusion and exclusion criteria, 8,577 PH patients were included in the cohort of patients with PH related to COPD for analysis. According to the definition of statin exposure, the final study population had 1,487 statin users and 7,090 statin non-users. The statin user group had a lower mortality related to PH than the non-user group (3.87 vs. 5.55 per 100 person-years, p < 0.001). The mortality rate for PH in the multivariate analysis (aHR = 0.78, 95% CI = 0.62–0.98, p = 0.046) was significantly lower for statin users than for non-users.

Conclusion: Statins seem to benefit patients with PH and COPD.

Introduction

Pulmonary hypertension (PH) is defined as a mean resting pulmonary arterial pressure exceeding 25 mm·Hg, as evaluated by right heart catheterization. In the sixth World Symposium on PH (WSPH) update from 2018, PH is divided into five groups according to clinical presentation, hemodynamic characteristics, pathophysiology, and therapeutic strategy (Simonneau et al., 2019). Patients with PH due to obstructive lung disease [e.g., chronic obstructive pulmonary disease (COPD)] or restrictive lung disease) are classified into group 3. Among older patients above the age of 65 years, group-2 PH is the most frequent diagnosis (17%), followed by group-3 PH (14%) (Cottin et al., 2005; Hurdman et al., 2013; Solidoro et al., 2015). In a Taiwan epidemiological survey, PH related to COPD was the most common form of PH (Chang et al., 2016).

Current medical treatments for PH, including supportive therapy and new specific drugs in PH related to COPD, are lacking, and the efficacy is still uncertain; only long-term oxygen therapy can improve symptoms.

Our previous study showed that the statin user group had a lower incidence rate of PH than the nonuser group (1.43 vs. 1.97 per 1,000 person-years) (Wu and Chen, 2020).

Therefore, we further performed a nationwide, population-based, retrospective cohort study to explore the association between statins and mortality among patients with COPD and PH.

Materials and Methods

Study Population

The Taiwan National health insurance program, which covers 99.6% of the 23 million people in Taiwan’s population and 93% of hospitals, clinics, and pharmacies, created one of the largest and most comprehensive medical population databases in the world. We used the full population database and multiple causes of death data from the National health insurance program. The study was conducted at the Center of Health and Welfare Data Science Centers. The medical record of the full population database was used from January 1, 2002, to December 31, 2017.

To assess whether the use of statins before PH diagnosis affects the study outcome, the study population was divided into statin users, defined as statin exposure within one year, and statin non-users, defined as no statin prescription within one year. This study was approved by the IRB of Kaohsiung Medical University Chung-Ho Memorial Hospital (KMUHIRB-EXEMPT(I)-20190032).

Inclusion Criteria

The inclusion criteria were 1) age between 40 and 90 years, 2) newly diagnosed COPD between January 1, 2002, and December 31, 2015 (according to International Classification of Diseases, Ninth Revision, Clinical Modification codes 490, 491, 492, and 496), 3) more than one inpatient diagnosis or more than two consecutive outpatient diagnoses of COPD within one year after the primary COPD diagnosis date, and 4) treatment with COPD medications indicated in outpatient claims for more than 28 days (COPD medication: long-acting beta 2 agonist, long-acting muscarinic antagonist, short-acting beta 2 agonist combined with oral methylxanthines or short-acting muscarinic antagonist combined with oral methylxanthines) (Wang et al., 2018), and 5) more than one inpatient diagnosis or emergency room admission for PH or more than two consecutive outpatient diagnoses of PH within one year. The date of the first diagnosis of PH was defined as the index date.

Exclusion Criteria

The exclusion criteria were 1) COPD with lung cancer, 2) COPD with other etiologies of PH, including human immunodeficiency virus, sleep apnea, pulmonary embolism, connective tissue disease, congenital heart defect, or portal hypertension, 4) death within 28 days after COPD diagnosis, 5) specific drug therapy for PH treatment (e.g., sildenafil, riociguat) before COPD diagnosis, 6) death within 1 year after the first PH diagnosis date, 7) lung transplantation within 1 year after the first PH diagnosis, of 8) a less than one-year observation time in the database.

Immortal Time Bias

To resolve the immortal time bias, we set a one-year confirmation period (Levesque et al., 2010). Patients who needed to start using statins and who had at least one prescription for statins within one year after the COPD diagnosis were defined as statin users. Conversely, patients who did not receive a prescription for statins within one year after the COPD diagnosis were defined as statin non-users.

Baseline Characteristics

Baseline characteristics included age, age group, sex, urbanization level, insurance premium, comorbidities, and co-medication. Comorbidities and co-medication were confirmed by medical records in the year following the PH diagnosis (i.e., one year).

To adjust for other confounders that could influence the outcome, we listed the PH treatment [e.g., oxygen therapy, digoxin, calcium channel blocker (CCB), warfarin, endothelin receptor antagonists, phosphodiesterase type 5 inhibitors, prostacyclin analogs] in the baseline characteristics.

Severity of Pulmonary Hypertension

The severity of PH is usually expressed by the WHO functional class, mean pulmonary artery pressure, or exercise capacity test, such as the six-minute walk distance. However, we could not collect those data from our database. This study used COPD severity within one year after the PH diagnosis, with oxygen therapy and heart failure (HF) used to represent PH severity. In the clinic, if patients have more severe COPD with concomitant HF, they are considered to have more severe PH. These three factors were included in the calculation of the propensity score to identify patients with the same PH severity.

Propensity-Score Matching

We used propensity-score matching to reduce selection bias. This matching method calculated a propensity score for each patient by means of dependent variables. If patients had the same propensity score, they had equal probability of developing events that we wanted to observe. The baseline characteristics of the original cohort were used to perform propensity-score matching by multivariate logistic regression. We adopted 1:1 matching to generate the statin user group and statin non-user group after follow-up.

Primary Outcome

The primary outcome was PH-related mortality, which was a composite outcome that included deaths caused by PH, deaths caused by HF, and deaths caused by respiratory exacerbation (RE). Deaths were identified using the “Multiple Cause of Death Data.” The “Multiple Cause of Death Data” contains patients’ dates of death, hospitals of death, and causes of death.

Secondary Outcome

The all-cause mortality, deaths caused by PH, deaths caused by HF, and deaths caused by RE were secondary outcomes. The cause of death was identified by International Classification of Diseases, Ten Revision, Clinical Modification code according to the “Multiple Cause 1 to 20 (MC1-MC20)” columns from Multiple Cause of Death Data.

Follow-Up Time

We set a five-year observation duration to reduce the influence of switching to different statins or refusal to take statins. The follow-up period started from one year after the index date to five years after the index date. The observation time ended at death or censor occurrence. The patients in both groups who were alive or underwent lung transplantation at the end of the observation period were censored.

Subgroup Analysis

We conducted a subgroup analysis to assess whether the association variables could affect the effectiveness of statins. Only the primary outcomes were included in the subgroup analysis.

Different Kinds of Statins

We grouped patients treated with seven kinds of statins according into the group of most commonly used statins. The seven statins were simvastatin, lovastatin, pravastatin, fluvastatin, atorvastatin, rosuvastatin, and pitavastatin.

Cumulative Defined Daily Dose

The cDDD was the total statin exposure over the whole follow-up period. The cDDD was calculated over the entire five-year observation duration. The cDDD was divided into six levels in the user group compared to the non-user group: non-user, cDDD <28, cDDD ≥28 but <90, cDDD ≥90 but <180, cDDD ≥180 but <365, cDDD ≥365 but <730, and cDDD ≥730.

Duration of Statin Use (years)

The duration of statin use was calculated by year and divided into six levels in the user group compared to the non-user group: non-user, <0.5 year, ≥0.5 year but <1 year, ≥1 year but <2 years, ≥2 years but <3 years, and ≥3 years.

Intensity of Statin Use (Cumulative Defined Daily Dose/Month)

The intensity of statin use was calculated by dividing cDDD by the duration of statin use over the whole five-year observation period. Then, statin use was divided into four categories to determine whether the intensity could affect the protective effect of statins: non-user, intensity <10, intensity ≥10 ≤ but <20, intensity ≥20.

Sensitivity Analysis

Different Observation Durations

The observation times for the sensitivity analysis were 1 year, 3, 5, 7, 9 years, and the end of the database to detect the effects of different observation durations.

Statistical Analysis

For baseline characteristics, continuous variables are presented as the means (standard deviations), and categorical variables are presented as percentages. Continuous variables were analyzed by Student’s t-test or analysis of variance, and categorical variables were analyzed by Fisher’s exact test or chi-square tests if numbers were less than 30.

The mortality of PH was assessed using a time-dependent Cox proportional hazard model to estimate the hazard ratio between the two groups. The multivariate model was adjusted for demographic characteristics, sex, comorbidity, co-medication, and COPD severity. Because of potential interactions and collinearity between characteristics, stepwise selection was performed to select the important factors to construct the regression model. All the above analyses were performed using SAS 9.4. software. The statistical significance was determined with two-tailed analyses and α = 0.05.

Results

Study Population

A total of 643,131 COPD patients were included in the study, and only 12,308 patients developed PH during follow-up. Based on the inclusion and exclusion criteria, 8,577 PH patients were included in the PH related to the COPD cohort for analysis. According to the definition of statin exposure, the final study population had 1,487 statin users and 7,090 statin non-users. Figure 1 showed the study flow chart.

FIGURE 1

FIGURE 1. Result of flow chart in the study.

Baseline Characteristics

The PH cohort contained 1,487 statin users and 7,090 non-users. The mean age of the patients in the user group (72.52 years) was lower than that of the non-user group (74.10 years, p < 0.001). The age distributions were similar between groups, with most patients over 60 years old and the majority being male. The insurance premium showed the same distribution, with more patients paying below 22,800 Taiwan Dollars and a significant difference between the two groups (p = <0.001).

Because of the indications for statins, patients with statins had a significantly higher prevalence of dyslipidemia (p < 0.001), coronary artery disease (p < 0.001), and ischemic stroke (p < 0.001). The proportions of most comorbidities were significantly higher in the user group than in the non-user group. There was no significant difference in chronic liver disease, arrhythmia, interstitial pulmonary diseases, asthma, malignancy, hemorrhagic stroke or left ventricular hypertrophy between groups.

Co-medication, except with digoxin, had the same trend for comorbidity. More patients in the non-user group used digoxin (18.16% vs. 21.28%, p = 0.007), In addition, there were no significant differences in oxygen therapy or specific drug therapy between groups. Regarding medication for COPD, more than half of the patients used methylxanthines. Regarding the severity of COPD, there were significant differences in both severe and moderate exacerbations (p < 0.001) between the two groups. The distributions of severe and moderate exacerbations showed that most patients had no exacerbation in the firth year after COPD was diagnosed, and the second highest number of patients had over two exacerbations. The detailed baseline characteristics of the PH patients before matching are presented in Table 1.

TABLE 1

TABLE 1. Baseline characteristics of PH patients before and after matching, stratified according to statin use.

To identify similar patients in the control group, the study population was subjected to 1-to-1 propensity-score matching, and age group, sex, insurance premium, urbanization, hypertension, diabetes, HF, asthma, oxygen therapy, specific drug therapy and COPD severity were used to calculate the propensity score. Before matching, the propensity scores of the user and non-user groups were significantly different (0.77 ± 0.12 vs. 0.84 ± 0.10, p < 0.001). After matching, the propensity scores in the two groups were similar, with no significant differences (0.66 ± 0.04 vs. 0.66 ± 0.04, p = 1.000).

After matching, there were still significant differences in the distributions of comorbidity and co-medication between the two groups. The Poisson regression model and Cox proportional hazard model were used to adjust the imbalanced characteristics in the subsequent analysis. The detailed baseline characteristics of the PH cohort after 1-to-1 propensity-score matching are shown in Table 1.

Mortality of Patients with Pulmonary Hypertension Related to Chronic Obstructive Pulmonary Disease

Regarding the primary outcome, during the five-year observation duration, 462 (34.87%) statin users and 647 (48.83%) statin non-users died, among whom 145 (10.94%) statin users and 210 (15.84%) statin non-users died of causes related to PH. Regarding the primary outcome, the total follow-up times of the user and non-user groups were 3,751.18 and 3,784.92 person-years, respectively. The statin user group had a lower mortality related to PH than the non-user group (3.87 vs. 5.55 per 100 person-years, p < 0.001). Regarding the secondary outcomes, most patients with mortality related to PH died due to RE (rates of mortality related to PH of the user and non-user groups: 86.21% and 79.52%), and few patients died due to PH or HF.

The Cox proportional hazard model analysis of primary outcome and secondary outcome of PH patients is summarized in Table 2. The univariate analysis (crude HR = 0.70, 95% CI = 0.56–0.86, p < 0.001) and multivariate analysis (aHR = 0.78, 95% CI = 0.62–0.98, p = 0.046) of mortality related to PH both showed significantly lower mortality for statin users than non-users. The above results suggested that statins could reduce the risk of mortality related to PH by approximately 22–30% and that they were effective in patients with PH related to COPD. Regarding secondary outcomes, there was significantly lower risk for deaths caused by PH or HF (aHR = 0.41, 95% CI = 0.23–0.74, p = 0.003) and all-cause mortality (aHR = 0.76, 95% CI = 0.67–0.87, p < 0.001).

TABLE 2

TABLE 2. Cox proportional hazard model analysis of mortality of PH patients, stratified according to statin use.

In the multivariate analysis, all baseline characteristics were included in the adjusted model. When statin users were over 70 years old, were male, had lower insurance premiums, had comorbidities including chronic kidney disease, interstitial pulmonary disease and HF, and were on oxygen therapy, they had a significantly higher risk of mortality than nonusers. Regarding COPD severity, patients with moderate exacerbation had a significantly higher risk than those with no moderate exacerbation. Furthermore, there was a trend in severe exacerbation in that patients with more severe exacerbations had a higher risk of PH (0, 1, and ≥2 times, aHR = 1.00, 1.25, and 2.55, respectively, p < 0.001). Table 3 showed multivariate Cox proportional hazard model analysis for variables related to all-cause mortality.

TABLE 3

TABLE 3. Multivariate Cox proportional hazard model analysis for variables related to all-cause mortality.

Stepwise multiple regression analysis was then conducted to identify the main factors affecting the mortality rate of patient with PH. Only age group, sex, insurance premium, chronic kidney disease, interstitial pulmonary disease, HF, oxygen therapy and COPD severity were significant risk factors. Those significant factors were used to construct the adjusted model for the subsequent subgroup analysis.

Subgroup Analysis

Different Kinds of Statins

The most commonly used statin was atorvastatin, with 579 (44%) statin users taking this statin. The least commonly used statin was pitavastatin. The Cox proportional hazard model analysis of all-cause mortality in different kinds of statins is summarized in Table 4. The reference for each statin was the matched patients in the statin non-user group.

TABLE 4

TABLE 4. Subgroup analysis of the risk of mortality-related PH for different statins stratified according to statin use.

All of the statins showed a trend toward low risk of mortality (crude HR < 1.00). Patients using lovastatin (aHR = 0.36, 95% CI = 0.14–0.97, p = 0.044) and rosuvastatin (aHR = 0.51, 95% CI = 0.31–0.82, p = 0.005) had a significantly lower risk of mortality. Among all statins, lovastatin had the lowest mortality rate (aHR = 0.36, 95% CI = 0.14–0.97, p = 0.044).

Cumulative Defined Daily Doses

The cDDD was calculated as total statin exposure standardized by DDD over the whole follow-up period. There were six levels of classification for cDDD, and a multivariate Cox proportional hazard model analysis was performed. The statin non-user group was a reference in the analysis. Forty percent of statin users had a cDDD of more than 730. Table 5 shows that the mortality risk of the statin user group was significantly lower than that of the nonuser group until the cDDD was over 365 (aHR = 0.60, 95% CI = 0.40–0.90, p = 0.013). Additionally, patients with a cDDD above 365 had dose-dependent response. The higher the cDDD was, the lower the risk of mortality would be (HR from 1.21 to 0.45, p trend=<0.001).

TABLE 5

TABLE 5. Subgroup analysis of the risk of PH-related mortality, stratified according to the classification of cDDD, duration of statin use, and intensity.

Duration of Statin Use (years)

The duration of statin use was calculated by year and divided into six levels. As shown in Table 5, multivariate Cox proportional hazard model analysis revealed that patients who had a longer duration of statin use had a lower risk of mortality. When patients used statins for longer than 1 year, they started to have a significantly lower risk of mortality (aHR = 0.56, 95% CI = 0.35–0.89, p = 0.015). This result suggested that statin use was effective against PH when the duration of use was longer than one year.

Intensity of Statin Use (Cumulative Defined Daily Doses/Month)

The intensity of statins was calculated by dividing cDDD by the duration of statin use over the entire five-year observation period. Then, the patients were divided into four levels, and multivariate Cox proportional hazard model analysis was performed. The normal DDD per month of statins is 30. Table 5 shows that patients who had a higher statin DDD per month had a lower risk of mortality (aHR = 1.46 to 0.66, p trend=<0.001). Patients whose statin DDD was above 10 per month started to have a significantly lower risk of mortality (aHR = 0.64, 95% CI = 0.49–0.84, p = 0.001). This result suggested that statins had a protective effect against PH until the patients the DDD exceeded 10 per month. Table 5. Subgroup analysis of the risk of PH-related mortality, stratified according to the classification of cDDD, duration of statin use, and intensity.

Sensitivity Analysis

Different Observation Durations

As presented in Table 6, except for that of the one-year observation duration, the aHR of each observation period was similar to the aHR originally observed for five years. After the multivariate Cox proportional hazard model analysis was performed, the 3-year, 5-year, 7-year, 9-year, and end-of-study timepoints still had significantly lower aHRs than the statin non-users. Table 6 showed Sensitivity analysis of risk of mortality during different follow-up times.

TABLE 6

TABLE 6. Sensitivity analysis of risk of mortality during different follow-up times, stratified according to statin use.

Discussion

In our study, statins were effective in improving PH outcomes in COPD. Regarding the primary outcome, the five-year cumulative mortality was 10.94% for statin users and 15.84% for statin non-users. The PH patients with the statin user group had a lower mortality rate than those in than the non-user group (3.87 vs. 5.55 per 100 person-years, p < 0.001) After adjustments for risk factors in the Cox proportional hazard model analysis, the results maintained a significantly lower risk of PH mortality. These results suggest that statins may decrease the risk of mortality from PH by 22% and show a benefit for PH patients.

Reverse Causality

To resolve the reverse causality, we set a six-month drug washout period according to the pharmacokinetic and pharmacodynamic properties of statins (Liu et al., 2016a; Liu et al., 2016b). The half-lives of statins are between 1.2–19 h. It is generally believed that four to five half-lives (i.e., 1–5 days for statins) are needed to excrete most of the drugs from the body (Schachter, 2005). A study involving a mouse model showed that in smooth muscle cells and vascular endothelial cells, discontinuation of atorvastatin (half-life: 14 h) for 2 days caused a 90% decrease in nitric oxide (NO) production. (Laufs et al., 2000). Another study of 21 patients showed reduced endothelial NO synthase activity and NO levels after withdrawal of fluvastatin (half-life: 1.2 h) for 24 h (Westphal et al., 2008; Pineda and Cubeddu, 2011).

According to the above, we think that the six-month washout period was reasonable and long enough to ensure that statins did not influence the COPD diagnosis. Therefore, patients with any statin exposure during the washout period were excluded.

Chronic Obstructive Pulmonary Disease Severity

Because of the lack of forced expiratory volume in one second (FEV1) and exam results in the NHI database, it was difficult to confirm COPD severity in patients. We used the assessment of exacerbation risk in the Global Initiative for Chronic Obstructive Lung disease guidelines to define the severity of COPD. The assessment tool is divided into high-risk and low-risk groups according to the number of exacerbations over the previous year. Severe exacerbations are exacerbations leading to emergency room or hospital admission, and moderate exacerbations are exacerbations were treated with SABA plus antibiotics or oral corticosteroids. If patients had ≥2 moderate exacerbations or ≥1 severe exacerbation over the previous year, they were considered have a high exacerbation risk. If patients had no exacerbations or one moderate exacerbation, they were considered to have a low exacerbation risk (Sun and Ku, 2008). The severity of COPD was confirmed one year after COPD diagnosis.

Pathogenesis of Pulmonary Hypertension in Chronic Obstructive Pulmonary Disease

Because of the pathogenesis of PH involves proliferation of lung vasculature, vasoconstriction, and thrombosis, the anti-proliferative effects and antithrombotic effects of statins as a RhoA/Rho-kinase inhibitions seem to be benefic on PH (Seeger et al., 2000; Sun and Ku, 2008).

Although statins improve systemic inflammatory, pulmonary vascular proliferation, and block the RhoA/Rho-kinase signaling pathway in vitro studies, the efficacy of statins in clinical trials is still unclear. Each statin seems to have a different contribution to protective effect and effectiveness. Compare to pharmacokinetic and pharmacodynamics property of each statin, they were no consistent in bioavailability, protein binding, lipophilicity, half-life, and other metabolism properties. Besides, we could not find studies about comparing the effects of anti-proliferative and anti-inflammatory or the potency of RhoA/Rho-kinase inhibition. Many randomized controlled trials support the benefit of statins in patients with PH related to COPD. In 2009, Lee et al. performed the first randomized controlled trial to explore the effectiveness of statins in this population. After using pravastatin 40 mg one per day for six months, the statin users showed significantly improved exercise tolerance and decreased PH and dyspnea during exercise compared to the placebo group (p < 0.001). Recent studies of simvastatin, rosuvastatin, and atorvastatin showed that they consistently and significantly improved 6 min walking distance results and directly reduced pulmonary artery pressure (Shahmiri et al., 2009; Chogtu et al., 2016; Arian et al., 2017). According to the comprehensive prognostic evaluation and risk assessment in PH, exercise tolerance and pulmonary artery pressures can predict one-year mortality. Therefore, we can presume that if statins can improve these assessment items, statins may be effective for reducing the mortality of PH.

Regarding the secondary outcome, the rates of all-cause mortality (12.34 vs. 17.09 per 100 person-years, p < 0.001) and deaths caused by PH or HF (0.53 vs. 1.14 per 100 person-years, p = 0.005) were a significantly lower in the statin user group than in the non-user group. Although most patients died due to RE, the findings was only a trend of lower risk, with no significant difference. The distribution of causes of death was similar to the most common symptoms that caused patient hospitalization in the Chen et al. study. The majority of patients were hospitalized for dyspnea (93.5%), cough (89.8%), and sputum (79.7%) rather than symptoms caused by PH or HF, such as edema (8.2%) (Chen et al., 2016).

Statin and Prognosis

In previous randomized controlled trials, although statins showed some benefits for PH, the assessments of lung function and heart function, including cardiac output, right ventricular size, FEV1, and lung capacity, did not show significant differences compared to the placebo groups. Only the findings of Chogtu et al.’s study, which showed a significant reduction in COPD exacerbations, are consistent with those of this study (Chogtu et al., 2016). Death can have different causes, so it is difficult to distinguish which statin has the greatest efficacy in reducing the incidence of a specific cause of death. The secondary outcomes only suggest that statins may have an association with a lower risk of RE-related mortality.

In our study, the most commonly used statin was atorvastatin, with nearly 40% of users taking this statin, and the least used statin was pitavastatin, because pitavastatin was the most recent statin to be included in the National health insurance program. Most statins had a trend toward reducing the incidence of PH incidence and the risk of mortality; however, patients in part I who used only pravastatin and atorvastatin had a significantly lower incidence, and patients in part II using lovastatin and rosuvastatin had a significantly lower risk of mortality. The results revealed that the statins seem to have a differing levels of protection and efficacy. When the pharmacokinetic and pharmacodynamic properties of each statin are compared, it becomes apparent that the statins are not the same in terms of bioavailability, protein binding, lipophilicity, half-life, or other metabolic properties. In addition, we could not find studies comparing the antiproliferative and anti-inflammatory effects or the potency of RhoA/Rho-kinase inhibition by these statins.

The probable reason for the inconsistent statin effect between incidence and mortality is that some patients switched to another statin during the observation duration. Therefore, it was difficult to distinguish which statins contributed to the effect. Another reason is the definition of statin grouping. The study divided patients into each statin group according to the statin with the highest cDDD during the observation duration. Perhaps the first statin that the patient was exposed to after the index date or the most recent statin treatment relative to the events of the patient have more influence on the effect. Until now, randomized controlled trials have been performed for only pravastatin, simvastatin, atorvastatin, and rosuvastatin to explore their effects on PH. More direct or indirect comparative studies of statins are needed to support these results.

Dose-dependent Effects

In the subgroup analysis, the results of two parts showed a dose-dependent effect and time-dependent effect. When patients had a higher cDDD and longer duration of statin use, they had a lower risk of PH and mortality related to PH.

To understand the association between cDDD and the duration of statin use, a subgroup analysis of intensity was performed. The intensity indicates the cDDD of statins per month. The normal cDDD per month is 30. However, the DDD is an average maintenance dose per day, so it may represent a moderate intensity of statin use according to their lipid-lowering effects (Wierzbicki et al., 2003; Grundy et al., 2018). Patients with a statin DDD above 20 per month started to have a significantly lower risk of PH (aHR = 0.54, 95% CI = 0.41–0.71), and those with a DDD above 10 per month started to have a significantly lower risk of mortality related to PH (aHR = 0.64, 95% CI = 0.49–0.84, p = 0.001).

The dose-dependent and time-dependent effects of statins are common to many studies in different populations, making the observed effect more robust (Liu et al., 2016a; Liu et al., 2016b; Langballe et al., 2018). In summary, the results suggest that the statins have a protective effect against PH and are effective against PH when patients use them for more than 1 year and at a DDD of more than 10 to 20 per month.

Strengths and Limitations

To our knowledge, this study is the first to explore the protective effect of statins against PH in COPD patients. Moreover, it is the first population-based cohort study to explore the protective effect and efficacy of statins in COPD patients. There are two strengths in our study. First, the NHI database is one of the largest and most comprehensive medical population databases in the world. The data can be extracted from 1995 to 2017 and cover 99.6% of the 23 million people in the Taiwan population. The study had a large sample size, even after propensity-score matching was conducted. There were 41,163 COPD patients and 1,325 patients with PH related to COPD in the final analysis of the study. In contrast to that of other studies, the sample size of patients with PH related to COPD is only 40 to 60 patients (Lee et al., 2009; Wilkins et al., 2010; Zeng et al., 2012; Moosavi et al., 2013; Chogtu et al., 2016; Arian et al., 2017). In general, the larger sample size may provide enough power for statistical analysis. Another strength is the observation period. The study used the database from 2002 to 2017 and set a five-year observation duration. The long observation period is very important in the study. Our study could evaluate the long-term mortality of PH patients, in contrast to other RCT studies, which had follow-up periods of only six months or less. To confirm that the observation time was appropriate, different observation times were tested in the sensitivity analysis. Moreover, the study used event-driven endpoints including the morbidity and mortality of PH.

In addition, RCTs include only more stable PH patients and relatively healthier patients. Chogtu et al. included patients with COPD who were stable for at least 2 weeks and excluded patients with asthma and cardiac disorders such as arrhythmias and unstable angina (Arian et al., 2017). Our study assessed whether COPD severity influences the protective effect and efficacy of statins. Furthermore, the results were robust across several different definitions of statin drug exposure and observation durations. However, we could not find a consistent benefit in terms of protective effect and efficacy between different kinds of statins. Further studies with different statins and accurate statin exposure control may be required.

This study has several limitations. The disease diagnosis and severity assessments of the morbid conditions may not accurately reflect the real situation in this study that could have influenced the outcome of the measurements. The absence of information relating to other factors such as cigarette smoking cannot be underestimated. The study population and outcomes were all defined based on the International Classification of Diseases, Ninth Revision, Clinical Modification or International Classification of Diseases, Ten Revision, Clinical Modification in the admissions record and not on clinical diagnosis because the National health insurance database lacks examination results and lab data. Although we used multiple definitions, such as combined treatment or examination, to improve accuracy, uncertain error is still expected. The limitation of the National health insurance database also leads to the study not having access to smoking status, lifestyle, lung function, or PH severity. Therefore, adjustments could not be made for these important factors or potential confounders, which may influence the outcomes.

Conclusion

Our study suggests that statins have protective effects and efficacy against PH related to COPD and can reduce the mortality of PH. Moreover, the benefits of statins are dose dependent and time dependent.

Age over 60 years, male, low income, HF, and COPD severity are risk factors for mortality, and chronic kidney disease and interstitial pulmonary disease are risk factors for mortality among PH patients.

Data Availability Statement

Data are available from the National Health Insurance Research Database (NHIRD) published by the Taiwan National Health Insurance (NHI) Bureau. Due to legal restrictions imposed by the government of Taiwan in relation to the “Personal Information Protection Act,” data cannot be made publicly available. Requests to access the datasets should be directed to, jk2975525@hotmail.com.

Ethics Statement

This study was approved by the IRB of Kaohsiung Medical University Chung-Ho Memorial Hospital (KMUHIRB-EXEMPT(I)-20190032).

Author Contributions

C-YC and W-TW designed the study, collected the data, interpreted the analyses. Y-LW and K-ML wrote the draft of the manuscript. C-YC and W-TW performed the statistical analyses. C-YC, W-TW, Y-LW, and K-ML reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Kaohsiung Medical University (KMU-Q106011) and Kaohsiung Medical University Hospital (KMUH95-5N16). We thank the Center for Medical Informatics and Statistics of Kaohsiung Medical University for providing administrative and funding support.

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.

Abbreviations

ACEI, angiotensin-converting enzyme inhibitor; ARB, angiotensin receptor blocker; ASCVD, atherosclerotic cardiovascular disease; CCB, calcium channel blocker; cDDD, cumulative defined daily doses; COPD, chronic obstructive pulmonary disease; HF, heart failure; HR, hazard ratio; LABA, long-acting β2-aginist; LAMA, long-acting muscarinic antagonists; PH, pulmonary hypertension; PY, person-year; SABA, short-acting β2-aginist; SAMA, short-acting muscarinic antagonists; TIA, transient ischemic attack; TWD, Taiwan Dollars.

References

Arian, A., Moghadam, S. G., Kazemi, T., and Hajihosseini, M. (2017). The effects of statins on pulmonary artery pressure in patients with chronic obstructive pulmonary disease: a randomized controlled trial. J. Res. Pharm. Pract. 6, 27–30. doi:10.4103/2279-042X.200985

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, W. T., Weng, S. F., Hsu, C. H., Shih, J. Y., Wang, J. J., Wu, C. Y., et al. (2016). Prognostic factors in patients with pulmonary hypertension-a nationwide cohort study. J. Am. Heart Assoc. 5, e003579. doi:10.1161/JAHA.116.003579

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Liu, C., Lu, W., Li, M., Hadadi, C., Wang, E. W., et al. (2016). Clinical characteristics and risk factors of pulmonary hypertension associated with chronic respiratory diseases: a retrospective study. J. Thorac. Dis. 8, 350–358. doi:10.21037/jtd.2016.02.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Chogtu, B., Kuriachan, S., Magazine, R., Shetty, K. R., Kamath, A., George, M. M., et al. (2016). A prospective, randomized study: evaluation of the effect of rosuvastatin in patients with chronic obstructive pulmonary disease and pulmonary hypertension. Indian J. Pharmacol. 48, 503–508. doi:10.4103/0253-7613.190721

PubMed Abstract | CrossRef Full Text | Google Scholar

Cottin, V., Nunes, H., Brillet, P. Y., Delaval, P., Devouassoux, G., Tillie-Leblond, I., et al. (2005). Combined pulmonary fibrosis and emphysema: a distinct underrecognised entity. Eur. Respir. J. 26, 586–593. doi:10.1183/09031936.05.00021005

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundy, S. M., Stone, N. J., Bailey, A. L., Beam, C., Birtcher, K. K., Blumenthal, R. S., et al. (2018). 2018 AHA/ACC/AACVPR/AAPA/ABC/ACPM/ADA/AGS/APhA/ASPC/NLA/PCNA guideline on the management of blood cholesterol. Circulation 139, e1082–e1143. doi:10.1161/CIR.0000000000000625

PubMed Abstract | CrossRef Full Text | Google Scholar

Hurdman, J., Condliffe, R., Elliot, C. A., Swift, A., Rajaram, S., Davies, C., et al. (2013). Pulmonary hypertension in COPD: results from the ASPIRE registry. Eur. Respir. J. 41, 1292–1301. doi:10.1183/09031936.00079512

PubMed Abstract | CrossRef Full Text | Google Scholar

Langballe, R., Cronin-Fenton, D., Dehlendorff, C., Jensen, M. B., Ejlertsen, B., Andersson, M., et al. (2018). Statin use and risk of contralateral breast cancer: a nationwide cohort study. Br. J. Canc. 119, 1297–1305. doi:10.1038/s41416-018-0252-1

CrossRef Full Text | Google Scholar

Laufs, U., Endres, M., Custodis, F., Gertz, K., Nickenig, G., Liao, J. K., et al. (2000). Suppression of endothelial nitric oxide production after withdrawal of statin treatment is mediated by negative feedback regulation of rho GTPase gene transcription. Circulation 102, 3104–3110. doi:10.1161/01.cir.102.25.3104

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, T. M., Chen, C. C., Shen, H. N., and Chang, N. C. (2009). Effects of pravastatin on functional capacity in patients with chronic obstructive pulmonary disease and pulmonary hypertension. Clin. Sci. 116, 497–505. doi:10.1042/CS20080241

PubMed Abstract | CrossRef Full Text | Google Scholar

Levesque, L. E., Hanley, J. A., Kezouh, A., and Suissa, S. (2010). Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ 340, b5087. doi:10.1136/bmj.b5087

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. C., Yang, T. Y., Hsu, Y. P., Hao, W. R., Kao, P. F., Sung, L. C., et al. (2016a). Statins dose-dependently exert a chemopreventive effect against lung cancer in COPD patients: a population-based cohort study. Oncotarget 7, 59618–59629. doi:10.18632/oncotarget.11162

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. C., Hao, W. R., Hsu, Y. P., Sung, L. C., Kao, P. F., Lin, C. F., et al. (2016b). Statins dose-dependently exert a significant chemopreventive effect on colon cancer in patients with chronic obstructive pulmonary disease: a population-based cohort study. Oncotarget 7, 65270–65283. doi:10.18632/oncotarget.11263

PubMed Abstract | CrossRef Full Text | Google Scholar

Moosavi, S. A. J., Raji, H., Faghankhani, M., Yazdani, R., and Esmaeili, M. (2013). Evaluation of the effects of atorvastatin on the treatment of secondary pulmonary hypertension due to chronic obstructive pulmonary diseases: a randomized controlled trial. Iran. Red Crescent Med. J. 15, 649–654. doi:10.5812/ircmj.8267

PubMed Abstract | CrossRef Full Text | Google Scholar

Pineda, A., and Cubeddu, L. X. (2011). Statin rebound or withdrawal syndrome: does it exist? Curr. Atherosclerosis Rep. 13, 23–30. doi:10.1007/s11883-010-0148-x

CrossRef Full Text | Google Scholar

Schachter, M. (2005). Chemical, pharmacokinetic and pharmacodynamic properties of statins: an update. Fundam. Clin. Pharmacol. 19, 117–125. doi:10.1111/j.1472-8206.2004.00299.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Seeger, H., Mueck, A. O., and Lippert, T. H. (2000). Fluvastatin increases prostacyclin and decreases endothelin production by human umbilical vein endothelial cells. Int. J. Clin. Pharm. Ther. 38, 270–272. doi:10.5414/cpp38270

CrossRef Full Text | Google Scholar

Shahmiri, S. S., Adeli, S. H., Mousavi, S. A., Mostafapour, E., Fereshtehnejad, S. M., Namegh, B., et al. (2009). Evaluation of the simvastatin effects in the treatment of pulmonary hypertension (PH). Lijec. Vjesn. 131 (Suppl. 6), 12–15.

CrossRef Full Text | Google Scholar

Simonneau, G., Montani, D., Celermajer, D. S., Denton, C. P., Gatzoulis, M. A., Krowka, M., et al. (2019). Haemodynamic definitions and updated clinical classification of pulmonary hypertension. Eur. Respir. J. 53, 1801913. doi:10.1183/13993003.01913-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Solidoro, P., Patrucco, F., Bonato, R., Boffini, M., Libertucci, D., Ricci, D., et al. (2015). Pulmonary hypertension in chronic obstructive pulmonary disease and pulmonary fibrosis: prevalence and hemodynamic differences in lung transplant recipients at transplant center’s referral time. Transplant. Proc. 47, 2161–2165. doi:10.1016/j.transproceed.2015.01.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., and Ku, D. D. (2008). Rosuvastatin provides pleiotropic protection against pulmonary hypertension, right ventricular hypertrophy, and coronary endothelial dysfunction in rats. Am. J. Physiol. Heart Circ. Physiol. 294, H801–H809. doi:10.1152/ajpheart.01112.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M. T., Liou, J. T., Lin, C. W., Tsai, C. L., Wang, Y. H., Hsu, Y. J., et al. (2018). Association of cardiovascular risk with inhaled long-acting bronchodilators in patients with chronic obstructive pulmonary disease: a nested case-control study. JAMA Intern. Med. 178, 229–238. doi:10.1001/jamainternmed.2017.7720

PubMed Abstract | CrossRef Full Text | Google Scholar

Westphal, S., Abletshauser, C., and Luley, C. (2008). Fluvastatin treatment and withdrawal: effects on endothelial function. Angiology 59, 613–618. doi:10.1177/0003319708316005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wierzbicki, A. S., Poston, R., and Ferro, A. (2003). The lipid and non-lipid effects of statins. Pharmacol. Ther. 99, 95–112. doi:10.1016/s0163-7258(03)00055-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, M. R., Ali, O., Bradlow, W., Wharton, J., Taegtmeyer, A., Rhodes, C. J., et al. (2010). Simvastatin as a treatment for pulmonary hypertension trial. Am. J. Respir. Crit. Care Med. 181, 1106–1113. doi:10.1164/rccm.200911-1699OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W.-T., and Chen, C.-Y. (2020). Protective effect of statins on pulmonary hypertension in chronic obstructive pulmonary disease patients: a nationwide retrospective, matched cohort study. Sci. Rep. 10, 3104. doi:10.1038/s41598-020-59828-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, W. J., Xiong, C. M., Zhao, L., Shan, G. L., Liu, Z. H., Xue, F., et al. (2012). Atorvastatin in pulmonary arterial hypertension (APATH) study. Eur. Respir. J. 40, 67–74. doi:10.1183/09031936.00149011

PubMed Abstract | CrossRef Full Text | Google Scholar