Ambulatory blood pressure in patients with rheumatoid arthritis: association with immune activation

OBJECTIVE

The prevalence of hypertension, a major cardiovascular risk factor, is increased in patients with rheumatoid arthritis (RA) and may be driven by immune activation. The purpose of this study was to determine if ambulatory 24-hour blood pressure (BP) is elevated in RA versus control participants and associated with immune activation.

METHODS

We conducted a cross-sectional study of 46 patients with RA and 23 control participants. Participants wore an ambulatory BP monitor which obtained diurnal BP every 15-30 minutes and nocturnal BP every 30 minutes. Inflammatory mediators in plasma were measured using an inflammation proteomics panel. Differences in BP measurements were assessed by Mann-Whitney U, and association with inflammatory mediators was assessed by Spearman correlation.

RESULTS

Patients with RA and control participants had similar office BP, but ambulatory systolic BP measurements: 24-hour (median RA =121 mmHg versus control =116 mmHg, P=0.007), diurnal (RA =128 mmHg versus control =120 mmHg, P=0.003), and nocturnal (RA =112 mmHg versus control =103 mmHg, P=0.002) were higher in patients with RA. Patients with RA also had higher nocturnal diastolic BP (RA =63 mmHg versus control =57 mmHg, P=0.02), but other diastolic BP measurements were similar. Nocturnal BP dipping was less in RA (12%) compared to control participants (16%, P=0.02). In RA, higher 24-hour and nocturnal systolic BPs and less nocturnal dipping were strongly correlated with a wide range of inflammatory mediators.

CONCLUSION

Despite similar office measurements, 24-hour and nocturnal systolic BP measurements were higher in RA than control participants and were strongly associated with inflammation.

Serious Bleeding in Patients With Atrial Fibrillation Using Diltiazem With Apixaban or Rivaroxaban

Importance

Diltiazem, a commonly prescribed ventricular rate-control medication for patients with atrial fibrillation, inhibits apixaban and rivaroxaban elimination, possibly causing overanticoagulation.

Objective

To compare serious bleeding risk for new users of apixaban or rivaroxaban with atrial fibrillation treated with diltiazem or metoprolol.

Design, Setting, and Participants

This retrospective cohort study included Medicare beneficiaries aged 65 years or older with atrial fibrillation who initiated apixaban or rivaroxaban use and also began treatment with diltiazem or metoprolol between January 1, 2012, and November 29, 2020. Patients were followed up to 365 days through November 30, 2020. Data were analyzed from August 2023 to February 2024.

Exposures

Diltiazem and metoprolol.

Main Outcomes and Measures

The primary outcome was a composite of bleeding-related hospitalization and death with recent evidence of bleeding. Secondary outcomes were ischemic stroke or systemic embolism, major ischemic or hemorrhagic events (ischemic stroke, systemic embolism, intracranial or fatal extracranial bleeding, or death with recent evidence of bleeding), and death without recent evidence of bleeding. Hazard ratios (HRs) and rate differences (RDs) were adjusted for covariate differences with overlap weighting.

Results

The study included 204 155 US Medicare beneficiaries, of whom 53 275 received diltiazem and 150 880 received metoprolol. Study patients (mean [SD] age, 76.9 [7.0] years; 52.7% female) had 90 927 person-years (PY) of follow-up (median, 120 [IQR, 59-281] days). Patients receiving diltiazem treatment had increased risk for the primary outcome (RD, 10.6 [95% CI, 7.0-14.2] per 1000 PY; HR, 1.21 [95% CI, 1.13-1.29]) and its components of bleeding-related hospitalization (RD, 8.2 [95% CI, 5.1-11.4] per 1000 PY; HR, 1.22 [95% CI, 1.13-1.31]) and death with recent evidence of bleeding (RD, 2.4 [95% CI, 0.6-4.2] per 1000 PY; HR, 1.19 [95% CI, 1.05-1.34]) compared with patients receiving metoprolol. Risk for the primary outcome with initial diltiazem doses exceeding 120 mg/d (RD, 15.1 [95% CI, 10.2-20.1] per 1000 PY; HR, 1.29 [95% CI, 1.19-1.39]) was greater than that for lower doses (RD, 6.7 [95% CI, 2.0-11.4] per 1000 PY; HR, 1.13 [95% CI, 1.04-1.24]). For doses exceeding 120 mg/d, the risk of major ischemic or hemorrhagic events was increased (HR, 1.14 [95% CI, 1.02-1.27]). Neither dose group had significant changes in the risk for ischemic stroke or systemic embolism or death without recent evidence of bleeding. When patients receiving high- and low-dose diltiazem treatment were directly compared, the HR for the primary outcome was 1.14 (95% CI, 1.02-1.26).

Conclusions and Relevance

In Medicare patients with atrial fibrillation receiving apixaban or rivaroxaban, diltiazem was associated with greater risk of serious bleeding than metoprolol, particularly for diltiazem doses exceeding 120 mg/d.

Clinical associations with a polygenic predisposition to benign lower white blood cell counts

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30-0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69-0.88], p = 4.0 × 10-5) or immunosuppressant (n = 354, HR = 0.61 [0.38-0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44-0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease

SIGNIFICANCE STATEMENT

African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease.

BACKGROUND

African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene.

METHODS

We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models.

RESULTS

In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants.

CONCLUSIONS

APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

Concurrent Gabapentin and Opioid Use and Risk of Mortality in Medicare Recipients with Non-Cancer Pain

Gabapentin is prescribed for pain and is perceived as safe generally. However, gabapentin can cause respiratory depression, exacerbated by concomitant central nervous system depressants (e.g., opioids), a concern for vulnerable populations. We compared mortality rates among new users of either gabapentin or duloxetine with or without concurrent opioids in the 20% Medicare sample. We conducted a new-user design retrospective cohort study, in Medicare enrollees ages 65-89 years with noncancer chronic pain and no severe illness who filled prescriptions between 2015 and 2018 for gabapentin (n = 233,060) or duloxetine (n = 34,009). Daily opioid doses, estimated in morphine milligram equivalents (MMEs), were classified into none, low (0 < MME < 50), and high (≥ 50 MME), based on Centers for Disease Control and Prevention (CDC) recommendations. The outcomes were all-cause mortality (primary) and out-of-hospital mortality (secondary). We used inverse probability of treatment weighting to adjust for differences between gabapentin and duloxetine users. During 116,707 person-years of follow-up, 1,379 patients died. All-cause mortality rate in gabapentin users was 12.16 per 1,000 person-years vs. 9.94 per 1,000 in duloxetine users. Risks were similar for users with no concurrent opioids (adjusted hazard ratio (aHR) = 1.03, 95% confidence interval (CI): 0.80-1.31) or low-dose daily opioids (aHR = 1.06, 95% CI: 0.63-1.76). However, gabapentin users receiving concurrent high-dose daily opioids had an increased rate of all-cause mortality compared with duloxetine users on high-dose opioids (aHR = 2.03, 95% CI: 1.19-3.46). Out-of-hospital mortality yielded similar results. In this retrospective cohort study of Medicare beneficiaries, concurrent use of high-dose opioids and gabapentin was associated with a higher all-cause mortality risk than that for concurrent use of high-dose opioids and duloxetine.

Association between APOL1 risk variants and the occurrence of sepsis in Black patients hospitalized with infections: a retrospective cohort study

Background

Two risk variants in the apolipoprotein L1 gene (APOL1) have been associated with increased susceptibility to sepsis in Black patients. However, it remains unclear whether APOL1 high-risk genotypes are associated with occurrence of either sepsis or sepsis-related phenotypes in patients hospitalized with infections, independent of their association with pre-existing severe renal disease.

Methods

A retrospective cohort study of 2242 Black patients hospitalized with infections. We assessed whether carriage of APOL1 high-risk genotypes was associated with the risk of sepsis and sepsis-related phenotypes in patients hospitalized with infections. The primary outcome was sepsis; secondary outcomes were short-term mortality, and organ failure related to sepsis.

Results

Of 2242 Black patients hospitalized with infections, 565 developed sepsis. Patients with high-risk APOL1 genotypes had a significantly increased risk of sepsis (odds ratio [OR]=1.29 [95% CI, 1.00-1.67; p=0.047]); however, this association was not significant after adjustment for pre-existing severe renal disease (OR = 1.14 [95% CI, 0.88-1.48; p=0.33]), nor after exclusion of those patients with pre-existing severe renal disease (OR = 0.99 [95% CI, 0.70-1.39; p=0.95]). APOL1 high-risk genotypes were significantly associated with the renal dysfunction component of the Sepsis-3 criteria (OR = 1.64 [95% CI, 1.21-2.22; p=0.001]), but not with other sepsis-related organ dysfunction or short-term mortality. The association between high-risk APOL1 genotypes and sepsis-related renal dysfunction was markedly attenuated by adjusting for pre-existing severe renal disease (OR = 1.36 [95% CI, 1.00-1.86; p=0.05]) and was nullified after exclusion of patients with pre-existing severe renal disease (OR = 1.16 [95% CI, 0.74-1.81; p=0.52]).

Conclusions

APOL1 high-risk genotypes were associated with an increased risk of sepsis; however, this increased risk was attributable predominantly to pre-existing severe renal disease.

Funding

This study was supported by R01GM120523 (QF), R01HL163854 (QF), R35GM131770 (CMS), HL133786 (WQW), and Vanderbilt Faculty Research Scholar Fund (QF). The dataset(s) used for the analyses described were obtained from Vanderbilt University Medical Center's BioVU which is supported by institutional funding, the 1S10RR025141-01 instrumentation award, and by the CTSA grant UL1TR0004from NCATS/NIH. Additional funding provided by the NIH through grants P50GM115305 and U19HL065962. The authors wish to acknowledge the expert technical support of the VANTAGE and VANGARD core facilities, supported in part by the Vanderbilt-Ingram Cancer Center (P30 CA068485) and Vanderbilt Vision Center (P30 EY08126). The funders had no role in design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Association between APOL1 risk variants and the occurrence of sepsis in Black patients hospitalized with infections: a retrospective cohort study

Background

Two risk variants in the apolipoprotein L1 gene ( APOL1 ) have been associated with increased susceptibility to sepsis in Black patients. However, it remains unclear whether APOL1 high-risk genotypes are associated with occurrence of either sepsis or sepsis-related phenotypes in patients hospitalized with infections, independent of their association with pre-existing severe renal disease.

Methods

A retrospective cohort study of 2,242 Black patients hospitalized with infections. We assessed whether carriage of APOL1 high-risk genotypes was associated with the risk of sepsis and sepsis-related phenotypes in patients hospitalized with infections. The primary outcome was sepsis; secondary outcomes were short-term mortality and organ failure related to sepsis.

Results

Of 2,242 Black patients hospitalized with infections, 565 developed sepsis. Patients with high-risk APOL1 genotypes had a significantly increased risk of sepsis (odds ratio [OR]=1.29 [95% CI, 1.00-1.67; p=0.047]); however, this association was not significant after adjustment for pre-existing severe renal disease (OR=1.14 [95% CI, 0.88-1.48; p=0.33]), nor after exclusion of those patients with pre-existing severe renal disease (OR=0.99 [95% CI, 0.70-1.39; p=0.95]. APOL1 high-risk genotypes were significantly associated with the renal dysfunction component of the Sepsis-3 criteria (OR=1.64 [95% CI, 1.21-2.22; p=0.001], but not with other sepsis-related organ dysfunction or short-term mortality. The association between high-risk APOL1 genotypes and sepsis-related renal dysfunction was markedly attenuated by adjusting for pre-existing severe renal disease (OR=1.36 [95% CI, 1.00-1.86; p=0.05]) and was nullified after exclusion of patients with pre-existing severe renal disease (OR=1.16 [95% CI, 0.74-1.81; p=0.52]).

Conclusion

APOL1 high-risk genotypes were associated with an increased risk of sepsis; however, this increased risk was attributable predominantly to pre-existing severe renal disease.

Funding

This study was supported by R01GM120523 (Q.F.), R01HL163854 (Q.F.), R35GM131770 (C.M.S.), HL133786 (W.Q.W.), and Vanderbilt Faculty Research Scholar Fund (Q.F.). The dataset(s) used for the analyses described were obtained from Vanderbilt University Medical Center's BioVU which is supported by institutional funding, the 1S10RR025141-01 instrumentation award, and by the CTSA grant UL1TR0004from NCATS/NIH. Additional funding provided by the NIH through grants P50GM115305 and U19HL065962. The authors wish to acknowledge the expert technical support of the VANTAGE and VANGARD core facilities, supported in part by the Vanderbilt-Ingram Cancer Center (P30 CA068485) and Vanderbilt Vision Center (P30 EY08126).The funders had no role in design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Genome-Wide Association Study of CKD Progression

SIGNIFICANCE STATEMENT

Rapid progression of CKD is associated with poor clinical outcomes. Most previous studies looking for genetic factors associated with low eGFR have used cross-sectional data. The authors conducted a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD, focusing on longitudinal data. They identified three loci (two of them novel) associated with longitudinal eGFR decline. In addition to the known UMOD/PDILT locus, variants within BICC1 were associated with significant differences in longitudinal eGFR slope. Variants within HEATR4 also were associated with differences in eGFR decline, but only among Black/African American individuals without diabetes. These findings help characterize molecular mechanisms of eGFR decline in CKD and may inform new therapeutic approaches for progressive kidney disease.

BACKGROUND

Rapid progression of CKD is associated with poor clinical outcomes. Despite extensive study of the genetics of cross-sectional eGFR, only a few loci associated with eGFR decline over time have been identified.

METHODS

We performed a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD-defined by two outpatient eGFR measurements of <60 ml/min per 1.73 m 2 , obtained 90-365 days apart-from the Million Veteran Program and Vanderbilt University Medical Center's DNA biobank. The primary outcome was the annualized relative slope in outpatient eGFR. Analyses were stratified by ethnicity and diabetes status and meta-analyzed thereafter.

RESULTS

In cross-ancestry meta-analysis, the strongest association was rs77924615, near UMOD / PDILT ; each copy of the G allele was associated with a 0.30%/yr faster eGFR decline ( P = 4.9×10 -27 ). We also observed an association within BICC1 (rs11592748), where every additional minor allele was associated with a 0.13%/yr slower eGFR decline ( P = 5.6×10 -9 ). Among participants without diabetes, the strongest association was the UMOD/PDILT variant rs36060036, associated with a 0.27%/yr faster eGFR decline per copy of the C allele ( P = 1.9×10 -17 ). Among Black participants, a significantly faster eGFR decline was associated with variant rs16996674 near APOL1 (R 2 =0.29 with the G1 high-risk genotype); among Black participants with diabetes, lead variant rs11624911 near HEATR4 also was associated with a significantly faster eGFR decline. We also nominally replicated loci with known associations with eGFR decline, near PRKAG2, FGF5, and C15ORF54.

CONCLUSIONS

Three loci were significantly associated with longitudinal eGFR change at genome-wide significance. These findings help characterize molecular mechanisms of eGFR decline and may contribute to the development of new therapeutic approaches for progressive CKD.

Clinical consequences of a polygenic predisposition to benign lower white blood cell counts: Consequences of benign WBC count genetics

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is undefined. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio=0.55 per standard deviation increase in PGSWBC [95%CI, 0.30 - 0.94], p=0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n=1,724, hazard ratio [HR]=0.78 [0.69 - 0.88], p=4.0×10-5) or immunosuppressant (n=354, HR=0.61 [0.38 - 0.99], p=0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n=1,466, HR=0.62 [0.44 - 0.87], p=0.006). Collectively, these findings suggest that a WBC count polygenic score identifies individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Clinical diagnoses associated with a positive antinuclear antibody test in patients with and without autoimmune disease

BACKGROUND

Antinuclear antibodies (ANA) are antibodies present in several autoimmune disorders. However, a large proportion of the general population (20%) also have a positive test; very few of these individuals will develop an autoimmune disease, and the clinical impact of a positive ANA in them is not known. Thus, we test the hypothesis that ANA + test reflects a state of immune dysregulation that alters risk for some clinical disorders in individuals without an autoimmune disease.

METHODS

We performed high throughput association analyses in a case-control study using real world data from the de-identified electronic health record (EHR) system from Vanderbilt University Medical Center. The study population included individuals with an ANA titer ≥ 1:80 at any time (ANA +) and those with negative results (ANA-). The cohort was stratified into sub-cohorts of individuals with and without an autoimmune disease. A phenome-wide association study (PheWAS) adjusted by sex, year of birth, race, and length of follow-up was performed in the study cohort and in the sub-cohorts. As secondary analyses, only clinical diagnoses after ANA testing were included in the analyses.

RESULTS

The cohort included 70,043 individuals: 49,546 without and 20,497 with an autoimmune disease, 26,579 were ANA + and 43,464 ANA-. In the study cohort and the sub-cohort with autoimmune disease, ANA + was associated (P ≤ 5 × 10-5) with 88 and 136 clinical diagnoses respectively, including lupus (OR ≥ 5.4, P ≤ 7.8 × 10-202) and other autoimmune diseases and complications. In the sub-cohort without autoimmune diseases, ANA + was associated with increased risk of Raynaud's syndrome (OR ≥ 2.1) and alveolar/perialveolar-related pneumopathies (OR ≥ 1.4) and decreased risk of hepatitis C, tobacco use disorders, mood disorders, convulsions, fever of unknown origin, and substance abuse disorders (OR ≤ 0.8). Analyses including only diagnoses after ANA testing yielded similar results.

CONCLUSION

A positive ANA test, in addition to known associations with autoimmune diseases, Raynaud's phenomenon, and idiopathic fibrosing alveolitis related disorders, is associated with decreased prevalence of several non-autoimmune diseases.

Risk for Bleeding-Related Hospitalizations During Use of Amiodarone With Apixaban or Rivaroxaban in Patients With Atrial Fibrillation : A Retrospective Cohort Study

BACKGROUND

Amiodarone, the most effective antiarrhythmic drug in atrial fibrillation, inhibits apixaban and rivaroxaban elimination, thus possibly increasing anticoagulant-related risk for bleeding.

OBJECTIVE

For patients receiving apixaban or rivaroxaban, to compare risk for bleeding-related hospitalizations during treatment with amiodarone versus flecainide or sotalol, antiarrhythmic drugs that do not inhibit these anticoagulants' elimination.

DESIGN

Retrospective cohort study.

SETTING

U.S. Medicare beneficiaries aged 65 years or older.

PATIENTS

Patients with atrial fibrillation began anticoagulant use between 1 January 2012 and 30 November 2018 and subsequently initiated treatment with study antiarrhythmic drugs.

MEASUREMENTS

Time to event for bleeding-related hospitalizations (primary outcome) and ischemic stroke, systemic embolism, and death with or without recent (past 30 days) evidence of bleeding (secondary outcomes), adjusted with propensity score overlap weighting.

RESULTS

There were 91 590 patients (mean age, 76.3 years; 52.5% female) initiating use of study anticoagulants and antiarrhythmic drugs, 54 977 with amiodarone and 36 613 with flecainide or sotalol. Risk for bleeding-related hospitalizations increased with amiodarone use (rate difference [RD], 17.5 events [95% CI, 12.0 to 23.0 events] per 1000 person-years; hazard ratio [HR], 1.44 [CI, 1.27 to 1.63]). Incidence of ischemic stroke or systemic embolism did not increase (RD, -2.1 events [CI, -4.7 to 0.4 events] per 1000 person-years; HR, 0.80 [CI, 0.62 to 1.03]). The risk for death with recent evidence of bleeding (RD, 9.1 events [CI, 5.8 to 12.3 events] per 1000 person-years; HR, 1.66 [CI, 1.35 to 2.03]) was greater than that for other deaths (RD, 5.6 events [CI, 0.5 to 10.6 events] per 1000 person-years; HR, 1.15 [CI, 1.00 to 1.31]) (HR comparison: P = 0.003). The increased incidence of bleeding-related hospitalizations for rivaroxaban (RD, 28.0 events [CI, 18.4 to 37.6 events] per 1000 person-years) was greater than that for apixaban (RD, 9.1 events [CI, 2.8 to 15.3 events] per 1000 person-years) (P = 0.001).

LIMITATION

Possible residual confounding.

CONCLUSION

In this retrospective cohort study, patients aged 65 years or older with atrial fibrillation treated with amiodarone during apixaban or rivaroxaban use had greater risk for bleeding-related hospitalizations than those treated with flecainide or sotalol.

PRIMARY FUNDING SOURCE

National Heart, Lung, and Blood Institute.

Duloxetine, Gabapentin, and the Risk for Acute Myocardial Infarction, Stroke, and Out-of-Hospital Death in Medicare Beneficiaries With Non-Cancer Pain

OBJECTIVE

Duloxetine is a serotonin-norepinephrine reuptake inhibitor prescribed for musculoskeletal and other forms of chronic pain. Its dual pharmacologic properties have the potential to either raise or lower cardiovascular risk: adrenergic activity may increase the risk for acute myocardial infarction (AMI) and stroke, but antiplatelet activity may decrease risk. Gabapentin is another nonopioid medication used to treat pain, which is not thought to have adrenergic/antiplatelet effects. With the current emphasis on the use of nonopioid medications to treat patients with chronic pain, assessing cardiovascular risks associated with these medications among high-risk patients is important.

MATERIALS AND METHODS

We conducted a retrospective cohort study among a 20% sample of Medicare enrollees, aged 65 to 89, with chronic pain who were new users between 2015 and 2018 of either duloxetine (n = 34,009) or gabapentin (n = 233,060). We excluded individuals with cancer or other life-threatening conditions at study drug initiation. The primary outcome was a composite of AMI, stroke, and out-of-hospital mortality. We adjusted for comorbidity differences with time-dependent inverse probability of treatment weighting.

RESULTS

During 115,668 person-years of follow-up, 2361 patients had the composite primary outcome; the rate among new users of duloxetine was 16.7/1000 person-years compared with new users of gabapentin (21.1/1000 person-years), adjusted hazard ratio = 0.98 (95% CI: 0.83, 1.16). Results were similar for the individual components of the composite outcome as well as in analyses stratified by demographic and clinical characteristics.

DISCUSSION

In summary, cohort Medicare patients with non-cancer pain beginning treatment with duloxetine had rates of AMI, stroke, and out-of-hospital mortality comparable to those who initiated gabapentin.

APOL1 and the risk of adverse renal outcomes in patients of African ancestry with systemic lupus erythematosus

BACKGROUND

Systemic lupus erythematosus (SLE) disproportionately affects individuals of African ancestry (AA) compared to European ancestry (EA). In the general population, high risk (HR) variants in the apolipoprotein L1 (APOL1) gene increase the risk of renal and hypertensive disorders in individuals of AA. Since SLE is characterized by an interferon signature and APOL1 expression is driven by interferon, we examined the hypothesis that APOL1 HR genotypes predominantly drive higher rates of renal and hypertensive-related comorbidities observed in SLE patients of AA versus those of EA.

METHODS

We performed a retrospective cohort study in patients with SLE of EA and AA using a genetic biobank linked to de-identified electronic health records. APOL1 HR genotypes were defined as G1/G1, G2/G2, or G1/G2 and low risk (LR) genotypes as 1 or 0 copies of the G1 and G2 alleles. To identify renal and hypertensive-related disorders that differed in prevalence by ancestry, we used a phenome-wide association approach. We then used logistic regression to compare the prevalence of renal and hypertensive-related disorders in EA and AA patients, both including and excluding patients with the APOL1 HR genotype. In a sensitivity analysis, we examined the association of end stage renal disease secondary to lupus nephritis (LN-related ESRD) with ancestry and the APOL1 genotype.

RESULTS

We studied 784 patients with SLE; 195 (24.9%) were of AA, of whom 27 (13.8%) had APOL1 HR genotypes. Eighteen renal and hypertensive-related phenotypes were more common in AA than EA patients (p-value ≤ 1.4E-4). All phenotypes remained significantly different after exclusion of patients with APOL1 HR genotypes, and most point odds ratios (ORs) decreased only slightly. Even among ORs with the greatest decrease, risk for AA patients without the APOL1 HR genotype remained significantly elevated compared to EA patients. In the sensitivity analysis, LN-related ESRD was more prevalent in SLE patients of AA versus EA and AA patients with the APOL1 HR genotype versus LR (p-value < .05 for both).

CONCLUSION

The higher prevalence of renal and hypertensive disorders in SLE patients of AA compared to those of EA is not fully explained by the presence of APOL1 high risk variants.

The relationship between high density lipoprotein cholesterol and sepsis: A clinical and genetic approach

Sepsis accounts for one in three hospital deaths. Higher concentrations of high-density lipoprotein cholesterol (HDL-C) are associated with apparent protection from sepsis, suggesting a potential therapeutic role for HDL-C or drugs, such as cholesteryl ester transport protein (CETP) inhibitors that increase HDL-C. However, these beneficial clinical associations might be due to confounding; genetic approaches can address this possibility. We identified 73,406 White adults admitted to Vanderbilt University Medical Center with infection; 11,612 had HDL-C levels, and 12,377 had genotype information from which we constructed polygenic risk scores (PRS) for HDL-C and the effect of CETP on HDL-C. We tested the associations between predictors (measured HDL-C, HDL-C PRS, CETP PRS, and rs1800777) and outcomes: sepsis, septic shock, respiratory failure, and in-hospital death. In unadjusted analyses, lower measured HDL-C concentrations were significantly associated with increased risk of sepsis (p = 2.4 × 10-23 ), septic shock (p = 4.1 × 10-12 ), respiratory failure (p = 2.8 × 10-8 ), and in-hospital death (p = 1.0 × 10-8 ). After adjustment (age, sex, electronic health record length, comorbidity score, LDL-C, triglycerides, and body mass index), these associations were markedly attenuated: sepsis (p = 2.6 × 10-3 ), septic shock (p = 8.1 × 10-3 ), respiratory failure (p = 0.11), and in-hospital death (p = 4.5 × 10-3 ). HDL-C PRS, CETP PRS, and rs1800777 significantly predicted HDL-C (p < 2 × 10-16 ), but none were associated with sepsis outcomes. Concordant findings were observed in 13,254 Black patients hospitalized with infections. Lower measured HDL-C levels were significantly associated with increased risk of sepsis and related outcomes in patients with infection, but a causal relationship is unlikely because no association was found between the HDL-C PRS or the CETP PRS and the risk of adverse sepsis outcomes.

Genetic Determinants of IL-6 Levels and Risk of ESKD

Genetically predicted IL-6 levels are associated with risk of ESKD. Therapeutic modulation of IL-6 could potentially reduce the risk of ESKD.

Association Between Genetically Predicted Expression of TPMT and Azathioprine Adverse Events

Polymorphisms thiopurine-S-methyltransferase (TPMT) and nudix hydrolase 15 (NUDT15) can increase the risk of azathioprine myelotoxicity, but little is known about other genetic factors that increase risk for azathioprine-associated side effects. PrediXcan is a gene-based association method that estimates the expression of individuals' genes and examines their correlation to specified phenotypes. As proof of concept for using PrediXcan as a tool to define the association between genetic factors and azathioprine side effects, we aimed to determine whether the genetically predicted expression of TPMT or NUDT15 was associated with leukopenia or other known side effects. In a retrospective cohort of 1364 new users of azathioprine with EHR-reported White race, we used PrediXcan to impute expression in liver tissue, tested its association with pre-specified phecodes representing known side effects (e.g., skin cancer), and completed chart review to confirm cases. Among confirmed cases, patients in the lowest tertile (i.e., lowest predicted) of TPMT expression had significantly higher odds of developing leukopenia (OR=3.30, 95%CI: 1.07-10.20, p=0.04) versus those in the highest tertile; no other side effects were significant. The results suggest that this methodology could be deployed on a larger scale to uncover associations between genetic factors and drug side effects for more personalized care.

Intermuscular adipose tissue in patients with systemic lupus erythematosus

OBJECTIVE

Patients with SLE frequently have debilitating fatigue and reduced physical activity. Intermuscular adipose tissue (IMAT) accumulation is associated with reduced physical exercise capacity. We hypothesised that IMAT is increased in patients with SLE and associated with increased fatigue, reduced physical activity and increased inflammation.

METHODS

In a cross-sectional study, 23 patients with SLE and 28 control participants were evaluated. IMAT was measured in the calf muscles using sequential T ₁-weighted MRI. Patient-reported physical activity and fatigue were measured and a multiplex proteomic assay was used to measure markers and mediators of inflammation.

RESULTS

IMAT accumulation (percentage of IMAT area to muscle area) was significantly higher in SLE versus control participants (7.92%, 4.51%-13.39% vs 2.65%, 1.15%-4.61%, median, IQR, p<0.001) and remained significant after adjustment for age, sex, race and body mass index (p<0.001). In patients with SLE, IMAT accumulation did not differ significantly among corticosteroid users and non-users (p=0.48). In the study cohort (patients and controls), IMAT was positively correlated with self-reported fatigue score (rho=0.52, p<0.001) and inversely correlated with self-reported walking distance (rho=-0.60, p<0.001). Several markers of inflammation were associated with IMAT accumulation in patients with SLE, and gene ontology analysis showed significant enrichment for pathways associated with macrophage migration and activation in relation to IMAT.

CONCLUSION

Patients with SLE have greater IMAT accumulation than controls in the calf muscles. Increased IMAT is associated with greater fatigue and lower physical activity. Future studies should evaluate the effectiveness of interventions that improve muscle quality to alleviate fatigue in patients with SLE.

Ancestry, ACKR1 and leucopenia in patients with systemic lupus erythematosus

OBJECTIVE

SLE is more prevalent in populations of African (AA) than European ancestry (EA) and leucopenia is common. A homozygous variant in ACKR1 (rs2814778-CC) is associated with lower white cell counts; the variant is common in AA but not EA populations. We hypothesised that in SLE: (1) leucopenia is more frequent in patients of AA than EA, and (2) the ACKR1-CC genotype accounts for the higher frequency of leucopenia in AA patients.

METHODS

We performed a retrospective cohort study in patients with SLE at a tertiary care system. Ancestry was defined by genetic principal components. We compared the rate of leucopenia, thrombocytopenia and anaemia between (a) EA and AA patients, and (b) ACKR1-CT/TT and CC genotype in AA patients.

RESULTS

The cohort included 574 patients of EA and 190 of AA; ACKR1-CC genotype was common in AA (70%) but not EA (0%) patients. Rates of leucopenia for ancestry and genotype were AA 60.0% vs EA 36.8 % (p=1.9E-08); CC 67.7% vs CT/TT 42.1% (p=9.8E-04). The rate of leucopenia did not differ by ancestry comparing EA patients versus AA with CT/TT genotype (p=0.59). Thrombocytopenia (22.2% vs 13.2%, p=0.004) and anaemia (88.4% vs 66.2%, p=3.7E-09) were more frequent in AA patients but were not associated with ACKR1 genotype (p=0.82 and p=0.84, respectively).

CONCLUSIONS

SLE of AA had higher rates of anaemia, leucopenia, and thrombocytopenia than those of EA; only the difference in leucopenia was explained by ACKR1-CC genotype. This genotype could affect clinical practice.

Genome-wide association analyses of common infections in a large practice-based biobank

INTRODUCTION

Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections.

METHODS

We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes.

RESULTS

We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10^-3) and rs2808290-C (p = 9.6 × 10^-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10^-8): rs113235453-G for otitis media (p = 3.04 × 10^-8), and rs10422015-T for candidiasis (p = 3.11 × 10^-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10^-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10^-7 and p = 1.95 × 10^-6); and AAMDC for hepatitis B (p = 1.51 × 10^-6).

CONCLUSION

We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

Race, Genotype, and Azathioprine Discontinuation : A Cohort Study

BACKGROUND

Thiopurines are an important class of immunosuppressants despite their risk for hematopoietic toxicity and narrow therapeutic indices. Benign neutropenia related to an ACKR1 variant (rs2814778-CC) is common among persons of African ancestries.

OBJECTIVE

To test whether rs2814778-CC was associated with azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine dosing.

DESIGN

Retrospective cohort study.

SETTING

Two tertiary care centers.

PATIENTS

Thiopurine users with White or Black race.

MEASUREMENTS

Azathioprine discontinuation attributed to hematopoietic toxicity. Secondary outcomes included weight-adjusted final dose, leukocyte count, and change in leukocyte count.

RESULTS

The rate of azathioprine discontinuation attributed to hematopoietic toxicity was 3.92 per 100 person-years among patients with the CC genotype (n = 101) and 1.34 per 100 person-years among those with the TT or TC genotype (n = 1365) (hazard ratio [HR] from competing-risk model, 2.92 [95% CI, 1.57 to 5.41]). The risk remained significant after adjustment for race (HR, 2.61 [CI, 1.01 to 6.71]). The risk associated with race alone (HR, 2.13 [CI, 1.21 to 3.75]) was abrogated by adjustment for genotype (HR, 1.13 [CI, 0.48 to 2.69]). Lower last leukocyte count and lower dosing were significant among patients with the CC genotype. Lower dosing was validated in an external cohort of 94 children of African ancestries prescribed the thiopurine 6-mercaptopurine (6-MP) for acute lymphoblastic leukemia. The CC genotype was independently associated with lower 6-MP dose intensity relative to the target daily dose of 75 mg/m2 (median, 0.83 [IQR, 0.70 to 0.94] for the CC genotype vs. 0.94 [IQR, 0.72 to 1.13] for the TT or TC genotype; P = 0.013).

LIMITATIONS

Unmeasured confounding; data limited to tertiary centers.

CONCLUSION

Patients with the CC genotype had higher risk for azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine doses. Genotype was associated with those risks, even after adjustment for race.

PRIMARY FUNDING SOURCE

National Institutes of Health.

Phenome-Wide Association Study of UMOD Gene Variants and Differential Associations With Clinical Outcomes Across Populations in the Million Veteran Program a Multiethnic Biobank

Introduction

Common variants in the UMOD gene are considered an evolutionary adaptation against urinary tract infections (UTIs) and have been implicated in kidney stone formation, chronic kidney disease (CKD), and hypertension. However, differences in UMOD variant-phenotype associations across population groups are unclear.

Methods

We tested associations between UMOD/PDILT variants and up to 1528 clinical diagnosis codes mapped to phenotype groups in the Million Veteran Program (MVP), using published phenome-wide association study (PheWAS) methodology. Associations were tested using logistic regression adjusted for age, sex, and 10 principal components of ancestry. Bonferroni correction for multiple comparisons was applied.

Results

Among 648,593 veterans, mean (SD) age was 62 (14) years; 9% were female, 19% Black, and 8% Hispanic. In White patients, the rs4293393 UMOD risk variant associated with increased uromodulin was associated with increased odds of CKD (odds ratio [OR]: 1.22, 95% CI: 1.20-1.24, P = 5.90 × 10-111), end-stage kidney disease (OR: 1.17, 95% CI: 1.11-1.24, P = 2.40 × 10-09), and hypertension (OR: 1.03, 95% CI: 1.05-1.05, P = 2.11 × 10-06) and significantly lower odds of UTIs (OR: 0.94, 95% CI: 0.92-0.96, P = 1.21 × 10-10) and kidney calculus (OR: 0.85, 95% CI: 0.83-0.86, P = 4.27 × 10-69). Similar findings were observed across UMOD/PDILT variants. The rs77924615 PDILT variant had stronger associations with acute cystitis in White female (OR: 0.73, 95% CI: 0.59-0.91, P = 4.98 × 10-03) versus male (OR: 0.99, 95% CI: 0.89-1.11, P = 8.80 × 10-01) (P interaction = 0.01) patients. In Black patients, the rs77924615 PDILT variant was significantly associated with pyelonephritis (OR: 0.65, 95% CI: 0.54-0.79, P = 1.05 × 10-05), whereas associations with UMOD promoter variants were attenuated.

Conclusion

Robust associations were observed between UMOD/PDILT variants linked with increased uromodulin expression and lower odds of UTIs and calculus and increased odds of CKD and hypertension. However, these associations varied significantly across ancestry groups and sex.

Pharmacogenetic variants and risk of remdesivir-associated liver enzyme elevations in Million Veteran Program participants hospitalized with COVID-19

Remdesivir is the first US Food and Drug Administration (FDA)-approved drug for the treatment of coronavirus disease 2019 (COVID-19). We conducted a retrospective pharmacogenetic study to examine remdesivir-associated liver enzyme elevation among Million Veteran Program participants hospitalized with COVID-19 between March 15, 2020, and June 30, 2021. Pharmacogene phenotypes were assigned using Stargazer. Linear regression was performed on peak log-transformed enzyme values, stratified by population, adjusted for age, sex, baseline liver enzymes, comorbidities, and 10 population-specific principal components. Patients on remdesivir had higher peak alanine aminotransferase (ALT) values following treatment initiation compared with patients not receiving remdesivir. Remdesivir administration was associated with a 33% and 24% higher peak ALT in non-Hispanic White (NHW) and non-Hispanic Black (NHB) participants (p < 0.001), respectively. In a multivariable model, NHW CYP2C19 intermediate/poor metabolizers had a 9% increased peak ALT compared with NHW normal/rapid/ultrarapid metabolizers (p = 0.015); this association was not observed in NHB participants. In summary, remdesivir-associated ALT elevations appear to be multifactorial, and further studies are needed.

Predicted expression of genes involved in the thiopurine metabolic pathway and azathioprine discontinuation due to myelotoxicity

TPMT and NUDT15 variants explain less than 25% of azathioprine-associated myelotoxicity. There are 25 additional genes in the thiopurine pathway that could also contribute to azathioprine myelotoxicity. We hypothesized that among TPMT and NUDT15 normal metabolizers, a score combining the genetically predicted expression of other proteins in the thiopurine pathway would be associated with a higher risk for azathioprine discontinuation due to myelotoxicity. We conducted a retrospective cohort study of new users of azathioprine who were normal TPMT and NUDT15 metabolizers. In 1201 White patients receiving azathioprine for an inflammatory disease, we used relaxed Least Absolute Shrinkage and Selection Operator (LASSO) regression to select genes that built a score for discontinuing azathioprine due to myelotoxicity. The score incorporated the predicted expression of AOX1 and NME1. Patients in the highest score tertile had a higher risk of discontinuing azathioprine compared to those in the lowest tertile (hazard ratio [HR] = 2.15, 95% confidence interval [CI] = 1.11-4.19, p = 0.024). Results remained significant after adjusting for a propensity score, including sex, tertile of calendar year at initial dose, initial dose, age at baseline, indication, prior TPMT testing, and the first 10 principal components of the genetic data (HR = 2.11, 95% CI = 1.08-4.13, p = 0.030). We validated the results in a cohort (N = 517 non-White patients and those receiving azathioprine to prevent transplant rejection) that included all other patients receiving azathioprine (HR = 2.00, (95% CI = 1.09-3.65, p = 0.024). In conclusion, among patients who were TPMT and NUDT15 normal metabolizers, a score combining the predicted expression of AOX1 and NME1 was associated with an increased risk for discontinuing azathioprine due to myelotoxicity.

APOL1 Risk Variants, Acute Kidney Injury, and Death in Participants With African Ancestry Hospitalized With COVID-19 From the Million Veteran Program

IMPORTANCE

Coronavirus disease 2019 (COVID-19) confers significant risk of acute kidney injury (AKI). Patients with COVID-19 with AKI have high mortality rates.

OBJECTIVE

Individuals with African ancestry with 2 copies of apolipoprotein L1 (APOL1) variants G1 or G2 (high-risk group) have significantly increased rates of kidney disease. We tested the hypothesis that the APOL1 high-risk group is associated with a higher-risk of COVID-19-associated AKI and death.

DESIGN, SETTING, AND PARTICIPANTS

This retrospective cohort study included 990 participants with African ancestry enrolled in the Million Veteran Program who were hospitalized with COVID-19 between March 2020 and January 2021 with available genetic information.

EXPOSURES

The primary exposure was having 2 APOL1 risk variants (RV) (APOL1 high-risk group), compared with having 1 or 0 risk variants (APOL1 low-risk group).

MAIN OUTCOMES AND MEASURES

The primary outcome was AKI. The secondary outcomes were stages of AKI severity and death. Multivariable logistic regression analyses adjusted for preexisting comorbidities, medications, and inpatient AKI risk factors; 10 principal components of ancestry were performed to study these associations. We performed a subgroup analysis in individuals with normal kidney function prior to hospitalization (estimated glomerular filtration rate ≥60 mL/min/1.73 m2).

RESULTS

Of the 990 participants with African ancestry, 905 (91.4%) were male with a median (IQR) age of 68 (60-73) years. Overall, 392 (39.6%) patients developed AKI, 141 (14%) developed stages 2 or 3 AKI, 28 (3%) required dialysis, and 122 (12.3%) died. One hundred twenty-five (12.6%) of the participants were in the APOL1 high-risk group. Patients categorized as APOL1 high-risk group had significantly higher odds of AKI (adjusted odds ratio [OR], 1.95; 95% CI, 1.27-3.02; P = .002), higher AKI severity stages (OR, 2.03; 95% CI, 1.37-2.99; P < .001), and death (OR, 2.15; 95% CI, 1.22-3.72; P = .007). The association with AKI persisted in the subgroup with normal kidney function (OR, 1.93; 95% CI, 1.15-3.26; P = .01). Data analysis was conducted between February 2021 and April 2021.

CONCLUSIONS AND RELEVANCE

In this cohort study of veterans with African ancestry hospitalized with COVID-19 infection, APOL1 kidney risk variants were associated with higher odds of AKI, AKI severity, and death, even among individuals with prior normal kidney function.

TPMT and NUDT15 Variants Predict Discontinuation of Azathioprine for Myelotoxicity in Patients with Inflammatory Disease: Real-World Clinical Results

Azathioprine is used frequently to treat several inflammatory conditions. However, treatment is limited by adverse events-in particular, myelotoxicity. Thiopurine-S-methyltransferase (TPMT) and nudix hydrolase-15 (NUDT15) are enzymes involved in azathioprine metabolism; variants in the genes encoding these enzymes increase the risk for azathioprine myelotoxicity. The Clinical Pharmacogenetics Implementation Consortium (CPIC) has recommended dose adjustments based on the results of TPMT and NUDT15 genotyping. However, little is known about the importance of this genetic information in routine clinical care. We hypothesized that in patients with inflammatory diseases, TPMT and NUDT15 genotype data predict the risk of discontinuing azathioprine due to myelotoxicity. This was a retrospective cohort study in 1,403 new adult azathioprine users for the management of inflammatory conditions for whom we had genetic information and clinical data. Among patients who discontinued azathioprine, we adjudicated the reason(s). Genotyping was performed using the Illumina Infinium Expanded Multi-Ethnic Genotyping Array plus custom content. We used CPIC guidelines to determine TPMT and NUDT15 metabolizer status; patients were grouped as either: (i) poor/intermediate, or (ii) normal/indeterminate metabolizers. We classified 110 patients as poor/intermediate, and 1,293 patients as normal/indeterminate metabolizers. Poor/intermediate status was associated with a higher risk for azathioprine discontinuation due to myelotoxicity compared to normal/indeterminate metabolizers (hazard ratio (HR) = 2.90, 95% confidence interval (CI): 1.58-5.31, P = 0.001). This association remained significant after adjustment for race, age at initiation, sex, primary indication, and initial daily dose of azathioprine (adjusted HR (aHR) = 2.67, 95% CI: 1.44-4.94, P = 0.002). In conclusion, TPMT and NUDT15 metabolizer status predicts discontinuation due to myelotoxicity for patients taking azathioprine for inflammatory conditions.

Association of Rivaroxaban vs Apixaban With Major Ischemic or Hemorrhagic Events in Patients With Atrial Fibrillation

IMPORTANCE

The comparative effectiveness of rivaroxaban and apixaban, the most frequently prescribed oral anticoagulants for ischemic stroke prevention in patients with atrial fibrillation, is uncertain.

OBJECTIVE

To compare major ischemic and hemorrhagic outcomes in patients with atrial fibrillation treated with rivaroxaban or apixaban.

DESIGN, SETTING, AND PARTICIPANTS

Retrospective cohort study using computerized enrollment and claims files for US Medicare beneficiaries 65 years or older. Between January 1, 2013, and November 30, 2018, a total of 581 451 patients with atrial fibrillation began rivaroxaban or apixaban treatment and were followed up for 4 years, through November 30, 2018.

EXPOSURES

Rivaroxaban (n = 227 572) and apixaban (n = 353 879), either standard or reduced dose.

MAIN OUTCOMES AND MEASURES

The primary outcome was a composite of major ischemic (stroke/systemic embolism) and hemorrhagic (intracerebral hemorrhage/other intracranial bleeding/fatal extracranial bleeding) events. Secondary outcomes were nonfatal extracranial bleeding and total mortality (fatal ischemic/hemorrhagic event or other death during follow-up). Rates, hazard ratios (HRs), and rate differences (RDs) were adjusted for baseline differences in comorbidity with inverse probability of treatment weighting.

RESULTS

Study patients (mean age, 77.0 years; 291 966 [50.2%] women; 134 393 [23.1%] receiving reduced dose) had 474 605 person-years of follow-up (median [IQR] of 174 [62-397] days). The adjusted primary outcome rate for rivaroxaban was 16.1 per 1000 person-years vs 13.4 per 1000 person-years for apixaban (RD, 2.7 [95% CI, 1.9-3.5]; HR, 1.18 [95% CI, 1.12-1.24]). The rivaroxaban group had increased risk for both major ischemic events (8.6 vs 7.6 per 1000 person-years; RD, 1.1 [95% CI, 0.5-1.7]; HR, 1.12 [95% CI, 1.04-1.20]) and hemorrhagic events (7.5 vs 5.9 per 1000 person-years; RD, 1.6 [95% CI, 1.1-2.1]; HR, 1.26 [95% CI, 1.16-1.36]), including fatal extracranial bleeding (1.4 vs 1.0 per 1000 person-years; RD, 0.4 [95% CI, 0.2-0.7]; HR, 1.41 [95% CI, 1.18-1.70]). Patients receiving rivaroxaban had increased risk of nonfatal extracranial bleeding (39.7 vs 18.5 per 1000 person-years; RD, 21.1 [95% CI, 20.0-22.3]; HR, 2.07 [95% CI, 1.99-2.15]), fatal ischemic/hemorrhagic events (4.5 vs 3.3 per 1000 person-years; RD, 1.2 [95% CI, 0.8-1.6]; HR, 1.34 [95% CI, 1.21-1.48]), and total mortality (44.2 vs 41.0 per 1000 person-years; RD, 3.1 [95% CI, 1.8-4.5]; HR, 1.06 [95% CI, 1.02-1.09]). The risk of the primary outcome was increased for rivaroxaban in both those receiving the reduced dose (27.4 vs 21.0 per 1000 person-years; RD, 6.4 [95% CI, 4.1-8.7]; HR, 1.28 [95% CI, 1.16-1.40]) and the standard dose (13.2 vs 11.4 per 1000 person-years; RD, 1.8 [95% CI, 1.0-2.6]; HR, 1.13 [95% CI, 1.06-1.21]) groups.

CONCLUSIONS AND RELEVANCE

Among Medicare beneficiaries 65 years or older with atrial fibrillation, treatment with rivaroxaban compared with apixaban was associated with a significantly increased risk of major ischemic or hemorrhagic events.

Host Genetic Determinants Associated With Helicobacter pylori Eradication Treatment Failure: A Systematic Review and Meta-analysis

BACKGROUND & AIMS

Helicobacter pylori infects approximately 50% of individuals worldwide. Successful H pylori eradication is associated with reduced risk of gastric cancer and peptic ulcer disease, among other conditions. We hypothesized that host genetic determinants, especially those affecting gastric pH, might contribute to eradication therapy failure, particularly when treatment adherence and antibiotic susceptibility are confirmed. We aimed to conduct a meta-analysis of host genetic variants associated with H pylori eradication failure.

METHODS

We searched the literature for studies comparing posttreatment H pylori eradication failure vs success (outcome) according to host genetic polymorphisms (exposure). Reference groups were defined according to genotypes (or corresponding phenotypes) hypothesized to be associated with successful eradication. We pooled estimates using a random-effects model and performed comprehensive sensitivity analyses.

RESULTS

We analyzed 57 studies from 11 countries; the vast majority analyzed CYP2C19 polymorphisms. Among individuals prescribed eradication regimens with proton pump inhibitors predominantly CYP2C19 metabolized, enhanced vs poor metabolizer phenotypes were associated with a 2.52-fold significantly higher likelihood of eradication failure and 4.44-fold significantly higher likelihood when treatment adherence and H pylori clarithromycin susceptibility (if relevant) were confirmed. There was no association between CYP2C19 variants and eradication failure if proton pump inhibitors less metabolized by or that bypass CYP2C19 metabolism were used. IL1B polymorphisms that are vs are not associated with less gastric acid suppression were associated with 1.72-fold significantly higher likelihood of eradication failure. There was no association between MDR1 polymorphisms and H pylori eradication failure. The certainty of evidence was moderate.

CONCLUSION

Based on meta-analysis, we identified host genetic polymorphisms significantly associated with H pylori eradication failure; host genetics might underlie eradication failure among treatment-adherent individuals with confirmed H pylori antibiotic susceptibility.

Mortality and concurrent use of opioids and hypnotics in older patients: A retrospective cohort study

BACKGROUND

Benzodiazepine hypnotics and the related nonbenzodiazepine hypnotics (z-drugs) are among the most frequently prescribed medications for older adults. Both can depress respiration, which could have fatal cardiorespiratory effects, particularly among patients with concurrent opioid use. Trazodone, frequently prescribed in low doses for insomnia, has minimal respiratory effects, and, consequently, may be a safer hypnotic for older patients. Thus, for patients beginning treatment with benzodiazepine hypnotics or z-drugs, we compared deaths during periods of current hypnotic use, without or with concurrent opioids, to those for comparable patients receiving trazodone in doses up to 100 mg.

METHODS AND FINDINGS

The retrospective cohort study in the United States included 400,924 Medicare beneficiaries 65 years of age or older without severe illness or evidence of substance use disorder initiating study hypnotic therapy from January 2014 through September 2015. Study endpoints were out-of-hospital (primary) and total mortality. Hazard ratios (HRs) were adjusted for demographic characteristics, psychiatric and neurologic disorders, cardiovascular and renal conditions, respiratory diseases, pain-related diagnoses and medications, measures of frailty, and medical care utilization in a time-dependent propensity score-stratified analysis. Patients without concurrent opioids had 32,388 person-years of current use, 260 (8.0/1,000 person-years) out-of-hospital and 418 (12.9/1,000) total deaths for benzodiazepines; 26,497 person-years,150 (5.7/1,000) out-of-hospital and 227 (8.6/1,000) total deaths for z-drugs; and 16,177 person-years,156 (9.6/1,000) out-of-hospital and 256 (15.8/1,000) total deaths for trazodone. Out-of-hospital and total mortality for benzodiazepines (respective HRs: 0.99 [95% confidence interval, 0.81 to 1.22, p = 0.954] and 0.95 [0.82 to 1.14, p = 0.513] and z-drugs (HRs: 0.96 [0.76 to 1.23], p = 0.767 and 0.87 [0.72 to 1.05], p = 0.153) did not differ significantly from that for trazodone. Patients with concurrent opioids had 4,278 person-years of current use, 90 (21.0/1,000) out-of-hospital and 127 (29.7/1,000) total deaths for benzodiazepines; 3,541 person-years, 40 (11.3/1,000) out-of-hospital and 64 (18.1/1,000) total deaths for z-drugs; and 2,347 person-years, 19 (8.1/1,000) out-of-hospital and 36 (15.3/1,000) total deaths for trazodone. Out-of-hospital and total mortality for benzodiazepines (HRs: 3.02 [1.83 to 4.97], p < 0.001 and 2.21 [1.52 to 3.20], p < 0.001) and z-drugs (HRs: 1.98 [1.14 to 3.44], p = 0.015 and 1.65 [1.09 to 2.49], p = 0.018) were significantly increased relative to trazodone; findings were similar with exclusion of overdose deaths or restriction to those with cardiovascular causes. Limitations included composition of the study cohort and potential confounding by unmeasured variables.

CONCLUSIONS

In US Medicare beneficiaries 65 years of age or older without concurrent opioids who initiated treatment with benzodiazepine hypnotics, z-drugs, or low-dose trazodone, study hypnotics were not associated with mortality. With concurrent opioids, benzodiazepines and z-drugs were associated with increased out-of-hospital and total mortality. These findings indicate that the dangers of benzodiazepine-opioid coadministration go beyond the documented association with overdose death and suggest that in combination with opioids, the z-drugs may be more hazardous than previously thought.

Pleiotropy of systemic lupus erythematosus risk alleles and cardiometabolic disorders: A phenome-wide association study and inverse-variance weighted meta-analysis

OBJECTIVES

To test the hypothesis that genetic predisposition to systemic lupus erythematosus (SLE) increases the risk of cardiometabolic disorders.

METHODS

Using 41 single nucleotide polymorphisms (SNPs) associated with SLE, we calculated a weighted genetic risk score (wGRS) for SLE. In a large biobank we tested the association between this wGRS and 9 cardiometabolic phenotypes previously associated with SLE: atrial fibrillation, ischemic stroke, coronary artery disease, type 1 and type 2 diabetes, obesity, chronic kidney disease, hypertension, and hypercholesterolemia. Additionally, we performed a phenome-wide association analysis (pheWAS) to discover novel clinical associations with a genetic predisposition to SLE. Findings were replicated in the Electronic Medical Records and Genomics (eMERGE) Network. To further define the association between SLE-related risk alleles and the selected cardiometabolic phenotypes, we performed an inverse variance weighted regression (IVWR) meta-analysis.

RESULTS

The wGRS for SLE was calculated in 74,759 individuals of European ancestry. Among the pre-selected phenotypes, the wGRS was significantly associated with type 1 diabetes (OR [95%CI] =1.11 [1.06, 1.17], P-value = 1.05x10-5). In the PheWAS, the wGRS was associated with several autoimmune phenotypes, kidney disorders, and skin neoplasm; but only the associations with autoimmune phenotypes were replicated. In the IVWR meta-analysis, SLE-related risk alleles were nominally associated with type 1 diabetes (P = 0.048) but the associations were heterogeneous and did not meet the adjusted significance threshold.

CONCLUSION

A weighted GRS for SLE was associated with an increased risk of several autoimmune-related phenotypes including type I diabetes but not with cardiometabolic disorders.

Association of Apparent Treatment-Resistant Hypertension With Differential Risk of End-Stage Kidney Disease Across Racial Groups in the Million Veteran Program

[Figure: see text].

High-throughput framework for genetic analyses of adverse drug reactions using electronic health records

Understanding the contribution of genetic variation to drug response can improve the delivery of precision medicine. However, genome-wide association studies (GWAS) for drug response are uncommon and are often hindered by small sample sizes. We present a high-throughput framework to efficiently identify eligible patients for genetic studies of adverse drug reactions (ADRs) using “drug allergy” labels from electronic health records (EHRs). As a proof-of-concept, we conducted GWAS for ADRs to 14 common drug/drug groups with 81,739 individuals from Vanderbilt University Medical Center’s BioVU DNA Biobank. We identified 7 genetic loci associated with ADRs at P < 5 × 10⁻⁸, including known genetic associations such as CYP2D6 and OPRM1 for CYP2D6-metabolized opioid ADR. Additional expression quantitative trait loci and phenome-wide association analyses added evidence to the observed associations. Our high-throughput framework is both scalable and portable, enabling impactful pharmacogenomic research to improve precision medicine.

Adverse drug reactions are a considerable burden on the healthcare system. Genetic studies can improve our understanding of the pathophysiological mechanisms of adverse drug reactions but have been hindered by small sample sizes. Drug responses are less often recorded than physiological traits and common diseases. Here, we present a high-throughput framework to efficiently identify eligible patients for genetic studies of adverse drug reactions from electronic health records. We validated our approach by conducting genome-wide association studies for adverse reactions to 14 common drug/drug groups with 81,739 individuals from Vanderbilt University Medical Centre’s BioVU DNA Biobank, identifying 7 genetic loci associated with adverse drug reactions. Our high-throughput framework can enable impactful pharmacogenomic research to help develop clinical guidelines for the delivery of the right drug to the right person.

Ambulatory blood pressure in patients with systemic lupus erythematosus: Association with markers of immune activation

Objectives: Ambulatory blood pressure monitoring measures 24-hour blood pressure, night-time blood pressure, and impaired dipping of nocturnal blood pressure, parameters that better predict cardiovascular risk than standard office blood pressure measurements. Systemic lupus erythematosus is characterized by immune system hyperactivity, elevated cardiovascular risk and high prevalence of hypertension; however, little is known about ambulatory blood pressure in lupus patients and its relationship to immune activation. Methods: We studied 26 patients with lupus and 26 control subjects. We obtained ambulatory 24-hour blood pressure measurements and report plasma concentrations of 77 markers of immune activation using a multiplex immunoassay and assessed their association with blood pressure measurements. Results: Despite similar office blood pressure measurements in patients with lupus and controls, lupus patients had higher 24-hour systolic [median (interquartile range) 129 (113 - 140) vs. 116 (111 - 121) mmHg, p = 0.03] and diastolic blood pressure [80 (69 - 86) vs. 72 (64 - 75) mmHg, p = 0.006] as well as less nocturnal dipping [7.8% (5.1 - 14.2%) vs. 12.0% (8.1 20.0%)] p = 0.03], compared to controls. In patients with lupus, markers of the innate (monocyte chemotactic protein-3) and adaptive immune systems [CUB domain-containing protein-1 and Interleukin-15 receptor subunit-α,] were associated with nocturnal blood pressure measurements and attenuated nocturnal dipping. In conclusion, 24-hour systolic and diastolic blood pressure was higher and nocturnal blood pressure dipping was attenuated in patients with lupus compared to control subjects. Conclusion: In patients with SLE, nocturnal blood pressure and attenuated nocturnal blood pressure dipping were significantly associated with several innate and adaptive immune system biomarkers.

Pleiotropy in the Genetic Predisposition to Rheumatoid Arthritis: A Phenome-Wide Association Study and Inverse Variance-Weighted Meta-Analysis

OBJECTIVE

This study was undertaken to investigate the hypothesis that a genetic predisposition toward rheumatoid arthritis (RA) increases the risk of 10 cardiometabolic and autoimmune disorders previously associated with RA in epidemiologic studies, and to define new genetic pleiotropy present in RA.

METHODS

Two approaches were used to test our hypothesis. First, we constructed a weighted genetic risk score (wGRS) and then examined its association with 10 prespecified disorders. Additionally, a phenome-wide association study (PheWAS) was carried out to identify potential new associations. Second, inverse variance-weighted regression (IVWR) meta-analysis was used to characterize the association between genetic susceptibility to RA and the prespecified disorders, with the results expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs).

RESULTS

The wGRS for RA was significantly associated with type 1 diabetes mellitus (DM) (OR 1.10 [95% CI 1.04-1.16]; P = 9.82 × 10^-4 ) and multiple sclerosis (OR 0.82 [95% CI 0.77-0.88]; P = 1.73 × 10^-8 ), but not with other cardiometabolic phenotypes. In the PheWAS, wGRS was also associated with an increased risk of several autoimmune phenotypes including RA, thyroiditis, and systemic sclerosis, and with a decreased risk of demyelinating disorders. In the IVWR meta-analyses, RA was significantly associated with an increased risk of type 1 DM (P = 1.15 × 10^-14 ), with evidence of horizontal pleiotropy (Mendelian Randomization-Egger intercept estimate P = 0.001) likely driven by rs2476601, a PTPN22 variant. The association between type 1 DM and RA remained significant (P = 9.53 × 10^-9 ) after excluding rs2476601, with no evidence of horizontal pleiotropy (intercept estimate P = 0.939). RA was also significantly associated with type 2 DM and C-reactive protein levels. These associations were driven by variation in the major histocompatibility complex region.

CONCLUSION

This study presents evidence of pleiotropy between the genetic predisposition to RA and associated phenotypes found in other autoimmune and cardiometabolic disorders, including type 1 DM.

Increased Incidence of Resistant Hypertension in Patients With Systemic Lupus Erythematosus: A Retrospective Cohort Study

OBJECTIVE

To compare the risk of resistant hypertension (RHTN) in patients with systemic lupus erythematosus (SLE) and in controls without SLE, and to define factors associated with RHTN in patients with SLE.

METHODS

We studied 1,044 patients with SLE and 5,241 control subjects using de-identified electronic health records from a tertiary care center. SLE was defined as ≥4 International Classification of Diseases, Ninth Revision codes for SLE and antinuclear antibody titer ≥1:160. RHTN was defined as uncontrolled blood pressure on 3 antihypertensive medications or requiring 4 or more antihypertensives to attain control. First, we compared the risk of RHTN between groups. Second, we examined the association between RHTN and all-cause mortality in patients with SLE.

RESULTS

RHTN was nearly twice as prevalent in patients with SLE compared to control subjects (10.2% and 5.3%, respectively), with an incidence rate of 10.2 versus 6.1 cases per 1,000 person-years of observation (hazard ratio [HR] 1.72 [95% confidence interval 1.28-2.30]; P < 0.001, adjusted for age, sex, race, baseline end-stage renal disease [ESRD], creatinine, and calendar year). In patients with SLE, we found associations between RHTN and black race, lower renal function, hypercholesterolemia, and increased inflammatory markers. RHTN was associated with a significantly higher mortality risk (HR 2.91, P = 0.0005) after adjustment for age, sex, race, calendar year, creatinine, baseline ESRD, and number of visits.

CONCLUSION

Patients with SLE have a higher risk of RHTN compared to frequency-matched controls, independent of multiple covariates. RHTN is an important comorbidity for clinicians to recognize in SLE, because it is associated with a higher risk of mortality.

Precision medicine for rheumatologists: lessons from the pharmacogenomics of azathioprine

Precision medicine aims to personalize treatment for both effectiveness and safety. As a critical component of this emerging initiative, pharmacogenomics seeks to guide drug treatment based on genetics. In this review article, we give an overview of pharmacogenomics in the setting of an immunosuppressant frequently prescribed by rheumatologists, azathioprine. Azathioprine has a narrow therapeutic index and a high risk of adverse events. By applying candidate gene analysis and unbiased approaches, researchers have identified multiple variants associated with an increased risk for adverse events associated with azathioprine, particularly bone marrow suppression. Variants in two genes, TPMT and NUDT15, are widely recognized, leading drug regulatory agencies and professional organizations to adopt recommendations for testing before initiation of azathioprine therapy. As more gene-drug interactions are discovered, our field will continue to face the challenge of balancing benefits and costs associated with genetic testing. However, novel approaches in genomics and the integration of clinical and genetic factors into risk scores offer unprecedented opportunities for the application of pharmacogenomics in routine practice. Key Points • Pharmacogenomics can help us understand how individuals' genetics may impact their response to medications. • Azathioprine is a success story for the clinical implementation of pharmacogenomics, particularly the effects of TPMT and NUDT15 variants on myelosuppression. • As our knowledge advances, testing and dosing recommendations will continue to evolve, with our field striving to balance costs and benefits to patients. • As we aim toward the goals of precision medicine, future research may integrate increasingly individualized traits-including clinical and genetic characteristics-to predict the safety and efficacy of particular medications for individual patients.

Heart Failure Risk Associated With Rheumatoid Arthritis-Related Chronic Inflammation

Background

Inflammation may contribute to incident heart failure (HF). Rheumatoid arthritis (RA), a prototypic inflammatory condition, may serve as a model for understanding inflammation‐related HF risk.

Methods and Results

Using the Vanderbilt University Medical Center electronic health record, we retrospectively identified 9889 patients with RA and 9889 control patients without autoimmune disease matched for age, sex, and race. Prevalent HF at entry into the electronic health record or preceding RA diagnosis was excluded. Incident HF was ascertained using International Classification of Diseases, Ninth Revision (ICD‐9), codes and medications. Over 177 566 person‐years of follow‐up, patients with RA were at 21% greater risk of HF (95% CI, 3–42%) independent of traditional cardiovascular risk factors. Among patients with RA, higher CRP (C‐reactive protein) was associated with greater HF risk (P<0.001), while the anti‐inflammatory drug methotrexate was associated with ≈25% lower HF risk (P=0.021). In a second cohort (n=115) of prospectively enrolled patients with and without RA, we performed proteomics and cardiac magnetic resonance imaging to discover circulating markers of inflammation associated with cardiac structure and function. Artemin levels were higher in patients with RA compared with controls (P=0.009), and higher artemin levels were associated with worse ventricular end‐systolic elastance and ventricular‐vascular coupling ratio (P=0.044 and P=0.031, respectively).

Conclusions

RA, a prototypic chronic inflammatory condition, is associated with increased risk of HF. Among patients with RA, higher levels of CRP were associated with greater HF risk, while methotrexate was associated with lower risk.

Tissue sodium content in patients with systemic lupus erythematosus: association with disease activity and markers of inflammation

OBJECTIVES

Sodium (Na⁺) is stored in the skin and muscle and plays an important role in immune regulation. In animal models, increased tissue Na⁺ is associated with activation of the immune system, and high salt intake exacerbates autoimmune disease and worsens hypertension. However, there is no information about tissue Na⁺ and human autoimmune disease. We hypothesized that muscle and skin Na⁺ content is (a) higher in patients with systemic lupus erythematosus (SLE) than in control subjects, and (b) associated with blood pressure, disease activity, and inflammation markers (interleukin (IL)-6, IL-10 and IL-17 A) in SLE.

METHODS

Lower-leg skin and muscle Na⁺ content was measured in 23 patients with SLE and in 28 control subjects using ²³Na⁺ magnetic resonance imaging. Demographic and clinical information was collected from interviews and chart review, and blood pressure was measured. Disease activity was assessed using the SLE Disease Activity Index (SLEDAI). Plasma inflammation markers were measured by multiplex immunoassay.

RESULTS

Muscle Na⁺ content was higher in patients with SLE (18.8 (16.7-18.3) mmol/L) than in control subjects (15.8 (14.7-18.3) mmol/L; p < 0.001). Skin Na⁺ content was also higher in SLE patients than in controls, but this difference was not statistically significant. Among patients with SLE, muscle Na⁺ was associated with SLEDAI and higher concentrations of IL-10 after adjusting for age, race, and sex. Skin Na⁺ was significantly associated with systolic blood pressure, but this was attenuated after covariate adjustment.

CONCLUSION

Patients with SLE had higher muscle Na⁺ content than control subjects. In patients with SLE, higher muscle Na⁺ content was associated with higher disease activity and IL-10 concentrations.

Individual short-acting opioids and the risk of opioid-related adverse events in adolescents

PURPOSE

Hydrocodone, codeine, oxycodone, and tramadol are frequently prescribed to adolescents for moderate pain related to minor trauma or dental, surgical, or medical procedures. Pharmacokinetic and pharmacodynamic differences between these opioids could affect their relative safety. We aimed to compare occurrence of opioid-related adverse events in adolescents without cancer or other severe conditions taking hydrocodone, codeine, oxycodone, and tramadol.

METHODS

Retrospective cohort study of 201 940 Tennessee Medicaid enrollees 12 to 17 years of age without cancer, other severe conditions, or evidence of substance abuse with 529 731 filled prescriptions for study opioids. Adverse events were defined as an emergency department visit, hospital admission, or death related to opioid use, confirmed by medical record review. Serious events had opioid-related escalation of care, hospitalization, or death. Propensity-score adjusted hazard ratios (HRs) were calculated with hydrocodone as the reference category.

RESULTS

The incidence of opioid-related adverse events per 10 000 person-years of opioid exposure was 97.5 for hydrocodone (127 events/13 026 person-years), 91.2 for codeine (58/6,359), 229.7 for oxycodone (43/1,872), and 317.7 for tramadol (47/1479). The HRs for tramadol in comparison with hydrocodone for all and serious events were 2.98 (2.03-4.39) and 2.94 (1.81-4.75), respectively. Increased risk for tramadol was consistently present when the adverse events were restricted to those with neurologic-respiratory depression/other symptoms of possible overdose.

CONCLUSION

In adolescents without cancer or other severe conditions prescribed short-acting opioids, the incidence of both all opioid-related adverse events and more serious events with opioid-related escalation of care, hospitalization, or death was consistently greater for tramadol than for hydrocodone.

A Genetic Approach to the Association Between PCSK9 and Sepsis

IMPORTANCE

Whether the PCSK9 gene is associated with the progress from infection to sepsis is unknown to date.

OBJECTIVE

To test the associations between PCSK9 genetic variants, a PCSK9 genetic risk score (GRS), or genetically estimated PCSK9 expression levels and the risk of sepsis among patients admitted to a hospital with infection.

DESIGN, SETTING, AND PARTICIPANTS

This retrospective cohort study used deidentified electronic health records to identify patients admitted to Vanderbilt University Medical Center, Nashville, Tennessee, with infection. Patients were white adults, had a code indicating infection from the International Classification of Diseases, Ninth Revision, Clinical Modification, or the International Statistical Classification of Diseases, Tenth Revision, Clinical Modification, and received an antibiotic within 1 day of hospital admission (N = 61 502). Data were collected from January 1, 1993, through December 31, 2017, and analyzed from April 1, 2018, to March 16, 2019.

EXPOSURES

Four known PCSK9 functional variants, a GRS for PCSK9, and genetically estimated PCSK9 expression.

MAIN OUTCOMES AND MEASURES

The primary outcome was sepsis; secondary outcomes included cardiovascular failure and in-hospital death.

RESULTS

Of patients with infection, genotype information was available in 10 922 white patients for PCSK9 functional variants (5628 men [51.5%]; mean [SD] age, 60.1 [15.7] years), including 7624 patients with PCSK9 GRS and 6033 patients with estimated PCSK9 expression. Of these, 3391 developed sepsis, 835 developed cardiovascular failure, and 366 died during hospitalization. None of the 4 functional PCSK9 variants were significantly associated with sepsis, cardiovascular failure, or in-hospital death, with or without adjustment for (1) age and sex or (2) age, sex, and Charlson-Deyo comorbidities (in model adjusted for age, sex, and comorbidities, odds ratios for any loss-of function variant were 0.96 [95% CI, 0.88-1.04] for sepsis, 1.05 [95% CI, 0.90-1.22] for cardiovascular failure, and 0.89 [95% CI, 0.72-1.11] for death). Similarly, neither the PCSK9 GRS nor genetically estimated PCSK9 expression were significantly associated with sepsis, cardiovascular failure, or in-hospital death in any of the analysis models. For GRS, in the full model adjusted for age, sex, and comorbidities, the odds ratios were 1.01 for sepsis (95% CI, 0.96-1.06; P = .70), 1.03 for cardiovascular failure (95% CI, 0.95-1.12; P = .48), and 1.05 for in-hospital death (95% CI, 0.92-1.19; P = .50). For genetically estimated PCSK9 expression, in the full model adjusted for age, sex, and comorbidities, the odds ratios were 1.01 for sepsis (95% CI, 0.95-1.06; P = .86), 0.96 for cardiovascular failure (95% CI, 0.88-1.05; P = .41), and 0.99 for in-hospital death (95% CI, 0.87-1.14; P = .94).

CONCLUSIONS AND RELEVANCE

In this study, PCSK9 genetic variants were not significantly associated with risk of sepsis or the outcomes of sepsis in patients hospitalized with infection.

Exploring the Lipid Paradox Theory in Rheumatoid Arthritis: Associations of Low Circulating Low-Density Lipoprotein Concentration With Subclinical Coronary Atherosclerosis

OBJECTIVE

Rheumatoid arthritis (RA) patients with the lowest circulating low-density lipoprotein (LDL) concentrations are at heightened risk of cardiovascular events. However, the atherosclerosis burden within this subgroup is unknown.

METHODS

RA patients pooled from 4 cohort studies of cardiovascular disease (CVD; n = 546) were compared with non-RA controls from the Multi-Ethnic Study of Atherosclerosis (n = 5,279). Those taking lipid-lowering medications were excluded. Differences in cardiac computed tomography-derived Agatston coronary artery calcium (CAC) scores between the RA and control groups were compared across strata of LDL concentration.

RESULTS

Among those with low LDL concentrations (<70 mg/dl), mean adjusted CAC scores were >4-fold higher for RA patients than for controls (18.6 versus 4.6 Agatston units, respectively; P < 0.001), a difference significantly greater than that in any other LDL concentration stratum except LDL concentration ≥160 mg/dl. Similarly, 32% of the RA patients with low LDL concentration had a CAC score of ≥100 Agatston units compared with only 7% of controls in the same LDL concentration stratum (odds ratio 5.97; P < 0.001), a difference significantly greater than that in all of the other LDL concentration strata. Low LDL concentration was most strongly associated with higher CAC score among RA patients who were white, had ever smoked, or were not obese. Other than a higher frequency of current smokers, RA patients with low LDL concentrations did not have more CVD risk factors or higher measures of RA disease activity or severity than RA patients with higher LDL concentrations.

CONCLUSION

RA patients with low LDL concentration may represent a group for whom heightened screening and prevention of atherosclerotic CVD is appropriate.

Increased blood pressure visit-to-visit variability in patients with systemic lupus erythematosus: association with inflammation and comorbidity burden

BACKGROUND

Blood pressure visit-to-visit variability is a novel risk factor for deleterious long-term cardiac and renal outcomes in the general population. We hypothesized that patients with systemic lupus erythematosus (SLE) have greater blood pressure visit-to-visit variability than control subjects and that blood pressure visit-to-visit variability is associated with a higher comorbidity burden.

METHODS

We studied 899 patients with SLE and 4172 matched controls using de-identified electronic health records from an academic medical center. We compared blood pressure visit-to-visit variability measures in patients with SLE and control subjects and examined the association between blood pressure visit-to-visit variability and patients' characteristics.

RESULTS

Patients with SLE had higher systolic blood pressure visit-to-visit variability 9.7% (7.8-11.8%) than the control group 9.2% (7.4-11.2%), P < 0.001 by coefficient of variation. Additional measures of systolic blood pressure visit-to-visit variability (i.e. standard deviation, average real variation, successive variation and maximum measure-to-measure change) were also significantly higher in patients with SLE than in control subjects. In patients with SLE, blood pressure visit-to-visit variability correlated significantly with age, creatinine, CRP, triglyceride concentrations and the Charlson comorbidity score (all P < 0.05). Hydroxychloroquine use was associated with reduced blood pressure visit-to-visit variability (P < 0.001), whereas the use of antihypertensives, cyclophosphamide, mycophenolate mofetil and corticosteroids was associated with increased blood pressure visit-to-visit variability (P < 0.05).

CONCLUSION

Patients with SLE had higher blood pressure visit-to-visit variability than controls, and this increased blood pressure visit-to-visit variability was associated with greater Charlson comorbidity scores, several clinical characteristics and immunosuppressant medications. In particular, hydroxychloroquine prescription was associated with lower blood pressure visit-to-visit variability.

Real-world electronic health record identifies antimalarial underprescribing in patients with lupus nephritis

Antimalarials (AMs) reduce disease activity and improve survival in patients with systemic lupus erythematosus (SLE), but studies have reported low AM prescribing frequencies. Using a real-world electronic health record cohort, we examined if patient or provider characteristics impacted AM prescribing. We identified 977 SLE cases, 94% of whom were ever prescribed an AM. Older patients and patients with SLE nephritis were less likely to be on AMs. Current age (odds ratio = 0.97, p < 0.01) and nephritis (odds ratio = 0.16, p < 0.01) were both significantly associated with ever AM use after adjustment for sex and race. Of the 244 SLE nephritis cases, only 63% were currently on AMs. SLE nephritis subjects who were currently prescribed AMs were more likely to be followed by a rheumatologist than a nephrologist and less likely to have undergone dialysis or renal transplant (both p < 0.001). Non-current versus current SLE nephritis AM users had higher serum creatinine (p < 0.001), higher urine protein (p = 0.05), and lower hemoglobin levels (p < 0.01). As AMs reduce disease damage and improve survival in patients with SLE, our results demonstrate an opportunity to target future efforts to improve prescribing rates among multi-specialty providers.

Association Between Low-Density Lipoprotein Cholesterol Levels and Risk for Sepsis Among Patients Admitted to the Hospital With Infection

IMPORTANCE

Whether low levels of low-density lipoprotein cholesterol (LDL-C) are associated with increased risk of sepsis and poorer outcomes is unknown.

OBJECTIVE

To examine the association between LDL-C levels and risk of sepsis among patients admitted to the hospital with infection.

DESIGN, SETTING, AND PARTICIPANTS

Cohort study in which deidentified electronic health records were used to define a cohort of patients admitted to Vanderbilt University Medical Center, Nashville, Tennessee, with infection. Patients were white adults, had a code indicating infection from the International Classification of Diseases, Ninth Revision, Clinical Modification, and received an antibiotic within 1 day of hospital admission (N = 61 502). Data were collected from January 1, 1993, through December 31, 2017, and analyzed from January 24 through October 31, 2018.

INTERVENTIONS

Clinically measured LDL-C levels (excluding measurements <1 year before hospital admission and those associated with acute illness) and a genetic risk score (GRS).

MAIN OUTCOMES AND MEASURES

The primary outcome was sepsis; secondary outcomes included admission to an intensive care unit (ICU) and in-hospital death.

RESULTS

Among the 3961 patients with clinically measured LDL-C levels (57.8% women; mean [SD] age, 64.1 [15.9] years) and the 7804 with a GRS for LDL-C (54.0% men; mean [SD] age, 59.8 [15.2] years), lower measured LDL-C levels were significantly associated with increased risk of sepsis (odds ratio [OR], 0.86; 95% CI, 0.79-0.94; P = .001) and ICU admission (OR, 0.85; 95% CI, 0.76-0.96; P = .008), but not in-hospital mortality (OR, 0.80; 95% CI, 0.63-1.00; P = .06); however, none of these associations were statistically significant after adjustment for age, sex, and comorbidity variables (OR for risk of sepsis, 0.96 [95% CI, 0.88-1.06]; OR for ICU admission, 0.94 [95% CI, 0.83-1.06]; OR for in-hospital death, 0.97 [95% CI, 0.76-1.22]; P > .05 for all). The LDL-C GRS correlated with measured LDL-C levels (r = 0.24; P < 2.2 × 10-16) but was not significantly associated with any of the outcomes.

CONCLUSIONS AND RELEVANCE

Results of this study suggest that lower measured LDL-C levels were significantly associated with increased risk of sepsis and admission to ICU in patients admitted to the hospital with infection; however, this association was due to comorbidities because both clinical models adjusted for confounders, and the genetic model showed no increased risk. Levels of LDL-C do not appear to directly alter the risk of sepsis or poor outcomes in patients hospitalized with infection.

Association of Oral Anticoagulants and Proton Pump Inhibitor Cotherapy With Hospitalization for Upper Gastrointestinal Tract Bleeding

IMPORTANCE

Anticoagulant choice and proton pump inhibitor (PPI) cotherapy could affect the risk of upper gastrointestinal tract bleeding, a frequent and potentially serious complication of oral anticoagulant treatment.

OBJECTIVES

To compare the incidence of hospitalization for upper gastrointestinal tract bleeding in patients using individual anticoagulants with and without PPI cotherapy, and to determine variation according to underlying gastrointestinal bleeding risk.

DESIGN, SETTING, AND PARTICIPANTS

Retrospective cohort study in Medicare beneficiaries between January 1, 2011, and September 30, 2015.

EXPOSURES

Apixaban, dabigatran, rivaroxaban, or warfarin with or without PPI cotherapy.

MAIN OUTCOMES AND MEASURES

Hospitalizations for upper gastrointestinal tract bleeding: adjusted incidence and risk difference (RD) per 10 000 person-years of anticoagulant treatment, incidence rate ratios (IRRs).

RESULTS

There were 1 643 123 patients with 1 713 183 new episodes of oral anticoagulant treatment included in the cohort (mean [SD] age, 76.4 [2.4] years, 651 427 person-years of follow-up [56.1%] were for women, and the indication was atrial fibrillation for 870 330 person-years [74.9%]). During 754 389 treatment person-years without PPI cotherapy, the adjusted incidence of hospitalization for upper gastrointestinal tract bleeding (n = 7119) was 115 per 10 000 person-years (95% CI, 112-118). The incidence for rivaroxaban (n = 1278) was 144 per 10 000 person-years (95% CI, 136-152), which was significantly greater than the incidence of hospitalizations for apixaban (n = 279; 73 per 10 000 person-years; IRR, 1.97 [95% CI, 1.73-2.25]; RD, 70.9 [95% CI, 59.1-82.7]), dabigatran (n = 629; 120 per 10 000 person-years; IRR, 1.19 [95% CI, 1.08-1.32]; RD, 23.4 [95% CI, 10.6-36.2]), and warfarin (n = 4933; 113 per 10 000 person-years; IRR, 1.27 [95% CI, 1.19-1.35]; RD, 30.4 [95% CI, 20.3-40.6]). The incidence for apixaban was significantly lower than that for dabigatran (IRR, 0.61 [95% CI, 0.52-0.70]; RD, -47.5 [95% CI,-60.6 to -34.3]) and warfarin (IRR, 0.64 [95% CI, 0.57-0.73]; RD, -40.5 [95% CI, -50.0 to -31.0]). When anticoagulant treatment with PPI cotherapy (264 447 person-years; 76 per 10 000 person-years) was compared with treatment without PPI cotherapy, risk of upper gastrointestinal tract bleeding hospitalizations (n = 2245) was lower overall (IRR, 0.66 [95% CI, 0.62-0.69]) and for apixaban (IRR, 0.66 [95% CI, 0.52-0.85]; RD, -24 [95% CI, -38 to -11]), dabigatran (IRR, 0.49 [95% CI, 0.41-0.59]; RD, -61.1 [95% CI, -74.8 to -47.4]), rivaroxaban (IRR, 0.75 [95% CI, 0.68-0.84]; RD, -35.5 [95% CI, -48.6 to -22.4]), and warfarin (IRR, 0.65 [95% CI, 0.62-0.69]; RD, -39.3 [95% CI, -44.5 to -34.2]).

CONCLUSIONS AND RELEVANCE

Among patients initiating oral anticoagulant treatment, incidence of hospitalization for upper gastrointestinal tract bleeding was the highest in patients prescribed rivaroxaban, and the lowest for patients prescribed apixaban. For each anticoagulant, the incidence of hospitalization for upper gastrointestinal tract bleeding was lower among patients who were receiving PPI cotherapy. These findings may inform assessment of risks and benefits when choosing anticoagulant agents.

Co-Prescription of Strong CYP1A2 Inhibitors and the Risk of Tizanidine-Associated Hypotension: A Retrospective Cohort Study

Tizanidine, a widely used muscle relaxant that can lower blood pressure, is metabolized by the cytochrome P450 1A2 (CYP1A2). We studied 1,626 patients prescribed tizanidine and 5,012 prescribed cyclobenzaprine concurrently with a strong CYP1A2 inhibitor. The primary outcome was severe hypotension, defined as systolic blood pressure (SBP) ≤ 70 mmHg during periods of drug co-exposure. Severe hypotension occurred more often in the tizanidine group (2.03%; n = 33) than the cyclobenzaprine group (1.28%; n = 64); odds ratio (OR) = 1.60; P = 0.029. This difference remained statistically significant after adjustment for a log-transformed propensity score that included age, sex, race, Charlson's comorbidity index, and concurrent use of antihypertensive medications (OR = 1.57; P = 0.049). A sensitivity analysis that defined hypotension as SBP < 90 mmHg also yielded higher rates of hypotension among patients prescribed tizanidine. In conclusion, CYP1A2 inhibition increases the risk of hypotensive episodes associated with the use of tizanidine in routine clinical practice.