Larson, Nicholas B.

Objectives Sex as a biological variable affects response to opioids. However, few reports describe the prevalence of specific adverse reactions to commonly prescribed opioids in men and women separately. A large cohort was used to investigate sex differences in type and occurrence of adverse reactions associated with use of codeine, tramadol, oxycodone and hydrocodone. Design Retrospective cohort study. Setting Participants in the Right Drug, Right Dose, Right Time (RIGHT) Study. Participants The medical records of 8457 participants in the RIGHT Study who received an opioid prescription between 1 January 2004 and 31 December 2017 were reviewed 61% women, 94% white, median age (Q1–Q3)=58 (47–66). Primary and secondary outcome measures Adverse reactions including gastrointestinal, skin, psychiatric and nervous system issues were collected from the allergy section of each patient’s medical record. Sex differences in the risk of adverse reactions due to prescribed opioids were modelled using logistic regression adjusted for age, body mass index, race and ethnicity. Results From 8457 participants (of which 449 (5.3%) reported adverse reactions), more women (6.5%) than men (3.4%) reported adverse reactions to at least one opioid (OR (95% CI)=2.3 (1.8 to 2.8), p<0.001). Women were more likely to report adverse reactions to tramadol (OR (95% CI)=2.8 (1.8 to 4.4), p<0.001) and oxycodone (OR (95% CI)=2.2 (1.7 to 2.9), p<0.001). Women were more likely to report gastrointestinal (OR (95% CI)=3.1 (2.3 to 4.3), p<0.001), skin (OR (95% CI)=2.1 (1.4 to 3.3), p=0.001) and nervous system issues (OR (95% CI)=2.3 (1.3 to 4.2), p=0.004). Conclusions These findings support the importance of sex as a biological variable to be factored into pain management studies.

Background: Prevention strategies for Alzheimer’s disease and Alzheimer’s disease related dementias (AD/ADRDs) are urgently needed for reducing incidence. Intra-patient variability in lipid levels is a potentially modifiable risk factor for incident AD/ADRD. Although regular lipid measurements are a part of common clinical practice and longitudinal data routinely available in electronic health records (EHR), research examining this association between AD/ADRD and lipid variability across multiple lipid types remains scarce. Methods: All residents living in Olmsted County, MN on 1/1/2006 age ≥60 years without an AD/ADRD diagnosis were identified using Centers for Medicare & Medicaid Services diagnostic codes. Persons with ≥3 lipid measurements (total cholesterol or triglycerides) in the 5 years prior to index date were retained. Lipid variability was quantified using variability independent of the mean (VIM). Models were adjusted for traditional risk factors. Associations between lipid variability quintiles and incident AD/ADRD were assessed using Cox proportional hazards regression. Multiple imputation was used for missing covariates. Participants were followed through 2018 for incident AD/ADRD. Results: The final analysis included data on 11,551 participants with total cholesterol and 11,502 participants with triglycerides. Participants had a mean age of 71 (range 60-99) years, and were primarily white (96%). Females (54%) were also slightly more frequent than males. Median follow up was 12.9 years (range: 0-13). Participants in the highest quintile of variability for total cholesterol and triglycerides had a 20% increased risk of incident AD/ADRD. Similar results were found in the subset with complete covariate data. Conclusion: In a large EHR derived cohort, persons in the highest quintile of lipid variation had an increased risk of incident AD/ADRD. Further studies to identify the mechanisms behind this risk factor and replication of these results across a more diverse population are needed.

Background: Larger within-patient variability of lipid levels has been associated with an increased risk of cardiovascular disease (CVD). However, measures of lipid variability are not currently used clinically. We investigated the feasibility of calculating lipid variability within a large electronic health record (EHR)-based population cohort and assessed associations with incident CVD. Methods: We identified all individuals ≥40 years of age who resided in Olmsted County, MN on 1/1/2006 (index date) without prior CVD. CVD was defined as myocardial infarction, coronary artery bypass graft surgery, percutaneous coronary intervention or stroke. Patients with ≥3 measurements of total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides during the 5 years before the index date were retained in the analyses. Lipid variability was calculated using variability independent of the mean (VIM). Patients were followed through 9/30/2017 for incident CVD (including CVD death). Cox regression was used to investigate the association between quintiles of lipid VIMs and incident CVD. Results: We identified 18,642 individuals (mean age 60; 55% female) who were free of CVD at baseline and VIM calculated for at least one lipid measurement. After adjustment, those in the highest VIM quintiles of total cholesterol had a 25% increased risk of CVD (Q5 vs. Q1 HR: 1.25, 95% CI: 1.08-1.45; Table). We observed similar results for LDL-C (Q5 vs. Q1 HR: 1.20, 95% CI: 1.04-1.39) and HDL-C (Q5 vs. Q1 HR: 1.25, 95% CI: 1.09-1.43). There was no association between triglyceride variability quintiles and CVD risk. Conclusion: In a large EHR-based population cohort, high variability in total cholesterol, HDL-C and LDL-C was associated with an increased risk of CVD, independently of traditional risk factors, suggesting it may be a target for intervention. Lipid variability can be calculated in the EHR environment but more research is needed to determine its clinical utility.

Background: Hepatocyte growth factor (HGF) is a cytokine released in response to vascular injury and a novel biomarker of cardiovascular disease (CVD) risk. However, the relationship between ideal cardiovascular health (CVH) and HGF is unknown. We examined whether ideal CVH is associated with lower HGF levels in a multi-ethnic cohort of adults free from clinical CVD at baseline. Methods: We analyzed data from the MESA study of 6,490 men and women aged 45-84 years. The independent variable was the CVH score derived from 7 metrics (smoking, body mass index, physical activity, diet, total cholesterol, blood pressure and blood glucose). Each metric was scored 0 points (poor), 1 point (intermediate) and 2 points (ideal). The total CVH score ranged from 0-14. An inadequate score was 0-8, average, 9-10 and optimal, 11-14. The number of ideal metrics was also counted. The dependent variable was logarithmically transformed HGF. We examined the association between the CVH score and HGF using linear regression models adjusted for age, sex, race/ethnicity, education, income, health insurance and study site. Results: The mean (SD) age of participants was 62 (10) years. Fifty-three percent were women. Participants with optimal CVH scores had the lowest HGF concentration [Median (IQR): 807 (678-962) pg/mL] compared to those with average [870 (740-1,036)] and inadequate scores [969 (821-1,159)]. A one-unit increment in the CVH score was significantly associated with a 3% lower HGF concentration (Table). Average and optimal CVH scores were also significantly associated with 8% and 12% lower HGF concentrations, respectively, compared to inadequate scores. Additionally, a greater number of ideal metrics was associated with lower HGF concentrations. Interactions by age, sex and race/ethnicity were not significant. Conclusion: In this ethnically diverse cohort, optimal CVH was significantly associated with lower HGF levels. Interventions aimed at promoting ideal CVH may reduce vascular injury as indicated by lower serum HGF levels.

Background: Serum levels of vascular cell adhesion molecule-1 (VCAM-1) are reflective of endothelial activation, a pathologic process that is associated with subclinical cardiac dysfunction. While VCAM-1 has been implicated in the pathogenesis of heart failure (HF) with preserved ejection fraction (HFpEF), the prospective association of VCAM-1 with clinically overt HF is unclear. Methods: In the Multi-Ethnic Study of Atherosclerosis, we evaluated the association of VCAM-1 at Exam 2 (2002-2004) with incident HF across ejection fraction (EF) categories (HFpEF and HF with reduced EF [HFrEF]) after adjustment for cardiovascular risk factors. Incident HF was adjudicated as first hospitalization for symptomatic HF, requiring specific clinical and/or imaging criteria. Results: Of 2,298 participants (mean age: 63.0 years, female: 53%), those with higher VCAM-1 were more likely white race, had higher blood pressure, and lower renal function. Over a median of 16.0 years, there were 102 HF events (HFpEF = 65; HFrEF = 37) ( Figure ). After covariate adjustment, VCAM-1 was independently associated with incident HF ( Table ). Upon evaluation of HF subtypes, VCAM-1 was associated with incident HFpEF, and risk effect estimates were consistent for incident HFrEF. The association of VCAM-1 with incident HF was consistent across the spectrum of age, sex, and BMI. Conclusion: In a multiethnic cohort, VCAM-1 was independently associated with incident HF over long-term follow up. These findings suggest a potential role for endothelial activation in driving clinical HF. Lifestyle and pharmacologic therapies that decrease endothelial activation may prevent the progression to clinical HF.

Introduction: Hepatocyte Growth Factor (HGF) is a mesenchymal cytokine linked to incident heart failure (HF), with recent data from our group showing a strong and independent association with HF with preserved ejection fraction (HFpEF). Cardiac MRI (cMRI) allows for precise analysis of morphologic changes in left ventricular (LV) structure. Increasing LV mass and concentric remodeling (defined by an increasing mass:volume ratio) are imaging markers of HFpEF risk. Whether HGF is associated with adverse LV remodeling over time is unknown. Hypothesis: Higher HGF will be associated with increasing LV mass, decreasing LV volume and increasing mass:volume ratio over 10 yrs. Methods: We studied 4762 participants of the MESA cohort, free of cardiovascular disease (CVD) and HF at baseline, who completed both HGF measurement and cMRI at baseline. Participants with LV EF<50% were excluded. Of these, 2855 completed a 2 nd cMRI at 10 yrs. We examined the cross-sectional and longitudinal associations of HGF and LV parameters using multivariable-adjusted linear mixed effect models. Results: The mean (SD) for age was 61 (10) yrs. Median (IQR) for HGF level was 888 pg/mL (745-1066); 53% women. At baseline, the 3 rd HGF tertile, compared to the 1 st , was associated with greater mass:volume ratio [relative difference 1.66 (0.43, 2.89)] and lower LV end diastolic volume [-1.87 mL (-3.45, -0.28)], after adjustment for CVD risk factors and NT- proBNP (model 2) [Table] . In longitudinal analysis, the 3 rd HGF tertile was also associated with increasing mass:volume ratio [difference in 10-yr change: 4.79 (2.73, 6.85)] and decreasing LV end diastolic volume [-4.97 (-7.10, -2.85)]. Conclusions: In a community cohort, higher HGF levels were independently associated with a concentric LV remodeling pattern of increasing mass:volume ratio and decreasing LV end diastolic volume over 10 yrs. This association may suggest an intermediate phenotype explaining the association of HGF with HFpEF risk.

Background Serum levels of vascular cell adhesion molecule‐1 (VCAM‐1) are reflective of endothelial activation. Although VCAM‐1 has been implicated in the pathogenesis of heart failure with preserved ejection fraction (HFpEF), the prospective association of VCAM‐1 with development of clinically overt heart failure (HF) across ejection fraction categories is unclear. Methods and Results In MESA (the Multi‐Ethnic Study of Atherosclerosis), we evaluated the association of VCAM‐1 at examination 2 (2002–2004) with incident HF (HFpEF and HF with reduced ejection fraction) after adjustment for cardiovascular risk factors. Incident HF was independently adjudicated as first hospitalization for symptomatic HF. Among 2297 participants (mean age, 63 years; women, 53%), those with higher VCAM‐1 were more likely to be White race, had higher blood pressure, and had lower kidney function. Over a median of 14.4 years, there were 102 HF events (HFpEF=65; HF with reduced ejection fraction=37). After covariate adjustment, each doubling of VCAM‐1 was associated with incident HF (hazard ratio [HR], 1.94; 95% CI, 1.17–3.23; P =0.01). This association appeared stronger among current/former smokers compared with never smokers. On evaluation of HF subtypes, VCAM‐1 was associated with incident HFpEF (HR, 1.97; 95% CI, 1.04–3.72; P =0.04) but not with incident HF with reduced ejection fraction, although risk estimates were consistent (HR, 1.82; 95% CI, 0.79–4.21; P =0.16). Conclusions In a multiethnic cohort, VCAM‐1 was significantly associated with incident HF over long‐term follow‐up. These findings suggest a potential role for endothelial activation in driving clinical HF, and specifically HFpEF. Therapies that decrease endothelial activation may prevent the progression from cardiovascular risk factors to clinical HF.

Abstract Background: Large genome-wide association studies (GWAS) have identified hundreds of single-nucleotide polymorphisms (SNPs) associated with increased risk of prostate cancer (PRCA). The majority of these risk loci are intergenic and the underlying causal variants are presumed to confer regulatory effects on gene expression. While follow-up eQTL studies have begun to explore putative regulatory targets of PRCA risk SNPs, these have largely focused on protein-coding genes. Recent evidence demonstrates that pri-miRNA transcripts are susceptible to the same cis-acting regulatory effects as other long RNAs, and PRCA risk SNPs may confer systemic impacts on the prostate transcriptome through miRNA expression dysregulation. However, no miRNA transcriptome-wide association studies (TWAS) of PRCA risk have been conducted to date, and the role of miRNA dysregulation in PRCA risk is poorly understood. Methods: We performed a comprehensive miRNA TWAS of PRCA risk using small RNA sequencing and genome-wide genotyping data from N = 441 normal prostate tissue samples from Mayo Clinic along with N = 411 prostate adenocarcinoma (PRAD) tumor samples from the Cancer Genome Atlas (TCGA). Both miRNA expression datasets were uniformly processed separately using the CAP-miRSeq miRNA workflow and trimmed mean of M (TMM) normalization. Expressed mature miRNAs were identified as ≥ 5 reads in &gt;25% of samples. PEER was used to remove known and latent factors on miRNA expression levels, and subsequent analyses were performed on residual values. Genetically regulated expression (GReX) prediction models were trained for all expressed miRNAs with evidence of cis-heritability (P&lt;0.01) in the respective normal and tumor datasets using the FUSION TWAS software. TWAS for PRCA risk was performed with both sets of GReX models using single-SNP summary statistics from the recent PRACTICAL consortium PRCA case-control OncoArray GWAS meta-analysis. Significant TWAS associations were declared at a false discovery rate of 0.05. Results: A total of 615 and 571 distinct expressed miRNAs were identified in the normal and tumor tissue datasets, respectively (overlap: 480). Among these, 98 (16%) normal tissue miRNAs demonstrated significant cis-heritability (median cis-h2 = 0.09, range: 0.04-0.79) for GReX model training. Similar results were obtained from TCGA tumor samples, with 63 (11%) miRNA GReX models successfully trained (median cis-h2 = 0.11, range: 0.05-0.60). Using normal tissue GReX models, we identified three significant TWAS miRNA associations with PRCA risk: over-expression of mir-941 family miRNAs (PTWAS = 3.1E-04) and reduced expression of miR-3617-5p (PTWAS = 1.7E-04) and hsa-miR-16-2-3p (PTWAS = 9.1E-04). The TCGA tumor TWAS also identified a significant association with miR-941 overexpression (PTWAS = 9.7E-04). All associated miRNAs have previously reported oncogenic or tumor-suppressive associations in other cancer types. Conclusions: We identified significant associations with PRCA genetic risk and miRNA expression dysregulation using tumor and normal prostate tissue data. Future work will examine downstream regulatory effects of miRNA dysregulation as well as miRNA-mediated risk mechanisms via competing endogenous RNA relationships. Citation Format: Nicholas B. Larson, Shannon K. McDonnell, Zachary Fogarty, Stephen N. Thibodeau. An miRNA transcriptome-wide association study of prostate cancer risk [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 54.

Abstract Copy number alterations and loss of heterozygosity (LOH) are common in ovarian tumors, including those on the X chromosome. In females, one of the two copies of the X chromosome is transcriptionally silenced during embryonic development in the epigenetic process of X chromosome inactivation. Specifically, the long non-coding RNA XIST is expressed from the inactive X (Xi) and triggers chromosome-wide DNA methylation. In breast and ovarian cancer cell line studies, Xi deletion with concomitant duplication of the active X (Xa) has been observed. We sought to examine which copy of the X chromosome, Xi or Xa, was gained or lost in a set of primary ovarian tumors, using levels of XIST expression as a marker of X chromosome inactivation. To determine which copy of the X chromosome harbors the alterations, we compared local copy number and LOH in Xq13.2 to the levels of XIST expression in a set of N=134 primary tumors from Mayo Clinic ovarian cancer patients. Associations with XIST expression were examined using ANOVA, with copy number of two and absence of LOH as the reference level. Associations between expression of XIST and expression of other genes were examined with a mixed effects model, to account for the repeated observations across patients, after adjusting for local copy number and LOH. Associations with XIST expression and overall survival were performed with Cox regression. Additional analyses were performed in a subset of 24 patient tumors with X chromosome-wide alterations to simplify interpretation, including a visual examination of chromosome-wide DNA methylation patterns. Analyses were replicated in N=106 ovarian tumors from TCGA with X chromosome-wide alterations. Across N=134 ovarian tumors from Mayo Clinic, abnormal copy number (not equal to two) was associated with reduced XIST expression (p=0.0016). Levels of XIST expression were similar for tumors with one copy and tumors with three or more copies. Similarly, levels of XIST expression were also lower for tumors with LOH (p=5.2E-5). Reduced XIST expression was also associated with poorer survival, after adjusting for age (p=0.02). Results were consistent in the subset of tumors with whole chromosome alterations (pCN=0.0013, pLOH=0.0019). For these tumors, when LOH was present, XIST expression was low and the DNA methylation distribution is bi-modal, suggesting that the inactive X was lost and the active X was gained. Furthermore, after adjusting for copy number and LOH, XIST expression was weakly positively correlated with the median expression of other genes (R2=0.04, p=0.026). In N=106 TCGA tumors with whole X chromosome alterations, XIST expression was also reduced in tumors with abnormal copy number (p=0.006) and with LOH (p=1.5E-10). Taken together, because XIST is expressed from Xi and not Xa, associations of reduced XIST expression with abnormal copy number and presence of LOH suggest that Xi is frequently lost and Xa is gained in ovarian tumors, consistent with cell line studies. These results should be further validated in larger sample sets and examined for association with clinical outcome. Citation Format: Yilin Song, Bryan M. McCauley, Melissa C. Larson, Sebastian M. Armasu, Kate Lawrenson, Ellen L. Goode, Nicholas B. Larson, Stacey J. Winham. Evidence of decreased X chromosome inactivation in primary ovarian tumors [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 2452.