Post-approval Safety Surveillance Study of Golimumab in the Treatment of Rheumatic Disease Using a United States Healthcare Claims Database.

BACKGROUND AND
OBJECTIVE: Golimumab is a fully human anti-tumor necrosis factor monoclonal antibody approved for the treatment of rheumatoid arthritis (RA), psoriatic arthritis (PsA), and ankylosing spondylitis (AS). This study estimated rates of prespecified outcomes in patients with RA, PsA or AS initiating golimumab versus matched patients initiating non-biologic systemic (NBS) medications.
METHODS: Patients enrolled in a US health plan with rheumatic disease who initiated a study medication were accrued between April 2009 and November 2014. Golimumab initiators were matched by propensity score to NBS initiators in a 1:4 ratio. Outcomes were identified through September 2015. As-treated, as-matched, and nested case-control (NCC) analyses were conducted in the matched cohorts. Sensitivity analyses evaluated the impact of residual confounding and nondifferential misclassification of exposure and outcomes.
RESULTS: Risks of outcomes were similar between golimumab and NBS initiators. In the as-treated analysis, the rate ratio (RR) for depression was elevated during current golimumab use versus golimumab non-use in the NBS cohort [RR 1.45, 95% confidence interval (CI) 1.31-1.61]. This finding was not replicated in as-matched (RR 1.08, 95% CI 0.97-1.19) or NCC (odds ratio 1.01, 95% CI 0.78-1.31) analyses, which focused on incident cases. Sensitivity analyses suggest that depression was sensitive to misclassification, and the RR changed from greater than to less than one across a plausible range of specificity.
CONCLUSIONS: This study suggests that there is no association between exposure to golimumab and an increased risk of prespecified outcomes. Increased depression risk in the as-treated analysis was not replicated in other analyses and may be associated with residual imbalance in baseline history or severity of depression.

Key Points

Introduction

Golimumab (Janssen Biotech, Inc.) is a fully-human tumor necrosis factor (TNF)-alpha inhibitor indicated for the treatment of rheumatoid arthritis (RA), psoriatic arthritis (PsA), ankylosing spondylitis (AS), and ulcerative colitis in the USA, with additional indications for non-radiographic axial spondyloarthritis and polyarticular juvenile idiopathic arthritis in the European Union. Golimumab was first approved in the USA in 2009 for once-monthly subcutaneous administration in RA, PsA, and AS and in 2013 for intravenous administration in RA.

RA is characterized by chronic, progressive inflammation of the joints that may lead to deformity, disability, and in some cases, premature death [1, 2]. The overall prevalence of RA ranges between 0.3 and 1% among adults worldwide and is more common with increased age and among women [3, 4]. AS and PsA are related chronic inflammatory arthritides that are distinct from and less prevalent than RA. The goals of rheumatic disease management are to decrease pain, to prevent or control joint damage, and to prevent loss of function [4]. Current treatments include nonsteroidal anti-inflammatory drugs, corticosteroids, and disease-modifying anti-rheumatic drugs (DMARDs) [5-7]. Tumor necrosis factor inhibitors (TNFi) were the first biologic DMARDs (bDMARDs) approved for the treatment of these conditions and have revolutionized pharmacologic treatment by preventing disease progression [8-10].

Because of inherent immunological abnormalities associated with rheumatic disease and the immunosuppressive effects of therapeutic agents, patients may be at an increased risk of infection, hematologic malignancies, and other autoimmune diseases [11-13]. The incidence of infections requiring hospitalization or parenteral antibiotics in RA patients prior to the introduction of bDMARDs is estimated to be in the range of 0.02–0.12 per person-year, with infection-related deaths 6–9 times higher in RA patients compared to the general population [13, 14]. Risk factors of infection in RA patients include age, extra-articular manifestations, leukopenia, comorbidities, use of corticosteroids, skin breakdown, and joint surgery [14, 15].

Case reports and post-marketing data have documented the occurrence of serious and opportunistic infections in patients treated with bDMARDs [13, 16–23], which should not be used in patients with active infections [14]. However, not all studies have found an increased risk associated with their use [13]. As per current prescribing guidelines, use of bDMARDs is contraindicated in patients with severe active infections, and patients receiving these agents should be monitored for signs and symptoms of infection [14]. In addition to infections, malignancy, autoimmune conditions, hepatotoxicity, hepatitis B reactivation, and heart failure have been observed with use of TNFi [14].

Clinical trials of patients with RA, PsA, and AS have shown that golimumab is effective and generally well tolerated [24-27]. Due to the rarity of some outcomes and potential differences between patients enrolled in clinical trials and those under standard care, comprehensive assessments of risks associated with golimumab, a newer TNFi, among large populations reflecting routine practice are warranted. In this study, we aimed to provide a comprehensive overview of the safety of golimumab using claims data from a large US health insurer along with validation of outcomes using linked medical records and National Death Index (NDI) records. This study question is important as there is a lack of agent-specific safety data. We sought to estimate the risks of prespecified outcomes in a cohort of patients with RA, PsA, or AS initiating golimumab compared to a cohort of similar patients initiating non-biologic systemic (NBS) treatments, which is considered the standard of care and preferred comparator for biologic treatment [1, 2, 9].

Patients and Methods

Data Source

This prospective, observational, cohort study was conducted in the Optum Research Database (ORD), a proprietary research database containing claims and enrollment data dating back to 1993 for approximately 14 million members of a large US health plan with both medical and pharmacy benefit coverage and an average of 2.5 years of enrollment. The medical and pharmacy claims for these individuals form a longitudinal record of reimbursed medical services, irrespective of treatment site, along with detailed information on drug dispensing. The database is routinely updated and maintained. Access to a subset of patients’ medical charts (i.e. approximately 35% of all patients within the ORD) allows for confirmation of outcomes identified through claims. The underlying insured population is geographically diverse across the USA, representing approximately 4% of its population.

This study was approved by the New England Institutional Review Board and Privacy Board.

Patient Population

Patients who initiated golimumab or NBS treatments (methotrexate, azathioprine, cyclosporine, penicillamine, hydroxychloroquine, sulfasalazine, apremilast, tofacitinib, leflunomide, or gold compounds) were accrued from 24 April 2009 (date of first availability of golimumab in the USA) until 30 November 2014. The date of cohort entry (i.e. index date) is the date of the first claim of the cohort-defining drug with no previous claim of that specific drug in the prior 6 months (baseline period). Eligible patients included those of any age with complete medical coverage and pharmacy benefits who were continuously enrolled for at least 6 months prior (baseline period) to the index date.

The operational definition of an initiator included patients who were naïve to antirheumatic treatments, those who switched to the cohort-defining therapy from another antirheumatic treatment, and those who started the cohort-defining therapy as add-on treatment to existing therapy. For example, golimumab initiators included patients who began therapy with golimumab as the first prescription treatment for rheumatic disease observed in the claims database, and those who switched from or added to another prescription rheumatic disease treatment regimen. The cohort for NBS treatment initiators may include, for example, methotrexate initiators without any other antirheumatic treatment in baseline, and methotrexate initiators without prior methotrexate use but with prior use of another antirheumatic (including other NBS and biologic) treatment during baseline.

Initiators were required to have at least one claim for RA, PsA, or AS [International Classification of Diseases (ICD), 9th Revision diagnosis codes 714.xx, 696.0x, and 720.0x, respectively] in the 6 months before or up to 3 months following the index date.

Propensity Score Modeling and Matching

For the prediction model discriminating golimumab initiators from patients in the NBS cohort, using claims and membership data, a set of potential a priori confounders in addition to the 100 most frequently occurring diagnoses, procedures, and drug dispensing were included in the logistic regression model for the propensity score calculation. Age, sex, region, number of days enrolled before index date, calendar time, cost (total and associated with rheumatic disease), and healthcare utilization variables were forced into the model. Additional predictors were retained based on a stepwise regression technique with an entry criterion p-value of ≤ 0.10 and an exit criterion p-value of ≥ 0.30.

The golimumab cohort was matched by propensity score to NBS initiators using a standard greedy matching algorithm [28-30]. In order to achieve adequate power in the context of a medication with a large pool of comparators, the final cohort included up to 4 matched comparators for each golimumab patient, with the comparison cohort member selected at random from all potential comparators whose propensity scores were within the variable caliper applied by the matching algorithm that allows a maximum difference of 0.10 to that of the golimumab patients.

Outcome Identification and Follow-Up

Prespecified outcomes of interest were identified by a claim for an outpatient physician visit or inpatient hospitalization associated with the specific disease, procedure, or drug code as listed on claims for medical services. These outcomes include:

serious infection (i.e. requiring hospitalization or parenteral antibiotic therapy)

tuberculosis (TB)/non-TB mycobacterial infection

malignancy

lymphoma

systemic hypersensitivity

congestive heart failure (CHF)

hepatotoxicity

new onset hypertension (identified among those without a baseline claim for hypertension)

autoimmune disease

hematologic reaction

depression

mortality.

See Electronic Supplementary Materials (ESM) Appendix 1 for codes that were used to identify these outcomes.

Outcome identification began the day after cohort entry through the earliest of: (1) health plan disenrollment; (2) death; or (3) 30 September 2015. For systemic hypersensitivity and mortality, follow-up began on the index date to identify those outcomes that occurred immediately after exposure. Patients were censored for subsequent occurrences of the same outcome but were eligible for others. Person-time was calculated from the day after index date (or on the index date for systemic hypersensitivity and mortality) through end of follow-up.

Medical Record Review and Refinement of Outcome Definitions

Medical record review was undertaken to confirm potential outcomes identified in the claims data between 24 April 2009 and 31 May 2014. Medical records were adjudicated by a panel of 3 clinicians to determine whether each potential case met the clinical definition of the outcome. Using the results from adjudication, a positive predictive value (PPV) was calculated for each outcome as the number of confirmed cases divided by the number of charts adjudicated, excluding those cases with insufficient information to determine case status.

Outcomes with high PPVs (> 75%) were included in analyses using the original claims definition. In order to improve the PPV and increase the likelihood of identifying true cases, the claims definitions for outcomes with low-to-moderate PPVs (≤ 75%) were revised by requiring the codes (diagnosis, procedure, drug) used to identify each outcome to appear in the principal (primary) position on an inpatient or outpatient claim. PPVs were recalculated for the revised definition of each outcome as the number of confirmed cases with codes in the primary position on the healthcare claim divided by the total number of confirmed cases and non-cases with codes in the primary position and were included in analyses using the revised definition.

National Death Index

The NDI is considered a gold standard mortality assessment in studies based on health insurance claims data, and a NDI Plus search (including fact, date, and causes of death) was conducted on all golimumab initiators and matched comparators who disenrolled from the health plan during the study period with no evidence of re-enrollment.

After restricting to patients considered to have a high probability of being a true match based on an algorithm developed by the NDI, deaths occurring within 60 days of health plan disenrollment were selected as cases. Since exposure information is only available during enrolled time and deaths occurring after disenrollment may have been exposed to other medications, a 60-day window was chosen as a trade-off between completeness of mortality identification and uncertainty regarding medication exposure. Cause-specific mortality outcomes were identified from the NDI data using ICD-10 coded causes of death.

The mortality outcome used in statistical analyses included those deaths identified in the claims data through end of follow-up with an NDI search used to validate vital status for these patients.

Statistical Analyses

All analyses used SAS (Statistical Analysis Software™) version 9.2 (SAS Institute).

For non-cancer outcomes, the primary analysis was the as-treated analysis, where incidence rates (IR) per 1000 person-years and incidence rate ratios (IRR) were estimated during periods of current drug use and non-use among the matched golimumab and matched NBS cohorts. Drug exposure was characterized based on claims records during follow-up, including date of dispensing and days supplied or recommended number of days between injections/infusions. A 60-day “grace period” was added to account for medication non-adherence and the uncertainty surrounding the duration of biologic effect for these drugs.

Current use: dispensing/injection/infusion date + days supplied or recommended number of days between injection/infusion + 60 days. With each new dispensing/injection/infusion, a patient continues or reinitiates on-therapy current status

Non-use: days during which patients are not using the study drug(s) of interest, i.e. neither current nor prior use of the drug(s) in question during the follow-up. Person-time attributed to non-use of golimumab came from the matched NBS cohort before any golimumab was dispensed during follow-up; and vice versa, person-time attributed to non-use of NBS came from the matched golimumab cohort before any NBS was dispensed during follow-up.

Multivariable Poisson regression was used to estimate the IRR of non-cancer outcomes during periods of current drug use compared with non-use. Each model was stratified by matching ratio (number of NBS matches for each golimumab initiator) and adjusted for age and sex.

For cancer outcomes, the primary analysis was the as-matched analysis where IRs per 1000 person-years were calculated based on the cohort-initiating drug. In the as-matched analysis, Cox proportional hazards regression was used to compare golimumab to NBS initiators in the matched cohorts. Models for all outcomes were stratified by matching ratio and adjusted for propensity score. To avoid overfitting or lack of model convergence, the propensity score was used for adjustment since the number of covariates that could be included in the Cox proportional hazards models was limited.

To address the different potential mechanisms of cancer occurrence, the incidence of malignancy and lymphoma was estimated using a cumulative dose measure (considered a secondary analysis for cancer outcomes) in the matched cohorts. With each claim for the dispensing of a study drug during follow-up, cumulative dose was updated, and patients were categorized into mutually exclusive categories of dose (i.e. 1–5, 6–10, and ≥ 11 dispensings). If a patient had a cancer that occurred during follow-up, it was counted within the dose category it occurred, and the patient was no longer at risk for that outcome. Similar to the as-treated analysis, the Poisson regression models were stratified by matching ratio and adjusted for age and sex.

As an additional assessment of medication exposure in relation to occurrence of outcomes, a nested case–control (NCC) analysis was performed among propensity score-matched cohorts. Cases were defined as the first outcome identified during follow-up with no baseline claims for the outcome. Each risk set included all individuals from the propensity score-matched quintuplet of the larger drug cohorts who were still at risk of the event (controls). After excluding those controls with baseline claims for a specific outcome, for each case, up to 2 controls enrolled in the health plan were randomly selected among the cohort members in the risk sets defined by the case (matching preserved), and a case/control date was assigned as the outcome date of the matched case.

Exposure to study medications was assessed among cases and controls during the 6 months prior to the case/control date, and age, sex, region, and calendar time were assessed on this date. Cases and controls were compared with respect to study drug exposure (golimumab and NBS), with multivariable-adjusted odds ratios (OR) generated from conditional logistic regression models (conditioned on the matched set) and adjusted for age, sex, region, and calendar time.

Sensitivity Analyses

Three sensitivity analyses were conducted. First, as-treated and as-matched analyses were repeated using a narrowed definition of depression where initiators with a baseline claim of a diagnosis of depression or dispensing of an antidepressant were excluded from each matched cohort. During follow-up, the first occurrence of a diagnosis for depression or antidepressant claim was identified. If both a diagnosis and dispensing occurred within 14 days of one another, then the date of the later diagnosis or dispensing was assigned as the outcome date.

Second, residual confounding in the propensity score model by covariates not included in the propensity score model or directly ascertainable in claims data (e.g. disease severity, duration) was assessed. A standard epidemiological spreadsheet was used to assess the robustness of the results in the presence of an unmeasured confounder [31]. The results of the primary analyses and ranges of association between an unmeasured confounder and disease outcome were entered into the spreadsheet and used to calculate the difference from the observed result if the unmeasured confounder was accounted for. In addition to using the adjusted IRR from the primary analyses, the value of the observed lower 95% CI was also plotted to determine the range in which the 95% CI would cross the null. The prevalence of the unmeasured confounder was fixed at 40%, which provided close to maximum potential confounding in many situations while the prevalence of the exposure was calculated as the number of person-years during current use divided by the number of person-years during non-use. The resulting figures include the range of unmeasured confounder characteristics assuming varying levels of strength of the confounder-disease association and strength of association with current golimumab use that would result in a change in the point estimates. Detailed assumptions for conducting this analysis are described in ESM Appendix 2.

For the last sensitivity analysis, another standard epidemiological spreadsheet was used to adjust primary analysis point estimates by exploring different pairs of sensitivity and specificity values in order to assess different degrees of nondifferential misclassification [32]. To assess the potential effect of nondifferential exposure misclassification, sensitivity and specificity were set to 99, 95, and 90%, respectively, to assess the effect on primary as-treated analysis point estimates in the golimumab cohort matched to NBS. Similarly, to determine the potential effect of nondifferential outcome misclassification, sensitivity was assumed to be 99% while specificity was defined as a transformation of the PPV based on contingency table analysis. The original PPV was used for malignancy, lymphoma, and autoimmune disease while the revised PPV was used for the remaining outcomes. A PPV of 50% was used to estimate misclassification bias for mortality. The specificity was varied by ± 5% to assess the effect on primary analysis point estimates across a plausible range of alternate specificity values. Detailed assumptions for conducting this analysis are described in ESM Appendix 3.

Results

Cohort Characteristics

From 24 April 2009 through 30 November 2014, 1515 golimumab (97.0% dispensed subcutaneous formulation) and 48,975 NBS initiators were identified, and 1337 golimumab initiators (88.3% of those eligible) were matched to 4227 NBS initiators (3.2:1 matching ratio).

Prior to matching, differences were present between the golimumab and NBS initiators with respect to baseline characteristics (Table 1). Afterwards, balance improved on most of these variables (Table 2), with absolute standardized differences generally less than 0.1. A higher proportion of matched golimumab initiators than NBS initiators had ≥ 4 different treatments for rheumatic disease (33.7% vs 28.5%, standardized difference = 0.11). Both disease-specific and total healthcare utilization costs were also higher in the golimumab group (standardized differences = 0.11).

Table 1

Distribution of select baseline characteristicsa for the golimumab and non-biologic cohorts before propensity score matching, identified 24 April 2009 through 30 November 2014, Optum Research Database

Baseline characteristics	Golimumab		Non-biologic		Standardized difference
	N = 1515		N = 48,975
Age [years, mean (median), SD]	48.1 (50.0)	11.5	49.9 (51.0)	13.1	0.15
Age group (years, N%)
0–5	0	0.0	53	0.1	0.00
6–11	0	0.0	201	0.4	0.00
12–17	1	0.1	523	1.1	0.13
18–29	114	7.5	2803	5.7	0.07
30–39	234	15.4	6447	13.2	0.06
40–49	403	26.6	11,469	23.4	0.07
50–59	489	32.3	16,055	32.8	0.01
60–64	215	14.2	7024	14.3	0.00
65 or older	59	3.9	4400	9.0	0.21
Sex (N%)
Female	981	64.8	36,455	74.4	0.21
Male	534	35.2	12,520	25.6	0.21
History of cohort-defining disease (N%)
Rheumatoid arthritis	1034	68.3	44,237	90.3	0.56
Psoriatic arthritis	460	30.4	6111	12.5	0.45
Ankylosing spondylitis	252	16.6	1568	3.2	0.46
History of outcomes (N%)
Serious infection	55	3.6	2383	4.9	0.06
TB/non-TB mycobacterial infection	8	0.5	73	0.1	0.07
Hepatotoxicity	53	3.5	1731	3.5	0.00
Systemic hypersensitivity	19	1.3	993	2.0	0.06
Malignancy	28	1.8	1877	3.8	0.12
Lymphoma	1	0.1	193	0.4	0.06
Autoimmune disease	14	0.9	2791	5.7	0.27
Hypertension	583	38.5	19,729	40.3	0.04
Congestive heart failure	12	0.8	610	1.2	0.04
Hematologic reaction	15	1.0	901	1.8	0.07
Depression	462	30.5	13,928	28.4	0.05
History of comorbidities and other diagnoses (N%)
Comorbidities
Cardiovascular events	62	4.1	2051	4.2	0.01
Chronic obstructive pulmonary disease	29	1.9	1573	3.2	0.08
Diabetes mellitus	141	9.3	5086	10.4	0.04
Renal insufficiency/failure	21	1.4	1359	2.8	0.10
Stroke	6	0.4	637	1.3	0.10
Asthma	60	4.0	2934	6.0	0.09
Hepatic insufficiency/failure	66	4.4	2350	4.8	0.02
Other diagnoses
Smoking (tobacco use disorder)	44	2.9	1899	3.9	0.06
Alcohol abuse	6	0.4	209	0.4	0.00
Hepatitis B infection	2	0.1	73	0.1	0.00
History of medications (N%)
Golimumab (excluding index drug)	0	0.0	143	0.3	0.00
Other anti-TNF biologics	668	44.1	4964	10.1	0.83
Adalimumab	286	18.9	1531	3.1	0.52
Etanercept	291	19.2	2251	4.6	0.46
Infliximab	95	6.3	1136	2.3	0.20
Certolizumab	7	0.5	101	0.2	0.05
Non-anti-TNF biologics	51	3.4	583	1.2	0.15
Tocilizumab	3	0.7	59	0.4	0.04
Abatacept	42	2.8	400	0.8	0.15
Anakinra	0	0.0	16	0.0	0.00
Rituximab	7	0.5	114	0.2	0.05
Non-biologics (excluding index drug)	687	45.3	7609	15.5	0.69
Apremilast	0	0.0	0	0.0	0.00
Tofacitinib	4	0.9	14	0.1	0.11
Leflunomide	116	7.7	728	1.5	0.30
Methotrexate	491	32.4	3316	6.8	0.68
Azathioprine	15	1.0	226	0.5	0.06
Cyclosporine	7	0.5	226	0.5	0.00
Penicillamine	1	0.1	5	0.0	0.05
Hydroxychloroquine	129	8.5	3254	6.6	0.07
Sulfasalazine	63	4.2	802	1.6	0.16
Gold compounds	1	0.1	20	0.0	0.05
Oral or IV corticosteroids	894	59.0	30,830	63.0	0.08
NSAIDs	662	43.7	24,645	50.3	0.13
Specified health services (N% unless otherwise noted)
Visits to rheumatologists
0	96	6.3	6677	13.6	0.25
1	67	4.4	7367	15.0	0.36
2	156	10.3	10,499	21.4	0.31
3 or more	1196	78.9	24,432	49.9	0.64
Number of different drugs dispensed
0	0	0.0	20	0.0	0.00
1	32	2.1	808	1.6	0.04
2	66	4.4	1890	3.9	0.03
3	75	5.0	3143	6.4	0.06
4 or more	1342	88.6	43,114	88.0	0.02
Number of different treatments for rheumatic disease
0	0	0.0	0	0.0	0.00
1	122	8.1	7015	14.3	0.20
2	405	26.7	18,838	38.5	0.25
3	455	30.0	18,256	37.3	0.16
4 or more	533	35.2	4866	9.9	0.64
Emergency room visits
0	1213	80.1	38,262	78.1	0.05
1	184	12.1	6233	12.7	0.02
2 or more	118	7.8	4480	9.1	0.05
Hospitalizations
0	1442	95.2	45,428	92.8	0.10
1	63	4.2	2940	6.0	0.08
2 or more	10	0.7	607	1.2	0.05
RA, PsA, or AS healthcare utilization costs [$, mean (median), SD]	8864.20 (6152.60)	8704.30	2646.60 (506.90)	7057.50	0.79
Total healthcare utilization costs [$, mean (median), SD]	12,971.00 (10,253.80)	12,016.40	8675.40 (3755.80)	17,260.40	0.29

SD standard deviation, TB tuberculosis, TNF tumor necrosis factor, IV intravenous, NSAID nonsteroidal anti-inflammatory drug, CCP cyclic citrullinated peptide, DNA deoxyribonucleic acid, RA rheumatoid arthritis, PsA psoriatic arthritis, AS ankylosing spondylitis

aAside from age and sex, which were identified on the index date, characteristics were identified during the 6-month period before and including the index date (baseline period) except where otherwise indicated

Table 2

Distribution of select baseline characteristicsa for the golimumab and non-biologic cohorts after propensity score matching, identified 24 April 2009 through 30 November 2014, Optum Research Database

Baseline characteristics	Matched initiators				Standardized difference	Unmatched golimumab initiators
	Golimumab		Non-biologic
	N = 1337		N = 4227			N = 178
Age [years, mean (median), SD]	48.4 (50.0)	11.5	48.8 (50.0)	11.8	0.03	46.1 (47.0)	11.2
Age group (years, N%)
0–5	0	0.0	0	0.0	0.00	0	0.0
6–11	0	0.0	0	0.0	0.00	0	0.0
12–17	1	0.1	3	0.1	0.00	0	0.0
18–29	98	7.3	321	7.6	0.01	16	9.0
30–39	203	15.2	592	14.0	0.03	31	17.4
40–49	348	26.0	1105	26.1	0.00	55	30.9
50–59	436	32.6	1421	33.6	0.02	53	29.8
60–64	193	14.4	585	13.8	0.02	22	12.4
65 or older	58	4.3	200	4.7	0.02	1	0.6
Sex (N%)
Female	888	66.4	2814	66.6	0.00	93	52.2
Male	449	33.6	1413	33.4	0.00	85	47.8
History of cohort-defining disease (N%)
Rheumatoid arthritis	941	70.4	3139	74.3	0.09	93	52.2
Psoriatic arthritis	395	29.5	1214	28.7	0.02	65	36.5
Ankylosing spondylitis	189	14.1	498	11.8	0.07	63	35.4
History of outcomes (N%)
Serious infection	49	3.7	144	3.4	0.02	6	3.4
TB/non-TB mycobacterial infection	6	0.4	8	0.2	0.04	2	1.1
Hepatotoxicity	48	3.6	132	3.1	0.03	5	2.8
Systemic hypersensitivity	17	1.3	65	1.5	0.02	2	1.1
Malignancy	26	1.9	122	2.9	0.07	2	1.1
Lymphoma	1	0.1	6	0.1	0.00	0	0.0
Autoimmune disease	13	1.0	52	1.2	0.02	1	0.6
Hypertension	514	38.4	1558	36.9	0.03	69	38.8
Congestive heart failure	11	0.8	31	0.7	0.01	1	0.6
Hematologic reaction	11	0.8	49	1.2	0.04	4	2.2
Depression	411	30.7	1181	27.9	0.06	51	28.7
History of comorbidities and other diagnoses (N%)
Comorbidities
Cardiovascular events	58	4.3	144	3.4	0.05	4	2.2
Chronic obstructive pulmonary disease	28	2.1	93	2.2	0.01	1	0.6
Diabetes mellitus	125	9.3	393	9.3	0.00	16	9.0
Renal insufficiency/failure	20	1.5	88	2.1	0.05	1	0.6
Stroke	6	0.4	33	0.8	0.05	0	0.0
Asthma	50	3.7	183	4.3	0.03	10	5.6
Hepatic insufficiency/failure	61	4.6	181	4.3	0.02	5	2.8
Other diagnoses
Smoking (tobacco use disorder)	43	3.2	138	3.3	0.01	1	0.6
Alcohol abuse	6	0.4	18	0.4	0.00	0	0.0
Hepatitis B infection	2	0.1	10	0.2	0.03	0	0.0
History of medications (N%)
Golimumab (excluding index drug)	0	0.0	60	1.4	0.00	0	0.0
Other anti-TNF biologics	553	41.4	1593	37.7	0.08	115	64.6
Adalimumab	222	16.6	585	13.8	0.08	64	36.0
Etanercept	250	18.7	744	17.6	0.03	41	23.0
Infliximab	84	6.3	257	6.1	0.01	11	6.2
Certolizumab	5	0.4	24	0.6	0.03	2	1.1
Non-anti-TNF biologics	48	3.6	150	3.5	0.01	3	1.7
Tocilizumab	3	0.9	10	1.0	0.01	0	0.0
Abatacept	39	2.9	110	2.6	0.02	3	1.7
Anakinra	0	0.0	7	0.2	0.00	0	0.0
Rituximab	7	0.5	24	0.6	0.01	0	0.0
Non-biologics (excluding index drug)	584	43.7	1721	40.7	0.06	103	57.9
Apremilast	0	0.0	0	0.0	0.00	0	0.0
Tofacitinib	2	0.6	4	0.4	0.03	2	1.8
Leflunomide	94	7.0	241	5.7	0.05	22	12.4
Methotrexate	411	30.7	1172	27.7	0.07	80	44.9
Azathioprine	13	1.0	35	0.8	0.02	2	1.1
Cyclosporine	6	0.4	22	0.5	0.02	1	0.6
Penicillamine	1	0.1	2	0.0	0.05	0	0.0
Hydroxychloroquine	117	8.8	346	8.2	0.02	12	6.7
Sulfasalazine	57	4.3	181	4.3	0.00	6	3.4
Gold compounds	1	0.1	6	0.1	0.00	0	0.0
Oral or IV corticosteroids	790	59.1	2348	55.5	0.07	104	58.4
NSAIDs	586	43.8	1800	42.6	0.02	76	42.7
Specified health services (N% unless otherwise noted)
Visits to rheumatologists
0	87	6.5	308	7.3	0.03	9	5.1
1	65	4.9	197	4.7	0.01	2	1.1
2	140	10.5	479	11.3	0.03	16	9.0
3 or more	1045	78.2	3243	76.7	0.04	151	84.8
Number of different drugs dispensed
0	0	0.0	0	0.0	0.00	0	0.0
1	31	2.3	116	2.7	0.03	1	0.6
2	56	4.2	193	4.6	0.02	10	5.6
3	73	5.5	245	5.8	0.01	2	1.1
4 or more	1177	88.0	3673	86.9	0.03	165	92.7
Number of different treatments for rheumatic disease
0	0	0.0	0	0.0	0.00	0	0.0
1	118	8.8	422	10.0	0.04	4	2.2
2	359	26.9	1221	28.9	0.05	46	25.8
3	409	30.6	1378	32.6	0.04	46	25.8
4 or more	451	33.7	1206	28.5	0.11	82	46.1
Emergency room visits
0	1085	81.2	3438	81.3	0.00	128	71.9
1	157	11.7	495	11.7	0.00	27	15.2
2 or more	95	7.1	294	7.0	0.00	23	12.9
Hospitalizations
0	1271	95.1	4034	95.4	0.01	171	96.1
1	57	4.3	158	3.7	0.03	6	3.4
2 or more	9	0.7	35	0.8	0.01	1	0.6
RA, PsA, or AS healthcare utilization costs [$, mean (median), SD]	8375.90 (5484.00)	8354.70	7228.60 (1799.10)	11,521.40	0.11	12,532.40 (10,666.80)	10,289.60
Total healthcare utilization costs [$, mean (median), SD]	12,452.40 (9675.70)	11,790.20	11,072.60 (7893.50)	14,043.90	0.11	16,866.10 (14,459.90)	12,983.20

SD standard deviation, TB tuberculosis, TNF tumor necrosis factor, IV intravenous, NSAID nonsteroidal anti-inflammatory drug, CCP cyclic citrullinated peptide, DNA deoxyribonucleic acid, RA rheumatoid arthritis, PsA psoriatic arthritis, AS ankylosing spondylitis

aAside from age and sex, which were identified on the index date, characteristics were identified during the 6-month period before and including the index date (baseline period) except where otherwise indicated

Outcome Chart Abstraction and Adjudication Results

Medical records were obtained and abstracted for 777 of 917 subjects with potential study outcomes (84.7%). Of these, 384 were assessed to have definite evidence supporting the diagnosis of the outcome while 364 did not contain supporting evidence. Determination of case status was not possible for 29 charts given the information available in the medical record.

The PPVs for malignancy, lymphoma, and autoimmune disease were > 75% (Table 3) so the original claims definitions for these outcomes were used for analyses. Claims definitions were refined for the remaining outcomes with moderate (30–75%) or low (< 30%) PPVs to require that the codes used to identify them occur in the primary position. Revisions resulted in increases in the PPV for outcomes other than serious infection, systemic hypersensitivity, and new onset hypertension, whose PPVs remained largely unchanged (Table 3).

Table 3

Results of chart abstraction and adjudication including calculation of PPVs, outcomes identified 24 April 2009 through 31 May 2014, Optum Research Database

Outcome of interest	Number of charts requested	Number of charts adjudicated			PPVa (%)	95% CI	Revised algorithm
		Confirmed case	Confirmed non-case	Insufficient information			Confirmed case with codes in the primary position	Confirmed non-case with codes in the primary position	PPVb (%)	95% CI
Serious infection	92	53	29	1	64.6	54.3–75.0	52	29	64.2	53.8–74.6
TB/non-TB mycobacterial infection	44	12	24	1	33.3	17.9–48.7	11	13	45.8	25.9–65.8
Hepatotoxicity	93	18	51	7	26.1	15.7–36.4	8	13	38.1	17.3–58.9
Systemic hypersensitivity	94	35	43	2	44.9	33.8–55.9	18	26	40.9	26.4–55.4
Malignancyc	45	27	7	3	79.4	65.8–93.0	24	3	88.9	77.0–100.0
Lymphomac	39	26	8	2	76.5	62.2–90.7	23	7	76.7	61.5–91.8
Autoimmune diseasec	47	35	4	0	89.7	80.2–99.3	22	3	88.0	75.3–100.0
New onset hypertension	93	19	51	6	27.1	16.7–37.6	8	22	26.7	10.8–42.5
Congestive heart failure	91	40	42	1	48.8	38.0–59.6	25	22	53.2	38.9–67.5
Hematologic reaction	91	54	21	1	72.0	61.8–82.2	21	5	80.8	65.6–95.9
Depression	96	31	42	0	42.5	31.1–53.8	5	2	71.4	38.0–100.0

PPV positive predictive value, CI confidence interval, TB tuberculosis

aPPV calculated as number of confirmed cases divided by total number of charts adjudicated excluding those cases with insufficient information

bPPV calculated as number of confirmed cases with codes in the primary position divided by total number of confirmed cases and non-cases with codes in the primary position

cBased on one round of medical chart abstraction and adjudication

NDI Results

Table 4 lists the number of NDI-identified death matches occurring within 60 days of health plan disenrollment for the matched golimumab and NBS cohort and the unmatched golimumab cohort by ICD-10 category. Person-time was calculated from the index date to the death date provided by the NDI or from the index date to the earliest of end of enrollment or data cut-off date (30 September 2015) plus 60 days, and crude IRs were calculated in order to compare the cohorts.

Table 4

NDI ICD-10 Causes of death in matched golimumab and non-biologic cohorts and in the unmatched golimumab cohort, deaths identified 24 April 2009 through 30 September 2015, Optum Research Database

ICD-10 diagnosis code category	ICD-10 diagnosis category description	Golimumab	Non-biologic	Unmatched golimumab
		N = 1337	N = 4227	N = 178
		na = 4	na = 12	na = 0
		Person-yearsb = 3269.75	Person-yearsb = 10,012.41	Person-yearsb = 382.99
A00-B99	Certain infectious and parasitic diseases	0	0	0
C00-D49	Neoplasms	1	7	0
D50-D89	Diseases of the blood and blood-forming organs and certain disorders involving the immune mechanism	0	0	0
E00-E90	Endocrine, nutritional and metabolic diseases	1	0	0
F00-F99	Mental, behavioral and neurodevelopmental disorders	0	0	0
G00-G99	Diseases of the nervous system	0	0	0
H00-H59	Diseases of the eye and adnexa	0	0	0
H60-H95	Diseases of the ear and mastoid process	0	0	0
I00-I99	Diseases of the circulatory system	1	0	0
J00-J99	Diseases of the respiratory system	0	3	0
K00-K95	Diseases of the digestive system	0	1	0
L00-L99	Diseases of the skin and subcutaneous tissue	0	0	0
M00-M99	Diseases of the musculoskeletal system and connective tissue	0	0	0
N00-N99	Diseases of the genitourinary system	0	0	0
O00-O99	Pregnancy, childbirth, and the puerperium	0	0	0
P00-P96	Certain conditions originating in the perinatal period	0	0	0
Q00-Q99	Congenital malformations, deformations, and chromosomal abnormalities	0	0	0
R00-R99	Symptoms, signs, and abnormal clinical and laboratory findings, not elsewhere classified	0	0	0
S00-T98	Injury, poisoning, and certain other consequences of external causes	0	0	0
V00-Y99	External causes of morbidity	0	0	0
Z00-Z99	Factors influencing health status and contact with health services	1	1	0
U00-U85	Codes for special purposes	0	0	0

NDI National Death Index; ICD-10 International Classification of Disease, tenth revision; TNF tumor necrosis factor

aNumber with a confirmed death in the National Death Index

bPerson-time was calculated from the index date to the death date provided by the NDI or from the index date to the earliest of end of enrollment or data cut-off date (30 September 2015) plus 60 days

Four deaths were identified in the golimumab cohort matched to NBS for a crude IR of 1.22 per 1000 person-years (4 deaths/3270 person-years), and no deaths were identified in the unmatched golimumab cohort (383 person-years).

Outcome Analyses

In as-treated analyses, most findings were consistent with no increase in risk of the outcomes in the golimumab cohort (Table 5). Compared to golimumab non-use in the matched NBS cohort, the risk for depression was increased during current golimumab use (RR 1.45, 95% CI 1.31–1.61). There was no association between cumulative dose of golimumab and elevated risk of malignancy and lymphoma (Table 6).

Table 5

Incidence rates and incidence rate ratios of outcomes, as-treated analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Cohort	Number of events	Person-years	IR	95% CI	IRRa	95% CI
Serious infection	Golimumab
	Current	40	1880.32	21.27	(15.20–28.97)	0.91	(0.65–1.27)
	Non-use	207	8553.85	24.20	(21.02–27.73)	Reference
	Non-biologic
	Current	172	7226.47	23.80	(20.38–27.64)	0.86	(0.58–1.27)
	Non-use	30	1180.65	25.41	(17.14–36.27)	Reference
TB/non-TB primary mycobacterial infection	Golimumab
	Current	4	1914.40	2.09	(0.57–5.35)	1.17	(0.39–3.51)
	Non-use	16	8858.00	1.81	(1.03–2.93)	Reference
	Non-Biologic
	Current	14	7439.40	1.88	(1.03–3.16)	0.76	(0.22–2.68)
	Non-use	3	1208.41	2.48	(0.51–7.26)	Reference
Hepatotoxicity	Golimumab
	Current	29	1888.58	15.36	(10.28–22.05)	1.02	(0.68–1.52)
	Non-use	131	8628.52	15.18	(12.69–18.02)	Reference
	Non-biologic
	Current	94	7293.46	12.89	(10.42–15.77)	0.56	(0.36–0.86)
	Non-use	26	1168.97	22.24	(14.53–32.59)	Reference
Systemic hypersensitivity	Golimumab
	Current	25	1894.01	13.20	(8.54–19.49)	0.78	(0.51–1.19)
	Non-use	150	8644.64	17.35	(14.69–20.36)	Reference
	Non-biologic
	Current	123	7305.92	16.84	(13.99–20.09)	1.30	(0.75–2.27)
	Non-use	14	1196.45	11.70	(6.40–19.63)	Reference
Autoimmune disease	Golimumab
	Current	19	1900.07	10.00	(6.02–15.62)	0.71	(0.44–1.15)
	Non-use	124	8689.47	14.27	(11.87–17.01)	Reference
	Non-biologic
	Current	128	7282.66	17.58	(14.66–20.90)	3.00	(1.40–6.43)
	Non-use	7	1200.27	5.83	(2.34–12.02)	Reference
New onset hypertension	Golimumab
	Current	111	1799.66	61.68	(50.74–74.28)	1.09	(0.88–1.34)
	Non-use	457	8002.38	57.11	(51.99–62.59)	Reference
	Non-biologic
	Current	398	6809.98	58.44	(52.84–64.48)	0.99	(0.76–1.29)
	Non-use	65	1122.87	57.89	(44.68–73.78)	Reference
Congestive heart failure	Golimumab
	Current	10	1912.76	5.23	(2.51–9.61)	1.11	(0.56–2.22)
	Non-use	43	8830.38	4.87	(3.52–6.56)	Reference
	Non-Biologic
	Current	41	7408.97	5.53	(3.97–7.51)	0.83	(0.37–1.86)
	Non-use	7	1202.01	5.82	(2.34–12.00)	Reference
Hematologic reaction	Golimumab
	Current	10	1911.92	5.23	(2.51–9.62)	0.83	(0.42–1.62)
	Non-use	57	8776.83	6.49	(4.92–8.41)	Reference
	Non-biologic
	Current	55	7383.52	7.45	(5.61–9.70)	2.02	(0.73–5.59)
	Non-use	4	1205.97	3.32	(0.90–8.49)	Reference
Depression	Golimumab
	Current	479	1376.97	347.86	(317.41–380.46)	1.45	(1.31–1.61)
	Non-use	1478	5991.55	246.68	(234.26–259.59)	Reference
	Non-biologic
	Current	1551	5054.21	306.87	(291.79–322.53)	0.98	(0.85–1.12)
	Non-use	244	842.81	289.51	(254.32–328.21)	Reference
Mortality	Golimumab
	Current	2	1924.49	1.04	(0.13–3.75)	0.27	(0.06–1.10)
	Non-use	37	8894.58	4.16	(2.93–5.73)	Reference
	Non-biologic
	Current	21	7472.01	2.81	(1.74–4.30)	1.37	(0.32–5.89)
	Non-use	2	1213.31	1.65	(0.20–5.95)	Reference

IR incidence rate, representing number of events per 1000 person-years; CI confidence interval; IRR incidence rate ratio; TB tuberculosis

aAll models were stratified by matching ratio and adjusted for age and sex

Table 6

Incidence rates and incidence rate ratios of malignancy and lymphoma outcomes, cumulative dose analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Number of study drug dispensings	Number of events	Person-years	IR	95% CI	IRRa	95% CI
Malignancy	Golimumab
	1–5	48	1794.95	26.74	(19.72–35.46)	0.65	(0.48–0.89)
	6–10	13	677.22	19.20	(10.22–32.83)	0.38	(0.22–0.68)
	11+	16	1073.69	14.90	(8.52–24.20)	0.43	(0.26–0.71)
	None	280	8196.08	34.16	(30.28–38.41)	Reference
	Non-biologic
	1–5	203	4952.10	40.99	(35.55–47.04)	1.57	(1.00–2.45)
	6–10	46	1973.85	23.30	(17.06–31.09)	0.75	(0.45–1.26)
	11+	86	3753.89	22.91	(18.32–28.29)	0.99	(0.62–1.59)
	None	22	1062.11	20.71	(12.98–31.36)	Reference
Lymphoma	Golimumab
	1–5	3	1852.47	1.62	(0.33–4.73)	0.74	(0.21–2.66)
	6–10	0	699.22	0.00	(0.00–4.28)	0.00	(0.00–0.00)
	11+	1	1126.31	0.89	(0.02–4.95)	0.49	(0.06–3.70)
	None	16	8637.11	1.85	(1.06–3.01)	Reference
	Non-biologic
	1–5	9	5146.32	1.75	(0.80–3.32)	1.29	(0.16–10.26)
	6–10	3	2059.52	1.46	(0.30–4.26)	0.78	(0.08–7.66)
	11+	7	4007.64	1.75	(0.70–3.60)	1.71	(0.21–13.97)
	None	1	1101.63	0.91	(0.02–5.06)	Reference

IR incidence rate, representing number of events per 1000 person-years; CI confidence interval; IRR incidence rate ratio

aAll models are stratified by matching ratio and adjusted for age and sex

As-matched analyses in Table 7 showed a decreased risk for malignancy and mortality in golimumab versus NBS initiators (RR 0.66, 95% CI 0.50–0.87 and RR 0.41, 95% CI 0.17–0.98, respectively). No elevated rate ratios were seen for any of the outcomes in the as-matched analyses, including depression (RR 1.08, 95% CI 0.97–1.19).

Table 7

Incidence rates and incidence rate ratios of outcomes, as-matched analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Cohort	Number of events	Person-years	IR	95% CI	IRRa	95% CI
Serious infection	Golimumab	75	2925.17	25.64	(20.17–32.14)	1.08	(0.83–1.42)
	Non-biologic	211	8966.87	23.53	(20.46–26.93)	Reference
TB/non-TB primary mycobacterial infection	Golimumab	5	3036.86	1.65	(0.53–3.84)	0.80	(0.28–2.27)
	Non-biologic	16	9280.67	1.72	(0.99–2.80)	Reference
Hepatotoxicity	Golimumab	46	2951.06	15.59	(11.41–20.79)	1.10	(0.78–1.54)
	Non-biologic	135	9040.88	14.93	(12.52–17.67)	Reference
Systemic hypersensitivity	Golimumab	42	2978.04	14.10	(10.16–19.06)	0.83	(0.59–1.17)
	Non-biologic	155	9051.82	17.12	(14.53–20.04)	Reference
Malignancy	Golimumab	64	2926.78	21.87	(16.84–27.92)	0.66	(0.50–0.87)
	Non-biologic	293	8815.17	33.24	(29.54–37.27)	Reference
Lymphoma	Golimumab	4	3038.02	1.32	(0.36–3.37)	0.81	(0.27–2.45)
	Non-biologic	16	9277.08	1.72	(0.99–2.80)	Reference
Autoimmune disease	Golimumab	28	2989.27	9.37	(6.22–13.54)	0.70	(0.46–1.07)
	Non-biologic	125	9109.37	13.72	(11.42–16.35)	Reference
New onset hypertension	Golimumab	149	2793.28	53.34	(45.12–62.63)	0.92	(0.76–1.11)
	Non-biologic	475	8379.38	56.69	(51.70–62.02)	Reference
Congestive heart failure	Golimumab	17	3004.47	5.66	(3.30–9.06)	1.16	(0.66–2.06)
	Non-biologic	45	9250.58	4.86	(3.55–6.51)	Reference
Hematologic reaction	Golimumab	14	3018.13	4.64	(2.54–7.78)	0.73	(0.40–1.32)
	Non-biologic	60	9191.20	6.53	(4.98–8.40)	Reference
Depression	Golimumab	516	1968.58	262.12	(239.99–285.74)	1.08	(0.97–1.19)
	Non-biologic	1,521	6253.87	243.21	(231.14–255.75)	Reference
Mortality	Golimumab	6	3050.69	1.97	(0.72–4.28)	0.41	(0.17–0.98)
	Non-biologic	38	9318.32	4.08	(2.89–5.60)	Reference

IR incidence rate, representing number of events per 1000 person-years; CI confidence interval; IRR incidence rate ratio; TB tuberculosis

aModels are stratified by matching ratio and adjusted for propensity score

In NCC analyses (Table 8), no associations were observed between golimumab exposure and most outcomes. Golimumab use was associated with a decreased risk of systemic hypersensitivity and mortality (OR 0.57, 95% CI 0.37–0.90 and OR 0.23, 95% CI 0.08–0.67, respectively). No association between golimumab exposure and depression was seen (OR 1.01, 95% CI 0.78–1.31).

Table 8

Odds ratios of outcomes comparing treatment exposures, nested case–control analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Treatment exposure	Case	Control	ORa	95% CI
Serious infection	Golimumab	58	101	0.87	(0.62–1.21)
	Non-biologic	889	1310	0.98	(0.84–1.14)
TB/non-TB primary mycobacterial infection	Golimumab	3	4	1.48	(0.28–7.79)
	Non-biologic	53	67	1.60	(0.74–3.44)
Hepatotoxicity	Golimumab	28	43	0.94	(0.58–1.55)
	Non-biologic	426	692	0.86	(0.69–1.06)
Systemic hypersensitivity	Golimumab	28	73	0.57	(0.37–0.90)
	Non-biologic	627	938	0.96	(0.81–1.16)
Malignancy	Golimumab	52	84	1.06	(0.73–1.54)
	Non-biologic	851	1200	1.17	(0.99–1.39)
Lymphoma	Golimumab	2	7	0.62	(0.11–3.50)
	Non-biologic	56	77	1.16	(0.61–2.21)
Autoimmune disease	Golimumab	19	32	1.03	(0.56–1.89)
	Non-biologic	506	654	1.91	(1.48–2.45)
New onset hypertension	Golimumab	105	117	1.28	(0.96–1.69)
	Non-biologic	1411	1782	1.20	(1.06–1.37)
Congestive heart failure	Golimumab	13	28	0.87	(0.42–1.81)
	Non-biologic	228	384	0.64	(0.46–0.88)
Hematologic reaction	Golimumab	13	22	1.03	(0.50–2.13)
	Non-biologic	240	366	1.03	(0.75–1.40)
Depression	Golimumab	109	146	1.01	(0.78–1.31)
	Non-biologic	1700	2377	0.95	(0.85–1.06)
Mortality	Golimumab	6	19	0.23	(0.08–0.67)
	Non-biologic	171	288	0.69	(0.47–1.01)

OR odds ratio, representing number of events per 1000 person-years; CI confidence interval; TB tuberculosis

aAll models were stratified by matched set and adjusted for age, sex, region, and calendar time

Results using the restricted definition of incident cases (as opposed to prevalent cases) of depression had corresponding 95% CIs that suggested chance differences (Tables 9 and 10). The results of this sensitivity as-treated analysis reversed those seen in the primary as-treated analysis; the risk of incident depression included the null during current golimumab use compared to golimumab non-use (RR 0.71, 95% CI 0.30–1.69) (Table 9).

Table 9

Sensitivity analysis of incidence rates and incidence rate ratios of depression, as-treated analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Cohort	Number of events	Person-years	IR	95% CI	IRRa	95% CI
Depression	Golimumab
	Current	6	1916.99	3.13	(1.15–6.81)	0.71	(0.30–1.69)
	Non-use	38	8815.21	4.31	(3.05–5.92)	Reference
	Non-biologic
	Current	27	7419.66	3.64	(2.40–5.29)	0.50	(0.23–1.06)
	Non-use	9	1208.41	7.45	(3.41–14.14)	Reference

IR incidence rate, representing number of events per 1000 person-years; CI confidence interval; IRR incidence rate ratio

aAll models were stratified by matching ratio and adjusted for age and sex

Table 10

Sensitivity analysis of incidence rates and incidence rate ratios of depression, as-matched analysis, matched treatment cohorts, identified 24 April 2009 through 30 November 2014 with follow-up through 30 September 2015, Optum Research Database

Outcome of interest	Cohort	Number of events	Person-years	IR	95% CI	IRRa	95% CI
Depression	Golimumab	16	3029.41	5.28	(3.02–8.58)	1.27	(0.70–2.30)
	Non-biologic	40	9236.82	4.33	(3.09–5.90)	Reference

IR incidence rate, representing number of events per 1000 person-years; CI confidence interval; IRR incidence rate ratio

aModels are stratified by matching ratio and adjusted for propensity score

The assessment of the effect of residual confounding (i.e. as a result of an unmeasured or unavailable confounding variable omitted from the propensity score model) indicates that strong risk factors for the outcome that are also not balanced among the matched golimumab and NBS cohorts would be required in order to change the results of the analyses substantially. Most known independent risk factors were already included in the propensity score model; therefore, any unmeasured confounder of sufficient strength to alter findings would also have to be independent of the confounders that were included in the propensity score. See ESM Appendix 2 for outcome-specific results and figures.

In the assessment of nondifferential exposure or outcome misclassification, estimates in the golimumab cohort were only slightly altered, suggesting no change in inferences related to the point estimates. For depression, the observed as-treated result of an increased risk in current golimumab use versus non-use changed substantially when corrected for potential misclassification, suggesting that this result was sensitive to misclassification. See ESM Appendix 3 for outcome-specific results.

Discussion

This post-approval observational study assessed the risk of several safety outcomes in a population of golimumab and NBS users over a 6-year period. The as-treated analysis showed no elevation in risk for most of the outcomes associated with use of golimumab compared to NBS, with the exception of depression during current golimumab use. In the cumulative drug exposure analysis for malignancy and lymphoma, no evidence of a dose–response effect for golimumab was observed. In the as-matched analysis that does not account for changes in exposure during follow-up, there were no elevated rates in the golimumab cohort versus the NBS cohort; similarly, no elevated rates were seen in the NCC analysis of golimumab.

A study strength is the use of automated medical claims that reflect routine care within a large, US insurance database, providing generalizability to patients aged < 65 with commercial health insurance. The comprehensive nature of the database means that any billable medical service would be recorded, including rare events. However, there are limitations inherent to using claims databases. Patients were identified based on claims for medication dispensings, but it is not known whether the medication was taken as prescribed. Actual use of the drug is inferred from the dispensings, which may result in exposure misclassification. Sensitivity analyses were conducted to examine the impact of a range of assumed misclassification of exposure on the effect estimates and indicated that nondifferential misclassification of exposure likely biased the effect estimates toward the null or did not change the inference of the point estimates. Also, diagnosis codes used for outcome identification do not necessarily indicate true disease. Accordingly, some degree of outcome misclassification may be present [33]. Medical record confirmation was used to revise the outcome claims definitions applied to the full population, which improved the specificity of outcome definitions.

The PPVs associated with each study outcome reflect the accuracy of the claims data and allow an interpretation of the RR with the knowledge that cases of outcomes with high to moderate PPVs are likely to be real cases of disease. For such outcomes, the results presented would be close to those that might be observed from a study that included confirmation of all study outcomes [34-37]. Among those outcomes with lower PPVs, the observed RRs in this study might diverge from that of RRs based on confirmed cases, and this divergence would in most cases tend to bias the effect estimate towards the null, potentially obscuring either a positive or negative treatment effect. The original and revised PPVs for depression in this study (42.5 and 71.4%, respectively) were consistent with a range of PPVs identified in a review that used various database algorithms (31.5–98.8%) [38], and with a study that focused on incident depression (48.6%) [33]. The results of the sensitivity analyses indicate that outcome misclassification of the extent indicated by the PPVs would not change effect estimates substantially, with the exception of depression. Additional limitations include limited follow-up time and lack of information on disease duration and disease activity.

Propensity score matching created cohorts that were similar with respect to available baseline data on underlying risk factors related to outcomes, claims proxies for disease activity and severity, healthcare utilization, and other factors that may be considered by clinicians to inform treatment selection. The labelled indications for golimumab do not specify a treatment sequence (such as requiring prior failure of NBS medications), although preferential prescribing (channeling) of sicker patients to biologic agents is plausible and reflected by the inability to find eligible matches for some golimumab patients (n = 178). However, this study accounted for an extensive list of important confounders, and comparison with NBS provides useful evidence for the benefit-risk balance in real-world practice not only for golimumab but also for other biologics when the majority of patients are treated by NBS. A limitation to this approach is the potential for confounding by unmeasured factors which may not be balanced by matching on the variables included in the propensity score model. Of note, a study drawing on medical chart data demonstrated that propensity score matching can lead to balance in potential confounders not directly measurable in claims (e.g. smoking status, BMI) but correlated with covariates included in the propensity score model [30]. A sensitivity analysis was conducted to address residual confounding associated with propensity score modeling and matching, but there are limitations associated with this analysis: (1) it was constrained to one binary confounder, which may not be informative if several confounders were unmeasured and the joint effect was unknown, and (2) it did not provide an assessment of the magnitude of existing residual confounding.

Given the purpose of the study to assess the safety profile of golimumab in the real-world and the value obtained by having a larger patient population and more potential outcomes, combining the different rheumatologic indications was warranted. Specifically, regarding the inclusion of ICD-9 diagnosis codes for juvenile idiopathic arthritis (714.3x), before matching, 5 patients in the golimumab cohort (5/1515 = 0.3%) and 557 patients in the NBS cohort (557/48,975 = 1.1%) had a claim for 714.3x. Among these 562 pre-matched initiators, zero golimumab and 3 NBS patients also had a claim for rheumatoid arthritis (ICD-9 714.0) during the baseline period, indicating little overlap in coding of these 2 conditions. After matching, the numbers with a claim for ICD-9 714.3 × dropped to 4 golimumab (4/1337 = 0.3%) and 16 non-biologic initiators (16/4227 = 0.4%), with zero golimumab and 2 non-biologics also identified with a claim for ICD-9 714.0 during baseline. The small numbers and good balance obtained on this variable means that these patients are unlikely to have any influence on study results.

Several types of analyses were undertaken in order to understand the association between golimumab use and the prespecified study outcomes. Based on the exposure definition used in this study, determining the temporal relationship between golimumab exposure and these outcomes, particularly chronic outcomes (i.e. malignancy and lymphoma) is difficult. Protopathic bias, a form of reverse causation in which treatment is initiated in response to symptoms of undiagnosed disease, cannot be ruled out. Protopathic symptoms of study outcomes may be mistaken for rheumatic disease activity leading to initiation of golimumab or NBS. The resulting association between the study drug and outcome could, in fact, be due to reverse causation, in which the underlying, undiagnosed disease leads to use of the study drug [39].

The current study design did not specifically account for the possibility of golimumab to be used in a manner that might produce protopathic bias. In order to address potential protopathic bias in the analysis of malignancy and lymphoma, a sensitivity analysis was conducted where person-time in the as-matched analysis that is particularly susceptible to protopathic bias (the first 6 months of exposure time) was removed, and the effect estimates recalculated [40]. Results from this sensitivity analysis were consistent with the primary findings, suggesting no strong protopathic bias effect (data not shown).

An elevated rate of depression during current golimumab use was observed in the as-treated analysis, but not in the as-matched or NCC analyses. Sensitivity analyses indicated depression is sensitive to misclassification, so this result should be interpreted cautiously. Furthermore, as patients with a baseline medical history of depression were excluded from the NCC and sensitivity analyses, the increased risk seen in the as-treated analysis may be related to an imbalance in baseline history of depression. Due to the design of the study, which focused on baseline balance of patient characteristics, the effect of changes in severity or treatment for depression after initiation of golimumab were not analyzed.

Assessment of depression in insurance claims databases presents challenges (e.g. patients underreporting symptoms due to stigma of mental illness, failure of physicians to recognize or treat depression). The definition of depression used here was comprehensive (i.e. any inpatient or outpatient claim with a diagnosis code in the primary position or a dispensing for an antidepressant) and captured more patients at the expense of generating false positive diagnoses. When a more stringent definition of depression that required a hospitalized diagnosis in the primary position was used, 3 cases were identified in the matched golimumab cohort among 1969 person-years (crude IR = 1.52 per 1000 person-years) and 19 in the NBS cohort among 6254 person-years (crude IR = 3.04 per 1000 person-years).

In chronic inflammatory diseases, comorbid depression is common. Among RA patients, prevalence of depression ranges from 14–48% [41-43] and is associated with higher levels of pain and disability, lower health-related quality of life, and increased mortality [44]. Evidence suggests a role of pro-inflammatory cytokines including interleukin (IL)-6, TNF-alpha, and IL-1 in promoting depression, and indeed, patients with severe RA have a higher risk of depression [45-49]. Specifically, TNF-alpha has been shown to be associated with depression and reducing its levels may actually reverse symptoms. In a randomized clinical trial investigating the effect of etanercept on fatigue and depression in patients with plaque psoriasis, improvements in depression were noted [50]. Smaller studies in patients with RA and inflammatory bowel disease confirm that therapy with TNFi could have a beneficial role in the treatment of depressive symptoms in patients with inflammatory conditions [51-53].

Conclusion

In summary, this study involving active surveillance of prespecified outcomes identified no specific safety concerns in patients with rheumatic disease treated with golimumab versus a NBS medication. The results are consistent with golimumab’s overall safety profile and generally comparable with observations from other studies in patients with treated rheumatic disease.

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 132 kb)

This study suggests that there is no association between exposure to golimumab and an increased risk of prespecified outcomes.

The results of this study are consistent with golimumab’s overall safety profile and generally comparable with observations from other studies in patients with treated rheumatic disease.