Bielinski, Suzette J.

Introduction: Phenome-wide association studies (PheWAS) using electronic medical record (EMR)-linked biobanks have been used not only to identify and replicate known associations of genetic variants with disease phenotypes but have also resulted in the discovery of potentially novel genotype-phenotype relationships. The natriuretic peptide (NP) system plays an important role in a broad range of disease processes including cardiovascular and inflammatory diseases. We hypothesized that performing a PheWAS using previously known functional genetic variants of the NP system may result in novel disease associations that could provide mechanistic insights in an unbiased manner. Methods: We scanned for associations between 9 single-nucleotide polymorphisms (SNPs) in the NP system and 27 EMR-derived chronic disease phenotypes in 3,025 individuals participating in a case-control study of peripheral arterial disease. The EMR phenotypes were identified using two or more ICD-9-CM diagnosis codes based on the AHRQ Clinical Classifications Software (CCS). The relationship of SNPs and phenotypes were modeled using logistic regression adjusting for gender. Results: We identified rs5065, a SNP located in the stop codon of exon 3 of the NPPA gene, to be the strongest associated SNP with rheumatoid arthritis (RA) (OR=0.78, p=0.0008, q-value=0.11). The SNP leads to the extension of atrial natriuretic peptide (ANP) by 2 additional arginines at the C terminus. Cardiovascular disease is known to be the leading cause of death in patients with RA and ANP plays an important immunomodulatory role by inhibiting inducible nitric oxide synthase, reducing TNF-α production and attenuating prostaglandin E2 production in macrophages. Circulating NPs have been used to screen for occult cardiac disease and are associated with mortality in RA. This study demonstrates for the first time the importance of the relationship between genetic variation in the NP system and RA. Conclusions: PheWAS was successfully used as a tool to identify a novel association of functional genetic variation in the NPPA gene with RA. The observation is hypothesis generating and further replication studies are required to determine the role of rs5065 in cardiovascular outcomes of RA.

Background: Cellular adhesion molecules (CAM) have a central role in the accumulation of circulating leukocytes at sites of vascular injury, infection and/ inflammation, and have been associated with the development of atherosclerotic plaque and coronary artery disease in mature adults. Objective: To test the hypothesis that higher overall circulating CAM levels in young adults predict cardiovascular disease (CVD) events over the next 18 years. Method: We measured several circulating CAM molecules (ICAM-1, P-selectin, E-selectin and VCAM) in the Coronary Artery Risk Development in Young Adults (CARDIA) study at exam year 7 (1992-93, black and white men and women, CVD-free, mean age 32, range 25-37 years, n=2428) and monitored incident CVD events (n=70, including coronary heart disease, stroke, and heart failure, adjudicated based on medical records) through exam year 25, mean age of 50 years. We ranked each CAM in quintiles (coded 0-4) and summed the ranks across CAMs into an index to examine the association with incident CVD with Cox regression models. Results: Unadjusted cumulative CVD event rates were 1.6% (sum of CAM quartile ranks 0-8: 22 events in 1353 participants), 2.3% (sum of ranks 9-12: 29/813), and 7.3% (sum of ranks 13-16: 19/262). In proportional hazards regression analysis adjusted for year 7 age, sex, race, clinic, education, smoking, diet, physical activity, body mass index, blood pressure, blood lipids, and blood glucose, sum of ranks 13-16 were associated with a higher risk of CVD compared to the referent (rank sum 0-8) (See Table). Conclusion: High levels of circulating CAMs at an early stage of adulthood (mean age 32, range 25-37 years) were associated with an increased risk of incident CVD events. CAMs may be an early biomarker for development of subclinical CVD, even in CVD-free young adults with low atherosclerosis burden and decades prior to the development of clinical CVD.

Several loci have been discovered and replicated for urine albumin-creatinine ratio (ACR), a measure of kidney function, in European populations; however, other ethnic populations such as African-Americans have been less studied for this trait. We examined the association of genetic variants on the Metabochip with natural log transformed ACR (ln(ACR)) in 4,629 African-Americans from the Population Architecture using Genomics and Epidemiology (PAGE) Study which includes participants from the Atherosclerosis Risk in Communities (ARIC), Cardiovascular Health Study (CHS), Coronary Artery Risk Development in Young Adults (CARDIA), and Epidemiologic Architecture for Genes Linked to the Environment (EAGLE-BioVU) studies. We then replicated our findings in 1,533 African-American participants from the Multi-Ethnic Study of Atherosclerosis (MESA). Linear regression models were used to examine associations, assuming additive effects of the SNPs. Within each study, models were adjusted for age, field center (where applicable), and principal components of ancestry, and stratified by sex. Fixed effects inverse variance weighted meta-analysis was used to combine results from contributing studies. Array-wide statistical significance was concluded for SNP associations of p-value less than 2.8 x 10-7. Mean ln(ACR) ranged from 1.4 to 2.9 mg/g, average age ranged from 44 to 77 years, and percent female participants ranged from 59-65% in the contributing studies. Four SNPs were significantly associated array-wide with ACR in the discovery sample: rs17399841 in GRHL2 (β=-0.67, MAF=0.01, p=2.7x10-9), rs153914 near SPATA9/RHOBTB3 (β=-0.59, MAF=0.01, p=1.0x10-7), rs17439845 in LEF1 (β=-0.49, MAF=0.02, p=1.3x10-7), and rs11836243 near ZNF641/ANP32D (β=0.23, MAF=0.12, p=1.8x10-7). Of these four loci, only the association of rs17439845 with ACR replicated in MESA (β=-0.36, p=0.038), with a combined discovery-replication beta of -0.46 and p value of 1.8x10-8. Combined discovery-replication for rs17399841, rs153914, and rs11836243 no longer achieved array-wide significance (p-values 9.84x10-7, 1.98x10-5, and 1.32x10-4, respectively). LEF1 is a mediator in the Wnt signaling pathway, and has been previously implicated in systemic lupus erythematosus and osteoporosis, as well as previously associated with waist circumference, echocardiographic traits, bone mineral density, and blood count traits. Continuing studies in diverse populations are warranted to identify a set of loci important in contributing to the variation in ACR.

Introduction: Proton pump inhibitors (PPIs) are used by an estimated 29 million Americans. PPIs increase the levels of asymmetrical dimethylarginine, a known risk factor for cardiovascular disease (CVD). Data from a select population of patients with CVD suggest that PPI use is associated with an increased risk of stroke, heart failure, and coronary heart disease. The impact of PPI use on incident CVD is largely unknown in the general population. Hypothesis: We hypothesized that PPI users have a higher risk of incident total CVD, coronary heart disease, stroke, and heart failure compared to nonusers. To demonstrate specificity of association, we additionally hypothesized that there is not an association between use of H 2 -blockers - another commonly used class of medications with similar indications as PPIs - and CVD. Methods: We used the Rochester Epidemiology Project’s medical records-linkage system to identify all residents of Olmsted County, MN on our baseline date of January 1, 2004 (N=140217). We excluded persons who did not grant permission for their records to be used for research, were <18 years old, had a history of CVD, had missing data for any variable included in our model, or had evidence of PPI use within the previous year.We followed our final cohort (N=58175) for up to 12 years. The administrative censoring date for CVD was 1/20/2014, for coronary heart disease was 8/3/2016, for stroke was 9/9/2016, and for heart failure was 1/20/2014. Time-varying PPI ever-use was ascertained using 1) natural language processing to capture unstructured text from the electronic health record, and 2) outpatient prescriptions. An incident CVD event was defined as the first occurrence of 1) validated heart failure, 2) validated coronary heart disease, or 3) stroke, defined using diagnostic codes only. As a secondary analysis, we calculated the association between time-varying H 2 -blocker ever-use and CVD among persons not using H 2 -blockers at baseline. Results: After adjustment for age, sex, race, education, hypertension, hyperlipidemia, diabetes, and body-mass-index, PPI use was associated with an approximately 50% higher risk of CVD (hazard ratio [95% CI]: 1.51 [1.37-1.67]; 2187 CVD events), stroke (hazard ratio [95% CI]: 1.49 [1.35-1.65]; 1928 stroke events), and heart failure (hazard ratio [95% CI]: 1.56 [1.23-1.97]; 353 heart failure events) compared to nonusers. Users of PPIs had a 35% greater risk of coronary heart disease than nonusers (95% CI: 1.13-1.61; 626 coronary heart disease events). Use of H 2 -blockers was also associated with a higher risk of CVD (adjusted hazard ratio [95% CI]: 1.23 [1.08-1.41]; 2331 CVD events). Conclusions: PPI use is associated with a higher risk of CVD, coronary heart disease, stroke and heart failure. Use of a drug with no known cardiac toxicity - H 2 -blockers - was also associated with a greater risk of CVD, warranting further study.

Introduction: Increased levels of hepatocyte growth factor (HGF), active in cell growth, motility, and morphogenesis, are associated with the presence of obesity, poor metabolic health, and cardiovascular disease. Hypothesis: We assessed the hypothesis that higher baseline levels of HGF will be associated with increased risk of diabetes. Methods: We examined the association between HGF and incident diabetes in MESA, including 5395 men and women 45-84 years of age at enrollment (2000-02). Fasting serum HGF was measured at baseline and on a subsample of participants at exam 2 (n = 1915). From 2000-11, incidence of diabetes was ascertained over 4 follow-up examinations, determined by new use of insulin or oral hypoglycemic medication or fasting glucose ≥ 126 mg/dL. Cox regression was used to estimate hazard ratios (HR) for incident diabetes according to 1 standard deviation unit (SDU) of HGF (1 SDU =256 pg/mL), before and after adjustment for age, sex, race/ethnicity, education, study center, smoking status, alcohol consumption, BMI, WC, fasting glucose and insulin, CRP, and IL-6 levels. Similarly, hazard ratios for incident diabetes were estimated according to change in HGF levels from exam 1 to exam 2 in the subsample. Results: At baseline, older age, male sex, current smoking, and higher body mass index (BMI), waist circumference (WC), fasting glucose and insulin, C-reactive protein (CRP) and interleukin-6 (IL-6) levels were all associated with higher levels of HGF, while greater education and physical activity were associated with lower serum HGF. Incidence of diabetes in this analytic sample was 12% (n cases = 670). Per 1 SDU increase in baseline HGF level, unadjusted risk for diabetes increased 1.46 fold (95% CI=1.37, 1.56). After adjustment, diabetes risk per 1 SDU increase in HGF was attenuated but remained significantly increased (HR=1.22; 95% CI=1.12, 1.32). No association was found between change in HGF level between exam 1 and exam 2 and incidence of diabetes. There was no evidence of effect modification by race/ethnicity for either analysis. Conclusion: In conclusion, in this ethnically diverse U.S. adult population, higher levels of serum HGF were independently associated with increased incidence of diabetes.

Background: Given the global burden of peripheral arterial disease (PAD) and the associated morbidity, continuing investigation of potential biomarkers as therapeutic targets for preventing PAD is warranted. The association of P-selectin, an endothelial cell adhesion molecule, with the ankle brachial index (ABI) and PAD is not well established. Thus, we examined the association of circulating P-selectin with prevalent and incident PAD, the ABI, and change in the ABI. Methods: The Multi-Ethnic Study of Atherosclerosis (MESA) is a prospective population-based cohort including 6814 non-Hispanic white, African American, Hispanic and Chinese men and women aged 45-84 at baseline (2000-02). For this analysis, 5700 participants with both P-selectin at Exam 2 (2002-04) and the ABI at Exam 3 (2004-05) were available, with change in ABI from Exam 3 to Exam 5 (2010-12) available on 4276 participants. P-selectin was modeled per standard deviation (SD), while prevalent and incident PAD were defined as ABI ≤ 0.90 after excluding those with ABI>1.4. Linear regression was used for ABI, growth curve models for change in ABI, and logistic regression for prevalent and incident PAD. Results: In models adjusted for demographics, lifestyle factors, comorbidities, blood pressure, lipids, fasting glucose, c-reactive protein (CRP), interleukin-6 (IL-6), D-dimer, and fibrinogen, each SD (13 ng/mL) higher P-selectin was significantly associated with 0.007 lower ABI (95% CI ((-0.011, -0.004)), p<0.001), and an average change in the ABI of -0.006 ((-0.010, -0.003, p<0.001). P-selectin was also significantly associated with a 1.14-fold greater odds of prevalent PAD ((1.00, 1.28), p=0.04), and a 1.30-fold greater odds of incident PAD ((1.10, 1.55), p<0.001). There were no significant interactions of p-selectin by race/ethnicity or gender. Addition of P-selectin to models containing traditional PAD risk factors as well as CRP, IL-6, D-dimer and fibrinogen did not improve the c-statistic or net reclassification for incident PAD. Conclusions: P-selectin appears to be strongly associated with the development of PAD. However, further research is needed in population-based studies to confirm prospective associations of p-selectin with incident PAD and change in the ABI.

PR interval (PR) prolongation is a predictor of atrial fibrillation, a common cardiac arrhythmia, and is an established risk factor for pacemaker implantation, heart failure, stroke, and all-cause mortality. Previous genome-wide association studies in European (EU) and African (AA) descent populations have identified multiple loci associated with PR. However, it is unclear whether these loci are relevant in other racial groups, including Hispanic/Latinos (HL). Extending studies to include admixed populations such as HL also provides the potential to narrow regions associated with PR and identify population specific variants. Therefore, we evaluated 1,135 SNPs from four previously identified PR loci ( SCN5A, SCN10A, CAV1-CAV2, and TBX5-TBX3 ), fine-mapped on the Illumina Metabochip for association with PR. Associations were examined using linear regression assuming an additive genetic model and adjusting for global ancestry and clinical covariates (age, sex, RR interval, body mass index, height, systolic blood pressure) in 11,800 HL participants from participating studies of the Population Architecture using Genomics and Epidemiology (PAGE) consortia. Three of the four loci ( SCN5A, SCN10A, and CAV1-CAV2 ) generalized to HL populations ( P <7.1E-4 based on the number of SNPs correlated at r 2 >0.2 with EU-identified index SNPs). At each of these loci, we identified SNPs with stronger evidence of association with PR than the EU-identified index SNP when examined in the HL population ( SCN5A : rs7374540; SCN10A : rs6801957; CAV-CAV2 : rs717957). Conditional analyses also identified two additional independent SNPs at each of the three generalized loci ( SCN5A : rs3796387, rs9861242; SCN10A : rs10428132, rs7433723; CAV1-CAV2 : rs3801995, rs3807994), representing potential secondary signals. Furthermore, linkage disequilibrium (LD) patterns in Hispanics enabled the narrowing of the regions flanking the EU-identified index SNPs. For example, at the CAV1-CAV2 locus, HL LD patterns indicated that seven SNPs were in LD with rs3807989, the EU-identified index SNP, compared to 17 SNPs that were in LD with rs3807989 when applying EU LD patterns, reducing the size of the region by 71 kilobases. Finally, we examined evidence for population specific signals outside of the specific regions identified in EU populations which were associated with PR in Hispanics by focusing on SNPs uncorrelated with the EU index signals. At all four loci, multiple SNPs had evidence of association independent of the signals examined based on EU results. In conclusion, the same genetic loci that influence PR in European descent populations are also important in Hispanic populations, although we detected evidence for population specific SNPs. These results emphasize the importance of examining genetic associations in diverse populations in order to narrow the genomic interval for future functional interrogation.

Incident hypertension varies across races/ethnicities. Mechanisms implicated in the pathogenesis are endothelial dysfunction with reduction of nitric oxide and alterations of the Renin-Angiotensin-Aldosterone System (RAAS). Stromal Derived Factor 1 (SDF1) and Transforming Growth Factor β1 (TGFβ1) are essential in these processes. SDF1 is a chemotactic protein that regulates migration of endothelial progenitor cells from the bone marrow to the arterial wall to repair endothelial damage; TGFβ1 promotes extra-cellular matrix deposition within arterial vessels, inhibits nitric oxide production and stimulates renin release. In order to investigate the association of circulating levels of these proteins with incident hypertension across different ethnicities, an ethnic stratified random sample of 2,880 participants from the Multi-Ethnic Study of Atherosclerosis (MESA) was used. Of those, 988 participants had available plasma samples and were without hypertension, defined as blood pressure ≥ 140/90 mmHg, or current use of antihypertensive medications, at baseline. During a median follow up of 8.5 years, 388 new cases of hypertension (101 in non-Hispanic whites, 94 in Chinese, 101 in African Americans and 92 in Hispanics) were observed. Time to hypertension incidence was estimated using Cox proportional hazards regression. Hazard Ratios were adjusted for age, sex, diabetes, LDL, HDL, triglycerides, BMI, smoking status and current alcohol use. The race/ethnicity by protein interaction p-values were 0.055 and 0.066 for TGFβ1 and SDF1 respectively, suggesting heterogeneity across races/ethnicities. In Chinese, a difference of one standard deviation of plasma TGFβ1 was associated with a 22% increase in incident hypertension (HR (95% CI) 1.22 (1.02, 1.47), p-value=0.03); and a difference of one standard deviation of plasma SDF1 was associated with a 31% decrease in incident hypertension (HR (95% CI) 0.69 (0.52, 0.91), p-value 0.008). No significant association was found between these proteins and incident hypertension among non-Hispanic whites, African or Hispanic Americans. In conclusion, soluble markers of endothelial integrity and of RAAS activation were associated with the incidence of hypertension among Chinese Americans. These results suggest heterogeneity in the pathophysiology of hypertension according to race/ethnicity.

L-selectin is an adhesion molecule constitutively expressed on leukocytes and mediating their interaction with the endothelium. As a member of the selectin family, it plays a crucial role in leukocyte adhesion to and migration through the endothelial wall during inflammation. A protective role of L-selectin in early atherosclerosis was demonstrated in animal models. Previous studies assessing the relationship of this protein with cardiovascular disease were inconsistent. Lower serum L-selectin levels were found in patients with coronary artery disease (CHD) and peripheral artery disease (PAD) compared to controls. In contrast, plasma levels of L-selectin were increased in patients with CHD and PAD. These inconsistencies might be due to an interaction of soluble L-selectin with clotting factors or anticoagulants present in plasma but not in serum. Also, the small number and Caucasian predominance of patients included in previous studies represent limitations that need to be addressed using a rigorous, population-based approach. Herein, we seek to investigate the association of serum L-selectin with subclinical outcomes in a large, ethnically heterogeneous population. In a randomly selected subset of 720 subjects for each of the four ethnicities included in the Multi-Ethnic Study of Atherosclerosis (MESA), serum samples were available for 2403 individuals (579 African Americans, 620 Caucasians, 600 Chinese and 604 Hispanics). Serum L-selectin levels and intima-media thickness (IMT) were measured at exam 2 and CAC was measured at exam 2 or 3 during the cohort follow up. Regression analysis and the Tobit model were used to assess the association of L-selectin levels with IMT and CAC respectively. The mean age was 63±10 years and 47% were men. The mean±SD protein levels were significantly different across ethnicities (867±192 ng/mL in African Americans, 956±215 ng/mL in Caucasians, 834±174 ng/mL in Chinese and 904±195 ng/mL in Hispanics; p<0.0001). In the univariate analysis, age was inversely correlated to the protein levels; in addition, levels were lower in males, regardless of ethnicity and in smokers, except among Hispanics. L-selectin was associated with modestly lower CAC only among Chinese (p=0.030), the association persisted after adjustment for traditional risk factors (p=0.046). No significant association was demontrated between the protein levels and IMT. In conclusion, in this large diverse population, we found that serum L-selectin was modestly associated to measures of subclinical atherosclerosis in Chinese but not in African, European or Hispanic Americans.

Pulse pressure (PP), defined as the difference between systolic and diastolic blood pressure, is considered a marker of large artery stiffness and has been shown to independently predict cardiovascular disease. Differences among ethnicities with regard to PP have been previously described, yet the underlying biological mechanisms are poorly understood. In addition, biomarkers of inflammation and matrix remodeling have been shown to be involved in the development of arterial stiffness. The aim of this study was to investigate the role of adhesion pathway components in the development of arterial stiffness measured by PP among four ethnicities. E-cadherin, HGF, L-selectin, P-selectin, 6-Ckine, MMP1, MMP2, RANTES, SDF1, SLPI, TGFβ1, TIMP2 and VCAM were measured in a random sample of 2402 participants (579 African Americans, 619 Caucasians, 600 Chinese and 604 Hispanics) in the Multi-Ethnic Study of Atherosclerosis (MESA) study at Exam 2. PP was measured concurrently and its association with each adhesion protein level was assessed via linear regression in each ethnicity. All analyses were adjusted for age, sex, and traditional cardiovascular risk factors. Mean PP was significantly different among ethnicities: 55.9±18.0 mmHg in African Americans, 52.1± 16.1 mmHg in Caucasians, 51.4±16.4 mmHg in Chinese, and 55.0±17.5 in Hispanics, p<0.001. MMP2 was positively associated with PP among African Americans (β=2.3 (0.6), p<0.001) and Caucasians (β=1.7 (0.5), p=0.001); TIMP2 among Caucasians (β=1.5 (0.5), p=0.002) and to a lesser extent among African Americans (β=1.2 (0.6) p=0.050) and Chinese (β=1.3 (0.6), p=0.034). SDF1 was associated to PP in African Americans (β=1.4 (0.6), p=0.017) and in Chinese (β=1.1 (0.5), p=0.042); finally an association existed between TGFβ1 and PP in C (β=1.1 (0.5), p=0.025). None of the studied proteins was associated with PP in Hispanics. In conclusion, components of the adhesion pathway were correlated to PP in this large multi-ethnic cohort, possibly giving new insight on the pathophysiology of arterial stiffness. In addition, heterogeneity by ethnicity suggests different mechanisms in the development of arterial stiffness.