Abstract
Objective
The heterogeneous spectrum of systemic lupus erythematosus (SLE) often presents with secondary complications such as cardiovascular disease (CVD), infections and neoplasms. Our study assessed whether the presence of SLE independently increases or reduces the disparities, accounting for the already higher risk of these outcomes among racial/ethnic minority groups without SLE.
Methods
We defined a cohort using electronic health records data (2005–2016) from a mixed-payer community-based outpatient setting in California serving patients of diverse racial/ethnic backgrounds. The eligible population included adult patients with SLE and matched non-SLE patients (≥18 years old). SLE was the primary exposure. The following outcomes were identified: pneumonia, other infections, CVD and neoplasms. For each racial/ethnic group, we calculated the proportion of incident co-morbidities by SLE exposure, followed by logistic regression for each outcome with SLE as the exposure. We evaluated interaction on the additive and multiplicative scales by calculating the relative excess risk due to interaction and estimating the cross-product term in each model.
Results
We identified 1036 SLE cases and 8875 controls. The incidence for all outcomes was higher among the SLE exposed. We found little difference in the odds of the outcomes associated with SLE across racial/ethnic groups, even after multivariable adjustment. This finding was consistent on the multiplicative and additive scales.
Conclusion
We demonstrated that SLE status does not independently confer substantial interaction or heterogeneity by race/ethnicity toward the risk of pneumonia, other infections, CVD or neoplasms. Further studies in larger datasets are necessary to validate this novel finding.
Introduction
Systemic lupus erythematosus (SLE) is an autoimmune disease whose clinical management is often complex and multifaceted. In the rheumatology literature, there have been some calls for SLE researchers and clinicians to shift their efforts from addressing the control of the primary disease to becoming more mindful of managing the constellation of secondary complications or outcomes in SLE such as cardiovascular disease (CVD),1–3 cerebrovascular manifestation 4 and infections. 5 In the United States, racial minorities (i.e. blacks, Hispanics, Native Americans and Asians) experience higher rates of these outcomes in the general population compared with whites6—health disparities that are also particularly pathognomonic of SLE. In addition to a much higher incidence and prevalence of SLE among female patients,7,8 consistent across multiple studies is the finding that racial and ethnic minorities experience certain clinical phenotypes.9–11 For example, a recent prevalence estimate indicates that women and non-whites comprise approximately 90% and 70% of prevalent SLE in the United States, respectively. 12
In general, racial and ethnic minorities with SLE also have more severe disease activity and worse outcomes.13–15 Within the past decade, important insights have emerged regarding the differential distribution of key co-morbidities or outcomes by race/ethnicity among patients with SLE. Medicaid data demonstrated that blacks and Native Americans with SLE had 13–37% increased risk of serious infections compared with whites with SLE. 16 Interestingly, a Hispanic and Asian paradox has been described within the Medicaid population, wherein these patients with SLE were found to have decreased mortality 17 and risk of myocardial infarction 18 compared with their white counterparts. Less well explored is the contribution of race/ethnicity to the risk of malignancy in SLE. While one study identified non-Hispanic whites as having the highest risk of neoplasms, 19 a subsequent larger study of patients in California found higher rates among Hispanic patients with SLE. 20
Prior studies have mainly focused on comparing secondary complications or outcomes by race/ethnicity among patients with SLE. However, we know that some of these same outcomes occur with increased frequency among minority populations without SLE compared with their white counterparts. A key question with relevance to the clinician emerges: does the presence of SLE independently increase or reduce the disparities in the frequency of these co-morbid conditions among racial/ethnic minorities compared with that already seen within the non-SLE population? The findings of this inquiry could improve the understanding of the heterogeneous spectrum of SLE by enabling better phenotyping of patients at risk of key outcomes or complications. In addition, understanding whether SLE contributes substantially to the existing white–racial minority gap in outcomes can guide efforts aimed at minimizing racial disparities in this realm. With this in mind, we assessed the independent contribution of SLE on numerous outcomes separately across different racial/ethnic groups in a large, well-defined racially diverse cohort of adults in Northern California. We explore the interaction between race/ethnicity and SLE on the risk of five outcomes and specifically assess whether the magnitude of health disparities in these outcomes is greater or lesser than would be expected for patients with SLE compared with controls. This study was funded as a supplement to the Cardiovascular Disease among Asians and Pacific Islanders (CASPER) study (funded by grant HL 126172 from NHLBI to Drs. Waitzfelder, Palaniappan, and Fortmann). 21
Methods
Setting and study population
Our study population was derived from a mixed-payer community-based outpatient care setting which serves approximately one million patients each year in Northern California. We defined a cohort of patients who were diagnosed with SLE and non-SLE comparators. The eligible population included adult patients with SLE and demographically matched non-SLE patients receiving care (i.e. those with ≥2 recorded visits) in the same setting between 2005 and 2016. Patients without SLE (referred to as controls throughout) were frequency-matched 10:1 to SLE subjects by calendar year of diagnosis (for cases) or first service date (for controls), age and sex. The earliest date where an SLE diagnosis is found in the electronic health records is the index date for each SLE case (aged ≥18 years on index date and did not die before 2005). The earliest date on or after 1 January 2005 when a subject meets all eligibility criteria (i.e. aged ≥18, in care) is the index date for controls. We used calendar year of index date (index year), sex and birth year as the predictors to calculate a propensity score and matched based on the value of the propensity score.22,23 In the first round of 1:10 matching, eight matched controls died before seven cases’ index date. Another two rounds of matching were conducted to find eight replacements for eight controls. We initially identified 1290 SLE cases who were matched by sex, age and index year to 12,900 controls. After excluding patients without race/ethnicity profiles, there were 1036 SLE cases and 8875 controls. The study was approved by the Sutter Health Institutional Review Board under the project ID number 976421-2.
Systemic lupus erythematosus
SLE was the primary exposure and was defined as having ≥2 clinic visits from 1 January 2005 onwards with an ICD9 code of 710.0 or an ICD10 code of M32.1, M32.8, M32.9. 16 We were underpowered to restrict to incident SLE, therefore have included prevalent and incident SLE.
Outcomes
We identified five separate co-morbidities/outcomes: pneumonia, other infections, CVD and neoplasms (all and malignancies, separately) using ICD codes for inpatient and outpatient visits (Supplementary Table 1). We included the records of patients with outcomes that were diagnosed after SLE diagnosis (for SLE cases) and after the index date (for controls). For each analysis, individuals with a history of the specific outcome in the medical record according to ICD codes and the problems list at baseline (SLE diagnosis or index date) were identified and then subsequently were excluded from the specific analysis to remove potential prevalent cases and focus on incident outcomes. Only the first event after SLE diagnosis or the index date was counted. Since there are multiple outcomes, an individual patient may have been counted under more than one outcome (e.g. CVD and infection). For each of the outcomes we noted the first time that they have that outcome—thus, the underlying population for each outcome was slightly different. For example, for CVD incidence, we counted instances among those who have never had a history of CVD (i.e. those at risk) in the codes or problem list before the start of follow-up, but these individuals could have had infections or other outcomes. Participants were followed until death, loss to follow-up (primarily leaving the health care system), or 31 December 2016.
Additional covariates
Additional covariates extracted from the electronic health records included: self-reported biological sex (male or female), race/ethnicity and baseline values of age, body mass index (BMI) in kg/m2 (categorized as: underweight ≤18.5, 18.5<normal weight ≤25, 25<overweight≤30, or obese>30), smoking status (categorized as: current, former, or never), and diagnoses of diabetes, hypertension and dyslipidemia. Race/ethnicity was categorized as non-Hispanic white, non-Hispanic black, Hispanic, Asian (including Asian Indian, Chinese, Filipino, Japanese, Korean, Vietnamese and other Asian). Potential confounders were different for each exposure-outcome model.
Statistical analysis
We compared demographic and clinical characteristics between SLE cases and controls using frequencies and proportions for categorical variables, and medians with interquartile ranges for continuous variables. These characteristics were compared using Pearson Chi-square test for categorical variables and Student’s t-test for continuous variables. We found that overall approximately 22% of BMI and smoking data were missing. This was relatively similar between SLE cases and controls for smoking but different for BMI (11% SLE and 23% non-SLE). We used multiple imputation to estimate missing values of smoking and BMI. 24 We calculated the proportion of newly diagnosed (incident) pneumonia, infections, CVD and neoplasms among SLE cases and controls by race/ethnicity. Separate logistic regression models estimated odds ratios (OR) and 95% confidence intervals (CI) for the association between SLE and outcome adjusting for potential confounding by sex, age, smoking status, BMI and history of diabetes by racial/ethnic group.
Evaluation of additive and multiplicative interactions
We evaluated interaction on additive and multiplicative scales using whites as the reference group because the study sample was not powered to examine multiple comparisons among race/ethnicity groups. We extracted three datasets comprising Asians and whites only (dataset 1), blacks and whites only (dataset 2), and Hispanics and whites only (dataset 3). First, we built five logistic regression models with SLE status (case versus control) as the exposure of interest for each outcome (pneumonia, other infections, CVD, neoplasms and malignant neoplasms) for each dataset, that is, we built a model for every outcome for each of the datasets. Thus, there were 15 models in total (see Table 4). In each model, we included potential confounders, as well as an interaction term for race/ethnicity and SLE status. For example, the regression equation modeling malignant neoplasms as the outcome comparing Asians with whites (using dataset 1) would be:
In another example, the regression equation modeling CVD as the outcome comparing blacks with whites (using dataset 2) would be:
Second, to detect additive interaction, we calculated the Relative Excess Risk due to Interaction (RERI), a measure that assesses whether the risk difference for each outcome between whites and each race/ethnic minority was different from zero.25,26 We used the following equation: RERI ≈ OR11 – OR10 – OR01 + 1 and estimated the corresponding 95% CIs. For the example of malignant neoplasms in Asians versus whites, the RERI was calculated as:
Third, multiplicative interaction was calculated using the cross-product term in each model. We estimated whether the expected OR was greater (for antagonistic interaction) or lesser (for synergistic interaction) than the observed OR for each outcome between whites and each race/ethnic minority group.25,26
Results
Baseline characteristics of study participants by SLE status in the PAMF population, 2005–2016
p < 0.05.
#Body mass index was categorized (kg/m2) as: underweight ≤18.5, 18.5<normal weight≤25, 25<overweight ≤30, or obese >30.
^Malignant neoplasms are a subset of all neoplasms
The incidence of major outcomes comparing SLE with non-SLE controls, 2005–2016*
NH Black: Non-Hispanic Black; NH White: Non-Hispanic White.
The denominators are not the same across outcomes because the analyses was restricted to those without a history of these outcomes
^Malignant neoplasms are a subset of all neoplasms.
Odds ratios of having outcomes comparing SLE cases to controls stratified by race/ethnicity
NH Black: Non-Hispanic Black; NH White: Non-Hispanic White.
non-SLE groups are reference categories.
Covariates: Gender, SLE status, Baseline age, Smoking status.
Covariates: Gender, SLE status, Baseline age, Smoking status.
Covariates: Gender, SLE status, Baseline age, Smoking status, BMI category, Diabetes.
Covariates: Gender, SLE status, Baseline age, Smoking status, BMI category.
Covariates: Gender, SLE status, Baseline age, Smoking status, BMI category.
Interaction between SLE status and race/ethnicity on the risk of outcomes
NH Black: Non-Hispanic Black; NH White: Non-Hispanic White.
RERI appear to be statistically significant.
To calculate interaction on the additive and multiplicative scales, we extracted three datasets comprising Asians and whites only (dataset 1), blacks and whites only (dataset 2), and Hispanics and whites only (dataset 3).
1
Dataset 1;
2
Dataset 2;
3
Dataset 3. In each model, we included potential confounders, as well as an interaction term for race/ethnicity and SLE status. For example, the regression equation modeling malignant neoplasms as the outcome comparing Asians with whites (using dataset 1) would be:
Discussion
Overall, we found that patients with SLE had more CVD, pneumonia, infections and neoplasms than the control population across all racial/ethnic groups in our study. There were also some differences in incidence of these outcomes across racial/ethnic groups within patients with SLE. However, when addressing our central clinical question of the independent association between SLE and these outcomes by racial/ethnic group, we found there was almost no evidence for substantial interaction or heterogeneity by race/ethnicity. This finding was generally consistent across both multiplicative and additive scales.
The reasoning behind the lack of an independent disease-mediated effect on the risk of outcomes by race/ethnicity in SLE is unknown but is an intriguing finding that warrants further study. We could not identify similar studies, limiting our ability to contextualize our findings in the literature. There are two ways to answer the questions that arise in these analyses. One way would be to explore whether the health disparities found in the outcomes in racial/ethnic minorities with SLE may be due to the disease uniformly activating pathways that already predispose racial minorities to outcomes. Another way would be to explore whether the heterogeneity in these outcomes has nothing to do with SLE pathogenesis but rather may be a consequence of underlying pathways present in the background risk within the general population. The findings of these inquiries would help anticipate the degree to which the white–racial minority gap may widen or narrow in the future because of shifting trends in SLE. Such knowledge may be useful in planning and targeting systems to deliver care to a growing population of patients with SLE.
There is an abundance of studies in the SLE literature comparing race/ethnic groups, with whites as the reference categories. Oftentimes, the regression models in these studies treat race/ethnicity as a confounder of the associations between the exposure(s) and outcome(s) of interest. For example, although the Medicaid studies found significant differences among the patients with SLE by race, they were not necessarily comparing patients with SLE with the general population or examining the effect of SLE on those outcomes for each racial group separately.8,16,18 The questions that arise in these settings are variations of the following: assuming that we have a population of patients with SLE, which racial/ethnic groups do worse? Our current work instead dissects this from a different angle and asks: given that we know that some of these outcomes are more common in certain racial/ethnic groups, are the racial differences still apparent when one factors out the increased incidence in the non-SLE minority populations? Thus, we suggest that future studies of racial disparities in SLE outcomes need to include general population comparators so that inferences can be made about the extent that SLE widens the gaps in existing inequalities.
We found higher prevalence of more “upstream” risk factors such as diabetes, hypertension and dyslipidemia among patients with SLE compared with controls, and these conditions are often precursors to CVD and neoplasms. When we adjusted for these risk factors in the regression models, the significant associations between SLE and the outcomes did not change. It is not clear how these traditional risk factors and those specific to SLE (e.g. autoantibodies, interferon signatures) impact CVD outcomes in racial minorities with SLE. In the general population, CVD risk algorithms such as the Framingham Risk Score and the American College of Cardiology/American Heart Association have been shown to validate well in racial minorities;27,28 however, these algorithms may underestimate CVD risk in patients with SLE. 29 The extent to which risk algorithms (and their constituents) can accurately predict CVD risk in individual race/ethnic groups with SLE is currently unknown and warrants research resources. The extent to which diabetes influences CVD and cancer among patients with SLE and how this association is patterned along race/ethnic groups is not well understood in rheumatology. In addition, one question we ask clinicians to ponder is whether the severe clinical phenotypes observed in racial minorities with SLE have correlates within the general population, and if the severe manifestations of these outcomes are due to the underlying predisposition of race/ethnic minorities for these conditions, and not SLE. There is increasing evidence that these risk factors are patterned differently for race/ethnic groups; for example, a 10 mm Hg higher systolic blood pressure level is associated with 24% increase in stroke risk for blacks and only an 8% increase in stroke risk for whites. 30
The increased incidence of key medical outcomes among patients with SLE in our study—across all racial/ethnic groups—reinforces the health challenges faced by these patients beyond their autoimmune condition. It also presents an opportunity for healthcare providers who care for these patients to enhance current efforts at disease prevention within this patient population. This is certainly true for rheumatologists, who are tasked with the chronic management of these patients’ underlying SLE. At the same time, primary care physicians are often the first point of contact prior to referral to a rheumatologist, and they remain essential to patients’ long-term health maintenance. Close partnership between the rheumatologist and primary care provider is thus critical, ensuring that appropriate screening and treatment of these outcomes occur in a timely manner that conforms to accepted guidelines and recommendations. Whether an updated screening and management approach toward these outcomes among patients with SLE is needed, in light of the increased risk within this patient population, is worthy of further discussion.
To our knowledge, this is one of a few studies to investigate heterogeneity by race/ethnicity in the association between SLE status and outcomes in a diverse US sample. Our study has a few limitations. We were unable to assess severity of these outcomes. Additional research is needed to understand whether there is interaction between race/ethnicity and SLE status on mortality and co-morbidity progression over time, and to understand the underlying biological mechanisms. Our findings may have limited generalizability due to the study sample being composed of individuals who are insured, likely skewing the population by socioeconomic status. In contrast, several recent epidemiologic studies assessing the impact of race/ethnicity on cardiovascular and infectious outcomes in SLE were drawn from the Medicaid population. Therefore, it is difficult to compare our findings with respect to these outcomes of interest in particular. Due to the study location in California, our sample may not represent the general population of patients with SLE in the US, as blacks comprise 43% of the national prevalence estimates, compared with 6% in our study. 12 However, our study was enriched with Hispanic and Asian patients (32% and 16%, respectively) while these groups comprise 16% and 13% of the national prevalent estimates. 12 We also cannot exclude the possibility that differences in the risk of outcomes by race/ethnicity exist at a population level and had small within-group sample sizes for some race/ethnic groups. We were underpowered to examine differences among more granular groupings of race/ethnicity. Future analyses will include Pacific Islanders and parse out Asian subgroups. We were also unable to account for person time, disease duration, clinical phenotypes or treatment consideration, due to the limited ability to ascertain when clinical manifestations first showed in the electronic health records. Thus, there may be some misclassification of prevalent outcomes being labeled as incident and we cannot predict whether this is differential or nondifferential.
The possibility of confounding by race/ethnicity is eliminated in the study design because we are interested in the heterogeneity and therefore are looking by groups. Stratified analyses such as these also implicitly account for confounding because they reduce the variability of that factor within strata. 31 We acknowledge that there may be some measurement error around the reporting of race/ethnicity. Studies of heterogeneity by race/ethnicity are often misleading, as race categories are social constructs that often do not reflect the genomic diversity in population subgroups. Thus, we cannot rule out that there is misclassification/measurement error of the race/ethnicity effect modifier and that this may impact the validity of our findings. The race/ethnicity categories reflected in the electronic health records are based on self-report, although it is possible that on rare occasions some are simply entered by clinic intake staff based on appearance. However, it is our understanding that this practice is relatively uncommon.
In summary, we demonstrated that while SLE increases the risk of secondary complications, it does not independently confer substantial heterogeneity by race/ethnicity toward the risk of pneumonia, other infections, CVD or neoplasms. Further studies in larger datasets are necessary to validate this novel finding. At the same time, our work confirms the higher risk of these important outcomes among patients with SLE overall—similar to prior studies—and also within each race/ethnicity category compared with controls. Future research should also include exploration of mental health (e.g. depression, schizophrenia, etc.), kidney disease and pregnancy/birth outcomes, as these conditions are associated with SLE and racial/ethnic minorities are at increased risk for these diseases. Our study brings further awareness toward the increased burden of medical outcomes among patients with SLE, and is a call to arms to rheumatologists and primary care providers alike to enhance preventative health efforts to reduce the incidence and severity of these outcomes within this at-risk patient population.
Supplemental Material
LUP884646 Supplemetal Material - Supplemental material for Does SLE widen or narrow race/ethnic disparities in the risk of five co-morbid conditions? Evidence from a community-based outpatient care system
Supplemental material, LUP884646 Supplemetal Material for Does SLE widen or narrow race/ethnic disparities in the risk of five co-morbid conditions? Evidence from a community-based outpatient care system by T Falasinnu, Y Chaichian, J Li, S Chung, B E Waitzfelder, S P Fortmann, L Palaniappan and J F Simard in Lupus
Footnotes
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded as a supplement to the Cardiovascular Disease among Asians and Pacific Islanders (CASPER) study (funded by grant HL 126172 from NHLBI to Drs. Waitzfelder, Palaniappan, and Fortmann). Dr. Simard is supported by a K01 award from NIAMS.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
