Association between bisphosphonate use and COVID-19 related outcomes

  1. Jeffrey Thompson
  2. Yidi Wang
  3. Tobias Dreischulte
  4. Olga Barreiro
  5. Rodrigo J Gonzalez
  6. Pavel Hanč
  7. Colette Matysiak
  8. Harold R Neely
  9. Marietta Rottenkolber
  10. Thomas Haskell
  11. Stefan Endres
  12. Ulrich H von Andrian  Is a corresponding author
  1. Cerner Enviza, United States
  2. Dept. of Immunology, Harvard Medical School, United States
  3. Institute of General Practice and Family Medicine, University Hospital of Ludwig Maximilians-University Munich, Germany
  4. Center of Integrated Protein Science Munich and Division of Clinical Pharmacology, University Hospital, LMU Munich, Germany, Germany

Abstract

Background:

Although there are several efficacious vaccines against COVID-19, vaccination rates in many regions around the world remain insufficient to prevent continued high disease burden and emergence of viral variants. Repurposing of existing therapeutics that prevent or mitigate severe COVID-19 could help to address these challenges. The objective of this study was to determine whether prior use of bisphosphonates is associated with reduced incidence and/or severity of COVID-19.

Methods:

A retrospective cohort study utilizing payer-complete health insurance claims data from 8,239,790 patients with continuous medical and prescription insurance January 1, 2019 to June 30, 2020 was performed. The primary exposure of interest was use of any bisphosphonate from January 1, 2019 to February 29, 2020. Bisphosphonate users were identified as patients having at least one bisphosphonate claim during this period, who were then 1:1 propensity score-matched to bisphosphonate non-users by age, gender, insurance type, primary-care-provider visit in 2019, and comorbidity burden. Main outcomes of interest included: (a) any testing for SARS-CoV-2 infection; (b) COVID-19 diagnosis; and (c) hospitalization with a COVID-19 diagnosis between March 1, 2020 and June 30, 2020. Multiple sensitivity analyses were also performed to assess core study outcomes amongst more restrictive matches between BP users/non-users, as well as assessing the relationship between BP-use and other respiratory infections (pneumonia, acute bronchitis) both during the same study period as well as before the COVID outbreak.

Results:

A total of 7,906,603 patients for whom continuous medical and prescription insurance information was available were selected. A total of 450,366 bisphosphonate users were identified and 1:1 propensity score-matched to bisphosphonate non-users. Bisphosphonate users had lower odds ratios (OR) of testing for SARS-CoV-2 infection (OR = 0.22; 95%CI:0.21–0.23; p<0.001), COVID-19 diagnosis (OR = 0.23; 95%CI:0.22–0.24; p<0.001), and COVID-19-related hospitalization (OR = 0.26; 95%CI:0.24–0.29; p<0.001). Sensitivity analyses yielded results consistent with the primary analysis. Bisphosphonate-use was also associated with decreased odds of acute bronchitis (OR = 0.23; 95%CI:0.22–0.23; p<0.001) or pneumonia (OR = 0.32; 95%CI:0.31–0.34; p<0.001) in 2019, suggesting that bisphosphonates may protect against respiratory infections by a variety of pathogens, including but not limited to SARS-CoV-2.

Conclusions:

Prior bisphosphonate-use was associated with dramatically reduced odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations. Prospective clinical trials will be required to establish a causal role for bisphosphonate-use in COVID-19-related outcomes.

Funding:

This study was supported by NIH grants, AR068383 and AI155865, a grant from MassCPR (to UHvA) and a CRI Irvington postdoctoral fellowship, CRI2453 (to PH).

Editor's evaluation

Using health insurance claims data, this valuable paper reports on a retrospective propensity score matched cohort study that was performed to quantify associations between bisphosphonate (BP) use and COVID-19-related outcomes (COVID-19 diagnosis, testing, and COVID-19 hospitalization). The evidence is solid showing that in primary and sensitivity analyses, BP use was consistently associated with lower odds for COVID-19, testing, and COVID-19 hospitalization. The study is of interest to a broad readership (clinicians, public health physicians, pharmacologists and epidemiologists).

https://doi.org/10.7554/eLife.79548.sa0

eLife digest

The COVID-19 pandemic challenged the world to rapidly develop strategies to combat the virus responsible for the disease. While several effective vaccines and new drugs have since become available, these therapies are not always easy to access and take time to generate and distribute. To address these challenges, researchers have tried to find ways to repurpose existing medications that are already commonly used and known to be safe.

One potential candidate are bisphosphonates, a family of drugs used to reduce bone loss in patients with osteoporosis. Bisphosphonates have been shown to boost the immune response to viral infections, and it has been observed that patients prescribed these drugs are less likely to develop or die from pneumonia. But whether bisphosphonates are effective against COVID-19 had not been fully explored.

To investigate, Thompson, Wang et al. analyzed insurance claims data from about 8 million patients between January 2019 and June 2020, including around 450,000 individuals that had filled a prescription for bisphosphonates. Patients prescribed bisphosphonates were then compared to non-users that were similar in terms of their gender, age, the type of health insurance they had, their access to healthcare, and other health comorbidities.

The study revealed that bisphosphonate users were around three to five times less likely to be tested for, diagnosed with, or hospitalized for COVID-19 during the first four months of the pandemic. They were also less commonly diagnosed with other respiratory infections in 2019, like bronchitis or pneumonia.

Although the results suggest that bisphosphonates provide some protection against COVID-19, they cannot directly prove it. Verifying that bisphosphonates can treat or prevent COVID-19 and/or other respiratory infections requires more studies that follow patients in real-time rather than studying previously collected data.

If such studies confirm the link, bisphosphonates could be a helpful tool to protect against COVID-19 or other virus outbreaks. The drugs are widely available, safe, and affordable, and therefore may provide an alternative for patients who cannot access other medications or vaccines.

Introduction

Throughout the COVID-19 pandemic, massive global efforts to repurpose existing drugs as potential therapeutic options for COVID-19 have been undertaken. Drug repurposing, whereby a drug already proven to be safe and effective in humans for another approved clinical indication is evaluated for novel clinical use, may allow for faster identification and deployment of therapeutic agents compared to traditional drug discovery pipelines. Using in silico and in vitro analyses, a growing list of drugs have been suggested to be potentially efficacious in treating COVID-19 by either direct or indirect antiviral actions (Sultana et al., 2020). Another potentially beneficial class of drugs may be agents that boost or modulate anti-viral immune responses to SARS-CoV-2 infection to reduce clinical symptoms and/or mitigate disease progression. Regardless of the mechanism of action, ultimately, randomized prospective clinical studies are needed to test the safety and efficacy of each candidate in treating or preventing COVID-19. Observational studies can help prioritize candidates for prospective clinical testing, by examining associations between the use of a candidate drug and the incidence or severity of disease in users compared to a matched group of non-users. Drugs with strong observational evidence for potential effectiveness against COVID-19 may then be considered for prospective trials (Sultana et al., 2020).

Here, we have investigated bisphosphonates (BPs), a class of small-molecule drugs that inhibit bone resorption by osteoclasts (Roelofs et al., 2010b). BPs are widely prescribed as either oral or intravenous formulations to treat osteoporosis, Paget disease, and malignancy-induced hypercalcemia. Additionally, BPs are used as adjuvant therapy for breast cancer (Dhesy-Thind et al., 2017). BPs are subdivided into two classes, nitrogen-containing (amino-BPs) and nitrogen-free BPs (non-amino-BPs; Russell et al., 2008). Both accumulate in bone but have distinct molecular mechanisms by which they kill osteoclasts to prevent bone resorption (Roelofs et al., 2010b).

Aside from depleting osteoclasts, clinical and experimental studies indicate that BPs exert a plethora of immunomodulatory effects, providing a rationale for exploring BPs as potential repurposed drug candidates for COVID-19 (Brufsky et al., 2020). Indeed, amino-BPs regulate the activation, expansion, and/or function of a major subset of human γδT cells (Poccia et al., 2006; Hewitt et al., 2005; Tu et al., 2011) as well as neutrophils (Favot et al., 2013), monocytes (Roelofs et al., 2010a), and macrophages (Rogers and Holen, 2011; Wolf et al., 2006); they can modulate the antigen-presentation capacity of dendritic cells (Xia et al., 2018); and in animal studies, both amino-BPs and non-amino-BPs exerted potent adjuvant-like activity to boost antibody and T cells responses to viral antigens (Tonti et al., 2013). Furthermore, observational studies have reported decreased in-hospital mortality for patients in the ICU (Lee et al., 2016), and reduced incidence of pneumoniae and pneumonia-related mortality in patients treated with amino-BPs versus controls (Sing et al., 2020). These immunological and clinical effects of BPs combine with several other characteristics that make BPs well-suited as repurposed drug candidates in the context of a pandemic: they are globally accessible as generics, affordable, straightforward to administer, and have known safety profiles in adult (Suresh et al., 2014) and paediatric populations (Sbrocchi et al., 2010; George et al., 2015).

In light of these considerations, we have analysed a database of health insurance claims in the U.S. to determine if prior BP-use is associated with a differential incidence and/or severity of COVID-19-related outcomes. Specifically, we assessed the relationship between use of BPs and COVID-19-related hospitalizations and COVID-19 diagnosis, as well as testing for SARS-CoV-2 infection (as a proxy for severe COVID-19 symptoms given the restricted access to testing during the initial surge). Outcomes were measured from March 1, 2020 to June 30, 2020, a period that roughly coincided with the first wave of COVID-19 in the U.S. and predated the advent of potential outcome modifiers, such as vaccines or other effective treatment options.

Methods

Study design

A retrospective cohort study was performed using health insurance claims data from January 1, 2019 to June 30, 2020 (study period) in order to assess the relationship between use of BPs and three COVID-19-related outcomes: (a) testing for SARS-CoV-2 infection; (b) COVID-19 diagnosis; and (c) hospitalization with a COVID-19 diagnosis, whereby COVID-19-related hospitalization was deemed the primary endpoint and COVID-19 diagnosis and testing were secondary endpoints. Primary and secondary endpoints were assessed during the observation period of March 1, 2020 to June 30, 2020, roughly corresponding to the first nation-wide surge of COVID-19 in the U.S. (Figure 1A). In the primary analysis, the risk of COVID-19-related outcomes was assessed among BP users compared to a matched sample of BP non-users with similar demographic and clinical characteristics.

Study periods, cohort selection, and analyses of BP use on COVID-19-related outcomes.

(A) Schematic overview of the study timeline. (B) Schematic flow diagram illustrating the identification of the study population and matched control populations for primary analysis and sensitivity analyses cohorts. BP: bisphosphonate; CA: California; CCI: Charlson comorbidity index; CI: confidence interval; COPD: chronic obstructive pulmonary disease; FL: Florida; IL: Illinois; NY: New York; OR: odds ratio; PCP: primary care physician; PS: propensity score; PSM: propensity score match.

Data source

Data used for this study included closed medical (inpatient and outpatient) and outpatient-pharmacy-dispensed claims between January 1, 2019 and June 30, 2020, from the Komodo Health payer-complete dataset (https://www.komodohealth.com). This dataset is derived from over 150 private insurers in the U.S. and includes patients with commercial, individual, state exchange-purchased, Medicare Advantage, and Medicaid managed-care insurance coverage. The dataset also provides information on insurance eligibility periods. Closed claims within this dataset represent those that had undergone insurance adjudication. In total, the Komodo Health payer-complete dataset includes health insurance claims data from over 140 million individuals in the U.S. from 2015 to 2020.

Cohort definition

All patients were required to have continuous medical and prescription insurance eligibility during the entire study period. Patients with missing information for age, gender, insurance type, or state/region were excluded.

Exposures of interest

The primary exposure of interest was the use of any amino- or non-amino BP medication. Exposure to BPs and all other medications of interest were assessed over a 14-month pre-observation period preceding the COVID-19 pandemic in the U.S. This long duration was chosen because of the extended bioavailability of BPs, which accumulate in bone where they are retained and slowly released for up to several years (Cremers et al., 2019). Patients were classified as BP users if they had any claim at any time during the pre-observation period for one of the following: alendronate, alendronic acid, etidronate, ibandronate, ibandronic acid, pamidronate, risedronate, and zoledronic acid (full details in Appendix 1).

Timing of BP dose

The effect of timing and formulation of BPs on COVID-19-related outcomes was more closely examined by varying the window between BP exposure and outcome measurement. The primary analysis BP user cohort, along with their propensity-score matched (see below for cohort matching) BP non-user cohort, were stratified as follows: two cohorts were used as the reference comparator with known BP-exposure during all or most of the pre-observation and the entire observation period, specifically (i) BP users who took oral alendronic acid (dosed daily or weekly) throughout the pre-observation period (i.e. at least one claim or drug-on-hand in each quarter in 2019 and in Jan/Feb. 2020) that also had a days-supply extending past June 30, 2020, and (ii) users of infusion zoledronic acid (dosed annually) with a claim in Q3 or Q4 2019; two cohorts with BP-exposure only during the pre-observation period, namely (iii) users of alendronic acid occurring during the first six months of 2019 with days-covered ending prior to June 30, 2019 and no other BP claims thereafter, and (iv) users of zoledronic acid in January or February 2019 with no other BP claims during the remainder of the study period; and, two cohorts with short-term BP exposure, specifically new users of (v) alendronic acid or (vi) zoledronic acid in February 2020, with no prior BP claims during the pre-observation period.

Covariates

As covariates, we considered factors that may influence either the use of BPs or potential modulators of primary or secondary study endpoints. These included: age; gender; insurance type (commercial, dual, Medicaid, Medicare); having had any primary care physician (PCP) visit in 2019; and comorbidity burden. The variable ‘PCP visit in 2019’ was used to control for prior healthcare-use behaviour and was assigned based on any physician office claim from January 1, 2019 to December 31, 2019 with one of the following provider types: family practice, general practice, geriatric medicine, internal medicine, and preventive medicine. Comorbidity score assignment was calculated following the Charlson Comorbidity Index (CCI) methodology (Quan et al., 2005), and was based on diagnosis codes present on any medical claim (inpatient or outpatient) occurring during the pre-observation period. The assigned CCI score was used as the comorbidity covariate for the primary cohort propensity score matching, but to better control for differences in comorbidity burden when assessing outcomes, all regression analyses involving the primary analysis cohort included the following individual comorbidity covariates in lieu of the aggregate CCI score: osteoporosis, cancer, chronic obstructive pulmonary disease (COPD), depression, dyslipidaemia, hypertension, obesity, type 2 diabetes, cardiovascular disease overall, sickle cell anemia, stroke, dementia, HIV/AIDS, chronic kidney disease/end-stage renal disease (CKD/ESRD), and liver disease (Appendix 1).

Cohort matching

For the primary analysis, BP users were propensity-score (PS) matched to BP non-users via a PS calculated using multiple variables, including age, gender, insurance type, CCI, and any PCP visit in 2019, to yield comparable populations by demographics and clinical characteristics (Figure 1B). To account for the differential geographic spread of COVID-19 across the U.S. during the observation period, matching was performed within each geographic region separately (Northeast, Midwest, South, West) and then combined. In addition to this within-region stratified match, a cohort build was also performed after restricting to patients from New York (NY) state only, since this state was the site of the largest outbreak in the initial COVID-19 surge in the U.S. All matching algorithms used a greedy-match propensity score technique (Parsons, 2001) to match BP users to non-users with a maximum permitted propensity-score difference of 0.015.

Definition of endpoints

Primary and secondary endpoints were assigned using inpatient and outpatient medical claims that occurred during the four-month observation period. The primary endpoint, COVID-19-related hospitalization, was assigned based on the presence of an International Classification of Diseases, Tenth Revision (ICD-10) code on any inpatient medical service claim indicating test-confirmed 2019 Novel Coronavirus (2019-nCoV) acute respiratory disease, specifically U07.1. The first secondary endpoint, SARS-CoV-2 testing, was assigned using Current Procedural Terminology (CPT) codes indicating a test for active infection, specifically 87635, 87636, and 87637. The second secondary endpoint, COVID-19-related diagnosis, was assigned based on any medical service claim with the ICD-10 diagnosis code U07.1.

Statistical analysis

Unadjusted analyses assessing the association between BP-use and COVID-19-related outcomes were performed for the primary analysis cohort using chi-square tests for categorical variables and calculation of the crude unadjusted odds ratio (OR) in the matched cohort groups overall, when stratified by region and in NY state alone, and when further stratified by age group and gender. Chi-square tests for categorical variables and t-tests for continuous variables were also performed to assess differences in demographic and clinical characteristics of BP users compared to BP non-users both pre-match and post-match to assess the success of the propensity-score match.

Multivariate logistic regression analyses, modelled separately to determine the adjusted OR for each COVID-19-related primary and secondary outcome while adjusting for demographic and clinical characteristics, were performed on the matched primary analysis cohort with all regions combined, when stratified by region, and in NY state alone. The primary exposure of interest was BP-use (yes/no) during the pre-observation period. Additional demographic/clinical characteristics also included as regression model covariates were: age group, gender, region (for all regions-combined analyses), insurance type, PCP visit in 2019, and the following comorbid conditions: osteoporosis, cancer, COPD, depression, dyslipidaemia, hypertension, obesity, type 2 diabetes, cardiovascular disease overall, sickle cell anaemia, stroke, dementia, HIV/AIDS, CKD/ESRD, and liver disease. Demographic characteristics used in the matching procedure were also included in the final outcome regressions to control for the impact of those characteristics on outcomes modelled.

All tests were two-tailed, and p-values of less than 0.05 were considered significant. All analyses were performed using SAS 9.4 (Cary, NC).

Sensitivity analyses

Multiple sensitivity analyses were performed to assess the reliability of the primary analysis results and/or to address potential unmeasured confounding (full details in Appendix 1).

  1. The first sensitivity analysis addressed potential confounding by indication (i.e. the possibility of the indication for BP use rather than BP use itself being responsible for differences in outcomes among BP users and non-users) by restricting the control group to an active comparator cohort of patients who had used non-BP anti-resorptive bone medications during the pre-observation period. Users of non-BP anti-resorptive bone medications, the smaller patient population, were then 1:1 matched to BP users, providing a sample where all patients had used bone health medications during the pre-observation period (‘Bone-Rx’ cohort) (Figure 1B). Cohort matching and regression modelling were performed following the same methodology employed for the primary analysis.

  2. The second sensitivity analysis further addressed potential baseline differences between users of BPs and users of non-BP anti-resorptive bone medications in terms of indication for treatment and risk of SARS-CoV-2 exposure. To homogenise indication for treatment, we restricted the ‘Bone-Rx’ cohort to females aged older than 50 years with an osteoporosis diagnosis (ICD-10: M80.x, M81.x, M82.x), which is the main (but not the only) indication for use of anti-resorptive bone medications. In order to homogenise risk of COVID-19 exposure, we additionally (a) restricted both groups to residents of New York, Illinois, Florida, and California (four states with a high incidence of COVID-19 cases during the observation period, with each representing a geographic region) (CDC, 2021a), and (b) matched within each state by insurance-type strata (i.e. BP non-users matched to BP users with Medicaid coverage residing in New York) to control for differences in socioeconomic characteristics. Non-BP anti-resorptive bone medication users were then matched to BP users by age, PCP visit in 2019, and the following select comorbid conditions that include those thought to impact COVID-19 severity: cancer, COPD, depression, dyslipidaemia, heart failure, hypertension, obesity, and type 2 diabetes (Rosenthal et al., 2020). In addition to assessing COVID-19-related outcomes, the matched cohorts that resulted from this analysis, older female patients from New York, Illinois, Florida, or California with a diagnosis of osteoporosis who were users of BP or non-BP anti-resorptive medications (‘Osteo-Dx-Rx’ cohort), were used for the third and fourth sensitivity analyses (see below).

  3. The third sensitivity analysis assessed the relationship between BP-use and exploratory positive control outcomes (anticipated to be impacted by the immunomodulatory pharmacological mechanism of BPs) occurring in 2019. For this analysis, the primary, ‘Bone-Rx’, and Osteo-Dx-Rx” cohorts were restricted to BP users who had any BP claim during the first half of 2019 and their previously-assigned BP non-user matched pair to assess the relationship between BP-use and medical services for other respiratory infectious diseases (acute bronchitis, pneumonia).

  4. The fourth sensitivity analysis addressed potential bias due to the 'healthy adherer' effect, whereby users of a preventive drug may have better disease outcomes due to their healthier behaviours rather than due to drug treatment itself (Ladova et al., 2014). Two strategies were employed to validate the findings from our primary analysis while controlling for the potential impact of healthy adherer effect-associated bias. First, we tested whether effects observed with exposure to BPs were similarly observed with exposure to other preventive drugs, namely statins, antihypertensives, antidiabetics, and antidepressants. Second, we assessed whether the association between BP-use and COVID-19-related outcomes was maintained among the matched user/non-user populations of these other preventive drugs, i.e. BP users were compared to BP non-users within, for example, the statin user population and separately within the matched statin non-user population.

Results

Study population

A total of 8,239,790 patients met the inclusion criterion of continuous medical and prescription insurance eligibility over the full study period, of which 333,107 were excluded due to missing demographic information, resulting in a total eligible sample of 7,906,603 patients (Figure 1B). Of this full population, 452,051 (5.7%) and 7,454,552 (94.3%) patients were classified as BP users and BP non-users, respectively. Within BP users, more than 99% were prescribed an amino-BP, with oral alendronic acid (75.4%), zoledronic acid infusion (11.5%), and oral ibandronic acid (8.4%) as the most prevalent formulations (Table 1).

Table 1
Most recent bisphosphonate claim among all users.
Drug (route)N%
Alendronate / alendronic acid (oral)340,81075.4%
Etidronate (oral)140.0%
Ibandronate / ibandronic acid (oral)37,9888.4%
Ibandronic acid (injection/infusion)11690.3%
Pamidronate (injection/infusion)11210.2%
Risedronate (oral)18,9914.2%
Zoledronic acid (injection/infusion)51,95811.5%

Prior to propensity-score matching, there were significant differences between BP users and non-users across all demographic and clinical characteristics. BP users were older (age >60: 82.7% vs 27.7%; p<0.001), predominantly female (91.0% vs 57.2%; p<0.001), with a higher comorbidity burden (mean CCI 0.95 vs 0.60; p<0.001), with a larger proportion of patients residing in the Western U.S. (21.1% vs 15.4%; p<0.001), covered by Medicare (43.3% vs 13.7%; p<0.001), and having visited a PCP in 2019 (63.8% versus 44.7%; p<0.001). Propensity-score matching yielded 450,366 BP users and 450,366 BP non-users with no significant differences across all characteristics used in matching (Table 2). Differences did exist, however, in the distribution of individual comorbid condition indicators that were used as covariates in the regression analysis, with the BP non-user cohort having a higher proportion of patients with COPD (10.2% vs 8.5%; p<0.001), cardiovascular disease (25.1% vs 18.7%; p<0.001), dyslipidemia (36.9% vs 34.6%; p<0.001), hypertension (46.4% vs 38.8%; p<0.001), obesity (10.3% vs 6.7%; p<0.001), and type 2 diabetes (22.9% vs 18.2%; p<0.001). Over 98% of all BP user/non-user matches for the primary analysis cohort were completed with differences in matched propensity scores <0.000001 (overall mean difference of 0.000004, max difference of 0.0147).

Table 2
Primary analysis cohort (all regions), patient characteristics pre/post match.
All Observations UnmatchedAll Observations Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients7,906,603100.00%7,454,55294.30%452,0515.70%900,732100.00%450,36650.00%450,36650.00%
Demographics
Age
 ≤201,840,05023.30%1,838,92224.70%1,1280.20%<0.0012,2530.30%1,1250.20%1,1280.30%1
 21-401,446,99918.30%1,443,90819.40%3,0910.70%6,1950.70%3,1040.70%3,0910.70%
 41-50925,30911.70%916,75812.30%8,5511.90%17,0961.90%8,5451.90%8,5511.90%
 51-601,250,19015.80%1,184,46915.90%65,72114.50%131,44514.60%65,72414.60%65,72114.60%
 61-701,181,26114.90%1,024,38313.70%156,87834.70%313,82234.80%156,94434.80%156,87834.80%
 71-80783,7759.90%642,0508.60%141,72531.40%280,80331.20%140,36631.20%140,43731.20%
 ≥81479,0196.10%404,0625.40%74,95716.60%149,11816.60%74,55816.60%74,56016.60%
Gender
 Female4,670,96059.10%4,263,52457.20%407,43690.10%<0.001811,49790.10%405,74690.10%405,75190.10%0.99
 Male3,235,64340.90%3,191,02842.80%44,6159.90%89,2359.90%44,6209.90%44,6159.90%
Region
 Midwest1,467,80218.60%1,391,83518.70%75,96716.80%<0.001151,80216.90%75,90116.90%75,90116.90%1
 Northeast2,152,56027.20%2,032,83227.30%119,72826.50%238,98826.50%119,49426.50%119,49426.50%
 South3,042,60438.50%2,881,71838.70%160,88635.60%319,40835.50%159,70435.50%159,70435.50%
 West1,243,63715.70%1,148,16715.40%95,47021.10%190,53421.20%95,26721.20%95,26721.20%
Insurance
 Commercial3,938,60349.80%3,791,54550.90%147,05832.50%<0.001294,07032.60%147,01232.60%147,05832.70%1
 Dual156,4972.00%125,0901.70%31,4076.90%59,9366.70%29,9806.70%29,9566.70%
 Medicaid2,594,50032.80%2,517,02033.80%77,48017.10%154,51917.20%77,27217.20%77,24717.20%
 Medicare1,217,00315.40%1,020,89713.70%196,10643.40%392,20743.50%196,10243.50%196,10543.50%
PCP Visit 2019
 No4,283,69754.20%4,119,83155.30%163,86636.20%<0.001327,38336.30%163,65936.30%163,72436.40%0.89
 Yes3,622,90645.80%3,334,72144.70%288,18563.80%573,34963.70%286,70763.70%286,64263.60%
Clinical Characteristics
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.621.380.61.350.951.76<0.0010.951.760.951.760.951.760.7
Regression Comorbidity Covariates
N%N%N%p-valueN%N%N%p-value
Osteoporosis267,0203.40%135,2311.80%131,78929.20%<0.001163,81418.20%32,3907.20%131,42429.20%<0.001
Cancer419,0835.30%366,7864.90%52,29711.60%<0.00194,14810.50%41,8619.30%52,28711.60%<0.001
CKD/ESRD361,4514.60%328,6334.40%32,8187.30%<0.00168,9997.70%36,1828.00%32,8177.30%<0.001
COPD466,0945.90%427,8505.70%38,2448.50%<0.00184,2349.40%45,99010.20%38,2448.50%<0.001
CVD1,084,03113.70%999,52613.40%84,50518.70%<0.001197,24321.90%112,93325.10%84,31018.70%<0.001
Dementia125,8111.60%113,7781.50%12,0332.70%<0.00124,9212.80%12,8892.90%12,0322.70%<0.001
Depression571,3037.20%531,3557.10%39,9488.80%<0.00186,2809.60%46,43110.30%39,8498.80%<0.001
Dyslipidemia1,532,25419.40%1,375,92018.50%156,33434.60%<0.001322,12535.80%166,36036.90%155,76534.60%<0.001
HIV/AIDS33,2290.40%31,7110.40%15180.30%<0.00128970.30%13790.30%1,5180.30%0.01
Hypertension1,899,06324.00%1,723,51923.10%175,54438.80%<0.001384,05942.60%209,18446.40%174,87538.80%<0.001
Liver Disease251,3313.20%231,6643.10%19,6674.40%<0.00138,6974.30%19,0314.20%19,6664.40%0.001
Obesity638,5068.10%608,0838.20%30,4236.70%<0.00176,8448.50%46,49810.30%30,3466.70%<0.001
Sickle Cell Anemia10,4990.10%10,2920.10%2070.00%<0.0014220.00%2150.00%2070.00%0.7
Stroke104,8591.30%97,0011.30%7,8581.70%<0.00119,3952.20%11,5692.60%7,8261.70%<0.001
Type 2 Diabetes978,23912.40%895,98312.00%82,25618.20%<0.001184,97820.50%103,03122.90%81,94718.20%<0.001

Similar profiles in pre-match versus post-match characteristics were seen when patients were stratified by region or restricted to NY-state (Appendix 2—tables 1–3, Appendix 2—table 4, Appendix 2—table 5). Demographic distributions, including differences between BP user versus BP non-user characteristics pre-match versus post-match characteristics were seens pre- and post-matching for all sensitivity analysis cohorts are detailed in Appendix 2.

BP use and COVID-19-related outcomes

Among the full matched cohort, BP users had significantly lower rates and unadjusted (crude) odds of testing (1.2% vs 5.1%; OR = 0.22; 95%CI:0.21–0.22; p<0.001), diagnosis (0.7% vs 2.9%; OR = 0.22; 95%CI:0.21–0.23; p<0.001), and hospitalization (0.2% vs 0.7%; OR = 0.24; 95%CI:0.22–0.26; p<0.001) as compared to BP non-users (Figure 2 and Appendix 3—figure 1). Consistent findings were seen when sub-stratifying the full matched cohort by age, gender, age*gender, within grouped regions, by individual region, and in NY-state alone (Appendix 2—tables 6–11).

Association of BP use and COVID-19-related outcomes incidence (left) and regression-adjusted results for odds (right) of SARS-CoV-2 testing (blue), COVID-19 diagnosis (purple), and COVID-19-related hospitalizations (red) of BP users compared with BP non-users in the all-regions combined primary analysis cohort (i) and when stratified by region/state into: Northeast (ii), Midwest (iii), South (iv), West (v), and New York state (vi).

For details see Figure 2—source data 1.

Figure 2—source data 1

COVID-19-related outcomes in the primary analysis cohort.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig2-data1-v2.xlsx

Multivariate regression analyses yielded similar results for all outcomes while additionally controlling for patient demographic and comorbidity characteristics. In the full matched cohort, BP users had lower adjusted odds of testing (OR = 0.22; 95%CI:0.21–0.23; p<0.001), diagnosis (OR = 0.23; 95%CI:0.22–0.24; p<0.001), and hospitalizations (OR = 0.26; 95%CI:0.24–0.29; p<0.001). These findings were robust when comparing BP users with BP non-users when stratified by geographic region or NY-state alone.

Timing of last BP exposure and COVID-19-related outcomes

The above results demonstrate that any BP exposure during the 14-months pre-observation period is associated with a marked reduction in each of the three COVID-19-related outcomes. To further investigate the relationship between COVID-19-related outcomes and the timing of BP exposure, we focused on the two most commonly prescribed BPs, alendronic acid (oral formulation dosed daily or weekly) and zoledronic acid (infusion dosed annually). For each BP type, COVID-19-related outcomes were assessed among users: (i-ii) with exposure or days covered (based on prescription frequency) during the pre-observation period and throughout the observation period; (iii-iv) with exposure or days covered ending prior to the observation period; and (v-vi) newly initiating therapy prior to the observation period (Figure 3A). Furthermore, all subgroups of BP users had decreased odds of COVID-19-related outcomes (Figure 3B) except for the odds of hospitalization among zoledronic acid users who were last dosed in January/February of 2019 (OR = 0.52; 95%CI:0.20–1.40; p=0.20) or newly initiated in February of 2020 (OR = 0.49; 95%CI:0.13–1.88; p=0.30).

Timing of BP use and COVID-19-related outcomes.

(A) Schematic of BP user sub-stratification by timing of exposure to alendronic acid or zoledronic acid prior to outcome assessment. Broken lines represent periods of active BP dosing. For zoledronic acid users, days covered was considered to extend 1 year past the dosing period based on dosing guidelines. (B) Incidence (left) and regression-adjusted results (right) for odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations of BP users compared with BP non-users in pre-specified subgroups. For further details see Figure 3—source data 1. CI: confidence interval; OR: odds ratio.

Figure 3—source data 1

Primary analysis cohort by timing of BP dosing, COVID-19-related outcomes.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig3-data1-v2.xlsx

Sensitivity analysis 1: COVID-19-related outcomes among all users of anti-resorptive medications (‘Bone-Rx’ cohort)

The first sensitivity analysis was performed to address potential confounding by indication. To validate our primary findings in more comparable cohorts, analysis was restricted to comparing BP users to patients using non-BP anti-resorptive bone medications during the pre-observation period. Compared to non-BP users of anti-resorptive medications, BP users had decreased odds of testing (OR = 0.31; 95%CI:0.28–0.33; p<0.001), diagnosis (OR = 0.35; 95%CI:0.31–0.38; p<0.001), and hospitalization (OR = 0.45; 95%CI:0.36–0.56; p<0.001) (Figure 4A and Appendix 3—figure 2). Furthermore, these findings were robust when assessed separately across every geographic region as well as NY state for all outcomes except hospitalizations when restricted to the Western U.S. (p=0.08; Appendix 2—table 12).

COVID-19-related outcomes among the Bone-RX and Osteo-Dx-Rx restricted cohorts.

Incidence and forest plots summarizing regression-adjusted odds ratios of SARS-CoV-2 testing (blue), COVID-19 diagnosis (purple), and COVID-19-related hospitalizations (red) in the (A) ‘Bone-Rx’ (see also Figure 4—source data 1) and (B) ‘Osteo-Dx-Rx’ sensitivity analysis cohorts (see also Figure 4—source data 2).

Figure 4—source data 1

Source data for Figure 4A: Bone-Rx cohort COVID-19-related outcomes.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig4-data1-v2.xlsx
Figure 4—source data 2

Source data for Figure 4B: Osteo-Dx-Rx cohort COVID-19-related outcomes.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig4-data2-v2.xlsx

Sensitivity analysis 2: COVID-19-related outcomes among users of anti-resorptive medications with a diagnosis of osteoporosis (‘Osteo-Dx-Rx’ cohort)

The second sensitivity analysis was performed to address the fact that, even after restricting the comparator cohort to users of anti-resorptive medications, differences may still exist between patient cohorts that could affect COVID-19-related outcomes, including different indications for anti-resorptive medication use and other uncontrolled patient characteristics. To address this, the association between BP use and COVID-19 related outcomes were examined in a cohort restricted to female patients over 50 years old, with a diagnosis of osteoporosis, using either a BP or a non-BP anti-resorptive bone medication, matched within insurance-type as a proxy for socioeconomic status, and selected from four states (NY, IL, FL, CA) with high incidences of COVID-19 cases during the observation period (CDC, 2021a; ‘Osteo-Dx-Rx’ cohort). In agreement with the results reported above, the decrease in odds of COVID-19-related outcomes in BP users remained robust for testing (OR = 0.28; 95%CI:0.23–0.35; p<0.001), diagnosis (OR = 0.40; 95%CI:0.32–0.49; p<0.001), and hospitalizations (OR = 0.45; 95%CI:0.26–0.75; p=0.003) (Figure 4B).

Sensitivity analysis 3: Association of BP-use with exploratory positive control outcomes

The third sensitivity analysis was performed to assess if there is an association between BP-use and incidence of other respiratory infections, which has been previously reported (Sing et al., 2020). Medical services for acute bronchitis or pneumonia were measured during the second half of 2019, prior to the advent of COVID-19, in the primary, ‘Bone-Rx’, and ‘Osteo-Dx-Rx’ cohorts. Regression modelling found that, among all cohort variations modelled, BP users had a decreased odds of any medical service related to acute bronchitis (point estimates of ORs ranged from 0.23 to 0.28) and pneumonia (point estimates of ORs ranged from 0.32 to 0.36) (Figure 5).

Exploratory outcomes among BP users versus BP non-users.

Incidence and adjusted odds ratios of other respiratory infections, in the primary, ‘Bone-Rx’, and ‘Osteo-Dx-Rx’ cohorts. For details, see Figure 5—source data 1. CI: confidence interval; OR: odds ratio.

Figure 5—source data 1

Positive control outcomes by primary, bone-Rx, and osteo-Dx-Rx cohorts.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig5-data1-v2.xlsx

Sensitivity analysis 4: Association of other preventive drugs with COVID-19-related outcomes

A potential pitfall in the interpretation of apparent effects of preventive medications on health outcomes is the so-called healthy adherer effect, whereby patients may have better outcomes due to their overall healthier behaviours and not due to active drug treatment itself (Ladova et al., 2014). To address this possibility of unmeasured confounding, a final sensitivity analysis was performed to evaluate the association between control exposures (i.e. use of other preventive medications such as statins, antihypertensives, antidiabetics, and antidepressants) and COVID-19-related outcomes (Figure 6A). In comparison to BPs, the impact of other preventive drug classes on COVID-19-related outcomes was much weaker overall (Figure 6B–E) and varied between geographic regions in terms of magnitude or direction (Appendix 2—tables 13–16). Furthermore, when assessing the impact of BP-use within matched user/non-user preventive drug cohorts (e.g. BP users compared to BP non-users among the matched statin user and statin non-user populations), we found BP-use to be consistently associated with lower odds of testing (point estimates of ORs ranged from 0.21 to 0.27), diagnosis (point estimates of ORs ranged from 0.22 to 0.30), and hospitalizations (point estimates of ORs ranged from 0.25 to 0.33) across all stratified preventive user/non-user cohorts (Figure 6B–E).

Association of other preventive drugs with COVID-19-related outcomes.

(A). Schematic illustrating the identification of study populations and matched controls for each drug class. (B–E) Incidence and adjusted odds ratios of SARS-CoV-2 testing (blue), COVID-19 diagnosis (purple), and COVID-19-related hospitalizations (red) in users and non-users of (B) statins (see also Figure 6—source data 1), (C) antihypertensive medications (see also Figure 6—source data 2), (D) non-insulin antidiabetic medications (see also Figure 6—source data 3), and (E) antidepressant medications (see also Figure 6—source data 4). For each class of preventive medications, further analysis was performed comparing BP users and BP non-users within matched cohorts of medication users (middle) and medication non-users (bottom). BP: bisphosphonate; CCI: Charlson comorbidity index; CI: confidence interval; COPD: chronic obstructive pulmonary disease; OR: odds ratio; PCP: primary care physician; PS: propensity score; PSM: propensity score match.

Figure 6—source data 1

Source data for Figure 6B: COVID-19-related outcomes by statin use overall & sub-stratified by BP use.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig6-data1-v2.xlsx
Figure 6—source data 2

Source data for Figure 6C: COVID-19-related outcomes by antihypertensive use overall & sub-stratified by BP use.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig6-data2-v2.xlsx
Figure 6—source data 3

Source data for Figure 6D: COVID-19-related outcomes by antidiabetic use overall & sub-stratified by BP use.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig6-data3-v2.xlsx
Figure 6—source data 4

Source data for Figure 6E: COVID-19-related outcomes by antidepressant use overall & sub-stratified by BP use.

https://cdn.elifesciences.org/articles/79548/elife-79548-fig6-data4-v2.xlsx

Discussion

This study examined the association between recent exposure to BPs and subsequent COVID-19-related outcomes during the initial outbreak of the COVID-19 pandemic in the U.S. Our findings demonstrate that amino-BP users experienced a three- to five-fold reduced incidence of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalization during this period. This dramatic difference in outcomes was consistently observed when comparing BP users to BP non-users in a propensity score-matched general population, when comparing to users of other anti-resorptive bone medications, when further restricting the latter cohort to female osteoporosis patients matched by comorbidities within state of residence and by insurance type, and when comparing BP users to BP non-users stratified by use of other preventive medications. Therefore, although there are confounding-related limitations inherent within retrospective studies, the consistency and strength of our observed associations when using various methods to control for unmeasured confounding support the contention that further prospective research should be performed to determine the true magnitude of the potential immunomodulatory effects of BP use.

Our findings are consistent with previous observational studies, prior to the advent of COVID-19, that had reported associations between BP use and reduced incidence of pneumonia and pneumonia-related mortality (Sing et al., 2020; Colón-Emeric et al., 2010; Reid et al., 2021). Accordingly, we observed in our population that BP use was associated with decreased odds of medical services for acute bronchitis and pneumonia during the second half of 2019. Taken together, these findings suggest that BPs may play a protective role in respiratory tract infections from a variety of causes, including SARS-CoV-2.

Other recent retrospective studies have explored, to some extent, associations of anti-resorptive medication use and COVID-19-related outcomes, albeit in much smaller patient populations than were analysed here. One study found no differences in the COVID-19-related risk of hospitalization (70.7% vs 72.7%, p = 0.16) and mortality (11.9% vs 12.8%, p = 0.386) among 1,997 female patients diagnosed with COVID-19 who received anti-osteoporosis medication as compared to propensity score-matched COVID-19 patients who were not receiving such drugs (Atmaca et al., 2022). This study did not examine the incidence of COVID-19 among BP users, but it raises the possibility that the subset of BP users who do develop sufficient pathology to be diagnosed with COVID-19 may have a similar clinical course as BP non-users. Another retrospective cohort study in Italy examining the association of oral amino-BP use and incidence of COVID-19-related hospitalization found no difference between BP users (12.32% (95% CI, 9.61–15.04)) and BP non-users (11.55% (95% CI, 8.91–14.20)) (Degli Esposti et al., 2021). However, the overall incidence of COVID-19 hospitalization in the primary cohort (151/126,370 patients, or 0.12%) of this study was markedly lower than in the present analysis (3,710/900,732 patients, or 0.41%). A third study examined the influence of various anti-osteoporosis drugs, including BPs, on the cumulative incidence of COVID-19 in 2,102 patients with non-inflammatory rheumatic conditions that were compared to population estimates in the same geographic region (Blanch-Rubió et al., 2020). In this analysis, users of non-BP anti-resorptive medications and zoledronate, but not users of oral BPs, had a lower incidence and relative risk of COVID-19 diagnosis and hospitalization. The observations with zoledronate are consistent with the findings reported here. However, we did not detect a significant impact of non-BP anti-resorptive medications in comparison to BPs, and we found a robust association between oral BP use and lower odds of COVID-19 diagnosis and related hospitalization. The reason for these discrepancies is unclear but could potentially reflect the large disparity in sample size between our study, which differed by more than three orders of magnitude. A fourth study, which used Israeli insurance data to perform an analysis involving two separate case-control matched cohorts to assess the risk of COVID-19 hospitalizations when stratified by recent medication use, also found that the odds COVID-19-related hospitalizations were lower among users of BPs, and ranged from an OR of 0.705 (95%CI: 0.522–0.935) to 0.567 (95%CI:0.400–0.789) (Israel et al., 2021).

The large size of our dataset allowed for a range of fully powered, stratified analyses to be performed to explore the robustness of our findings and to address unmeasured confounding factors and other sources of potential bias that can occur in retrospective studies using insurance claims data. Notwithstanding, a retrospective analysis of insurance claims data has inevitable limitations that should be considered. Specifically, there is the potential that key patient characteristics impacting outcomes could not be derived from claims data. For example, the interpretation of our findings depends, in part, on the assumption that BP users and non-users had a similar risk of SARS-CoV-2 infection during the observation period. However, our dataset does not allow us to restrict patient observations to those with known exposure to SARS-CoV-2. Therefore, to minimize potential differences in SARS-CoV-2 exposure between BP users and non-users in our primary study cohort, we implemented additional analytical strategies, including the sensitivity analyses, as well as matching BP users to BP non-users within geographical regions and specific states.

Despite these efforts, it is important to note that we have limited information to assess and match BP users to BP non-users by sociodemographic risk factors, such as socio-economic status and racial/ethnic minority status, that are associated with COVID-19 incidence and mortality (Karmakar et al., 2021; Rogers et al., 2020). Notably, Black/African-American and Hispanic patients have been shown to have significantly higher test positivity rates (Kaufman et al., 2021; Escobar et al., 2021; Jacobson et al., 2021; Rubin-Miller, 2020) and severity of disease at the time of testing (Rubin-Miller, 2020). Furthermore, Black/African American (Azar et al., 2020) and Hispanic patients were found to have a higher incidence of COVID-19 infection (Escobar et al., 2021; CDC, 2021b) and odds of COVID-19 related hospital admission even after adjustment for comorbidities (Nau et al., 2021), residence in a low-income area (Rubin-Miller, 2020), and insurance plan (Azar et al., 2020; Price-Haywood et al., 2020; Muñoz-Price et al., 2020). The greater COVID-19 burden in these groups is likely due to a combination of systemic health inequities as well as a disproportionate representation among essential workers (Selden and Berdahl, 2020; US Bureau of Labor Statistics, 2019), which could potentially increase their exposure risk to SARS-CoV-2. In addition, there are known variations in the prevalence of osteoporosis between different racial groups, which could potentially result in disproportionate frequencies of BP prescriptions (No authors listed, 2021). The potential confounding due to socio-economic status and differential prevalence of osteoporosis among racial/ethnic groups was addressed in our analysis of the ‘Osteo-Dx-Rx’ cohort where we compared BP users to non-users after restricting to female patients with a diagnosis of osteoporosis, all using anti-resorptive bone medications, and matched by insurance type (proportion of Medicaid and dual Medicare/Medicaid users) as a proxy for social-economic status (Figure 4B). Nevertheless, this strategy cannot rigorously rule out a potential under-representation of groups with higher sociodemographic risk factors among BP users that could have contributed to the observed decreased odds of COVID-19 related outcomes in our primary analyses.

The potential bias introduced by a putative differential racial/ethnic group composition of BP users versus BP non-users is at least partially addressed by a recent study of a large Californian cohort of female BP users (Black et al., 2020). Compared to the racial composition of California at-large (a proxy for BP non-users) (United States Census Bureau, 2019), BP users were predominantly Non-Hispanic White (36.5% in California versus 53.3% among BP users). The proportions of Black/African-Americans and Asians among BP users in that study were similar to those in California at-large, whereas Hispanic patients represented a smaller percentage (24%) of BP users as compared to Hispanics in the state’s general population (39.4%). Based on these findings and the reported differential case rates of COVID-19 infections among racial groups in California (Reitsma et al., 2021), we can estimate the race-adjusted incidence of COVID-19 in populations reflecting the composition of BP users and non-users (Black et al., 2020) to be 1.7% and 2.1%, respectively. By comparison, in our study the actual rate of COVID-19 diagnosis in the Western US was 2.5% for BP non-users versus 0.46% for BP users (Figure 2), indicating that the uneven representation of ethnic/racial groups cannot fully explain the observed differences in COVID-19 related outcomes. Moreover, we note that racial/ethnic minorities are also under-represented among statin users (Salami et al., 2017), but statin-users in our primary cohort had similar odds of COVID-19 hospitalization as statin non-users (Figure 6B). Similarly, Black/African-Americans and Hispanics have lower utilization rates of antidepressants (Chen and Rizzo, 2008) and Hispanics were also reported to be undertreated with antihypertensive medications (Gu et al., 2017). Our analysis of COVID-19-related outcomes among users and non-users of antihypertensives showed a modest decrease in COVID-19 diagnosis and minimal association with COVID-19-related hospitalization (Figure 6C). By contrast, users of antidepressants had uniformly lower odds for both endpoints (Figure 6E), which is consistent with other recent studies (Israel et al., 2021; Hoertel et al., 2021; Zimniak et al., 2021). However, regardless of the class of non-BP preventive drugs analysed, concomitant BP use was consistently associated with dramatically decreased odds of COVID-19 diagnosis and hospitalization as well as testing for SARS-CoV-2 (Figure 6B–E).

Furthermore, specifically looking at the rate of SARS-CoV-2 testing in California (Escobar et al., 2021; Jacobson et al., 2021) or nation-wide (Kaufman et al., 2021), the proportions of different racial and ethnic groups among tested patients were nearly identical to estimates for the state or national population. Thus, the observed association between BP use and reduced testing for SARS-CoV-2 infection in our nation-wide cohorts is unlikely to be explained by potential differences in racial composition between BP users and non-users. It also seems unlikely that exposure to BPs reduces the actual incidence of SARS-CoV-2 infections. More likely, we propose that immune-modulatory effects of BPs may enhance the anti-viral response of BP users to SARS-CoV-2 and mitigate the development of symptoms. Milder or absent symptoms may have caused infected BP users to be less likely to seek testing. Moreover, because there was a nationwide shortage of available tests for SARS-CoV-2 during the observation period, patients needed to present with sufficiently severe disease symptoms to be eligible for testing, so fewer test-seeking BP users may have qualified. Consequently, a larger proportion of uncaptured ’silent' infections among BP users could explain why fewer diagnoses and hospitalizations were observed in this group.

The scarceness of COVID-19 tests combined with the strain on healthcare systems during the observation period could potentially have resulted in a misclassification bias whereby some patients may have been falsely diagnosed and/or hospitalized with COVID-19 without having received a confirmatory test. However, this bias should equally affect BP users and BP non-users and bias our findings towards the null. Relatedly, limited hospital capacity during the observation period could have led to rationing of inpatient hospital beds based on severity of disease and likelihood to survive (Emanuel et al., 2020). However, matching by age and comorbidities should produce patient populations with similar characteristics used for rationing.

A further limitation of our study is the lack of information on the result of COVID-19 tests received by patients. Therefore, as discussed above, the incidence and odds of COVID-19 testing should not be viewed as a proxy for the rate of infection, but rather reflects the incidence of patients with severe enough symptoms or exposure to warrant testing. Another potential source of confounding is the possibility that some patients in our study were classified as BP non-users due to the absence of BP exposure during the pre-observation period but may have received a BP during the observation period. The potential misclassification of BP non-users, however, would bias towards the null hypothesis, and was only seen in 1.92% of the matched BP non-user population.

An additional limitation is potential censoring of patients who died during the observation period, resulting in truncated insurance eligibility and exclusion based on the continuous insurance eligibility requirement. However, modelling the impact of censoring by using death rates observed in BP users and non-users in the first six months of 2020 and attributing all deaths as COVID-19-related did not significantly alter the decreased odds of COVID-19 diagnosis in BP users (see Appendix 3).

Another limitation in the current study is related to a potential ‘double correction’ of patient characteristics that were included in both the propensity score matching procedure as well as the outcome regression modelling, which could lead to overfitting of the regression models and an overestimation of the measured treatment effect. Covariates were included in the regression models since these characteristics could have differential impacts on the outcomes themselves, and our results show that the adjusted ORs were in fact slightly larger (showing a decreased effect size) when compared to unadjusted ORs, which show the difference in effect sizes of the matched populations alone.

Furthermore, another potential limitation in both the primary and ‘Bone-Rx’ cohorts is imbalanced comorbidity burden in BP user and non-user cohorts post-match. Table 1 shows there is differential prevalence of most co-morbid diseases despite matched cumulative CCI score between BP user and BP non-user cohorts. However, this limitation is in part addressed given (1) these covariates were controlled for during our regression analyses on study outcomes, and (2) that the key study findings were also observed in the ‘Osteo-Dx-Rx’ cohort, which matched based on individual comorbidities.

Additionally, limitations may be present due to misclassification bias of study outcomes due to the specific procedure/diagnostic codes used as well as the potential for residual confounding occurring for patient characteristics related to study outcomes that are unable to be operationalized in claims data, which would impact all cohort comparisons. For SARS-CoV-2 testing, procedure codes were limited to those testing for active infection, and therefore observations could be missed if they were captured via antibody testing (CPT 86318, 86328). These codes were excluded a priori due to the focus on the symptomatic COVID-19 population. Furthermore, for the COVID-19 diagnosis and hospitalization outcomes, all events were identified using the ICD-10 code for lab-confirmed COVID-19 (U07.1), and therefore events with an associated diagnosis code for suspected COVID-19 (U07.2) were not included. This was done to have a more stringent algorithm when identifying COVID-19-related events, and any impact of events identified using U07.2 is considered minimal, as previous studies of the early COVID-19 outbreak have found that U07.1 alone has a positive predictive value of 94% (Kluberg et al., 2022), and for this study U07.1 captured 99.2%, 99.0%, and 97.5% of all COVID-19 patient-diagnoses for the primary, ‘Bone-Rx’, and ‘Osteo-Dx-Rx’ cohorts, respectively.

Another potential limitation of this study relates to the positivity assumption, which when building comparable treatment cohorts is violated when the comparator population does not have an indication for the exposure being modelled (Petersen et al., 2012). This limitation is present in the primary cohort comparisons between BP users and BP non-users, as well as in the sensitivity analyses involving other preventive medications. This limitation, however, is mitigated by the fact that the outcomes in this study are related to infectious disease and are not direct clinical outcomes of known treatment benefits of BPs. The fact that the clinical benefits being assessed – the impact of BPs on COVID-related outcomes – was essentially unknown clinically at the time of the study data minimizes the impact of violation of the positivity assumption. Furthermore, our sensitivity analyses involving the ‘Bone-Rx’ and ‘Osteo-Dx-Rx’ cohorts did not suffer this potential violation, and the results from those analyses support those from the primary analysis cohort comparisons.

Moreover, we note that the propensity score-matched BP users and BP non-users in the primary analysis cohort mainly consisted of older females. According to the CDC,~75% and 95% of US women between 60–69 and 70–79 suffer from either low bone mass or osteoporosis, respectively (https://www.cdc.gov/nchs/data/databriefs/db93.pdf). Essentially all women (and 70% of men) above age 80 suffer from these conditions, which often go undiagnosed. Women aged 60 and older represent ~75% of our study population (Table 1). Although bone density measurements are not available for non-BP users in the matched primary cohort, there is a high probability that the incidence of osteoporosis and/or low bone mass in these patients was similar to the national average. Thus, BP therapy would have been indicated for most non-BP users in the matched primary cohort, and arguably, for these patients the positivity assumption was not violated.

One large potential bias to consider when comparing BP users to BP non-users is the healthy adherer effect, whereby adherence to drug therapy is associated with overall healthier behavior (Dormuth et al., 2009; Curtis et al., 2011). During the COVID-19 pandemic, this could have potentially resulted in differences between BP users and non-users such as, for example, adherence to mask-wearing, hand washing, or social distancing. However, if this effect accounted for the observed association between BP use and COVID-19-related outcomes, one would expect that users of other preventive medications would show similar associations. However, as discussed above, other preventive drug classes had a variable directional impact on the odds of COVID-19-related events, and sub-analyses within each drug class identified a strong association between concomitant BP use and decreased COVID-19-related events (Figure 6B–E). These analyses were based on the assumption that the association of unmeasured confounders with other drugs is comparable in magnitude and quality as for BPs. Taken together, these results suggest the observed association between BP use and COVID-19-related outcomes cannot solely be attributed to general behaviors associated with the healthy adherer effect.

Notably, several observational studies have reported that the use of one of our comparator preventive drug classes, statins, is associated with a lower risk of mortality in hospitalized COVID-19 patients (Israel et al., 2021; Lohia et al., 2021; Zhang et al., 2020). Indeed, statins are currently being tested as an adjunct therapy for COVID-19 (NCT04380402). In our study population, statin use was associated with moderately decreased odds of SARS-CoV-2 testing and COVID-19 diagnosis, though at a much smaller magnitude than BPs, and was not consistently associated with reduced odds of COVID-19-related hospitalizations. Our analysis did not address the clinical course of hospitalized patients, so these results are not necessarily conflicting. However, we note that in our primary cohort, as many as 15.2% of statin users concomitantly used a BP. Indeed, within statin users, stratification by BP use revealed that the decreased odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations remained regardless of statin use. Future studies on disease outcomes of hospitalized COVID-19 patients with antecedent use of BPs and statins alone or in combination are needed to clarify the effects of each drug class.

The differential association of amino-BPs versus statins with COVID-19 related outcomes is somewhat unexpected because both target the same biochemical pathway, albeit at different enzymatic steps (Xia et al., 2018). Statins block HMG-CoA reductase, the first and key rate-limiting enzyme in the mevalonate pathway (Istvan and Deisenhofer, 2001). Amino-BPs, which account for >99% of BPs prescribed in our study, inhibit a downstream enzyme in the same metabolic pathway, farnesyl pyrophosphate synthase (FPPS), which converts geranyl pyrophosphate to farnesyl pyrophosphate (Kavanagh et al., 2006). FPPS blockade disrupts protein prenylation and interferes with cytoskeletal rearrangement, membrane ruffling and vesicular trafficking in osteoclasts, thus preventing bone resorption (Russell, 2007). However, the anti-osteolytic activity of BPs per se is unlikely to account for the observed association between BP use and decreased incidence of COVID-19 and, more broadly, respiratory tract infections, because patients treated with non-BP anti-resorptive bone health medications have higher odds of respiratory infections (Sing et al., 2020 and this study).

Another consequence of mevalonate pathway inhibition by both statins and amino-BPs is arrested endosomal maturation in antigen-presenting cells resulting in enhanced antigen presentation, T cell activation and humoral immunity (Xia et al., 2018). In addition to this adjuvant-like effect, FPPS blockade by amino-BPs causes the intracellular accumulation of the enzyme’s substrate, isopentyl diphosphate (IPP), in myeloid leukocytes, which then stimulate Vγ9Vδ2 T cells (Wang et al., 2011; Nada et al., 2017), a large population of migratory innate lymphocytes in humans that are thought to play an important role in host defense against infectious pathogens (Ribot et al., 2021), including SARS-CoV-16. Experiments in humanized mice that were challenged with influenza viruses have shown that amino-BP-induced expansion of Vγ9Vδ2 T cells markedly improves viral control and mitigates disease severity and mortality (Tu et al., 2011; Zheng et al., 2015). However, since statins act upstream of FPPS, they are expected to inhibit IPP synthesis and, hence, have been shown to counteract the stimulatory effect of amino-BPs on Vγ9Vδ2 T cells (Wang et al., 2011). However, statins and amino-BPs do not always antagonize each other. In vitro, concomitant statin and amino-BP use has been shown to be synergistic in inhibition of cancer cell growth, but mainly through downstream inhibition of geranylgeranyl transferases and subsequent protein prenylation by statins (Abdullah et al., 2017). The fact that the observed reduction in COVID-19-related outcomes in BP users was not altered by concomitant statin use implies that the apparent protective effects of amino-BPs may not rely solely on stimulation of Vγ9Vδ2 T cells. Indeed, in mice (in which BPs are not known to stimulate γδ T cells), BPs potently boost systemic and mucosal antiviral antibody and T cell responses (Tonti et al., 2013). This effect was also seen with non-nitrogenous BPs, which do not antagonize FPPS (Tonti et al., 2013). In the present study, the number of patients who used non-nitrogenous BPs was less than 20, and therefore too small to determine any impact on COVID-19-related outcomes. Nevertheless, in aggregate, these clinical and pre-clinical findings raise the possibility that BPs may exert (at least some) immuno-stimulatory effects by engaging an as yet unidentified additional pathway, regardless of their nitrogen content.

Irrespective of the precise molecular mechanism of action, BPs have been reported to exert a plethora of effects on additional immune cell populations in humans, including NK cells (Sarhan et al., 2017) and regulatory T cells (Liu et al., 2016). Moreover, studies of patients treated with amino-BPs found impaired chemotaxis and generation of reactive oxygen species by neutrophils (Kuiper et al., 2012; Chadwick et al., 2020), a population of inflammatory cells whose dysregulated recruitment and activation are strongly implicated in the pathogenesis of severe COVID-19 (Meizlish et al., 2021; Reusch et al., 2021). Thus, BPs may provide therapeutic benefits during infections with SARS-CoV-2 through modulation of both innate and adaptive immune responses. However, further studies to directly test these pleiotropic immuno-modulatory effects of BPs and to assess their relative contribution to the host response to SARS-CoV-2 infection are needed.

We conclude that, despite several caveats discussed above, the association between BP use and decreased odds of COVID-19-related endpoints was robust in analyses comparing BP users to BP non-users. Large differences were detected regardless of age, sex or geographic location that remained robust when using multiple approaches to address unmeasured confounding and/or potential sources of bias. These retrospective findings strongly suggest that BPs should be considered for prophylactic use in individuals at risk of SARS-CoV-2 infection. However, additional well-controlled prospective clinical studies will be needed to rigorously assess whether the observed reduction in COVID-19-related outcomes is directly caused by BPs and remains true in patient populations not commonly prescribed BPs.

A number of BPs are globally available as relatively affordable generics that are generally well tolerated and could be prescribed for off-label use. Rare, but severe adverse events that have been linked to BP use include osteonecrosis of the jaw (Migliorati et al., 2006) and atypical femur fractures (Saita et al., 2015), which are both associated with long-term BP therapy. In this context, it is important to consider the relationship between the timing of BP exposure and COVID-19-related outcomes. Remarkably, BP users of alendronic acid whose prescription ended more than eight months prior to the observation period, as well as users who initiated alendronic acid therapy immediately preceding the observation period, had similarly decreased odds of COVID-19-related outcomes (Figure 3B). A likely explanation for the observed long-term protection after transient BP use may be the well-documented retention of BPs in bone resulting in half-lives of several years (Cremers et al., 2019). Small amounts of stored BPs are continuously released, especially in regions of high bone turnover, which may result in persistent exposure of immune cells either systemically or preferentially in bone marrow, a site of active immune cell trafficking (Mazo et al., 2005; Zhao et al., 2012) where anti-viral immune responses can be initiated in response to respiratory infection (Hermesh et al., 2010). Thus, BP use at the time of infection may not be necessary for protection against COVID-19. Rather, our results suggest that prophylactic BP therapy may be sufficient to achieve a potentially rapid and sustained immune modulation resulting in profound mitigation of the incidence and/or severity of infections by SARS-CoV-2.

Appendix 1

Study Methods

Section 1: Variable Assignment

Outcomes

The following details the identification algorithms and associated codes that were used to identify outcomes of interest, including COVID-19-related as well as the exploratory outcomes that were assessed during sensitivity analyses.

Primary outcomes

SARS-CoV-2 testing

  • Any medical services claim with a procedure code indicating polymerase chain reaction (PCR) testing for active SARS-CoV-2 infection 3/1/2020-6/30/2020

  • Identified using HCPCS codes: 87635, 87636, 87637

COVID-19 diagnosis

  • Any medical services claim with a diagnosis code indicating COVID-19 3/1/2020-6/30/2020

  • Identified using ICD-10 code U07.1x

COVID-19-related hospitalization

  • Any medical services claim occurring in an inpatient setting with a diagnosis code indicating COVID-19 3/1/2020-6/30/2020

  • Identified using ICD-10 code U07.1x

Exploratory outcomes (study observation period)

Acute cholecystitis-related service

  • Any medical services claim occurring in an emergency room/inpatient setting with a diagnosis indicating acute cholecystitis 3/1/2020-6/30/2020

  • Identified using ICD-10 codes K81.0x

Acute pancreatitis-related service

  • Any medical services claim occurring in an emergency room/inpatient setting with a diagnosis indicating acute pancreatitis 3/1/2020-6/30/2020

  • Identified using ICD-10 codes K85.x

Exploratory outcomes (2019)

Acute cholecystitis-related service

  • Any medical services claim occurring in an emergency room/inpatient setting with a diagnosis indicating acute cholecystitis 7/1/2019-12/31/2019

  • Identified using ICD-10 codes K81.0x

Acute pancreatitis-related service

  • Any medical services claim occurring in an emergency room/inpatient setting with a diagnosis indicating acute pancreatitis 7/1/2019-12/31/2019

  • Identified using ICD-10 codes K85.x

Acute bronchitis-related service

  • Any medical services claim with a diagnosis indicating acute bronchitis 7/1/2019-12/31/2019

  • Identified using ICD-10 codes J20.x-J21.x

Acute pneumonia-related service

  • Any medical services claim with a diagnosis indicating acute bronchitis 7/1/2019-12/31/2019

  • Identified using ICD-10 codes J13.x-J18.x

Osteonecrosis

Osteonecrosis

  • Any medical services claim with a diagnosis indicating drug-induced osteonecrosis 1/1/2019-6/30/2020

  • Identified using ICD-10 codes M87.1x

Drug-exposure assignment

The following details the identification algorithms and associated inputs used for drug-exposure classification of study subjects into users/non-users of bisphosphonates, non-bisphosphonates osteoporosis medications, statins, antihypertensives, non-insulin antidiabetics, and antidepressants.

Bisphosphonates

  • Any outpatient prescription or in-office dispensing 1/1/2019-2/29/2020

  • Drugs included: alendronate, alendronic acid, etidronate, ibandronate, ibandronic acid, pamidronate, risedronate, and zoledronic acid

Non-BP anti-resorptive bone health medications

  • Any outpatient prescription or in-office dispensing 1/1/2019-2/29/2020

  • Drugs included: denosumab, calcitonin, raloxifene, romosozumab-aqqg, teriparatide, abaloparatide, or bazedoxifene

Statins

  • Any outpatient prescription 1/1/2019-2/29/2020

  • Drugs included: pravastatin, rosuvastatin, fluvastatin, atorvastatin, pitavastatin, or simvastatin

Antihypertensives

  • Any non-ophthalmic, non-injection, outpatient prescription claim for a beta-blocker, calcium channel blocker, or renin angiotensin system antagonist 1/1/2019-2/29/2020

  • Drugs included: acebutolol, atenolol, betaxolol, bisoprolol, carvedilol, labetalol, metoprolol, nadolol, nebivolol, penbutolol, pindolol, propranolol, timolol, amlodipine, diltiazem, felodipine, isradipine, nicardipine, nifedipine, nisoldipine, verapamil, aliskiren, azilsartan, benazepril, candesartan, captopril, enalapril, eprosartan, fosinopril, irbesartan, lisinopril, losartan, moexipril, olmesartan, perindopril, quinapril, ramipril, sacubitril, telmisartan, trandolapril, valsartan

Antidiabetics

  • Any outpatient prescription claim for a non-insulin antidiabetic medication 1/1/2019-2/29/2020

  • Drugs included: metformin, chlorpropamide, glimepiride, glipizide, glyburide, tolazamide, tolbutamide, pioglitazone, rosiglitazone, alogliptin, linagliptin, saxagliptin, sitagliptin, albiglutide, dulaglutide, exenatide, liraglutide, lixisenatide, semaglutide, nateglinide, repaglinide, canagliflozin, dapagliflozin, empagliflozin, ertugliflozin

Antidepressants

  • Any outpatient prescription claim for a selective serotonin reuptake inhibitor, norepinephrine-dopamine reuptake inhibitor, serotonin-norepinephrine reuptake inhibitor, tricyclic, tetracyclic, modified cyclic, or MAO inhibitor medication 1/1/2019-2/29/2020

  • Drugs included: amoxapine, bupropion, citalopram, clomipramine, desipramine, desvenlafaxine, doxepin, duloxetine, escitalopram, esketamine, fluoxetine, fluvoxamine, imipramine, isocarboxazid, levomilnacipran, maprotiline, mirtazapine, nefazodone, nortriptyline, paroxetine, phenelzine, protriptyline, selegiline, sertraline, tranylcypromine, trazodone, trimipramine, venlafaxine, vilazodone, vortioxetine

Charlson comorbidity condition assignment

The following ICD-10 codes were used to assign the CCI condition-specific indicators that are used to calculate the overall CCI score. The time period used for identification of condition-specific indicators was the entire pre-observation period (1/1/2019-2/29/2020).

Myocardial infarction

  • ICD-10 codes: I21.x, I22.x, I25.2

Congestive heart failure

  • ICD-10 codes: I09.9, I11.0, I13.0, I13.2, I25.5, I42.0, I42.5 - I42.9, I43.x, I50.x, P29.0

Peripheral vascular disease

  • ICD-10 codes: I70.x, I71.x, I73.8, I73.9, I77.1, I79.0, I79.2, K55.1, K55.8, K55.9, Z95.8, Z95.9

Cerebrovascular disease

  • ICD-10 codes: G45.x, G46.x, H34.0, I60.x-I69.x

Dementia

  • ICD-10 codes: F00.x - F03.x, F05.1, G30.x, G31.1

Chronic pulmonary disease

  • ICD-10 codes: I27.8, I27.9, J40.x - J47.x, J60.x - J67.x, J68.4, J70.1, J70.3

Rheumatologic disease

  • ICD-10 codes: M05.x, M06.x, M31.5, M32.x - M34.x, M35.1, M35.3, M36.0

Peptic ulcer disease

  • ICD-10 codes: K25.x-K28.x

Mild liver disease

  • ICD-10 codes: B18.x, K70.0 - K70.3, K70.9, K71.3 - K71.5, K71.7, K73.x, K74.x, K76.0, K76.2 - K76.4, K76.8, K76.9, Z94.4

Diabetes without chronic complications

  • ICD-10 codes: E10.0, E10.1, E10.6, E10.8, E10.9, E11.0, E11.1, E11.6, E11.8, E11.9, E12.0, E12.1, E12.6, E12.8, E12.9, E13.0, E13.1, E13.6, E13.8, E13.9, E14.0, E14.1, E14.6, E14.8, E14.9

Diabetes with chronic complications

  • ICD-10 codes: E10.2 - E10.5, E10.7, E11.2 - E11.5, E11.7, E12.2 - E12.5, E12.7, E13.2 - E13.5, E13.7, E14.2 - E14.5, E14.7

Hemiplegia or paraplegia

  • ICD-10 codes: G04.1, G11.4, G80.1, G80.2, G81.x, G82.x, G83.0 - G83.4, G83.9

Renal disease

  • ICD-10 codes: I12.0, I13.1, N03.2 - N03.7, N05.2 - N05.7, N18.x, N19.x, N25.0, Z49.0 - Z49.2, Z94.0, Z99.2

Any tumor, leukemia, or lymphoma

  • ICD-10 codes: C00.x - C26.x, C30.x - C34.x, C37.x - C41.x, C43.x, C45.x - C58.x, C60.x - C76.x, C81.x - C85.x, C88.x, C90.x - C97.x

Moderate or severe liver disease

  • ICD-10 codes: I85.0, I85.9, I86.4, I98.2, K70.4, K71.1, K72.1, K72.9, K76.5, K76.6, K76.7

Metastatic solid tumor

  • ICD-10 codes: C77.x - C80.x

AIDS/HIV

  • ICD-10 codes: B20.x - B22.x, B24.x

Additional condition covariate assignment

The following details the ICD-10 diagnosis codes that were used to identify comorbid conditions. For all condition indicators classification was based on all medical claims occurring during the pre-observation period (1/1/2019-2/29/2020).

Osteoporosis: M80.x, M81.x, M82.x

  • Cardiovascular disease overall: I3x.x-I4x.x, I20.x-I28.x, I50.x-I52.x

  • Cancer: C0x.x - C9x.x

  • Chronic kidney disease (CKD)/ end-stage renal disease (ESRD): I12.0, I13.1, N03.2 - N03.7, N05.2 - N05.7, N18.x, N19.x, N25.0, Z49.0 - Z49.2, Z94.0, Z99.2

  • Chronic obstructive pulmonary disease (COPD): J43.x, J44.x

  • Dementia: F00.x - F03.x, F05.1, G30.x, G31.1

  • Depression: F32.x, F33.x

  • Dyslipidemia: E78.x

  • Heart failure: I50.x, I11.0xx, I13.0xx, I13.2xx

  • HIV/AIDS: B20.x - B22.x, B24.x

  • Hypertension: I10.x, I12.x, I11.9xx, I13.1xx

  • Liver disease: B18.x, K70.0 - K70.3, K70.9, K71.3 - K71.5, K71.7, K73.x, K74.x, K76.0, K76.2 - K76.4, K76.8, K76.9, Z94.4, I85.0, I85.9, I86.4, I98.2, K70.4, K71.1, K72.1, K72.9, K76.5, K76.6, K76.7

  • Obesity: E66.x

  • Sickle cell disease: D57.x

  • Stroke: I63.x

  • Type 2 diabetes: E11.x

Sensitivity Analysis (1): COVID-19-related outcomes in “Bone-Rx” cohort

Overview and rationale
  • The first sensitivity analysis was performed to validate the robustness of the primary findings by limiting all BP non-users to those who had used non-BP anti-resorptive bone health medications during the pre-observation period, thus yielding a more comparable comparator cohort that was also receiving bone health medication therapy.

  • The use of an active-comparator cohort was done to reduce the impact of unmeasured confounding that may have occurred in the primary analysis due to the use of the derived Charlson Comorbidity Index composite score as the only comorbidity matching covariate. Restriction of the patient population to users of any non-BP anti-resorptive bone health medication prior to propensity-score matching improves the probability of having drug user/non-user matches with more similar clinical characteristics.

  • This sensitivity analysis, further, also acted to increase the robustness and reliability of the matched user/non-user outcome comparisons since non-BP anti-resorptive bone health medication users represented the smaller portion of the total bone health medication-user population (“Bone-Rx” cohort) and therefore were matched to their best BP-user pair.

Analysis cohort definition(s)
  • Continuous medical and prescription insurance coverage 1/1/2019-6/30/2020

  • Patients with ≥1 claim for any anti-resorptive bone health medication 1/1/2019-2/29/2020

Exposures of interest
  • Patients were assigned into the BP user cohort if they had any claim 1/1/2019-2/29/2020 for one of the following: alendronate, alendronic acid, etidronate, ibandronate, ibandronic acid, pamidronate, risedronate, and zoledronic acid.

  • Patients were assigned into the non-BP any anti-resorptive bone health medication user cohort if: (1) they had any claim 1/1/2019-2/29/2020 for one of the following: denosumab, calcitonin, raloxifene, romosozumab-aqqg, teriparatide, abaloparatide, or bazedoxifene; and (2) they had no BP claims 1/1/2019-2/29/2020.

Outcomes
  • SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations

Cohort matching
  • Non-BP anti-resorptive bone health medication users were matched to BP users based on age, gender, insurance type, any PCP visit in 2019, and comorbidity score. Matching was performed within each region separately (northeast, midwest, south, west) and then combined as well as in NY-state alone.

Statistical analyses
  • Same as was performed for the primary analysis cohort.

Sensitivity Analysis (2): COVID-19-related outcomes in “Osteo-Dx-Rx” cohort

Overview and rationale
  • The second sensitivity analysis was performed to further assess the robustness of the primary analysis findings by performing a highly restricted comparator cohort matching that included patients diagnosed and treated for osteoporosis (“Osteo-Dx-Rx” cohort).

  • The relationship between COVID-19-related outcomes and BP-exposure was modelled after restricting anti-resorptive bone health medication users to those most likely to use BPs and matching BP non-users to BP users based on the presence of comorbid diagnoses within insurance type in four states with early COVID-19 spread representing each to further reduce confounding related to differences in demographic/clinical characteristics amongst BP users/non-users, confounding due to socioeconomic status (insurance type as proxy), and confounding due to differences in COVID-19-exposure risk based on geography.

Analysis cohort definition(s)
  • Continuous medical and prescription insurance coverage 1/1/2019-6/30/2020

  • Patients with ≥1 claim for any osteoporosis medication 1/1/2019-2/29/2020 who also met the following criteria: (i) female; (ii) age 51 or older; (iii) identified as residing in New York, Illinois, Florida, or California; and (iv) had ≥1 medical claim indicating a diagnosis of osteoporosis 1/1/2019-2/29/2020

Exposures of interest
  • Patients were assigned into the BP user cohort if they had any claim 1/1/2019-2/29/2020 for one of the following: alendronate, alendronic acid, etidronate, ibandronate, ibandronic acid, pamidronate, risedronate, and zoledronic acid.

  • Patients were assigned into the non-BP anti-resorptive bone health medication user cohort if: (1) they had any claim 1/1/2019-2/29/2020 for one of the following: denosumab, calcitonin, raloxifene, romosozumab-aqqg, teriparatide, abaloparatide, or bazedoxifene; and (2) they had no BP claims 1/1/2019-2/29/2020.

Outcomes
  • SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations

Cohort matching
  • Non- anti-resorptive bone health medication users were matched to BP users based on age, PCP visit in 2019, and the presence of the following comorbid conditions (assigned using ICD-10 codes on claims occurring 1/1/2019-2/29/2020): cancer, chronic obstructive pulmonary disease, depression, dyslipidaemia, heart failure, hypertension, obesity, and type 2 diabetes.

  • Matching was performed within each state when stratified by insurance type (commercial, dual, Medicaid, Medicare).

Statistical analyses

Multivariate logistic regression analyses, modelled separately for each COVID-19-related outcome of interest, were performed on the unmatched and matched samples after combining all patient observations. In addition to the key exposure variable (indicating BP user versus non-BP user), the regression model also included demographic/clinical covariate for age group, region, insurance type, PCP visit in 2019, and the following comorbid conditions: osteoporosis, cancer, chronic obstructive pulmonary disease, depression, dyslipidaemia, hypertension, obesity, type 2 diabetes, cardiovascular disease overall, sickle cell anemia, stroke, dementia, HIV/AIDS, chronic kidney disease/end-stage renal disease, and liver disease.

Sensitivity Analysis (3): Association of BP-use with exploratory negative control outcomes

Overview and rationale
  • The third sensitivity analysis was performed to assess the relationship between BP-use and outcomes not anticipated to be impacted by the pharmacological mechanism of BPs.

  • This was performed by modelling the relationship between BP-exposure and other outcomes occurring (1) during the study observation, and (2) during the second half of 2019 among BP users with claims during the first half of 2019 and their previously-assigned BP non-user matched pair, in the primary, “Bone-Rx”, and “Osteo-Dx-Rx” cohorts.

  • Outcomes modelled included any acute cholecystitis-related or acute pancreatitis-related inpatient/emergency-room (ER) service, used as exploratory outcomes not predicted to be modulated by BP exposure to assess the validity of the core COVID-19-related outcomes.

Analysis cohort definition(s)
  • Patients who were included in the primary analysis cohort for assessment of (1) outcomes occurring during the study observation period; for (2) outcomes assessed during the second half of 2019 the cohort was restricted to among BP users with claims during the first half of 2019 and their previously-assigned BP non-user matched pair.

  • Patients who met all eligibility criteria to be included in the ‘Bone-Rx’ cohort for assessment of (1) outcomes occurring during the study observation period; for (2) outcomes assessed during the second half of 2019 the cohort was restricted to among BP users with claims during the first half of 2019 and their previously-assigned BP non-user matched pair.

  • Patients who met all eligibility criteria to be included in the ‘Osteo-Dx-Rx’ cohort for assessment of (1) outcomes occurring during the study observation period; for (2) outcomes assessed during the second half of 2019 the cohort was restricted to among BP users with claims during the first half of 2019 and their previously-assigned BP non-user matched pair.

Exposures of interest
  • For the primary analysis cohort, the BP user / BP non-user assignment was the same as used in the core analyses.

  • For the “Bone-Rx” and “Osteo-Dx-Rx” cohorts, assignment was the same as used in those analyses stratifying medication users into BP users and non-BP medication users.

Outcomes
  • Any medical claim from an ER/inpatient setting with a diagnosis indicating acute cholecystitis (ICD-10 code K81.0x) occurring 3/1/2020-6/30/2020 (observation period)

  • Any medical claim from an ER/inpatient setting with a diagnosis indicating acute pancreatitis (ICD-10 code K85.x) occurring 3/1/2020-6/30/2020 (observation period)

  • Any medical claim from an ER/inpatient setting with a diagnosis indicating acute cholecystitis (ICD-10 code K81.0x) occurring 7/1/2019-12/31/2019 (2019)

  • Any medical claim from an ER/inpatient setting with a diagnosis indicating acute pancreatitis (ICD-10 code K85.x) occurring 7/1/2019-12/31/2019 (2019)

Cohort matching

NA; all cohorts previously matched.

Statistical analyses

Multivariate logistic regression analyses were performed using the same methodologies employed when assessing COVID-19 outcomes that were cohort-build-specific (i.e. followed previous approach detailed for each respective cohort build) to assess the odds of acute cholecystitis or acute pancreatitis.

Sensitivity Analysis (4): Association of BP-use with exploratory positive control outcomes in 2019

Overview and rationale
  • The fourth sensitivity analysis was performed to assess the relationship between BP-use and select outcomes occurring in 2019 to validate the theorized BP mechanism of action.

  • This was performed by modelling the relationship between BP-exposure in the first half of 2019 and other outcomes occurring during the second half of 2019 in the primary, “Bone-Rx”, and “Osteo-Dx-Rx” cohorts, specifically medical services for other infectious respiratory conditions (acute bronchitis, pneumonia), used to assess the validity of the relationship between BP-use and decreased respiratory infections.

Analysis cohort definition(s)

- The following criteria were applied to all three cohort build variations (primary analysis cohort, “Bone-Rx” cohort, “Osteo-Dx-Rx” cohort): (i) BP users were restricted to those with any BP claim 1/1/2019-6/30/2019, and the remaining previously-classified BP-user patients with their first BP-claim date occurring on/after 7/1/2019 were excluded; (ii) BP non-users were restricted to their BP-user matched-pair previously assigned.

Exposures of interest
  • In all cohort build variations, the previously-classified BP user cohorts were restricted to those with any BP-claim 1/1/2019-6/30/2019; all other previously-classified BP users were excluded.

Outcomes
  • Any medical claim with a diagnosis indicating acute bronchitis (ICD-10 code J20.x-J21.x) occurring 7/1/2019-12/31/2019

  • Any medical claim with a diagnosis indicating pneumonia (ICD-10 code J13.x-J18.x) occurring 7/1/2019-12/31/2019

Cohort matching
  • NA; all cohorts previously matched.

Statistical analyses
  • Multivariate logistic regression analyses were performed using the same methodologies employed when assessing COVID-19-related outcomes that were cohort-build-specific (i.e. followed previous approach detailed for each respective cohort build) to assess the odds of acute bronchitis, or pneumonia.

Sensitivity Analysis (5): Association between use of other drug classes and COVID-19-related outcomes

Overview and rationale
  • The fifth sensitivity analysis was performed to assess whether the observed protective effect of BPs may be associated with general healthier behaviours in patients using any medication rather than specifically BP use. To assess this unmeasured confounding due to the healthy adherer effect, which is a type of potential bias where patients may have better outcomes due to their heathier behaviours and not better outcomes related to active drug treatment itself, the first sensitivity analysis evaluated the association between use of other preventive medications (statin, antihypertensive, antidiabetic, antidepressant) and COVID-19-related outcomes were evaluated.

  • This was performed following the same techniques used in the primary cohort matching and analyses but when assigned drug exposure cohorts based on the use of statin, antihypertensive, antidiabetic, or antidepressant medications. The consistency of methods was done to permit direct comparison on the association between drug-use and COVID-19-related outcomes to assess whether the healthy adherer effect alone accounts for the decrease in the odds of COVID-19 outcomes when comparing BP users to non-users in the primary analysis. Evidence to support the contention that the HAE is a significant source of unmeasured confounding would necessitate that other drug classes display a similar statistically significant trend and/or magnitude when comparing drug users to non-users. Variability in directional impact, magnitude, and/or statistical significance would, conversely, suggest that the healthy adherer effect itself does not account for the differences seen when comparing BP users to BP non-users.

  • This sensitivity analysis, additionally, also employed a unique nested-matching technique wherein BP users were matched to BP non-users within the other-medication-class matched populations when stratified into the already matched but mutually exclusive user/non-user cohorts. This was performed to: (1) assess whether the decreased odds of COVID-19-realted outcomes in BP users compared to BP non-users was robust, even amongst cohorts displaying an increase in the odds of COVID-19-related outcomes; and (2) to assess whether the magnitude of decrease in odds of COVID-19-related outcomes amongst BP users compared to BP non-users seen in the primary analysis is impacted by use of other medication classes, including some that have also been identified as being associated with a reduced incidence and/or severity of COVID-19-related outcomes.

Analysis cohort definition(s)
  • Continuous medical and prescription insurance coverage 1/1/2019-6/30/2020 (all)

  • Patients with any claim for another drug class of interest (statin, antihypertensive, antidiabetic, antidepressant) medication 1/1/2019-2/29/2020 were classified users

  • Among the propensity-score matched drug user/non-user cohorts, a further stratification and propensity-score matching based on BP use 1/1/2019-2/29/2020 to yield the following: (i) drug user/BP user matched to drug user/BP non-user, (ii) drug non-user/BP user matched to drug non-user/BP non-user.

Exposures of interest
  • Patients were assigned into the statin user cohort if they had any claim 1/1/2019-2/29/2020 for one of the following: pravastatin, rosuvastatin, fluvastatin, atorvastatin, pitavastatin, or simvastatin

  • Patients were assigned into the antihypertensive user cohort if they had any non-ophthalmic, non-injection claim 1/1/2019-2/29/2020 for a beta blocker, calcium channel blocker, or renin-angiotensin system antagonist medication.

  • Patients were assigned into the antidiabetic user cohort if they had any claim 1/1/2019-2/29/2020 for one of the following non-insulin medications: metformin, chlorpropamide, glimepiride, glipizide, glyburide, tolazamide, tolbutamide, pioglitazone, rosiglitazone, alogliptin, linagliptin, saxagliptin, sitagliptin, albiglutide, dulaglutide, exenatide, liraglutide, lixisenatide, semaglutide, nateglinide, repaglinide, canagliflozin, dapagliflozin, empagliflozin, ertugliflozin

  • Patients were assigned into the antidepressant user cohort if they had any claim 1/1/2019-2/29/2020 for one of the following: amoxapine, bupropion, citalopram, clomipramine, desipramine, desvenlafaxine, doxepin, duloxetine, escitalopram, esketamine, fluoxetine, fluvoxamine, imipramine, isocarboxazid, levomilnacipran, maprotiline, mirtazapine, nefazodone, nortriptyline, paroxetine, phenelzine, protriptyline, selegiline, sertraline, tranylcypromine, trazodone, trimipramine, venlafaxine, vilazodone, vortioxetine

Outcomes
  • SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations

Cohort matching
  • For the larger drug-class analyses, matching was performed following the same methods used in the primary analysis: users were matched to non-users based on age, gender, insurance type, any PCP visit in 2019, and comorbidity score. Matching was performed within each region separately (northeast, midwest, south, west) and then combined, as well as in NY-state alone.

  • Following this matching procedure, a nested BP user to BP non-user propensity score match was then performed on the aforementioned matched populations (i.e. within the separate and already matched statin user and statin non-user populations). Matching was performed using the same list of demographic/clinical characteristics, and was also performed within each region separately (northeast, midwest, south, west) and then combined as well as in NY-state alone.

Statistical analyses
  • Same as was performed for the primary analysis cohort.

Appendix 2

Additional study results; cohort characteristics pre/post match

Primary analysis study population

Northeast region

A total of 2,152,560 patients identified as residing in the northeast were included in the unmatched primary analysis cohort comparisons, of which 119,728 (5.6%) and 2,032,832 (94.4%) were classified as BP users and BP non-users, respectively (Appendix 2—table 1). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (97.5% age ≥51 vs 49.8%; P<0.001), predominantly female (90.5% vs 57.4%; P<0.001), with higher comorbidity burden (mean CCI = 0.93 versus 0.65; P<0.001), insured by Medicare (46.5% vs 18.0%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (58.3% vs 42.8%; P<0.001). Propensity-score matching yielded 119,494 BP users and 119,494 BP non-users with no significant differences across examined characteristics. A total of 234 BP users from the northeast region in the unmatched primary analysis cohort were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Midwest region

A total of 1,467,802 patients identified as residing in the midwest were included in the unmatched primary analysis cohort comparisons, of which 75,967 (5.2%) and 1,391,835 (94.8%) were classified as BP users and BP non-users, respectively (Appendix 2—table 2). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (96.6% age ≥51 vs 44.0%; P<0.001), predominantly female (90.3% vs 57.1%; P<0.001), with higher comorbidity burden (mean CCI = 0.99 versus 0.56; P<0.001), insured by Medicare (43.6% vs 14.5%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (62.2% vs 51.0%; P<0.001). Propensity-score matching yielded 75,901 BP users and 75,901 BP non-users with no significant differences across examined characteristics. A total of 66 BP users from the midwest region in the unmatched primary analysis cohort were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

South region

A total of 3,042,604 patients identified as residing in the south were included in the unmatched primary analysis cohort comparisons, of which 160,886 (5.3%) and 2,881,718 (94.7%) were classified as BP users and BP non-users, respectively (Appendix 2—table 3). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (96.8% age ≥51 vs 39.2%; P<0.001), predominantly female (90.6% vs 57.4%; P<0.001), with higher comorbidity burden (mean CCI = 0.86 versus 0.55; P<0.001), insured by Medicare (41.0% vs 11.3%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (66.1% vs 49.2%; P<0.001). Propensity-score matching yielded 159,704 BP users and 159,704 BP non-users with no significant differences across examined characteristics. A total of 1,182 BP users from the south region in the unmatched primary analysis cohort were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

West region

A total of 1,243,637 patients identified as residing in the west were included in the unmatched primary analysis cohort comparisons, of which 95,470 (7.7%) and 1,148,167 (92.3%) were classified as BP users and BP non-users, respectively (Appendix 2—table 4). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (97.8% age ≥51 vs 43.5%; P<0.001), predominantly female (88.7% vs 56.4%; P<0.001), with higher comorbidity burden (mean CCI = 1.08 versus 0.66; P<0.001), insured by Medicare (43.5% vs 11.0%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (67.7% vs 45.3%; P<0.001). Propensity-score matching yielded 95,267 BP users and 95,267 BP non-users with no significant differences across examined characteristics. A total of 203 BP users from the west region in the unmatched primary analysis cohort were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

New York State

A total of 968,296 patients identified as residing in New York state were included in the unmatched primary analysis NY-state restricted cohort, of which 50,035 (5.2%) and 918,261 (94.8%) were classified as BP users and BP non-users, respectively (Appendix 2—table 5). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.1% age ≥51 vs 50.7%; P<0.001), predominantly female (90.9% vs 57.5%; P<0.001), with higher comorbidity burden (mean CCI = 0.95 versus 0.63; P<0.001), insured by Medicare (57.7% vs 19.5%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (62.7% vs 45.3%; P<0. 001). Propensity-score matching yielded 49,862 BP users and 49,862 BP non-users with no significant differences across examined characteristics. A total of 173 BP users from the unmatched New York state primary analysis cohort were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Bone-Rx analysis study population

All observations (all regions combined)

A total of 502,895 patients were included in the unmatched “Bone-Rx” analysis cohort comparisons, of which 452,051 (89.9%) and 50,844 (10.1%) were classified as BP users and BP non-users, respectively (Appendix 2—table 17). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were younger (47.9% age ≥71 vs 55.2%; P<0.001), predominantly female (90.1% vs 87.2%; P<0.001), with a lower comorbidity burden (mean CCI = 0.95 vs 1.99; P<0.001), with a larger proportion of patients residing in the west (21.1% versus 15.8%; P<0.001), a lower proportion covered by Medicare (43.4% vs 47.5%; P<0.001), and a lower proportion have had a primary-care physician (PCP) visit in 2019 (63.8% vs 64.3%; P=0.009). Propensity-score matching yielded 50,498 BP users and 50,498 BP non-users with no significant differences across examined characteristics. A total of 346 BP non-users from the unmatched “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

Northeast region

A total of 135,867 patients identified as residing in the northeast were included in the unmatched “Bone-Rx” analysis cohort comparisons, of which 119,728 (88.1%) and 16,139 (11.9%) were classified as BP users and BP non-users, respectively (Appendix 2—table 18). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics except for any PCP visit in 2019 (P=0.95). Compared to BP non-users, BP users were younger (48.1% age ≥71 vs 54.8%; P<0.001), predominantly female (90.5% vs 87.5%; P<0.001), with a lower comorbidity burden (mean CCI = 0.93 vs 1.97; P<0.001), and a lower proportion insured by Medicare (46.5% vs 54.0%; P<0.001). Propensity-score matching yielded 15,993 BP users and 15,993 BP non-users with no significant differences across examined characteristics. A total of 146 BP non-users from the northeast region in the unmatched “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

Midwest region

A total of 85,391 patients identified as residing in the midwest were included in the unmatched “Bone-Rx” analysis cohort comparisons, of which 75,967 (89.0%) and 9,424 (11.0%) were classified as BP users and BP non-users, respectively (Appendix 2—table 19). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were younger (43.0% age ≥71 vs 54.1%; P<0.001), predominantly female (90.3% versus 86.1%; P<0.001), with a lower comorbidity burden (mean CCI = 0.99 versus 2.12; P<0.001), had a lower proportion insured by Medicare (43.6% versus 51.9%; P<0.001), with a lower proportion having a primary-care physician (PCP) visit in 2019 (62.2% vs 64.7%; P<0.001). Propensity-score matching yielded 9,360 BP users and 9,360 BP non-users with no significant differences across examined characteristics. A total of 64 BP non-users from the midwest region in the unmatched “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

South region

A total of 178,118 patients identified as residing in the south were included in the unmatched “Bone-Rx” analysis cohort comparisons, of which 160,886 (90.3%) and 17,232 (9.7%) were classified as BP users and BP non-users, respectively (Appendix 2—table 20). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics except for any PCP visit in 2019 (P=0.45). Compared to BP non-users, BP users were younger (46.6% age ≥71 vs 53.3%; P<0.001), predominantly female (90.6% vs 88.1%; P<0.001), with a lower comorbidity burden (mean CCI = 0.86 vs 1.86; P<0.001), and a lower proportion insured by Medicare (41.0% vs 44.0%; P<0.001). Propensity-score matching yielded 17,140 BP users and 17,140 BP non-users with no significant differences across examined characteristics. A total of 92 BP non-users from the south region in the unmatched “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

West region

A total of 103,519 patients identified as residing in the west were included in the unmatched “Bone-Rx” analysis cohort comparisons, of which 95,470 (92.2%) and 8,049 (7.8%) were classified as BP users and BP non-users, respectively (Appendix 2—table 21). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were younger (54.1% age ≥71 vs 61.6%; P<0.001), predominantly female (88.7% vs 86.2%; P<0.001), with a lower comorbidity burden (mean CCI = 1.08 versus 2.17; P<0.001), insured by Medicare (43.5% vs 36.9%; P<0.001), with a lower proportion having a primary-care physician (PCP) visit in 2019 (67.7% vs 71.6%; P<0.001). Propensity-score matching yielded 8,005 BP users and 8,005 BP non-users with no significant differences across examined characteristics. A total of 44 BP non-users from the west region in the unmatched “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

New York State

A total of 57,397 patients identified as residing in New York state were included in the unmatched “Bone-Rx” analysis NY-state restricted cohort, of which 50,035 (87.2%) and 7,362 (12.8%) were classified as BP users and BP non-users, respectively (Appendix 2—table 22). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics except for any PCP visit in 2019 (P=0.35). Compared to BP non-users, BP users were younger (53.2% age ≥11 vs 54.5%; P<0.001), predominantly female (90.9% vs 89.5%; P<0.001), with a lower comorbidity burden (mean CCI = 0.95 vs 1.81; P<0.001), and a higher proportion insured by Medicaid (18.3% vs 13.8%; P<0.001). Propensity-score matching yielded 7,254 BP users and 7,254 BP non-users with no significant differences across examined characteristics. A total of 108 BP non-users from the unmatched New York state “Bone-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

Osteo-Dx-Rx analysis study population

A total of 60,043 female patients age ≥51 with a diagnosis of osteoporosis who resided in New York (NY), Illinois (IL), Florida (FL), or California (CA) were included in the unmatched “Osteo-Dx-Rx” analysis cohort comparison, of which 51,651 (86.0%) and 8,392 (14.0%) were classified as BP users and BP non-users, respectively (Appendix 2—table 23). Prior to propensity-score matching, which was performed within each state by insurance type, there were significant differences across all demographic and clinical characteristics except the proportion of patients with a diagnosis of dyslipidemia (P=0.08). Compared to BP non-users, BP users were younger (18.8% age ≥81 vs 26.0%; P<0.001), with a larger proportion of patients residing in CA (42.5% vs 30.5%; P<0.001), insured by Medicaid (23.1% versus 21.3%; P<0.001), have had a primary-care physician (PCP) visit in 2019 (77.4% vs 71.1%; P<0.001), had a higher proportion with a diagnosis of obesity (11.2% vs 9.6%; P<0.001), and had a lower proportion diagnosed with the following: cancer (11.8% vs 19.4%; P<0.001), COPD (10.1% vs 16.2%; P<0.001), heart failure (6.1% vs 10.7%; P<0.001), hypertension (58.0% vs 60.9%; P<0.001), type 2 diabetes (25.6% vs 26.9%; P<0.01), and depression (13.9% vs 15.2%; P<0.001). Propensity-score matching yielded 7,949 BP users and 7,949 BP non-users with no significant differences across examined characteristics. A total of 443 BP non-users from the unmatched “Osteo-Dx-Rx” analysis cohort were not assigned an applicable BP user pair during the matching procedure and were excluded from the matched BP non-user population.

Statin user/non-user analysis

Statin-use comparison: All observations (all regions combined)

A total of 7,906,603 patients were included in the unmatched analysis cohort comparison of statin-use, of which 1,503,395 (19.0%) and 6,403,208 (81.0%) were classified as statin users and statin non-users, respectively (Appendix 2—table 24). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to statin non-users, statin users were older (87.9% age ≥51 vs 37.1%; P<0.001), with a higher proportion of males (41.1% vs 40.9%; P<0.001), from the northeast (29.7% versus 26.6%; P<0.001), with higher comorbidity burden (mean CCI = 1.15 vs 0.49; P<0.001), insured by Medicare (32.7% vs 11.3%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (66.1% vs 44.1%; P<0.001). Propensity-score matching yielded 1,436,300 statin users and 1,436,300 statin non-users with no significant differences across age group, region, insurance type, and having had any PCP visit in 2019. The final matched population did, however, display statistically significant differences between statin users and statin non-users for gender (58.7% vs 58.4% male; P<0.001) and mean CCI (1.11 vs 1.12; P<0.001). These differences, however, are small in magnitude, and were statistically significant due to the underlying statistical power associated with the large sample size. A total of 67,095 statin users from the unmatched analysis cohort were not assigned an applicable statin non-user pair during the matching procedure and were excluded from the matched statin user population.

Statin-use comparison: New York State

A total of 968,296 patients identified as residing in New York state were included in the unmatched analysis cohort comparison of statin-use, of which 206,301 (21.3%) and 761,995 (78.7%) were classified as statin users and statin non-users, respectively (Appendix 2—table 25). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to statin non-users, statin users were older (90.3% age ≥51 vs 43.1%; P<0.001), with a higher proportion of males (42.0% vs 40.4%; P<0.001), with higher comorbidity burden (mean CCI = 0.94=0.94=0.94 1.17 vs 0.51; P<0.001), insured by Medicare (47.4% versus 14.5%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (64.0% vs 41.3%; P<0.001). Propensity-score matching yielded 185,536 statin users and 185,536 statin non-users with no significant differences across age group, gender, insurance type, and having had any PCP visit in 2019. The final matched population did, however, display statistically significant differences between statin users and statin non-users for mean CCI (1.06 vs 1.08; P<0.001). This difference, however, is small in magnitude, and was statistically significant due to the underlying statistical power associated with the large sample size. A total of 20,765 statin users from the unmatched analysis cohort were not assigned an applicable statin non-user pair during the matching procedure and were excluded from the matched statin user population.

BP-use comparison within statin users: All regions combined

Of the 1,436,300 statin users from the statin user/non-user propensity-score matching analysis, a total of 217,981 (15.2%) and 1,218,319 (84.8%) were classified as BP users and BP non-users, respectively (Appendix 2—table 26). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics except for any PCP visit in 2019 (P=0.27). Compared to BP non-users, BP users were older (98.9% age ≥51 vs 85.3%; P<0.001), with a higher proportion of females (90.1% vs 53.1%; P<0.001), from the west (21.7% vs 14.0%; P<0.001), with lower comorbidity burden (mean CCI = 0.95 vs 1.13; P<0.001), and insured by Medicare (50.8% vs 29.7%; P<0.001). Propensity-score matching yielded 213,480 BP users and 213,480 BP non-users with no significant differences across examined characteristics. A total of 4,501 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within statin users: New York State

Of the 185,536 statin users from the statin user/non-user propensity-score matching analysis on patients residing in New York state, a total of 23,863 (12.9%) and 161,673 (87.1%) were classified as BP users and BP non-users, respectively (Appendix 2—table 27). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics except for any PCP visit in 2019 (P=0.33). Compared to BP non-users, BP users were older (99.3% age ≥51 vs 87.7%; P<0.001), with a higher proportion of females (91.2% vs 53.3%; P<0.001), with lower comorbidity burden (mean CCI = 0.92 versus 1.08; P<0.001), and insured by Medicare (66.4% vs 41.9%; P<0.001). Propensity-score matching yielded 23,736 BP users and 23,736 BP non-users with no significant differences across examined characteristics. A total of 127 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within statin non-users: All regions combined

Of the 1,436,300 statin non-users from the statin user/non-user propensity-score matching analysis, a total of 124,843 (8.7%) and 1,311,457 (91.3%) were classified as BP users and BP non-users, respectively (Appendix 2—table 28). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.7% age ≥51 vs 86.3%; P<0.001), with a higher proportion of females (89.6% vs 55.5%; P<0.001), from the west (21.4% vs 14.6%; P<0.001), with lower comorbidity burden (mean CCI = 1.02 versus 1.13; P<0.001), insured by Medicare (45.8% vs 31.7%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (71.7% vs 63.9%; P<0.001). Propensity-score matching yielded 124,716 BP users and 124,716 BP non-users with no significant differences across examined characteristics. A total of 127 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within statin non-users: New York State

Of the 185,536 statin non-users from the statin user/non-user propensity-score matching analysis on patients residing in New York state, a total of 14,546 (7.8%) and 170,990 (92.2%) were classified as BP users and BP non-users, respectively (Appendix 2—table 29). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (99.2% age ≥51 vs 88.4%; P<0.001), with a higher proportion of females (90.6% vs 55.1%; P<0.001), with lower comorbidity burden (mean CCI = 0.95 vs 1.09; P<0.001), insured by Medicare (59.7% versus 43.7%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (70.5% vs 59.4%; P<0.001). Propensity-score matching yielded 14,521 BP users and 14,521 BP non-users with no significant differences across examined characteristics. A total of 25 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Antihypertensive user/non-user analysis

Antihypertensive-use comparison: All observations (all regions combined)

A total of 7,906,603 patients were included in the unmatched analysis cohort comparison of antihypertensive-use, of which 2,101,120 (26.6%) and 5,805,483 (73.4%) were classified as antihypertensive users and antihypertensive non-users, respectively (Appendix 2—table 30). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antihypertensive non-users, antihypertensive users were older (80.8% age ≥51 vs 34.4%; P<0.001), with a higher proportion of females (60.4% vs 58.6%; P<0.001), from the northeast (27.8% vs 27.0%; P<0.001), with higher comorbidity burden (mean CCI = 1.13 vs 0.43; P<0.001), insured by Medicare (29.5% vs 10.3%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (64.2% vs 39.2%; P<0.001). Propensity-score matching yielded 1,786,001 antihypertensive users and 1,786,001 antihypertensive non-users with no significant differences across age group, gender, region, insurance type, and having had any PCP visit in 2019. The final matched population did, however, display statistically significant difference between antihypertensive users and antihypertensive non-users for mean CCI (1.64 vs 1.66; P<0.05). This difference, however, is small in magnitude, and was statistically significant due to the underlying statistical power associated with the large sample size. A total of 315,119 antihypertensive users from the unmatched analysis cohort were not assigned an applicable antihypertensive non-user pair during the matching procedure and were excluded from the matched antihypertensive user population.

Antihypertensive-use comparison: New York State

A total of 968,296 patients identified as residing in New York state were included in the unmatched analysis cohort comparison of antihypertensive-use, of which 258,652 (26.7%) and 709,644 (73.3%) were classified as antihypertensive users and antihypertensive non-users, respectively (Appendix 2—table 31). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antihypertensive non-users, antihypertensive users were older (86.6% age ≥51 vs 40.9%; P<0.001), with a higher proportion of females (59.4% vs 59.2%; P=0.02), with higher comorbidity burden (mean CCI = 1.17 vs 0.46; P<0.001), insured by Medicare (45.9% vs 12.6%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (62.4% vs 40.3%; P<0.001). Propensity-score matching yielded 203,624 antihypertensive users and 203,624 antihypertensive non-users with no significant differences across examined characteristics. A total of 55,028 antihypertensive users from the unmatched analysis cohort were not assigned an applicable antihypertensive non-user pair during the matching procedure and were excluded from the matched antihypertensive user population.

BP-use comparison within antihypertensive users: All regions combined

Of the 1,786,001 antihypertensive users from the antihypertensive user/non-user propensity-score matching analysis, a total of 206,613 (11.6%) and 1,579,388 (88.4%) were classified as BP users and BP non-users, respectively (Appendix 2—table 32). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.2% age ≥51 vs 75.2%; P<0.001), with a higher proportion of females (89.7% vs 56.6%; P<0.001), from the west (22.0% vs 14.3%; P<0.001), with lower comorbidity burden (mean CCI = 0.94 versus 0.95; P=0.02), insured by Medicare (48.6% vs 24.4%; P<0.001), and have not had a primary-care physician (PCP) visit in 2019 (41.2% vs 40.1%; P<0.001). Propensity-score matching yielded 204,396 BP users and 204,396 BP non-users with no significant differences across examined characteristics. A total of 2,217 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antihypertensive users: New York State

Of the 203,624 antihypertensive users from the antihypertensive user/non-user propensity-score matching analysis on patients residing in New York state, a total of 21,213 (10.4%) and 182,411 (89.6%) were classified as BP users and BP non-users, respectively (Appendix 2—table 33). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.8% age ≥51 vs 81.4%; P<0.001), with a higher proportion of females (90.9% vs 55.5%; P<0.001), with lower comorbidity burden (mean CCI = 0.88 vs 0.95; P<0.001), insured by Medicare (64.1% vs 35.9%; P<0.001), and have not had a primary-care physician (PCP) visit in 2019 (53.4% vs 55.7%; P<0.001). Propensity-score matching yielded 21,126 BP users and 21,126 BP non-users with no significant differences across examined characteristics. A total of 87 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antihypertensive non-users: All regions combined

Of the 1,786,001 antihypertensive non-users from the antihypertensive user/non-user propensity-score matching analysis, a total of 136,016 (7.6%) and 1,649,985 (92.4%) were classified as BP users and BP non-users, respectively (Appendix 2—table 34). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (97.7% age ≥51 vs 76.3%; P<0.001), with a higher proportion of females (90.5% vs 58.0%; P<0.001), from the west (20.3% vs 14.8%; P<0.001), with lower comorbidity burden (mean CCI = 0.88 versus 0.96; P<0.001), insured by Medicare (40.7% vs 26.0%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (68.0% vs 59.0%; P<0.001). Propensity-score matching yielded 135,724 BP users and 135,724 BP non-users with no significant differences across examined characteristics. A total of 292 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antihypertensive non-users: New York State

Of the 203,624 antihypertensive non-users from the antihypertensive user/non-user propensity-score matching analysis on patients residing in New York state, a total of 14,051 (6.9%) and 189,573 (93.1%) were classified as BP users and BP non-users, respectively (Appendix 2—table 35). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.7% age ≥51 vs 82.1%; P<0.001), with a higher proportion of females (91.3% vs 56.8%; P<0.001), with lower comorbidity burden (mean CCI = 0.81 vs 0.96; P<0.001), insured by Medicare (54.9% vs 37.7%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (66.3% vs 54.7%; P<0.001). Propensity-score matching yielded 13,983 BP users and 13,983 BP non-users with no significant differences across examined characteristics. A total of 68 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Antidiabetic user/non-user analysis

Antidiabetic-use cComparison: All observations (all regions combined)

A total of 7,906,603 patients were included in the unmatched analysis cohort comparison of antidiabetic-use, of which 755,252 (9.6%) and 7,151,351 (90.4%) were classified as antidiabetic users and antidiabetic non-users, respectively (Appendix 2—table 36). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antidiabetic non-users, antidiabetic users were older (79.4% age ≥51 vs 43.3%; P<0.001), with a higher proportion of females (60.8% vs 58.9%; P<0.001), from the northeast (28.8% vs 27.1%; P<0.001), with higher comorbidity burden (mean CCI = 1.25 vs 0.55; P<0.001), insured by Medicare (26.2% vs 14.2%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (66.5% vs 43.6%; P<0.001). Propensity-score matching yielded 754,553 antidiabetic users and 754,553 antidiabetic non-users with no significant differences across examined characteristics. A total of 699 antidiabetic users from the unmatched analysis cohort were not assigned an applicable antidiabetic non-user pair during the matching procedure and were excluded from the matched antidiabetic user population.

Antidiabetic-use comparison: New York State

A total of 968,296 patients identified as residing in New York state were included in the unmatched analysis cohort comparison of antidiabetic-use, of which 105,117 (10.9%) and 863,179 (89.1%) were classified as antidiabetic users and antidiabetic non-users, respectively (Appendix 2—table 37). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antidiabetic non-users, antidiabetic users were older (83.8% age ≥51 vs 49.4%; P<0.001), with a higher proportion of males (42.2% vs 40.6%; P<0.001), with higher comorbidity burden (mean CCI = 1.34 vs 0.56; P<0.001), insured by Medicare (40.5% vs 19.2%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (64.6% vs 43.9%; P<0.001). Propensity-score matching yielded 104,691 antidiabetic users and 104,691 antidiabetic non-users with no significant differences across examined characteristics. A total of 426 antidiabetic users from the unmatched analysis cohort were not assigned an applicable antidiabetic non-user pair during the matching procedure and were excluded from the matched antidiabetic user population.

BP-use comparison within antidiabetic users: All regions combined

Of the 754,553 antidiabetic users from the antidiabetic user/non-user propensity-score matching analysis, a total of 80,529 (10.7%) and 674,024 (89.3%) were classified as BP users and BP non-users, respectively (Appendix 2—table 38). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.2% age ≥51 vs 75.2%; P<0.001), with a higher proportion of females (98.5% vs 77.1%; P<0.001), from the west (22.2% versus 14.2%; P<0.001), with a higher comorbidity burden (mean CCI = 1.32 versus 1.23; P<0.001), insured by Medicare (45.2% vs 24.0%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (69.5% vs 66.1%; P<0.001). Propensity-score matching yielded 79,500 BP users and 79,500 BP non-users with no significant differences across examined characteristics. A total of 1,029 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidiabetic users: New York State

Of the 104,691 antidiabetic users from the antidiabetic user/non-user propensity-score matching analysis on patients residing in New York state, a total of 9,529 (9.1%) and 95,162 (90.9%) were classified as BP users and BP non-users, respectively (Appendix 2—table 39). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (99.1% age ≥51 vs 82.2%; P<0.001), with a higher proportion of females (90.1% vs 54.5%; P<0.001), with a higher comorbidity burden (mean CCI = 1.46 vs 1.31; P<0.001), insured by Medicare (64.6% vs 38.2%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (66.3% vs 64.4%; P<0.001). Propensity-score matching yielded 9,456 BP users and 9,456 BP non-users with no significant differences across examined characteristics. A total of 73 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidiabetic non-users: All regions combined

Of the 754,553 antidiabetic non-users from the antidiabetic user/non-user propensity-score matching analysis, a total of 73,173 (9.7%) and 681,380 (90.3%) were classified as BP users and BP non-users, respectively (Appendix 2—table 40). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic characteristics, but no difference was seen in mean CCI (1.24 vs 1.24; P=0.92). Compared to BP non-users, BP users were older (98.0% age ≥51 vs 77.3%; P<0.001), with a higher proportion of females (88.9% vs 57.7%; P<0.001), from the west (20.1% vs 14.5%; P<0.001), insured by Medicare (40.0% vs 24.8%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (74.1% vs 65.7%; P<0.001). Propensity-score matching yielded 72,514 BP users and 72,514 BP non-users with no significant differences across examined characteristics. A total of 659 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidiabetic non-users: New York State

Of the 104,691 antidiabetic non-users from the antidiabetic user/non-user propensity-score matching analysis on patients residing in New York state, a total of 9,275 (8.9%) and 95,416 (91.1%) were classified as BP users and BP non-users, respectively (Appendix 2—table 41). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (99.0% age ≥51 vs 82.2%; P<0.001), with a higher proportion of females (89.2% vs 54.7%; P<0.001), with a higher comorbidity burden (mean CCI = 1.37 vs 1.32; P<0.01), insured by Medicare (57.7% vs 38.9%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (72.5% vs 63.8%; P<0.001). Propensity-score matching yielded 13,983 BP users and 13,983 BP non-users with no significant differences across examined characteristics. A total of 131 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Antidepressant user/non-user analysis

Antidepressant-use comparison: All observations (all regions combined)

A total of 7,906,603 patients were included in the unmatched analysis cohort comparison of antidepressant-use, of which 1,571,005 (19.9%) and 6,335,598 (80.1%) were classified as antidepressant users and antidepressant non-users, respectively (Appendix 2—table 42). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antidepressant non-users, antidepressant users were older (58.6% age ≥51 vs 43.8%; P<0.001), with a higher proportion of females (72.8% vs 55.7%; P<0.001), from the midwest (22.1% vs 17.7%; P<0.001), with higher comorbidity burden (mean CCI = 0.90 vs 0.55; P<0.001), insured by Medicare (18.5% vs 14.6%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (61.1% versus 42.0%; P<0.001). Propensity-score matching yielded 1,536,048 antidepressant users and 1,536,048 antidepressant non-users with no significant differences across examined characteristics. A total of 34,957 antidepressant users from the unmatched analysis cohort were not assigned an applicable antidepressant non-user pair during the matching procedure and were excluded from the matched antidepressant user population.

Antidepressant-use comparison: New York State

A total of 968,296 patients identified as residing in New York state were included in the unmatched analysis cohort comparison of antidepressant-use, of which 136,081 (14.1%) and 832,215 (85.9%) were classified as antidepressant users and antidepressant non-users, respectively (Appendix 2—table 43). Prior to propensity-score matching, there were significant differences across all demographic and clinical characteristics. Compared to antidepressant non-users, antidepressant users were older (66.3% age ≥51 vs 51.0%; P<0.001), with a higher proportion of females (71.2% vs 57.3%; P<0.001), with higher comorbidity burden (mean CCI = 0.98 vs 0.59; P<0.001), insured by Medicare (32.2% vs 19.8%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (60.7% vs 43.8%; P<0.001). Propensity-score matching yielded 135,516 antidepressant users and 135,516 antidepressant non-users with no significant differences across examined characteristics. A total of 565 antidepressant users from the unmatched analysis cohort were not assigned an applicable antidepressant non-user pair during the matching procedure and were excluded from the matched antidepressant user population.

BP-use comparison within antidepressant users: All regions combined

Of the 1,536,048 antidepressant users from the antidepressant user/non-user propensity-score matching analysis, a total of 145,109 (9.4%) and 1,390,939 (90.6%) were classified as BP users and BP non-users, respectively (Appendix 2—table 44). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (96.7% age ≥51 vs 54.4%; P<0.001), with a higher proportion of females (91.9% vs 70.2%; P<0.001), from the west (19.6% versus 13.9%; P<0.001), with a higher comorbidity burden (mean CCI = 1.09 versus 0.84; P<0.001), insured by Medicare (42.4% vs 16.2%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (64.6% vs 60.2%; P<0.001). Propensity-score matching yielded 144,282 BP users and 144,282 BP non-users with no significant differences across examined characteristics. A total of 827 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidepressant users: New York State

Of the 135,516 antidepressant users from the antidepressant user/non-user propensity-score matching analysis on patients residing in New York state, a total of 12,950 (9.6%) and 122,566 (90.4%) were classified as BP users and BP non-users, respectively (Appendix 2—table 45). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (97.8% age ≥51 vs 63.0%; P<0.001), with a higher proportion of females (92.6% vs 68.9%; P<0.001), with a higher comorbidity burden (mean CCI = 1.13 vs 0.95; P<0.001), insured by Medicare (60.8% vs 29.1%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (65.3% vs 60.1%; P<0.001). Propensity-score matching yielded 12,859 BP users and 12,859 BP non-users with no significant differences across examined characteristics. A total of 91 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidepressant non-users: All regions combined

Of the 1,536,048 antidepressant non-users from the antidepressant user/non-user propensity-score matching analysis, a total of 113,110 (7.4%) and 1,422,938 (92.6%) were classified as BP users and BP non-users, respectively (Appendix 2—table 46). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic characteristics. Compared to BP non-users, BP users were older (97.1% age ≥51 vs 55.4%; P<0.001), with a higher proportion of females (93.2% vs 70.6%; P<0.001), from the west (20.0% versus 14.0%; P<0.001), with a higher comorbidity burden (mean CCI = 1.06 versus 0.85; P<0.001), insured by Medicare (40.4% vs 17.0%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (71.2% vs 59.8%; P<0.001). Propensity-score matching yielded 112,402 BP users and 112,402 BP non-users with no significant differences across examined characteristics. A total of 708 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

BP-use comparison within antidepressant non-users: New York State

Of the 135,516 antidepressant non-users from the antidepressant user/non-user propensity-score matching analysis on patients residing in New York state, a total of 10,174 (7.5%) and 125,342 (92.5%) were classified as BP users and BP non-users, respectively (Appendix 2—table 47). Prior to propensity-score matching based on BP-use, there were significant differences across all demographic and clinical characteristics. Compared to BP non-users, BP users were older (98.4% age ≥51 vs 63.7%; P<0.001), with a higher proportion of females (93.6% vs 69.4%; P<0.001), with a higher comorbidity burden (mean CCI = 1.13 vs 0.95; P<0.01), insured by Medicare (60.0% vs 29.9%; P<0.001), and have had a primary-care physician (PCP) visit in 2019 (71.7% vs 59.7%; P<0.001). Propensity-score matching yielded 10,091 BP users and 10,091 BP non-users with no significant differences across examined characteristics. A total of 83 BP users were not assigned an applicable BP non-user pair during the matching procedure and were excluded from the matched BP user population.

Appendix 2—table 1
Primary Analysis Cohort (Region=Northeast), Patient Characteristics Pre/Post Match.
Region=Northeast UnmatchedRegion=Northeast Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients2,152,560100.00%2,032,83294.40%119,7285.60%238,988100.00%119,49450.00%119,49450.00%
Age
 ≤20363,63716.90%363,40117.90%2360.20%<0.0014740.20%2380.20%2360.20%1
 21-40397,37718.50%396,61319.50%7640.60%1,5280.60%7640.60%7640.60%
 41-50261,57012.20%259,52812.80%2,0421.70%4,0841.70%2,0421.70%2,0421.70%
 51-60372,23817.30%354,22817.40%18,01015.00%36,02015.10%18,01015.10%18,01015.10%
 61-70354,33116.50%313,23715.40%41,09434.30%82,23334.40%41,13934.40%41,09434.40%
 71-80252,71211.70%215,15110.60%37,56131.40%74,83131.30%37,39331.30%37,43831.30%
 ≥81150,6957.00%130,6746.40%20,02116.70%39,81816.70%19,90816.70%19,91016.70%
Gender
 Female1,275,61159.30%1,167,24157.40%108,37090.50%<0.001216,27390.50%108,13790.50%108,13690.50%0.99
 Male876,94940.70%865,59142.60%11,3589.50%22,7159.50%11,3579.50%11,3589.50%
Insurance
 Commercial1,050,79548.80%1,017,50250.10%33,29327.80%<0.00166,55227.80%33,25927.80%33,29327.90%0.99
 Dual47,7732.20%40,1682.00%7,6056.40%15,1146.30%7,5766.30%7,5386.30%
 Medicaid631,86329.40%608,64929.90%23,21419.40%46,09419.30%23,04719.30%23,04719.30%
 Medicare422,12919.60%366,51318.00%55,61646.50%111,22846.50%55,61246.50%55,61646.50%
PCP Visit 2019
 No1,212,39456.30%1,162,52757.20%49,86741.70%<0.00199,74141.70%49,87441.70%49,86741.70%0.98
 Yes940,16643.70%870,30542.80%69,86158.30%139,24758.30%69,62058.30%69,62758.30%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.671.420.651.40.931.71<0.0010.931.710.931.710.931.710.96
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 2
Primary Analysis Cohort (Region=Midwest), Patient Characteristics Pre/Post Match.
Region=Midwest UnmatchedRegion=Midwest Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients1,467,802100.0%1,391,83594.8%75,9675.2%151,802100.0%75,90150.0%75,90150.0%
Age
 ≤20310,02721.1%309,75922.3%2680.4%<0.0015370.4%2690.4%2680.4%1.00
 21-40287,23619.6%286,64320.6%5930.8%1,1880.8%5950.8%5930.8%
 41-50185,24012.6%183,55613.2%1,6842.2%3,3672.2%1,6832.2%1,6842.2%
 51-60246,23016.8%233,99216.8%12,23816.1%24,47816.1%12,24016.1%12,23816.1%
 61-70224,66815.3%196,17214.1%28,49637.5%56,99137.5%28,49537.5%28,49637.5%
 71-80130,5638.9%109,4427.9%21,12127.8%42,15327.8%21,07527.8%21,07827.8%
 ≥8183,8385.7%72,2715.2%11,56715.2%23,08815.2%11,54415.2%11,54415.2%
Gender
 Female863,15658.8%794,57857.1%68,57890.3%<0.001137,02890.3%68,51690.3%68,51290.3%0.97
 Male604,64641.2%597,25742.9%7,3899.7%14,7749.7%7,3859.7%7,3899.7%
Insurance
 Commercial885,65160.3%854,51861.4%31,13341.0%<0.00162,24341.0%31,11041.0%31,13341.0%1.00
 Dual28,1901.9%24,5841.8%3,6064.7%7,2114.8%3,6054.7%3,6064.8%
 Medicaid318,59621.7%310,47322.3%8,12310.7%16,13610.6%8,07910.6%8,05710.6%
 Medicare235,36516.0%202,26014.5%33,10543.6%66,21243.6%33,10743.6%33,10543.6%
PCP Visit 2019
 No711,30848.5%682,60149.0%28,70737.8%<0.00157,39837.8%28,69137.8%28,70737.8%0.93
 Yes756,49451.5%709,23451.0%47,26062.2%94,40462.2%47,21062.2%47,19462.2%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.591.370.561.340.991.86<0.0010.991.860.991.851.001.860.77
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 3
Primary Analysis Cohort (Region=South), Patient Characteristics Pre/Post Match.
Region=South UnmatchedRegion=South Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients3,042,604100.0%2,881,71894.7%160,8865.3%319,408100.0%159,70450.0%159,70450.0%
Age
 ≤20890,67729.3%890,20330.9%4740.3%<0.0019430.3%4690.3%4740.3%1.00
 21-40527,97117.4%526,79418.3%1,1770.7%2,3640.7%1,1870.7%1,1770.7%
 41-50338,26211.1%334,84111.6%3,4212.1%6,8392.1%3,4182.1%3,4212.1%
 51-60442,75714.6%417,66414.5%25,09315.6%50,18615.7%25,09315.7%25,09315.7%
 61-70409,85413.5%353,95812.3%55,89634.7%111,80035.0%55,90435.0%55,89635.0%
 71-80272,7619.0%222,1567.7%50,60531.5%99,22331.1%49,60531.1%49,61831.1%
 ≥81160,3225.3%136,1024.7%24,22015.1%48,05315.0%24,02815.0%24,02515.0%
Gender
 Female1,800,16659.2%1,654,35157.4%145,81590.6%<0.001289,26390.6%144,63090.6%144,63390.6%0.99
 Male1,242,43840.8%1,227,36742.6%15,0719.4%30,1459.4%15,0749.4%15,0719.4%
Insurance
 Commercial1,475,45648.5%1,416,16649.1%59,29036.9%<0.001118,58737.1%59,29737.1%59,29037.1%1.00
 Dual53,4741.8%39,4141.4%14,0608.7%25,7528.1%12,8748.1%12,8788.1%
 Medicaid1,121,60636.9%1,099,95738.2%21,64913.5%43,29913.6%21,65013.6%21,64913.6%
 Medicare392,06812.9%326,18111.3%65,88741.0%131,77041.3%65,88341.3%65,88741.3%
PCP Visit 2019
 No1,701,04055.9%1,646,57257.1%54,46833.9%<0.001108,60134.0%54,27534.0%54,32634.0%0.85
 Yes1,341,56444.1%1,235,14642.9%106,41866.1%210,80766.0%105,42966.0%105,37866.0%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.571.310.551.280.861.70<0.0010.861.700.861.700.861.710.84
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 4
Primary Analysis Cohort (Region=West), Patient Characteristics Pre/Post Match.
Region=West UnmatchedRegion=West Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients1,243,637100.0%1,148,16792.3%95,4707.7%190,534100.0%95,26750.0%95,26750.0%
Age
 ≤20275,70922.2%275,55924.0%1500.2%<0.0012990.2%1490.2%1500.2%1.00
 21-40234,41518.8%233,85820.4%5570.6%1,1150.6%5580.6%5570.6%
 41-50140,23711.3%138,83312.1%1,4041.5%2,8061.5%1,4021.5%1,4041.5%
 51-60188,96515.2%178,58515.6%10,38010.9%20,76110.9%10,38110.9%10,38010.9%
 61-70192,40815.5%161,01614.0%31,39232.9%62,79833.0%31,40633.0%31,39233.0%
 71-80127,73910.3%95,3018.3%32,43834.0%64,59633.9%32,29333.9%32,30333.9%
 ≥8184,1646.8%65,0155.7%19,14920.1%38,15920.0%19,07820.0%19,08120.0%
Gender
 Female732,02758.9%647,35456.4%84,67388.7%<0.001168,93388.7%84,46388.7%84,47088.7%0.96
 Male511,61041.1%500,81343.6%10,79711.3%21,60111.3%10,80411.3%10,79711.3%
Insurance
 Commercial526,70142.4%503,35943.8%23,34224.4%<0.00146,68824.5%23,34624.5%23,34224.5%1.00
 Dual27,0602.2%20,9241.8%6,1366.4%11,8596.2%5,9256.2%5,9346.2%
 Medicaid522,43542.0%497,94143.4%24,49425.7%48,99025.7%24,49625.7%24,49425.7%
 Medicare167,44113.5%125,94311.0%41,49843.5%82,99743.6%41,50043.6%41,49743.6%
PCP Visit 2019
 No658,95553.0%628,13154.7%30,82432.3%<0.00161,64332.4%30,81932.4%30,82432.4%0.98
 Yes584,68247.0%520,03645.3%64,64667.7%128,89167.6%64,44867.6%64,44367.6%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.691.460.661.421.081.84<0.0011.091.831.081.831.091.840.73
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 5
Primary Analysis Cohort (Region=New York State), Patient Characteristics Pre/Post Match.
Region=New York State UnmatchedRegion=New York State Matched
AllBP Non-usersBP Usersp-valueAllBP Non-usersBP Usersp-value
N%N%N%N%N%N%
All Patients968,296100.0%918,26194.8%50,0355.2%99,724100.0%49,86250.0%49,86250.0%
Age
≤20133,17813.8%133,12814.5%500.1%<0.0011020.1%520.1%500.1%1.00
21-40192,95919.9%192,73121.0%2280.5%4530.5%2250.5%2280.5%
41-50127,79413.2%127,13913.8%6551.3%1,3111.3%6561.3%6551.3%
51-60172,44417.8%166,08018.1%6,36412.7%12,73212.8%6,36812.8%6,36412.8%
61-70159,91216.5%143,77615.7%16,13632.2%32,26532.4%16,12932.3%16,13632.4%
71-80120,11712.4%102,65511.2%17,46234.9%34,69334.8%17,35234.8%17,34134.8%
≥8161,8926.4%52,7525.7%9,14018.3%18,16818.2%9,08018.2%9,08818.2%
Gender
Female573,61059.2%528,15257.5%45,45890.9%<0.00190,56790.8%45,28290.8%45,28590.8%0.97
Male394,68640.8%390,10942.5%4,5779.1%9,1579.2%4,5809.2%4,5779.2%
Insurance
Commercial500,91851.7%490,50353.4%10,41520.8%<0.00120,83020.9%10,41520.9%10,41520.9%1.00
Dual6,8140.7%5,2180.6%1,5963.2%3,1543.2%1,5813.2%1,5733.2%
Medicaid252,36626.1%243,19126.5%9,17518.3%18,04418.1%9,01918.1%9,02518.1%
Medicare208,19821.5%179,34919.5%28,84957.7%57,69657.9%28,84757.9%28,84957.9%
PCP Visit 2019
No521,28253.8%502,60954.7%18,67337.3%<0.00137,25337.4%18,61637.3%18,63737.4%0.89
Yes447,01446.2%415,65245.3%31,36262.7%62,47162.6%31,24662.7%31,22562.6%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI0.651.390.631.370.951.68<0.0010.951.680.951.670.951.680.93
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 6
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, All-Regions Combined.
Primary Analysis Cohort, All Regions Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients900,732100.0%28,1373.1%16,2891.8%3,7100.4%
BP User450,36650.0%5,1891.2%0.22<0.0013,0240.7%0.22<0.0017150.2%0.24<0.001
BP Non-user450,36650.0%22,9485.1%0.210.2213,2652.9%0.210.232,9950.7%0.220.26
By Age
Age ≤202,253100.0%673.0%140.6%20.1%
BP User1,12850.1%292.6%0.750.2620.2%0.160.00720.2%NANA
BP Non-user1,12549.9%383.4%0.461.23121.1%0.040.7400.0%NANA
Age 21-406,195100.0%3355.4%1151.9%130.2%
BP User3,09149.9%581.9%0.20<0.001150.5%0.15<0.00140.1%0.450.27
BP Non-user3,10450.1%2778.9%0.150.261003.2%0.080.2590.3%0.141.45
Age 41-5017,096100.0%8945.2%2701.6%540.3%
BP User8,55150.0%1882.2%0.25<0.001480.6%0.21<0.001140.2%0.35<0.001
BP Non-user8,54550.0%7068.3%0.210.292222.6%0.150.29400.5%0.190.64
Age 51-60131,445100.0%5,7654.4%2,3711.8%3970.3%
BP User65,72150.0%1,1041.7%0.22<0.0014560.7%0.23<0.001830.1%0.26<0.001
BP Non-user65,72450.0%4,6617.1%0.210.241,9152.9%0.210.263140.5%0.210.34
Age 61-70313,822100.0%10,4383.3%5,0291.6%1,0350.3%
BP User156,87850.0%1,8431.2%0.21<0.0019390.6%0.23<0.0011730.1%0.20<0.001
BP Non-user156,94450.0%8,5955.5%0.200.224,0902.6%0.210.248620.5%0.170.24
Age 71-80280,803100.0%7,1792.6%4,8271.7%1,2120.4%
BP User140,43750.0%1,3090.9%0.22<0.0018770.6%0.22<0.0012340.2%0.24<0.001
BP Non-user140,36650.0%5,8704.2%0.200.233,9502.8%0.200.239780.7%0.210.27
Age ≥81149,118100.0%3,4592.3%3,6632.5%9970.7%
BP User74,56050.0%6580.9%0.23<0.0016870.9%0.22<0.0012050.3%0.26<0.001
BP Non-user74,55850.0%2,8013.8%0.210.252,9764.0%0.210.247921.1%0.220.30
Female Patients811,497100.0%24,9363.1%14,3671.8%3,1270.4%
BP User405,75150.0%4,5191.1%0.21<0.0012,6670.7%0.22<0.0015930.1%0.23<0.001
BP Non-user405,74650.0%20,4175.0%0.210.2211,7002.9%0.210.232,5340.6%0.210.25
By Age
Age ≤20885100.0%262.9%70.8%10.1%
BP User44249.9%112.5%0.730.4310.2%0.170.1210.2%NANA
BP Non-user44350.1%153.4%0.331.6061.4%0.021.3800.0%NANA
Age 21-403,765100.0%2185.8%641.7%90.2%
BP User1,87949.9%402.1%0.21<0.001120.6%0.23<0.00130.2%0.500.51
BP Non-user1,88650.1%1789.4%0.150.30522.8%0.120.4360.3%0.132.01
Age 41-5013,542100.0%7305.4%2061.5%370.3%
BP User6,77450.0%1572.3%0.26<0.001430.6%0.26<0.001110.2%0.420.01
BP Non-user6,76850.0%5738.5%0.210.311632.4%0.180.36260.4%0.210.85
Age 51-60119,205100.0%5,2004.4%2,0931.8%3270.3%
BP User59,60250.0%9731.6%0.22<0.0013990.7%0.23<0.001640.1%0.24<0.001
BP Non-user59,60350.0%4,2277.1%0.200.231,6942.8%0.210.262630.4%0.180.32
Age 61-70290,276100.0%9,4743.3%4,5061.6%8850.3%
BP User145,13150.0%1,6391.1%0.20<0.0018510.6%0.23<0.0011440.1%0.19<0.001
BP Non-user145,14550.0%7,8355.4%0.190.213,6552.5%0.210.257410.5%0.160.23
Age 71-80253,094100.0%6,3042.5%4,2541.7%1,0260.4%
BP User126,55950.0%1,1400.9%0.21<0.0017690.6%0.22<0.0011930.2%0.23<0.001
BP Non-user126,53550.0%5,1644.1%0.200.233,4852.8%0.200.238330.7%0.200.27
Age ≥81130,730100.0%2,9842.3%3,2372.5%8420.6%
BP User65,36450.0%5590.9%0.22<0.0015920.9%0.22<0.0011770.3%0.26<0.001
BP Non-user65,36650.0%2,4253.7%0.200.252,6454.0%0.200.246651.0%0.220.31
Male Patients89,235100.0%3,2013.6%1,9222.2%5830.7%
BP User44,61550.0%6701.5%0.25<0.0013570.8%0.22<0.0011220.3%0.26<0.001
BP Non-user44,62050.0%2,5315.7%0.230.281,5653.5%0.200.254611.0%0.220.32
By Age
Age ≤201,368100.0%413.0%70.5%10.1%
BP User68650.1%182.6%0.770.4210.1%0.160.0710.1%NANA
BP Non-user68249.9%233.4%0.411.4460.9%0.021.3700.0%NANA
Age 21-402,430100.0%1174.8%512.1%40.2%
BP User1,21249.9%181.5%0.17<0.00130.2%0.06<0.00110.1%0.330.63
BP Non-user1,21850.1%998.1%0.100.28483.9%0.020.1930.2%0.033.22
Age 41-503,554100.0%1644.6%641.8%170.5%
BP User1,77750.0%311.7%0.22<0.00150.3%0.08<0.00130.2%0.210.01
BP Non-user1,77750.0%1337.5%0.150.33593.3%0.030.21140.8%0.060.74
Age 51-6012,240100.0%5654.6%2782.3%700.6%
BP User6,11950.0%1312.1%0.29<0.001570.9%0.25<0.001190.3%0.37<0.001
BP Non-user6,12150.0%4347.1%0.240.352213.6%0.190.34510.8%0.220.63
Age 61-7023,546100.0%9644.1%5232.2%1500.6%
BP User11,74749.9%2041.7%0.26<0.001880.7%0.20<0.001290.2%0.24<0.001
BP Non-user11,79950.1%7606.4%0.220.304353.7%0.160.251211.0%0.160.36
Age 71-8027,709100.0%8753.2%5732.1%1860.7%
BP User13,87850.1%1691.2%0.23<0.0011080.8%0.23<0.001410.3%0.28<0.001
BP Non-user13,83149.9%7065.1%0.190.274653.4%0.180.281451.0%0.200.40
Age ≥8118,388100.0%4752.6%4262.3%1550.8%
BP User9,19650.0%991.1%0.26<0.001951.0%0.28<0.001280.3%0.22<0.001
BP Non-user9,19250.0%3764.1%0.200.323313.6%0.220.351271.4%0.140.33
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 7
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=Northeast.
Region=Northeast Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients238,988100.0%8,8313.7%7,8203.3%1,5050.6%
BP User119,49450.0%1,6841.4%0.22<0.0011,5781.3%0.24<0.0013140.3%0.26<0.001
BP Non-user119,49450.0%7,1476.0%0.210.246,2425.2%0.230.261,1911.0%0.230.30
By Age
Age ≤20474100.0%143.0%71.5%20.4%
BP User23649.8%73.0%1.010.9920.8%0.400.4520.8%NANA
BP Non-user23850.2%72.9%0.352.9252.1%0.082.0700.0%NANA
Age 21-401,528100.0%936.1%553.6%50.3%
BP User76450.0%141.8%0.16<0.00170.9%0.14<0.00110.1%0.250.37
BP Non-user76450.0%7910.3%0.090.29486.3%0.060.3140.5%0.032.23
Age 41-504,084100.0%2345.7%1182.9%180.4%
BP User2,04250.0%532.6%0.27<0.001170.8%0.16<0.00160.3%0.500.16
BP Non-user2,04250.0%1818.9%0.200.371014.9%0.100.27120.6%0.191.33
Age 51-6036,020100.0%1,8635.2%1,1903.3%1600.4%
BP User18,01050.0%3532.0%0.22<0.0012371.3%0.24<0.001380.2%0.31<0.001
BP Non-user18,01050.0%1,5108.4%0.190.259535.3%0.210.281220.7%0.220.45
Age 61-7082,233100.0%3,2003.9%2,4242.9%4030.5%
BP User41,09450.0%5971.5%0.22<0.0015071.2%0.26<0.001790.2%0.24<0.001
BP Non-user41,13950.0%2,6036.3%0.200.241,9174.7%0.230.283240.8%0.190.31
Age 71-8074,831100.0%2,2663.0%2,3063.1%4930.7%
BP User37,43850.0%4421.2%0.23<0.0014751.3%0.25<0.001990.3%0.25<0.001
BP Non-user37,39350.0%1,8244.9%0.210.261,8314.9%0.230.283941.1%0.200.31
Age ≥8139,818100.0%1,1612.9%1,7204.3%4241.1%
BP User19,91050.0%2181.1%0.22<0.0013331.7%0.23<0.001890.4%0.26<0.001
BP Non-user19,90850.0%9434.7%0.190.261,3877.0%0.200.263351.7%0.210.33
Female Patients216,273100.0%7,8973.7%6,9413.2%1,2630.6%
BP User108,13650.0%1,4831.4%0.22<0.0011,3921.3%0.24<0.0012550.2%0.25<0.001
BP Non-user108,13750.0%6,4145.9%0.210.235,5495.1%0.230.261,0080.9%0.220.29
By Age
Age ≤20180100.0%42.2%31.7%10.6%
BP User9050.0%22.2%1.001.0011.1%0.491.0011.1%NANA
BP Non-user9050.0%22.2%0.147.2622.2%0.045.5500.0%NANA
Age 21-40864100.0%596.8%323.7%40.5%
BP User43149.9%102.3%0.19<0.00161.4%0.22<0.00110.2%0.330.62
BP Non-user43350.1%4911.3%0.090.37266.0%0.090.5430.7%0.033.22
Age 41-503,176100.0%1765.5%872.7%130.4%
BP User1,58850.0%402.5%0.28<0.001150.9%0.20<0.00150.3%0.620.40
BP Non-user1,58850.0%1368.6%0.190.40724.5%0.110.3580.5%0.201.91
Age 51-6032,612100.0%1,6905.2%1,0483.2%1250.4%
BP User16,30650.0%3101.9%0.21<0.0012061.3%0.24<0.001310.2%0.33<0.001
BP Non-user16,30650.0%1,3808.5%0.180.248425.2%0.200.27940.6%0.220.49
Age 61-7076,403100.0%2,9333.8%2,1812.9%3430.4%
BP User38,20050.0%5361.4%0.21<0.0014561.2%0.26<0.001630.2%0.22<0.001
BP Non-user38,20350.0%2,3976.3%0.190.231,7254.5%0.230.282800.7%0.170.29
Age 71-8067,857100.0%2,0213.0%2,0633.0%4160.6%
BP User33,93050.0%3931.2%0.23<0.0014131.2%0.24<0.001770.2%0.23<0.001
BP Non-user33,92750.0%1,6284.8%0.210.261,6504.9%0.220.273391.0%0.180.29
Age ≥8135,181100.0%1,0142.9%1,5274.3%3611.0%
BP User17,59150.0%1921.1%0.23<0.0012951.7%0.23<0.001770.4%0.27<0.001
BP Non-user17,59050.0%8224.7%0.190.261,2327.0%0.200.262841.6%0.210.34
Male Patients22,715100.0%9344.1%8793.9%2421.1%
BP User11,35850.0%2011.8%0.26<0.0011861.6%0.26<0.001590.5%0.32<0.001
BP Non-user11,35750.0%7336.5%0.220.316936.1%0.220.301831.6%0.240.43
By Age
Age ≤20294100.0%103.4%41.4%10.3%
BP User14649.7%53.4%1.010.9810.7%0.330.6210.7%NANA
BP Non-user14850.3%53.4%0.293.5832.0%0.033.2400.0%NANA
Age 21-40664100.0%345.1%233.5%10.2%
BP User33350.2%41.2%0.12<0.00110.3%0.04<0.00100.0%NANA
BP Non-user33149.8%309.1%0.040.35226.6%0.010.3210.3%NANA
Age 41-50908100.0%586.4%313.4%50.6%
BP User45450.0%132.9%0.27<0.00120.4%0.06<0.00110.2%0.250.37
BP Non-user45450.0%459.9%0.140.50296.4%0.020.2740.9%0.032.23
Age 51-603,408100.0%1735.1%1424.2%351.0%
BP User1,70450.0%432.5%0.31<0.001311.8%0.27<0.00170.4%0.25<0.001
BP Non-user1,70450.0%1307.6%0.220.451116.5%0.180.40281.6%0.110.57
Age 61-705,830100.0%2674.6%2434.2%601.0%
BP User2,89449.6%612.1%0.29<0.001511.8%0.26<0.001160.6%0.37<0.001
BP Non-user2,93650.4%2067.0%0.210.381926.5%0.190.35441.5%0.210.65
Age 71-806,974100.0%2453.5%2433.5%771.1%
BP User3,50850.3%491.4%0.24<0.001621.8%0.33<0.001220.6%0.39<0.001
BP Non-user3,46649.7%1965.7%0.170.321815.2%0.240.44551.6%0.240.64
Age ≥814,637100.0%1473.2%1934.2%631.4%
BP User2,31950.0%261.1%0.21<0.001381.6%0.23<0.001120.5%0.23<0.001
BP Non-user2,31850.0%1215.2%0.130.321556.7%0.160.33512.2%0.120.43
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 8
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=Midwest.
Region=Midwest Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients151,802100.0%4,4512.9%2,0991.4%6360.4%
BP User75,90150.0%8681.1%0.23<0.0013830.5%0.22<0.0011210.2%0.23<0.001
BP Non-user75,90150.0%3,5834.7%0.220.251,7162.3%0.200.255150.7%0.190.29
By Age
Age ≤20537100.0%152.8%20.4%00.0%
BP User26849.9%62.2%0.660.4400.0%NANA00.0%NANA
BP Non-user26950.1%93.3%0.231.8920.7%NANA00.0%NANA
Age 21-401,188100.0%625.2%171.4%10.1%
BP User59349.9%71.2%0.12<0.00120.3%0.130.00200.0%NANA
BP Non-user59550.1%559.2%0.050.26152.5%0.030.5710.2%NANA
Age 41-503,367100.0%1845.5%461.4%160.5%
BP User1,68450.0%362.1%0.23<0.001100.6%0.27<0.00120.1%0.140.002
BP Non-user1,68350.0%1488.8%0.160.33362.1%0.140.55140.8%0.030.62
Age 51-6024,478100.0%9513.9%2931.2%800.3%
BP User12,23850.0%1801.5%0.22<0.001520.4%0.21<0.001150.1%0.23<0.001
BP Non-user12,24050.0%7716.3%0.190.262412.0%0.160.29650.5%0.130.40
Age 61-7056,991100.0%1,7643.1%6711.2%1890.3%
BP User28,49650.0%3221.1%0.21<0.0011230.4%0.22<0.001350.1%0.23<0.001
BP Non-user28,49550.0%1,4425.1%0.190.245481.9%0.180.271540.5%0.160.33
Age 71-8042,153100.0%1,0092.4%5771.4%2000.5%
BP User21,07850.0%2091.0%0.25<0.001950.5%0.19<0.001370.2%0.23<0.001
BP Non-user21,07550.0%8003.8%0.220.304822.3%0.160.241630.8%0.160.32
Age ≥8123,088100.0%4662.0%4932.1%1500.6%
BP User11,54450.0%1080.9%0.30<0.0011010.9%0.25<0.001320.3%0.27<0.001
BP Non-user11,54450.0%3583.1%0.240.373923.4%0.200.311181.0%0.180.40
Female Patients137,028100.0%3,9452.9%1,8281.3%5430.4%
BP User68,51250.0%7621.1%0.23<0.0013330.5%0.22<0.0011030.2%0.23<0.001
BP Non-user68,51650.0%3,1834.6%0.210.251,4952.2%0.190.254400.6%0.190.29
By Age
Age ≤20226100.0%73.1%10.4%00.0%
BP User11350.0%32.7%0.741.0000.0%NANA00.0%NANA
BP Non-user11350.0%43.5%0.163.4010.9%NANA00.0%NANA
Age 21-40700100.0%344.9%71.0%00.0%
BP User34949.9%61.7%0.20<0.00110.3%0.170.1200.0%NANA
BP Non-user35150.1%288.0%0.080.4961.7%0.021.3800.0%NANA
Age 41-502,639100.0%1575.9%321.2%100.4%
BP User1,31950.0%312.4%0.23<0.00180.6%0.330.00510.1%0.110.02
BP Non-user1,32050.0%1269.5%0.150.34241.8%0.150.7490.7%0.010.87
Age 51-6022,101100.0%8563.9%2601.2%700.3%
BP User11,05050.0%1591.4%0.22<0.001470.4%0.22<0.001130.1%0.23<0.001
BP Non-user11,05150.0%6976.3%0.180.262131.9%0.160.30570.5%0.120.42
Age 61-7052,520100.0%1,5943.0%5911.1%1650.3%
BP User26,26050.0%2861.1%0.21<0.0011070.4%0.22<0.001290.1%0.21<0.001
BP Non-user26,26050.0%1,3085.0%0.180.244841.8%0.180.271360.5%0.140.32
Age 71-8038,367100.0%8772.3%5011.3%1720.4%
BP User19,18450.0%1800.9%0.25<0.001850.4%0.20<0.001330.2%0.24<0.001
BP Non-user19,18350.0%6973.6%0.210.304162.2%0.160.251390.7%0.160.35
Age ≥8120,475100.0%4202.1%4362.1%1260.6%
BP User10,23750.0%970.9%0.29<0.001850.8%0.24<0.001270.3%0.27<0.001
BP Non-user10,23850.0%3233.2%0.230.373513.4%0.190.30991.0%0.180.41
Male Patients14,774100.0%5063.4%2711.8%930.6%
BP User7,38950.0%1061.4%0.25<0.001500.7%0.22<0.001180.2%0.24<0.001
BP Non-user7,38550.0%4005.4%0.200.322213.0%0.160.30751.0%0.140.40
By Age
Age ≤20311100.0%82.6%10.3%00.0%
BP User15549.8%31.9%0.600.7200.0%NANA00.0%NANA
BP Non-user15650.2%53.2%0.142.5410.6%NANA00.0%NANA
Age 21-40488100.0%285.7%102.0%10.2%
BP User24450.0%10.4%0.03<0.00110.4%0.110.0200.0%NANA
BP Non-user24450.0%2711.1%0.000.2593.7%0.010.8510.4%NANA
Age 41-50728100.0%273.7%141.9%60.8%
BP User36550.1%51.4%0.22<0.00120.5%0.160.00710.3%0.200.12
BP Non-user36349.9%226.1%0.080.57123.3%0.040.7351.4%0.021.69
Age 51-602,377100.0%954.0%331.4%100.4%
BP User1,18850.0%211.8%0.27<0.00150.4%0.18<0.00120.2%0.250.11
BP Non-user1,18950.0%746.2%0.170.44282.4%0.070.4680.7%0.051.17
Age 61-704,471100.0%1703.8%801.8%240.5%
BP User2,23650.0%361.6%0.26<0.001160.7%0.24<0.00160.3%0.330.01
BP Non-user2,23550.0%1346.0%0.180.37642.9%0.140.42180.8%0.130.84
Age 71-803,786100.0%1323.5%762.0%280.7%
BP User1,89450.0%291.5%0.27<0.001100.5%0.15<0.00140.2%0.16<0.001
BP Non-user1,89250.0%1035.4%0.180.41663.5%0.080.29241.3%0.060.48
Age ≥812,613100.0%461.8%572.2%240.9%
BP User1,30750.0%110.8%0.31<0.001161.2%0.38<0.00150.4%0.260.004
BP Non-user1,30650.0%352.7%0.160.61413.1%0.210.69191.5%0.100.70
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 9
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=South.
Region=South Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients319,408100.0%8,4182.6%3,5351.1%8490.3%
BP User159,70450.0%1,5531.0%0.22<0.0016240.4%0.21<0.0011670.1%0.24<0.001
BP Non-user159,70450.0%6,8654.3%0.210.232,9111.8%0.190.236820.4%0.210.29
By Age
Age ≤20943100.0%293.1%40.4%00.0%
BP User47450.3%153.2%1.060.8700.0%NANA00.0%NANA
BP Non-user46949.7%143.0%0.512.2340.9%NANA00.0%NANA
Age 21-402,364100.0%1134.8%251.1%40.2%
BP User1,17749.8%201.7%0.20<0.00140.3%0.19<0.00120.2%1.011.00
BP Non-user1,18750.2%937.8%0.120.33211.8%0.060.5520.2%0.147.17
Age 41-506,839100.0%3294.8%731.1%100.1%
BP User3,42150.0%722.1%0.26<0.001180.5%0.32<0.00150.1%1.000.99
BP Non-user3,41850.0%2577.5%0.200.34551.6%0.190.5550.1%0.293.45
Age 51-6050,186100.0%1,9994.0%5841.2%1030.2%
BP User25,09350.0%3931.6%0.23<0.0011140.5%0.24<0.001230.1%0.29<0.001
BP Non-user25,09350.0%1,6066.4%0.210.264701.9%0.190.29800.3%0.180.46
Age 61-70111,800100.0%3,2462.9%1,1061.0%2470.2%
BP User55,89650.0%5831.0%0.21<0.0011910.3%0.21<0.001380.1%0.18<0.001
BP Non-user55,90450.0%2,6634.8%0.190.239151.6%0.180.242090.4%0.130.26
Age 71-8099,223100.0%1,9422.0%1,0291.0%2600.3%
BP User49,61850.0%3220.6%0.19<0.0011700.3%0.20<0.001550.1%0.27<0.001
BP Non-user49,60550.0%1,6203.3%0.170.228591.7%0.170.232050.4%0.200.36
Age ≥8148,053100.0%7601.6%7141.5%2250.5%
BP User24,02550.0%1480.6%0.24<0.0011270.5%0.21<0.001440.2%0.24<0.001
BP Non-user24,02850.0%6122.5%0.200.285872.4%0.180.261810.8%0.170.34
Female Patients289,263100.0%7,5192.6%3,1591.1%7450.3%
BP User144,63350.0%1,3650.9%0.21<0.0015620.4%0.21<0.0011430.1%0.24<0.001
BP Non-user144,63050.0%6,1544.3%0.200.232,5971.8%0.190.236020.4%0.200.28
By Age
Age ≤20372100.0%113.0%30.8%00.0%
BP User18549.7%63.2%1.220.7500.0%NANA00.0%NANA
BP Non-user18750.3%52.7%0.374.0731.6%NANA00.0%NANA
Age 21-401,543100.0%815.2%161.0%30.2%
BP User77049.9%141.8%0.20<0.00140.5%0.330.0820.3%2.010.62
BP Non-user77350.1%678.7%0.110.35121.6%0.111.0310.1%0.1822.22
Age 41-505,569100.0%2734.9%661.2%90.2%
BP User2,78750.0%652.3%0.30<0.001180.6%0.37<0.00150.2%1.251.00
BP Non-user2,78250.0%2087.5%0.220.39481.7%0.210.6440.1%0.334.65
Age 51-6046,012100.0%1,8194.0%5211.1%890.2%
BP User23,00750.0%3581.6%0.23<0.0011000.4%0.23<0.001160.1%0.22<0.001
BP Non-user23,00550.0%1,4616.4%0.210.264211.8%0.190.29730.3%0.130.38
Age 61-70103,825100.0%2,9482.8%1,0071.0%2180.2%
BP User51,91050.0%5171.0%0.20<0.0011770.3%0.21<0.001330.1%0.18<0.001
BP Non-user51,91550.0%2,4314.7%0.190.238301.6%0.180.251850.4%0.120.26
Age 71-8089,474100.0%1,7291.9%9151.0%2300.3%
BP User44,74250.0%2830.6%0.19<0.0011530.3%0.20<0.001470.1%0.26<0.001
BP Non-user44,73250.0%1,4463.2%0.170.227621.7%0.170.241830.4%0.190.35
Age ≥8142,468100.0%6581.5%6311.5%1960.5%
BP User21,23250.0%1220.6%0.22<0.0011100.5%0.21<0.001400.2%0.26<0.001
BP Non-user21,23650.0%5362.5%0.180.275212.5%0.170.251560.7%0.180.36
Male Patients30,145100.0%8993.0%3761.2%1040.3%
BP User15,07150.0%1881.2%0.26<0.001620.4%0.19<0.001240.2%0.30<0.001
BP Non-user15,07450.0%7114.7%0.220.303142.1%0.150.26800.5%0.190.47
By Age
Age ≤20571100.0%183.2%10.2%00.0%
BP User28950.6%93.1%0.980.9600.0%NANA00.0%NANA
BP Non-user28249.4%93.2%0.382.4910.4%NANA00.0%NANA
Age 21-40821100.0%323.9%91.1%10.1%
BP User40749.6%61.5%0.22<0.00100.0%NANA00.0%NANA
BP Non-user41450.4%266.3%0.090.5592.2%NANA10.2%NANA
Age 41-501,270100.0%564.4%70.6%10.1%
BP User63449.9%71.1%0.13<0.00100.0%NANA00.0%NANA
BP Non-user63650.1%497.7%0.060.3071.1%NANA10.2%NANA
Age 51-604,174100.0%1804.3%631.5%140.3%
BP User2,08650.0%351.7%0.23<0.001140.7%0.28<0.00170.3%1.000.99
BP Non-user2,08850.0%1456.9%0.160.33492.3%0.150.5170.3%0.352.86
Age 61-707,975100.0%2983.7%991.2%290.4%
BP User3,98650.0%661.7%0.27<0.001140.4%0.16<0.00150.1%0.21<0.001
BP Non-user3,98950.0%2325.8%0.210.36852.1%0.090.29240.6%0.080.54
Age 71-809,749100.0%2132.2%1141.2%300.3%
BP User4,87650.0%390.8%0.22<0.001170.3%0.17<0.00180.2%0.360.01
BP Non-user4,87350.0%1743.6%0.150.31972.0%0.100.29220.5%0.160.81
Age ≥815,585100.0%1021.8%831.5%290.5%
BP User2,79350.0%260.9%0.34<0.001170.6%0.25<0.00140.1%0.16<0.001
BP Non-user2,79250.0%762.7%0.210.53662.4%0.150.43250.9%0.060.46
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 10
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=West.
Region=West Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients190,534100.0%6,4373.4%2,8351.5%7200.4%
BP User95,26750.0%1,0841.1%0.19<0.0014390.5%0.18<0.0011130.1%0.19<0.001
BP Non-user95,26750.0%5,3535.6%0.180.212,3962.5%0.160.206070.6%0.150.23
By Age
Age ≤20299100.0%93.0%10.3%00.0%
BP User15050.2%10.7%0.120.0200.0%NANA00.0%NANA
BP Non-user14949.8%85.4%0.010.9610.7%NANA00.0%NANA
Age 21-401,115100.0%676.0%181.6%30.3%
BP User55750.0%173.1%0.32<0.00120.4%0.120.00110.2%0.501.00
BP Non-user55850.0%509.0%0.180.56162.9%0.030.5320.4%0.055.53
Age 41-502,806100.0%1475.2%331.2%100.4%
BP User1,40450.0%271.9%0.21<0.00130.2%0.10<0.00110.1%0.110.01
BP Non-user1,40250.0%1208.6%0.140.32302.1%0.030.3290.6%0.010.87
Age 51-6020,761100.0%9524.6%3041.5%540.3%
BP User10,38050.0%1781.7%0.22<0.001530.5%0.21<0.00170.1%0.15<0.001
BP Non-user10,38150.0%7747.5%0.180.262512.4%0.150.28470.5%0.070.33
Age 61-7062,798100.0%2,2283.5%8281.3%1960.3%
BP User31,39250.0%3411.1%0.17<0.0011180.4%0.16<0.001210.1%0.12<0.001
BP Non-user31,40650.0%1,8876.0%0.150.197102.3%0.130.201750.6%0.080.19
Age 71-8064,596100.0%1,9623.0%9151.4%2590.4%
BP User32,30350.0%3361.0%0.20<0.0011370.4%0.17<0.001430.1%0.20<0.001
BP Non-user32,29350.0%1,6265.0%0.180.227782.4%0.140.212160.7%0.140.27
Age ≥8138,159100.0%1,0722.8%7361.9%1980.5%
BP User19,08150.0%1841.0%0.20<0.0011260.7%0.20<0.001400.2%0.25<0.001
BP Non-user19,07850.0%8884.7%0.170.236103.2%0.170.241580.8%0.180.36
Female Patients168,933100.0%5,5753.3%2,4391.4%5760.3%
BP User84,47050.0%9091.1%0.19<0.0013800.4%0.18<0.001920.1%0.19<0.001
BP Non-user84,46350.0%4,6665.5%0.170.202,0592.4%0.160.204840.6%0.150.24
By Age
Age ≤20107100.0%43.7%00.0%00.0%
BP User5450.5%00.0%NANA00.0%NANA00.0%NANA
BP Non-user5349.5%47.5%NANA00.0%NANA00.0%NANA
Age 21-40658100.0%446.7%91.4%20.3%
BP User32950.0%103.0%0.27<0.00110.3%0.120.0400.0%NANA
BP Non-user32950.0%3410.3%0.130.5682.4%0.020.9820.6%NANA
Age 41-502,158100.0%1245.7%211.0%50.2%
BP User1,08050.0%211.9%0.19<0.00120.2%0.10<0.00100.0%NANA
BP Non-user1,07850.0%1039.6%0.120.30191.8%0.020.4550.5%NANA
Age 51-6018,480100.0%8354.5%2641.4%430.2%
BP User9,23950.0%1461.6%0.20<0.001460.5%0.21<0.00140.0%0.10<0.001
BP Non-user9,24150.0%6897.5%0.170.242182.4%0.150.29390.4%0.040.29
Age 61-7057,528100.0%1,9993.5%7271.3%1590.3%
BP User28,76150.0%3001.0%0.17<0.0011110.4%0.18<0.001190.1%0.14<0.001
BP Non-user28,76750.0%1,6995.9%0.150.196162.1%0.140.221400.5%0.080.22
Age 71-8057,396100.0%1,6772.9%7751.4%2080.4%
BP User28,70350.0%2841.0%0.20<0.0011180.4%0.18<0.001360.1%0.21<0.001
BP Non-user28,69350.0%1,3934.9%0.170.226572.3%0.140.211720.6%0.150.30
Age ≥8132,606100.0%8922.7%6432.0%1590.5%
BP User16,30450.0%1480.9%0.19<0.0011020.6%0.18<0.001330.2%0.26<0.001
BP Non-user16,30250.0%7444.6%0.160.235413.3%0.150.231260.8%0.180.38
Male Patients21,601100.0%8624.0%3961.8%1440.7%
BP User10,79750.0%1751.6%0.24<0.001590.5%0.17<0.001210.2%0.17<0.001
BP Non-user10,80450.0%6876.4%0.210.293373.1%0.130.231231.1%0.110.27
By Age
Age ≤20192100.0%52.6%10.5%00.0%
BP User9650.0%11.0%0.240.3700.0%NANA00.0%NANA
BP Non-user9650.0%44.2%0.032.2111.0%NANA00.0%NANA
Age 21-40457100.0%235.0%92.0%10.2%
BP User22849.9%73.1%0.420.0610.4%0.120.0410.4%NANA
BP Non-user22950.1%167.0%0.171.0583.5%0.020.9800.0%NANA
Age 41-50648100.0%233.5%121.9%50.8%
BP User32450.0%61.9%0.340.0210.3%0.090.00610.3%0.250.37
BP Non-user32450.0%175.2%0.130.88113.4%0.010.6941.2%0.032.23
Age 51-602,281100.0%1175.1%401.8%110.5%
BP User1,14150.0%322.8%0.36<0.00170.6%0.21<0.00130.3%0.370.15
BP Non-user1,14050.0%857.5%0.240.54332.9%0.090.4780.7%0.101.41
Age 61-705,270100.0%2294.3%1011.9%370.7%
BP User2,63149.9%411.6%0.21<0.00170.3%0.07<0.00120.1%0.06<0.001
BP Non-user2,63950.1%1887.1%0.150.29943.6%0.030.16351.3%0.010.24
Age 71-807,200100.0%2854.0%1401.9%510.7%
BP User3,60050.0%521.4%0.21<0.001190.5%0.15<0.00170.2%0.16<0.001
BP Non-user3,60050.0%2336.5%0.160.291213.4%0.090.25441.2%0.070.35
Age ≥815,553100.0%1803.2%931.7%390.7%
BP User2,77750.0%361.3%0.24<0.001240.9%0.34<0.00170.3%0.22<0.001
BP Non-user2,77650.0%1445.2%0.170.35692.5%0.210.55321.2%0.100.49
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 11
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=New York State.
Region=New York State Matched
AllSARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
N%N%ORp-valueN%ORp-valueN%ORp-value
LLULLLULLLUL
All Patients99,724100.0%3,5983.6%3,6073.6%6220.6%
BP User49,86250.0%7721.5%0.26<0.0018111.6%0.28<0.0011360.3%0.28<0.001
BP Non-user49,86250.0%2,8265.7%0.240.282,7965.6%0.260.304861.0%0.230.34
By Age
Age ≤20102100.0%43.9%22.0%11.0%
BP User5049.0%24.0%1.041.0012.0%1.041.0012.0%NANA
BP Non-user5251.0%23.8%0.147.6911.9%0.0617.1100.0%NANA
Age 21-40453100.0%214.6%153.3%10.2%
BP User22850.3%31.3%0.15<0.00120.9%0.140.00410.4%NANA
BP Non-user22549.7%188.0%0.040.53135.8%0.030.6500.0%NANA
Age 41-501,311100.0%775.9%362.7%40.3%
BP User65550.0%223.4%0.38<0.00181.2%0.28<0.00110.2%0.330.62
BP Non-user65650.0%558.4%0.230.63284.3%0.130.6130.5%0.033.21
Age 51-6012,732100.0%6885.4%5274.1%580.5%
BP User6,36450.0%1552.4%0.27<0.0011181.9%0.28<0.001170.3%0.410.002
BP Non-user6,36850.0%5338.4%0.230.334096.4%0.220.34410.6%0.230.73
Age 61-7032,265100.0%1,2944.0%1,1503.6%1410.4%
BP User16,13650.0%2771.7%0.26<0.0012671.7%0.29<0.001270.2%0.24<0.001
BP Non-user16,12950.0%1,0176.3%0.230.308835.5%0.250.331140.7%0.150.36
Age 71-8034,693100.0%9572.8%1,1963.4%2400.7%
BP User17,34150.0%2041.2%0.26<0.0012571.5%0.26<0.001450.3%0.23<0.001
BP Non-user17,35250.0%7534.3%0.220.319395.4%0.230.301951.1%0.170.32
Age ≥8118,168100.0%5573.1%6813.7%1771.0%
BP User9,08850.0%1091.2%0.23<0.0011581.7%0.29<0.001440.5%0.33<0.001
BP Non-user9,08050.0%4484.9%0.190.295235.8%0.240.351331.5%0.230.46
Female Patients90,567100.0%3,2553.6%3,2353.6%5370.6%
BP User45,28550.0%6871.5%0.26<0.0017261.6%0.28<0.0011080.2%0.25<0.001
BP Non-user45,28250.0%2,5685.7%0.240.282,5095.5%0.260.304290.9%0.200.31
By Age
Age ≤2033100.0%00.0%13.0%13.0%
BP User1648.5%00.0%NANA16.3%NANA16.3%NANA
BP Non-user1751.5%00.0%NANA00.0%NANA00.0%NANA
Age 21-40261100.0%166.1%83.1%10.4%
BP User13250.6%21.5%0.130.00221.5%0.320.1710.8%NANA
BP Non-user12949.4%1410.9%0.030.5764.7%0.061.5900.0%NANA
Age 41-501,032100.0%585.6%282.7%30.3%
BP User51650.0%183.5%0.430.00371.4%0.320.00700.0%NANA
BP Non-user51650.0%407.8%0.240.76214.1%0.140.7730.6%NANA
Age 51-6011,699100.0%6375.4%4824.1%470.4%
BP User5,84950.0%1382.4%0.26<0.0011101.9%0.28<0.001140.2%0.420.006
BP Non-user5,85050.0%4998.5%0.210.313726.4%0.230.35330.6%0.230.79
Age 61-7030,115100.0%1,2044.0%1,0703.6%1260.4%
BP User15,06050.0%2571.7%0.26<0.0012481.6%0.29<0.001230.2%0.22<0.001
BP Non-user15,05550.0%9476.3%0.220.308225.5%0.250.331030.7%0.140.35
Age 71-8031,385100.0%8582.7%1,0523.4%2080.7%
BP User15,68850.0%1761.1%0.25<0.0012211.4%0.26<0.001330.2%0.19<0.001
BP Non-user15,69750.0%6824.3%0.210.308315.3%0.220.301751.1%0.130.27
Age ≥8116,042100.0%4823.0%5943.7%1510.9%
BP User8,02450.0%961.2%0.24<0.0011371.7%0.29<0.001360.4%0.31<0.001
BP Non-user8,01850.0%3864.8%0.190.304575.7%0.240.351151.4%0.210.45
Male Patients9,157100.0%3433.7%3724.1%850.9%
BP User4,57750.0%851.9%0.32<0.001851.9%0.28<0.001280.6%0.490.002
BP Non-user4,58050.0%2585.6%0.250.412876.3%0.220.36571.2%0.310.77
By Age
Age ≤2069100.0%45.8%11.4%00.0%
BP User3449.3%25.9%1.031.0000.0%NANA00.0%NANA
BP Non-user3550.7%25.7%0.147.7712.9%NANA00.0%NANA
Age 21-40192100.0%52.6%73.6%00.0%
BP User9650.0%11.0%0.240.3700.0%NANA00.0%NANA
BP Non-user9650.0%44.2%0.032.2177.3%NANA00.0%NANA
Age 41-50279100.0%196.8%82.9%10.4%
BP User13949.8%42.9%0.250.0210.7%0.140.0710.7%NANA
BP Non-user14050.2%1510.7%0.080.7675.0%0.021.1300.0%NANA
Age 51-601,033100.0%514.9%454.4%111.1%
BP User51549.9%173.3%0.490.0281.6%0.21<0.00130.6%0.370.22
BP Non-user51850.1%346.6%0.270.88377.1%0.090.4481.5%0.101.42
Age 61-702,150100.0%904.2%803.7%150.7%
BP User1,07650.0%201.9%0.27<0.001191.8%0.30<0.00140.4%0.360.08
BP Non-user1,07450.0%706.5%0.160.45615.7%0.180.50111.0%0.111.14
Age 71-803,308100.0%993.0%1444.4%321.0%
BP User1,65350.0%281.7%0.38<0.001362.2%0.32<0.001120.7%0.600.16
BP Non-user1,65550.0%714.3%0.250.601086.5%0.220.47201.2%0.291.23
Age ≥812,126100.0%753.5%874.1%261.2%
BP User1,06450.0%131.2%0.20<0.001212.0%0.30<0.00180.8%0.440.05
BP Non-user1,06250.0%625.8%0.110.37666.2%0.180.50181.7%0.191.02
  1. BP: bisphosphonate; LL: lower 95% confidence interval level; NA: not applicable; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 12
Unadjusted COVID-19-Related Outcomes Stratified by Age, Sex, & Age by Sex; Matched Primary Analysis Cohort, Region=New York State.
SARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.220.210.22<0.0010.220.210.23<0.0010.240.220.26<0.001
Adjusted0.220.210.23<0.0010.230.220.24<0.0010.260.240.29<0.001
NortheastUnadjusted0.220.210.24<0.0010.240.230.26<0.0010.260.230.30<0.001
Adjusted0.230.210.24<0.0010.250.230.26<0.0010.290.260.33<0.001
MidwestUnadjusted0.230.220.25<0.0010.220.200.25<0.0010.230.190.29<0.001
Adjusted0.240.220.26<0.0010.240.220.27<0.0010.260.210.32<0.001
SouthUnadjusted0.220.210.23<0.0010.210.190.23<0.0010.240.210.29<0.001
Adjusted0.220.210.23<0.0010.220.200.24<0.0010.260.230.30<0.001
WestUnadjusted0.190.180.21<0.0010.180.160.20<0.0010.190.150.23<0.001
Adjusted0.200.180.21<0.0010.190.170.21<0.0010.200.160.25<0.001
New YorkUnadjusted0.260.240.28<0.0010.280.260.30<0.0010.280.230.34<0.001
Adjusted0.260.240.28<0.0010.280.260.31<0.0010.330.270.40<0.001
  1. LL: lower 95% confidence interval level; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 13
Statin Use Sensitivity Analysis, Unadjusted/Adjusted Odds Ratio for COVID-19-Related Outcomes, Stratified by Region and New York State.
SARS-CoV-2 TestCOVID-19 DiagnosisCOVID-19 Hospitalization
Statin Uses versus Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.900.890.91<0.0010.910.900.92<0.0011.121.091.15<0.001
Adjusted0.870.860.87<0.0010.790.780.81<0.0010.990.961.020.48
NortheastUnadjusted0.870.850.88<0.0010.880.860.90<0.0011.161.111.21<0.001
Adjusted0.850.840.87<0.0010.770.750.78<0.0011.030.981.070.22
MidwestUnadjusted0.970.950.990.021.101.071.14<0.0011.271.191.36<0.001
Adjusted0.920.900.94<0.0010.990.961.030.751.151.081.23<0.001
SouthUnadjusted0.900.880.91<0.0010.900.880.93<0.0011.000.951.060.90
Adjusted0.850.840.87<0.0010.800.780.83<0.0010.880.830.94<0.001
WestUnadjusted0.880.860.90<0.0010.830.800.86<0.0011.020.951.100.58
Adjusted0.860.830.88<0.0010.710.680.74<0.0010.870.800.94<0.001
New YorkUnadjusted0.910.890.93<0.0010.930.910.96<0.0011.211.141.29<0.001
Adjusted0.920.900.95<0.0010.790.770.82<0.0011.050.981.130.15
BP Users versus BP Non-users among Statin Users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.230.220.24<0.0010.260.250.28<0.0010.260.230.29<0.001
Adjusted0.230.220.24<0.0010.270.250.29<0.0010.280.250.32<0.001
NortheastUnadjusted0.250.230.27<0.0010.290.270.31<0.0010.280.240.34<0.001
Adjusted0.250.230.27<0.0010.290.270.32<0.0010.320.260.38<0.001
MidwestUnadjusted0.240.220.27<0.0010.220.190.25<0.0010.210.160.27<0.001
Adjusted0.250.230.29<0.0010.230.220.25<0.0010.220.170.30<0.001
SouthUnadjusted0.220.210.24<0.0010.260.230.29<0.0010.260.210.33<0.001
Adjusted0.220.200.24<0.0010.270.240.31<0.0010.280.220.36<0.001
WestUnadjusted0.200.180.22<0.0010.220.190.25<0.0010.250.200.33<0.001
Adjusted0.200.180.22<0.0010.230.200.27<0.0010.280.210.36<0.001
New YorkUnadjusted0.270.240.30<0.0010.310.280.35<0.0010.300.230.39<0.001
Adjusted0.280.250.32<0.0010.310.280.35<0.0010.330.250.44<0.001
BP Users versus BP Non-users among Statin Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.230.210.24<0.0010.210.190.23<0.0010.210.170.25<0.001
Adjusted0.240.220.25<0.0010.230.210.25<0.0010.250.210.30<0.001
NortheastUnadjusted0.250.220.27<0.0010.220.200.25<0.0010.240.190.31<0.001
Adjusted0.260.230.29<0.0010.250.220.28<0.0010.290.220.37<0.001
MidwestUnadjusted0.240.210.28<0.0010.220.180.27<0.0010.210.140.31<0.001
Adjusted0.240.200.28<0.0010.250.200.32<0.0010.260.170.39<0.001
SouthUnadjusted0.230.210.25<0.0010.190.150.22<0.0010.180.120.27<0.001
Adjusted0.240.210.27<0.0010.210.170.25<0.0010.220.150.33<0.001
WestUnadjusted0.190.170.22<0.0010.180.150.22<0.0010.160.110.25<0.001
Adjusted0.200.170.23<0.0010.190.180.21<0.0010.180.110.29<0.001
New YorkUnadjusted0.260.230.30<0.0010.260.220.30<0.0010.270.190.39<0.001
Adjusted0.260.220.31<0.0010.250.210.30<0.0010.350.230.52<0.001
  1. LL: lower 95% confidence interval level; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 14
Antihypertensive Use Sensitivity Analysis, Unadjusted/Adjusted Odds Ratio for COVID-19-Related Outcomes, Stratified by Region and New York State.
Odds of SARS-CoV-2 TestOdds of COVID-19 DiagnosisOdds of COVID-19 Hospitalization
Antihypertensive Users versus Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.910.900.92<0.0010.860.850.87<0.0011.131.101.17<0.001
Adjusted0.870.860.88<0.0010.750.740.76<0.0010.980.951.000.10
NortheastUnadjusted0.860.840.87<0.0010.830.820.85<0.0011.201.151.25<0.001
Adjusted0.820.810.83<0.0010.720.710.73<0.0011.040.991.080.10
MidwestUnadjusted1.000.981.020.981.061.031.10<0.0011.281.201.36<0.001
Adjusted0.940.910.96<0.0010.940.900.97<0.0011.111.041.190.002
SouthUnadjusted0.930.920.94<0.0010.880.860.90<0.0011.020.961.070.58
Adjusted0.880.870.89<0.0010.780.760.80<0.0010.890.840.94<0.001
WestUnadjusted0.900.880.92<0.0010.750.730.78<0.0010.990.921.060.83
Adjusted0.870.850.89<0.0010.650.620.67<0.0010.840.780.90<0.001
New YorkUnadjusted0.920.900.94<0.0010.900.870.92<0.0011.231.151.31<0.001
Adjusted0.900.870.92<0.0010.750.730.77<0.0011.010.951.090.70
BP Users versus BP Non-users among Antihypertensive Users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.230.220.24<0.0010.260.250.28<0.0010.260.230.29<0.001
Adjusted0.230.220.24<0.0010.260.250.28<0.0010.270.240.30<0.001
NortheastUnadjusted0.240.220.26<0.0010.280.260.31<0.0010.270.220.32<0.001
Adjusted0.230.210.26<0.0010.280.260.31<0.0010.290.240.34<0.001
MidwestUnadjusted0.260.230.29<0.0010.270.230.31<0.0010.270.210.35<0.001
Adjusted0.270.240.30<0.0010.280.260.30<0.0010.270.200.35<0.001
SouthUnadjusted0.230.210.25<0.0010.240.220.28<0.0010.260.200.32<0.001
Adjusted0.230.210.25<0.0010.240.210.28<0.0010.250.200.32<0.001
WestUnadjusted0.200.180.22<0.0010.210.180.25<0.0010.240.180.31<0.001
Adjusted0.200.180.22<0.0010.220.180.25<0.0010.240.180.33<0.001
New YorkUnadjusted0.260.230.29<0.0010.300.260.33<0.0010.290.220.38<0.001
Adjusted0.250.220.29<0.0010.300.260.34<0.0010.330.240.44<0.001
BP Users versus BP Non-users among Antihypertensive Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.210.200.22<0.0010.200.180.22<0.0010.210.170.25<0.001
Adjusted0.210.200.22<0.0010.220.200.24<0.0010.270.220.32<0.001
NortheastUnadjusted0.210.190.23<0.0010.220.190.24<0.0010.230.180.31<0.001
Adjusted0.220.200.25<0.0010.250.220.28<0.0010.300.220.40<0.001
MidwestUnadjusted0.220.190.25<0.0010.160.120.20<0.0010.200.130.31<0.001
Adjusted0.210.180.25<0.0010.180.140.23<0.0010.260.160.42<0.001
SouthUnadjusted0.200.180.22<0.0010.190.160.22<0.0010.220.150.32<0.001
Adjusted0.200.180.22<0.0010.210.170.25<0.0010.280.190.41<0.001
WestUnadjusted0.190.170.22<0.0010.180.150.22<0.0010.150.090.24<0.001
Adjusted0.200.170.22<0.0010.200.160.25<0.0010.190.110.32<0.001
New YorkUnadjusted0.260.230.31<0.0010.250.210.29<0.0010.230.150.36<0.001
Adjusted0.270.230.32<0.0010.260.220.31<0.0010.260.160.43<0.001
  1. LL: lower 95% confidence interval level; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 15
Antidiabetic Use Sensitivity Analysis, Unadjusted/Adjusted Odds Ratio for COVID-19-Related Outcomes, Stratified by Region and New York State.
Odds of SARS-CoV-2 TestOdds of COVID-19 DiagnosisOdds of COVID-19 Hospitalization
Antidiabetic Users versus Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.980.970.990.011.151.131.18<0.0011.501.451.56<0.001
Adjusted0.920.900.93<0.0010.880.860.90<0.0011.131.081.18<0.001
NortheastUnadjusted1.000.981.020.921.111.091.14<0.0011.551.471.64<0.001
Adjusted0.940.920.97<0.0010.840.810.86<0.0011.181.111.27<0.001
MidwestUnadjusted1.041.011.080.011.391.331.46<0.0011.611.471.76<0.001
Adjusted0.950.910.990.011.111.041.17<0.0011.301.171.44<0.001
SouthUnadjusted0.970.950.990.011.161.121.21<0.0011.391.291.50<0.001
Adjusted0.900.880.93<0.0010.910.870.95<0.0011.040.951.140.40
WestUnadjusted0.910.880.94<0.0011.071.011.120.011.431.301.58<0.001
Adjusted0.860.820.89<0.0010.800.750.85<0.0010.970.861.090.60
New YorkUnadjusted1.061.031.10<0.0011.151.111.19<0.0011.591.461.72<0.001
Adjusted1.061.021.100.0070.870.830.90<0.0011.181.071.300.001
BP Users versus BP Non-users among Antidiabetic Users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.260.240.28<0.0010.290.270.32<0.0010.280.240.33<0.001
Adjusted0.260.240.28<0.0010.290.270.32<0.0010.290.250.34<0.001
NortheastUnadjusted0.280.240.32<0.0010.320.280.35<0.0010.290.230.36<0.001
Adjusted0.280.240.32<0.0010.310.270.35<0.0010.300.240.39<0.001
MidwestUnadjusted0.270.220.33<0.0010.300.240.38<0.0010.280.190.41<0.001
Adjusted0.270.220.34<0.0010.320.260.41<0.0010.290.190.42<0.001
SouthUnadjusted0.290.260.33<0.0010.310.260.36<0.0010.350.260.47<0.001
Adjusted0.300.260.34<0.0010.300.250.36<0.0010.360.260.48<0.001
WestUnadjusted0.190.160.22<0.0010.200.170.25<0.0010.210.150.30<0.001
Adjusted0.190.160.23<0.0010.210.170.26<0.0010.220.150.31<0.001
New YorkUnadjusted0.330.270.40<0.0010.340.290.39<0.0010.350.260.49<0.001
Adjusted0.320.260.40<0.0010.320.280.36<0.0010.400.280.56<0.001
BP Users versus BP Non-users among Antidiabetic Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.240.230.26<0.0010.240.220.26<0.0010.240.200.29<0.001
Adjusted0.250.230.27<0.0010.250.230.28<0.0010.270.220.33<0.001
NortheastUnadjusted0.240.220.28<0.0010.260.220.29<0.0010.250.190.34<0.001
Adjusted0.250.220.29<0.0010.270.240.32<0.0010.280.200.39<0.001
MidwestUnadjusted0.270.220.32<0.0010.220.170.30<0.0010.260.160.42<0.001
Adjusted0.280.240.31<0.0010.230.170.31<0.0010.260.160.45<0.001
SouthUnadjusted0.240.210.27<0.0010.250.200.30<0.0010.290.200.43<0.001
Adjusted0.240.210.27<0.0010.240.210.28<0.0010.330.220.49<0.001
WestUnadjusted0.230.200.27<0.0010.180.140.24<0.0010.130.070.23<0.001
Adjusted0.230.200.28<0.0010.200.150.26<0.0010.150.080.28<0.001
New YorkUnadjusted0.300.250.37<0.0010.300.250.36<0.0010.220.140.36<0.001
Adjusted0.300.250.37<0.0010.310.250.37<0.0010.240.140.41<0.001
  1. LL: lower 95% confidence interval level; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 16
Antidepressant Use Sensitivity Analysis, Unadjusted/Adjusted Odds Ratio for COVID-19-Related Outcomes, Stratified by Region and New York State.
Odds of SARS-CoV-2 TestOdds of COVID-19 DiagnosisOdds of COVID-19 Hospitalization
Antidepressant Users versus Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted1.041.031.05<0.0010.710.700.72<0.0010.810.780.83<0.001
Adjusted1.000.991.010.610.650.640.66<0.0010.750.730.78<0.001
NortheastUnadjusted1.010.991.020.540.710.690.72<0.0010.840.800.88<0.001
Adjusted0.970.950.990.0010.650.630.66<0.0010.770.730.82<0.001
MidwestUnadjusted1.101.081.12<0.0010.750.720.78<0.0010.840.780.90<0.001
Adjusted1.051.031.07<0.0010.690.660.71<0.0010.780.730.84<0.001
SouthUnadjusted1.041.021.05<0.0010.680.660.70<0.0010.740.700.79<0.001
Adjusted0.990.981.010.490.640.620.66<0.0010.720.680.77<0.001
WestUnadjusted1.041.021.060.000.700.670.73<0.0010.770.700.84<0.001
Adjusted0.990.971.020.460.640.610.66<0.0010.700.640.77<0.001
New YorkUnadjusted1.000.971.030.860.770.740.80<0.0010.830.760.91<0.001
Adjusted0.980.951.010.270.720.700.75<0.0010.770.700.85<0.001
BP Users versus BP Non-users among Antidepressant Users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.270.260.28<0.0010.300.280.32<0.0010.310.270.36<0.001
Adjusted0.270.250.28<0.0010.300.280.32<0.0010.330.280.38<0.001
NortheastUnadjusted0.280.260.31<0.0010.330.300.37<0.0010.360.290.45<0.001
Adjusted0.280.250.30<0.0010.320.290.36<0.0010.370.290.47<0.001
MidwestUnadjusted0.300.270.34<0.0010.260.220.31<0.0010.250.180.34<0.001
Adjusted0.300.260.34<0.0010.270.220.33<0.0010.260.180.36<0.001
SouthUnadjusted0.260.240.29<0.0010.270.230.31<0.0010.320.240.41<0.001
Adjusted0.260.240.28<0.0010.270.230.32<0.0010.320.240.43<0.001
WestUnadjusted0.250.220.28<0.0010.270.220.32<0.0010.290.200.41<0.001
Adjusted0.240.210.27<0.0010.290.280.30<0.0010.330.230.48<0.001
New YorkUnadjusted0.300.260.34<0.0010.330.280.38<0.0010.240.160.36<0.001
Adjusted0.300.250.34<0.0010.310.270.37<0.0010.250.160.39<0.001
BP Users versus BP Non-users among Antidepressant Non-users
ORLLULp valueORLLULp valueORLLULp value
AllUnadjusted0.200.190.22<0.0010.220.200.24<0.0010.240.200.28<0.001
Adjusted0.210.190.22<0.0010.230.210.25<0.0010.270.220.32<0.001
NortheastUnadjusted0.210.190.24<0.0010.230.200.26<0.0010.250.190.32<0.001
Adjusted0.220.190.25<0.0010.240.220.25<0.0010.290.220.39<0.001
MidwestUnadjusted0.220.190.26<0.0010.230.180.28<0.0010.280.190.39<0.001
Adjusted0.210.180.25<0.0010.260.240.27<0.0010.320.220.47<0.001
SouthUnadjusted0.200.180.22<0.0010.210.180.25<0.0010.210.150.30<0.001
Adjusted0.200.180.23<0.0010.230.190.27<0.0010.220.160.32<0.001
WestUnadjusted0.180.160.21<0.0010.200.160.25<0.0010.200.130.30<0.001
Adjusted0.190.160.22<0.0010.200.200.21<0.0010.220.140.35<0.001
New YorkUnadjusted0.260.220.32<0.0010.270.230.32<0.0010.290.190.43<0.001
Adjusted0.260.230.30<0.0010.260.220.32<0.0010.350.220.54<0.001
  1. LL: lower 95% confidence interval level; OR: odds ratio; UL: upper 95% confidence interval level.

Appendix 2—table 17
Bone-Rx” Cohort (All Regions), Patient Characteristics Pre/Post Match.
"Bone-Rx" Cohort / All Observations Unmatched"Bone-Rx" Cohort / All Observations Matched
AllBP Non-userBP Userp-valueAllBP Non-userBP Userp-value
N%N%N%N%N%N%
All Patients502,895100.0%50,84410.1%452,05189.9%100,996100.0%50,49850.0%50,49850.0%
Age
≤201,1640.2%360.1%1,1280.2%<0.001670.1%360.1%310.1%0.97
21-403,5010.7%4100.8%3,0910.7%7900.8%4030.8%3870.8%
41-509,6311.9%1,0802.1%8,5511.9%2,1072.1%1,0692.1%1,0382.1%
51-6072,13914.3%6,41812.6%65,72114.5%12,77712.7%6,39512.7%6,38212.6%
61-70171,68734.1%14,80929.1%156,87834.7%29,50929.2%14,75129.2%14,75829.2%
71-80157,87731.4%16,15231.8%141,72531.4%32,12931.8%16,05531.8%16,07431.8%
≥8186,89617.3%11,93923.5%74,95716.6%23,61723.4%11,78923.3%11,82823.4%
Gender
Female451,79089.8%44,35487.2%407,43690.1%<0.00188,55287.7%44,23587.6%44,31787.8%0.43
Male51,10510.2%6,49012.8%44,6159.9%12,44412.3%6,26312.4%6,18112.2%
Region
Midwest85,39117.0%9,42418.5%75,96716.8%<0.00118,72018.5%9,36018.5%9,36018.5%1.00
Northeast135,86727.0%16,13931.7%119,72826.5%31,98631.7%15,99331.7%15,99331.7%
South178,11835.4%17,23233.9%160,88635.6%34,28033.9%17,14033.9%17,14033.9%
West103,51920.6%8,04915.8%95,47021.1%16,01015.9%8,00515.9%8,00515.9%
Insurance
Commercial164,15032.6%17,09233.6%147,05832.5%<0.00133,97733.6%16,96333.6%17,01433.7%0.91
Dual33,9696.8%2,5625.0%31,4076.9%5,0565.0%2,5475.0%2,5095.0%
Medicaid84,51416.8%7,03413.8%77,48017.1%13,92513.8%6,98613.8%6,93913.7%
Medicare220,26243.8%24,15647.5%196,10643.4%48,03847.6%24,00247.5%24,03647.6%
PCP Visit 2019
No181,99636.2%18,13035.7%163,86636.2%0.00935,94335.6%17,97935.6%17,96435.6%0.92
Yes320,89963.8%32,71464.3%288,18563.8%65,05364.4%32,51964.4%32,53464.4%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI1.051.911.992.710.951.76<0.0011.932.591.932.601.922.590.76
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 18
Bone-Rx” Cohort (Region=Northeast), Patient Characteristics Pre/Post Match.
"Bone-Rx" Cohort / Region=Northeast Unmatched"Bone-Rx" Cohort / Region=Northeast Matched
AllBP Non-userBP Userp-valueAllBP Non-userBP Userp-value
N%N%N%N%N%N%
All Patients135,867100.0%16,13911.9%119,72888.1%31,986100.0%15,99350.0%15,99350.0%
Age
≤202450.2%≤100.1%2360.2%<0.001150.0%≤100.1%≤100.0%0.99
21-408910.7%1270.8%7640.6%2500.8%1240.8%1260.8%
41-502,3401.7%2981.8%2,0421.7%5701.8%2901.8%2801.8%
51-6020,06914.8%2,05912.8%18,01015.0%4,08812.8%2,04912.8%2,03912.7%
61-7045,89633.8%4,80229.8%41,09434.3%9,52629.8%4,76729.8%4,75929.8%
71-8042,82831.5%5,26732.6%37,56131.4%10,46532.7%5,22632.7%5,23932.8%
≥8123,59817.4%3,57722.2%20,02116.7%7,07222.1%3,52822.1%3,54422.2%
Gender
Female122,48590.2%14,11587.5%108,37090.5%<0.00128,15788.0%14,06287.9%14,09588.1%0.57
Male13,3829.8%2,02412.5%11,3589.5%3,82912.0%1,93112.1%1,89811.9%
Insurance
Commercial37,81027.8%4,51728.0%33,29327.8%<0.0018,92727.9%4,45927.9%4,46827.9%0.99
Dual8,4346.2%8295.1%7,6056.4%1,6375.1%8245.2%8135.1%
Medicaid25,29618.6%2,08212.9%23,21419.4%4,12212.9%2,06712.9%2,05512.8%
Medicare64,32747.3%8,71154.0%55,61646.5%17,30054.1%8,64354.0%8,65754.1%
PCP Visit 2019
No56,59341.7%6,72641.7%49,86741.7%0.9513,30741.6%6,65441.6%6,65341.6%0.99
Yes79,27458.3%9,41358.3%69,86158.3%18,67958.4%9,33958.4%9,34058.4%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI1.061.891.972.700.931.71<0.0011.892.571.892.581.892.570.91
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 19
“Bone-Rx” Cohort (Region=Midwest), Patient Characteristics Pre/Post Match.
"Bone-Rx" Cohort / Region=Midwest Unmatched"Bone-Rx" Cohort / Region=Midwest Matched
AllBP Non-userBP Userp-valueAllBP Non-userBP Userp-value
N%N%N%N%N%N%
All Patients85,391100.0%9,42411.0%75,96789.0%18,720100.0%9,36050.0%9,36050.0%
Age
≤202740.3%≤100.1%2680.4%<0.001130.1%≤100.1%≤100.1%1.00
21-406720.8%790.8%5930.8%1540.8%780.8%760.8%
41-501,8862.2%2022.1%1,6842.2%3892.1%2002.1%1892.0%
51-6013,52215.8%1,28413.6%12,23816.1%2,55913.7%1,28013.7%1,27913.7%
61-7031,25636.6%2,76029.3%28,49637.5%5,51229.4%2,75429.4%2,75829.5%
71-8023,88728.0%2,76629.4%21,12127.8%5,49229.3%2,74829.4%2,74429.3%
≥8113,89416.3%2,32724.7%11,56715.2%4,60124.6%2,29424.5%2,30724.6%
Gender
Female76,69689.8%8,11886.1%68,57890.3%<0.00116,22386.7%8,10286.6%8,12186.8%0.68
Male8,69510.2%1,30613.9%7,3899.7%2,49713.3%1,25813.4%1,23913.2%
Insurance
Commercial34,49440.4%3,36135.7%31,13341.0%<0.0016,69935.8%3,34535.7%3,35435.8%0.96
Dual4,0424.7%4364.6%3,6064.7%8524.6%4294.6%4234.5%
Medicaid8,85610.4%7337.8%8,12310.7%1,4417.7%7297.8%7127.6%
Medicare37,99944.5%4,89451.9%33,10543.6%9,72852.0%4,85751.9%4,87152.0%
PCP Visit 2019
No32,03737.5%3,33035.3%28,70737.8%<0.0016,62835.4%3,31235.4%3,31635.4%0.95
Yes53,35462.5%6,09464.7%47,26062.2%12,09264.6%6,04864.6%6,04464.6%
Continuous Outcomes
meanSDmeanSDmeanSDp-valuemeanSDmeanSDmeanSDp-value
CCI1.122.022.122.830.991.86<0.0012.052.722.062.722.052.720.91
  1. BP: bisphosphonate; CCI: Charlson Comorbidity Index; PCP: primary care physician; SD: standard deviation.

Appendix 2—table 20
Bone-Rx” Cohort (Region=South), Patient Characteristics Pre/Post Match.
"Bone-Rx" Cohort / Region=South Unmatched"Bone-Rx" Cohort / Region=South Matched
AllBP Non-userBP Userp-valueAllBP Non-userBP Userp-value
N%N%N%N%N%N%
All Patients178,118100.0%17,2329.7%160,88690.3%34,280100.0%17,14050.0%17,14050.0%
Age
≤204900.3%160.1%4740.3%<0.001310.1%160.1%150.1%1.00
21-401,3130.7%1360.8%1,1770.7%2620.8%1340.8%1280.7%
41-503,8662.2%4452.6%3,4212.1%8842.6%4442.6%4402.6%
51-6027,38915.4%2,29613.3%25,09315.6%4,57413.3%2,29013.4%2,28413.3%
61-7061,03834.3%5,14229.8%55,89634.7%10,27130.0%5,12929.9%5,14230.0%
71-8056,12631.5%5,52132.0%50,60531.5%10,99032.1%5,49332.0%5,49732.1%
≥8127,89615.7%3,67621.3%24,22015.1%7,26821.2%3,63421.2%3,63421.2%
Gender
Female160,99490.4%15,17988.1%145,81590.6%<0.00130,32288.5%15,14988.4%15,17388.5%0.69
Male17,1249.6%2,05311.9%15,0719.4%3,95811.5%1,99111.6%1,96711.5%
Insurance
Commercial66,33237.2%7,04240.9%59,29036.9%<0.00114,05241.0%7,00740.9%7,04541.1%0.95
Dual14,8298.3%7694.5%14,0608.7%1,5234.4%7694.5%7544.4%
Medicaid23,49213.2%1,84310.7%21,64913.5%3,63910.6%1,82910.7%1,81010.6%
Medicare