Under-specification as the source of ambiguity and vagueness in narrative phenotype algorithm definitions.

INTRODUCTION: Currently, one of the commonly used methods for disseminating electronic health record (EHR)-based phenotype algorithms is providing a narrative description of the algorithm logic, often accompanied by flowcharts. A challenge with this mode of dissemination is the potential for under-specification in the algorithm definition, which leads to ambiguity and vagueness. METHODS: This study examines incidents of under-specification that occurred during the implementation of 34 narrative phenotyping algorithms in the electronic Medical Record and Genomics (eMERGE) network. We reviewed the online communication history between algorithm developers and implementers within the Phenotype Knowledge Base (PheKB) platform, where questions could be raised and answered regarding the intended implementation of a phenotype algorithm. RESULTS: We developed a taxonomy of under-specification categories via an iterative review process between two groups of annotators. Under-specifications that lead to ambiguity and vagueness were consistently found across narrative phenotype algorithms developed by all involved eMERGE sites. DISCUSSION AND CONCLUSION: Our findings highlight that under-specification is an impediment to the accuracy and efficiency of the implementation of current narrative phenotyping algorithms, and we propose approaches for mitigating these issues and improved methods for disseminating EHR phenotyping algorithms.

Geometric Features Associated with Middle Cerebral Artery Bifurcation Aneurysm Formation: A Matched Case-Control Study.

OBJECTIVES: The pathogenesis of intracranial aneurysms is multifactorial and includes genetic, environmental, and anatomic influences. We aimed to identify image-based morphological parameters that were associated with middle cerebral artery (MCA) bifurcation aneurysms. MATERIALS AND METHODS: We evaluated three-dimensional morphological parameters obtained from CT angiography (CTA) or digital subtraction angiography (DSA) from 317 patients with unilateral MCA bifurcation aneurysms diagnosed at the Brigham and Women's Hospital and Massachusetts General Hospital between 1990 and 2016. We chose the contralateral unaffected MCA bifurcation as the control group, in order to control for genetic and environmental risk factors. Diameters and angles of surrounding parent and daughter vessels of 634 MCAs were examined. RESULTS: Univariable and multivariable statistical analyses were performed to determine statistical significance. Sensitivity analyses with smaller (≤ 3 mm) aneurysms only and with angles excluded, were also performed. In a multivariable conditional logistic regression model we showed that smaller diameter size ratio (OR 0.0004, 95% CI 0.0001-0.15), larger daughter-daughter angles (OR 1.08, 95% CI 1.06-1.11) and larger parent-daughter angle ratios (OR 4.24, 95% CI 1.77-10.16) were significantly associated with MCA aneurysm presence after correcting for other variables. In order to account for possible changes to the vasculature by the aneurysm, a subgroup analysis of small aneurysms (≤ 3 mm) was performed and showed that the results were similar. CONCLUSIONS: Easily measurable morphological parameters of the surrounding vasculature of the MCA may provide objective metrics to assess MCA aneurysm formation risk in high-risk patients.

Use of Narrative Concepts in Electronic Health Records to Validate Associations Between Genetic Factors and Response to Treatment of Inflammatory Bowel Diseases.

Crohn's disease (CD) and ulcerative colitis (UC) are heterogeneous. With availability of therapeutic classes with distinct immunologic mechanisms of action, it has become imperative to identify markers that predict likelihood of response to each drug class. However, robust development of such tools has been challenging because of need for large prospective cohorts with systematic and careful assessment of treatment response using validated indices. Most hospitals in the United States use electronic health records (EHRs) that warehouse a large amount of narrative (free-text) and codified (administrative) data generated during routine clinical care. These data have been used to construct virtual disease cohorts for epidemiologic research as well as for defining genetic basis of disease states or discrete laboratory values.1-3 Whether EHR-based data can be used to validate genetic associations for more nuanced outcomes such as treatment response has not been examined previously.

Making work visible for electronic phenotype implementation: Lessons learned from the eMERGE network.

BACKGROUND: Implementation of phenotype algorithms requires phenotype engineers to interpret human-readable algorithms and translate the description (text and flowcharts) into computable phenotypes - a process that can be labor intensive and error prone. To address the critical need for reducing the implementation efforts, it is important to develop portable algorithms. METHODS: We conducted a retrospective analysis of phenotype algorithms developed in the Electronic Medical Records and Genomics (eMERGE) network and identified common customization tasks required for implementation. A novel scoring system was developed to quantify portability from three aspects: Knowledge conversion, clause Interpretation, and Programming (KIP). Tasks were grouped into twenty representative categories. Experienced phenotype engineers were asked to estimate the average time spent on each category and evaluate time saving enabled by a common data model (CDM), specifically the Observational Medical Outcomes Partnership (OMOP) model, for each category. RESULTS: A total of 485 distinct clauses (phenotype criteria) were identified from 55 phenotype algorithms, corresponding to 1153 customization tasks. In addition to 25 non-phenotype-specific tasks, 46 tasks are related to interpretation, 613 tasks are related to knowledge conversion, and 469 tasks are related to programming. A score between 0 and 2 (0 for easy, 1 for moderate, and 2 for difficult portability) is assigned for each aspect, yielding a total KIP score range of 0 to 6. The average clause-wise KIP score to reflect portability is 1.37 ±â€¯1.38. Specifically, the average knowledge (K) score is 0.64 ±â€¯0.66, interpretation (I) score is 0.33 ±â€¯0.55, and programming (P) score is 0.40 ±â€¯0.64. 5% of the categories can be completed within one hour (median). 70% of the categories take from days to months to complete. The OMOP model can assist with vocabulary mapping tasks. CONCLUSION: This study presents firsthand knowledge of the substantial implementation efforts in phenotyping and introduces a novel metric (KIP) to measure portability of phenotype algorithms for quantifying such efforts across the eMERGE Network. Phenotype developers are encouraged to analyze and optimize the portability in regards to knowledge, interpretation and programming. CDMs can be used to improve the portability for some 'knowledge-oriented' tasks.

Facilitating phenotype transfer using a common data model.

BACKGROUND: Implementing clinical phenotypes across a network is labor intensive and potentially error prone. Use of a common data model may facilitate the process. METHODS: Electronic Medical Records and Genomics (eMERGE) sites implemented the Observational Health Data Sciences and Informatics (OHDSI) Observational Medical Outcomes Partnership (OMOP) Common Data Model across their electronic health record (EHR)-linked DNA biobanks. Two previously implemented eMERGE phenotypes were converted to OMOP and implemented across the network. RESULTS: It was feasible to implement the common data model across sites, with laboratory data producing the greatest challenge due to local encoding. Sites were then able to execute the OMOP phenotype in less than one day, as opposed to weeks of effort to manually implement an eMERGE phenotype in their bespoke research EHR databases. Of the sites that could compare the current OMOP phenotype implementation with the original eMERGE phenotype implementation, specific agreement ranged from 100% to 43%, with disagreements due to the original phenotype, the OMOP phenotype, changes in data, and issues in the databases. Using the OMOP query as a standard comparison revealed differences in the original implementations despite starting from the same definitions, code lists, flowcharts, and pseudocode. CONCLUSION: Using a common data model can dramatically speed phenotype implementation at the cost of having to populate that data model, though this will produce a net benefit as the number of phenotype implementations increases. Inconsistencies among the implementations of the original queries point to a potential benefit of using a common data model so that actual phenotype code and logic can be shared, mitigating human error in reinterpretation of a narrative phenotype definition.

Elevated International Normalized Ratio Is Associated With Ruptured Aneurysms.

Background and Purpose- The effects of anticoagulation therapy and elevated international normalized ratio (INR) values on the risk of aneurysmal subarachnoid hemorrhage are unknown. We aimed to investigate the association between anticoagulation therapy, elevated INR values, and rupture of intracranial aneurysms. Methods- We conducted a case-control study of 4696 patients with 6403 intracranial aneurysms, including 1198 prospective patients, diagnosed at the Massachusetts General Hospital and the Brigham and Women's Hospital between 1990 and 2016 who were on no anticoagulant therapy or on warfarin for anticoagulation. Patients were divided into ruptured and nonruptured groups. Univariable and multivariable logistic regression analyses were performed to evaluate the association of anticoagulation therapy, INR values, and presentation with a ruptured intracranial aneurysm, taking into account the interaction between anticoagulant use and INR. Inverse probability weighting using propensity scores was used to minimize differences in baseline demographics characteristics. The marginal effects of anticoagulant use on rupture risk stratified by INR values were calculated. Results- In unweighted and weighted multivariable analyses, elevated INR values were significantly associated with rupture status among patients who were not anticoagulated (unweighted odds ratio, 22.78; 95% CI, 10.85-47.81 and weighted odds ratio, 28.16; 95% CI, 12.44-63.77). In anticoagulated patients, warfarin use interacts significantly with INR when INR ≥1.2 by decreasing the effects of INR on rupture risk. Conclusions- INR elevation is associated with intracranial aneurysm rupture, but the effects may be moderated by warfarin. INR values should, therefore, be taken into consideration when counseling patients with intracranial aneurysms.

Antihyperglycemic Agents Are Inversely Associated With Intracranial Aneurysm Rupture.

BACKGROUND AND PURPOSE: Previous studies have suggested a protective effect of diabetes mellitus on aneurysmal subarachnoid hemorrhage risk. However, reports are inconsistent, and objective measures of hyperglycemia in these studies are lacking. Our aim was to investigate the association between aneurysmal subarachnoid hemorrhage and antihyperglycemic agent use and glycated hemoglobin levels. METHODS: The medical records of 4701 patients with 6411 intracranial aneurysms, including 1201 prospective patients, diagnosed at the Massachusetts General Hospital and Brigham and Women's Hospital between 1990 and 2016 were reviewed and analyzed. Patients were separated into ruptured and nonruptured groups. Univariate and multivariate logistic regression analyses were performed to determine the association between aneurysmal subarachnoid hemorrhage and antihyperglycemic agents and glycated hemoglobin levels. Propensity score weighting was used to account for selection bias. RESULTS: In both unweighted and weighted multivariate analysis, antihyperglycemic agent use was inversely and significantly associated with ruptured aneurysms (unweighted odds ratio, 0.58; 95% confidence interval, 0.39-0.87; weighted odds ratio, 0.57; 95% confidence interval, 0.34-0.96). In contrast, glycated hemoglobin levels were not significantly associated with rupture status. CONCLUSIONS: Antihyperglycemic agent use rather than hyperglycemia is associated with decreased risk of aneurysmal subarachnoid hemorrhage, suggesting a possible protective effect of glucose-lowering agents in the pathogenesis of aneurysm rupture.

Low Serum Calcium and Magnesium Levels and Rupture of Intracranial Aneurysms.

BACKGROUND AND PURPOSE: Both low serum calcium and magnesium levels have been associated with the extent of bleeding in patients with intracerebral hemorrhage, suggesting hypocalcemia- and hypomagnesemia-induced coagulopathy as a possible underlying mechanism. We hypothesized that serum albumin-corrected total calcium and magnesium levels are associated with ruptured intracranial aneurysms. METHODS: The medical records of 4701 patients, including 1201 prospective patients, diagnosed at the Brigham and Women's Hospital and Massachusetts General Hospital between 1990 and 2016 were reviewed and analyzed. One thousand two hundred seventy-five patients had available serum calcium, magnesium, and albumin values within 1 day of diagnosis. Individuals were divided into cases with ruptured aneurysms and controls with unruptured aneurysms. Univariable and multivariable logistic regression analyses were performed to determine the association between serum albumin-corrected total calcium and magnesium levels and ruptured aneurysms. RESULTS: In multivariable analysis, both albumin-corrected calcium (odds ratio, 0.33; 95% confidence interval, 0.27-0.40) and magnesium (odds ratio, 0.40; 95% confidence interval, 0.28-0.55) were significantly and inversely associated with ruptured intracranial aneurysms. CONCLUSIONS: In this large case-control study, hypocalcemia and hypomagnesemia at diagnosis were significantly associated with ruptured aneurysms. Impaired hemostasis caused by hypocalcemia and hypomagnesemia may explain this association.

Lipid-Lowering Agents and High HDL (High-Density Lipoprotein) Are Inversely Associated With Intracranial Aneurysm Rupture.

BACKGROUND AND PURPOSE: Growing evidence from experimental animal models and clinical studies suggests the protective effect of statin use against rupture of intracranial aneurysms; however, results from large studies detailing the relationship between intracranial aneurysm rupture and total cholesterol, HDL (high-density lipoprotein), LDL (low-density lipoprotein), and lipid-lowering agent use are lacking. METHODS: The medical records of 4701 patients with 6411 intracranial aneurysms diagnosed at the Massachusetts General Hospital and the Brigham and Women's Hospital between 1990 and 2016 were reviewed and analyzed. Patients were separated into ruptured and nonruptured groups. Univariable and multivariable logistic regression analyses were performed to determine the effects of lipids (total cholesterol, LDL, and HDL) and lipid-lowering medications on intracranial aneurysm rupture risk. Propensity score weighting was used to account for differences in baseline characteristics of the cohorts. RESULTS: Lipid-lowering agent use was significantly inversely associated with rupture status (odds ratio, 0.58; 95% confidence interval, 0.47-0.71). In a subgroup analysis of complete cases that includes both lipid-lowering agent use and lipid values, higher HDL levels (odds ratio, 0.95; 95% confidence interval, 0.93-0.98) and lipid-lowering agent use (odds ratio, 0.41; 95% confidence interval, 0.23-0.73) were both significantly and inversely associated with rupture status, whereas total cholesterol and LDL levels were not significant. A monotonic exposure-response curve between HDL levels and risk of aneurysmal rupture was obtained. CONCLUSIONS: Higher HDL values and the use of lipid-lowering agents are significantly inversely associated with ruptured intracranial aneurysms.

Natural Language Processing Accurately Calculates Adenoma and Sessile Serrated Polyp Detection Rates.

BACKGROUND: ADR is a widely used colonoscopy quality indicator. Calculation of ADR is labor-intensive and cumbersome using current electronic medical databases. Natural language processing (NLP) is a method used to extract meaning from unstructured or free text data. AIMS: (1) To develop and validate an accurate automated process for calculation of adenoma detection rate (ADR) and serrated polyp detection rate (SDR) on data stored in widely used electronic health record systems, specifically Epic electronic health record system, Provation® endoscopy reporting system, and Sunquest PowerPath pathology reporting system. METHODS: Screening colonoscopies performed between June 2010 and August 2015 were identified using the Provation® reporting tool. An NLP pipeline was developed to identify adenomas and sessile serrated polyps (SSPs) on pathology reports corresponding to these colonoscopy reports. The pipeline was validated using a manual search. Precision, recall, and effectiveness of the natural language processing pipeline were calculated. ADR and SDR were then calculated. RESULTS: We identified 8032 screening colonoscopies that were linked to 3821 pathology reports (47.6%). The NLP pipeline had an accuracy of 100% for adenomas and 100% for SSPs. Mean total ADR was 29.3% (range 14.7-53.3%); mean male ADR was 35.7% (range 19.7-62.9%); and mean female ADR was 24.9% (range 9.1-51.0%). Mean total SDR was 4.0% (0-9.6%). CONCLUSIONS: We developed and validated an NLP pipeline that accurately and automatically calculates ADRs and SDRs using data stored in Epic, Provation® and Sunquest PowerPath. This NLP pipeline can be used to evaluate colonoscopy quality parameters at both individual and practice levels.

Statin Use Is Associated With Reduced Risk of Colorectal Cancer in Patients With Inflammatory Bowel Diseases.

BACKGROUND & AIMS: Inflammatory bowel diseases (IBDs) such as Crohn's disease and ulcerative colitis are associated with an increased risk of colorectal cancer (CRC). Chemopreventive strategies have produced weak or inconsistent results. Statins have been associated inversely with sporadic CRC. We examined their role as chemopreventive agents in patients with IBD. METHODS: We collected data from 11,001 patients with IBD receiving care at hospitals in the Greater Boston metropolitan area from 1998 through 2010. Diagnoses of CRC were determined using validated International Classification of Diseases, 9th revision, Clinical Modification codes. Statin use before diagnosis was assessed through analysis of electronic prescriptions. We performed multivariate logistic regression analyses, adjusting for potential confounders including primary sclerosing cholangitis, smoking, increased levels of inflammation markers, and CRC screening practices to identify an independent association between statin use and CRC. We performed sensitivity analyses using propensity score adjustment and variation in the definition of statin use. RESULTS: In our cohort, 1376 of the patients (12.5%) received 1 or more prescriptions for a statin. Patients using statins were more likely to be older, male, white, smokers, and have greater comorbidity than nonusers. Over a follow-up period of 9 years, 2% of statin users developed CRC compared with 3% of nonusers (age-adjusted odds ratio, 0.35; 95% confidence interval, 0.24-0.53). On multivariate analysis, statin use remained independently and inversely associated with CRC (odds ratio, 0.42; 95% confidence interval, 0.28-0.62). Our findings were robust on a variety of sensitivity and subgroup analyses. CONCLUSIONS: Statin use was associated inversely with the risk of CRC in a large IBD cohort. Prospective studies on the role of statins as chemopreventive agents are warranted.

Colonoscopy is associated with a reduced risk for colon cancer and mortality in patients with inflammatory bowel diseases.

BACKGROUND & AIMS: Crohn's disease and ulcerative colitis are associated with an increased risk of colorectal cancer (CRC). Surveillance colonoscopy is recommended at 2- to 3-year intervals beginning 8 years after diagnosis of inflammatory bowel disease (IBD). However, there have been no reports of whether colonoscopy examination reduces the risk for CRC in patients with IBD. METHODS: In a retrospective study, we analyzed data from 6823 patients with IBD (2764 with a recent colonoscopy, 4059 without a recent colonoscopy) seen and followed up for at least 3 years at 2 tertiary referral hospitals in Boston, Massachusetts. The primary outcome was diagnosis of CRC. We examined the proportion of patients undergoing a colonoscopy within 36 months before a diagnosis of CRC or at the end of the follow-up period, excluding colonoscopies performed within 6 months before a diagnosis of CRC, to avoid inclusion of prevalent cancers. Multivariate logistic regression was performed, adjusting for plausible confounders. RESULTS: A total of 154 patients developed CRC. The incidence of CRC among patients without a recent colonoscopy (2.7%) was significantly higher than among patients with a recent colonoscopy (1.6%) (odds ratio [OR], 0.56; 95% confidence interval [CI], 0.39-0.80). This difference persisted in multivariate analysis (OR, 0.65; 95% CI, 0.45-0.93) and was robust when adjusted for a range of assumptions in sensitivity analyses. Among patients with CRC, a colonoscopy within 6 to 36 months before diagnosis was associated with a reduced mortality rate (OR, 0.34; 95% CI, 0.12-0.95). CONCLUSIONS: Recent colonoscopy (within 36 months) is associated with a reduced incidence of CRC in patients with IBD, and lower mortality rates in those diagnosed with CRC.

Identification of subjects with polycystic ovary syndrome using electronic health records.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a heterogeneous disorder because of the variable criteria used for diagnosis. Therefore, International Classification of Diseases 9 (ICD-9) codes may not accurately capture the diagnostic criteria necessary for large scale PCOS identification. We hypothesized that use of electronic medical records text and data would more specifically capture PCOS subjects. METHODS: Subjects with PCOS were identified in the Partners Healthcare Research Patients Data Registry by searching for the term "polycystic ovary syndrome" using natural language processing (n = 24,930). A training subset of 199 identified charts was reviewed and categorized based on likelihood of a true Rotterdam PCOS diagnosis, i.e. two out of three of the following: irregular menstrual cycles, hyperandrogenism and/or polycystic ovary morphology. Data from the history, physical exam, laboratory and radiology results were codified and extracted from notes of definite PCOS subjects. Thirty-two terms were used to build an algorithm for identifying definite PCOS cases and applied to the rest of the dataset. The positive predictive value cutoff was set at 76.8 % to maximize the number of subjects available for study. A true positive predictive value for the algorithm was calculated after review of 100 charts from subjects identified as definite PCOS cases with at least two documented Rotterdam criteria. The positive predictive value was compared to that calculated using 200 charts identified using the ICD-9 code for PCOS (256.4; n = 13,670). In addition, a cohort of previously recruited PCOS subjects was submitted for algorithm validation. RESULTS: Chart review demonstrated that 64 % were confirmed as definitely PCOS using the algorithm, with a 9 % false positive rate. 66 % of subjects identified by ICD-9 code for PCOS could be confirmed as definitely PCOS, with an 8.5 % false positive rate. There was no significant difference in the positive predictive values using the two methods (p = 0.2). However, the number of charts that had insufficient confirmatory data was lower using the algorithm (5 % vs 11 %; p < 0.04). Of 477 subjects with PCOS recruited and examined individually and present in the database as patients, 451 were found within the algorithm dataset. CONCLUSIONS: Extraction of text parameters along with codified data improves the confidence in PCOS patient cohorts identified using the electronic medical record. However, the positive predictive value was not significantly different when using ICD-9 codes or the specific algorithm. Further studies are needed to determine the positive predictive value of the two methods in additional electronic medical record datasets.

Mode of childbirth and long-term outcomes in women with inflammatory bowel diseases.

INTRODUCTION: Inflammatory bowel diseases [IBD; Crohn's disease (CD), ulcerative colitis] often affect women in their reproductive years. Few studies have analyzed the impact of mode of childbirth on long-term IBD outcomes. METHODS: We used a multi-institutional IBD cohort to identify all women in the reproductive age-group with a diagnosis of IBD prior to pregnancy. We identified the occurrence of a new diagnosis code for perianal complications, IBD-related hospitalization and surgery, and initiation of medical therapy after either a vaginal delivery or caesarean section (CS). Cox proportional hazards models adjusting for potential confounders were used to estimate independent effect of mode of childbirth on IBD outcomes. RESULTS: Our cohort included 360 women with IBD (161 CS). Women in the CS group were likely to be older and more likely to have complicated disease behavior prior to pregnancy. During follow-up, there was no difference in the likelihood of IBD-related surgery (multivariate hazard ratio 1.75, 95 % confidence interval (CI) 0.40-7.75), IBD-related hospitalization (HR 1.39), initiation of immunomodulator therapy (HR 1.45), or anti-TNF therapy (HR 1.11). Among the 133 CD pregnancies with no prior perianal disease, we found no excess risk of subsequent new diagnosis perianal fistulae with vaginal delivery compared to CS (HR 0.19, 95 % CI 0.04-1.05). CONCLUSIONS: Mode of delivery did not influence natural history of IBD. In our cohort, vaginal delivery was not associated with increased risk of subsequent perianal disease in women with CD.

Improving the power of genetic association tests with imperfect phenotype derived from electronic medical records.

To reduce costs and improve clinical relevance of genetic studies, there has been increasing interest in performing such studies in hospital-based cohorts by linking phenotypes extracted from electronic medical records (EMRs) to genotypes assessed in routinely collected medical samples. A fundamental difficulty in implementing such studies is extracting accurate information about disease outcomes and important clinical covariates from large numbers of EMRs. Recently, numerous algorithms have been developed to infer phenotypes by combining information from multiple structured and unstructured variables extracted from EMRs. Although these algorithms are quite accurate, they typically do not provide perfect classification due to the difficulty in inferring meaning from the text. Some algorithms can produce for each patient a probability that the patient is a disease case. This probability can be thresholded to define case-control status, and this estimated case-control status has been used to replicate known genetic associations in EMR-based studies. However, using the estimated disease status in place of true disease status results in outcome misclassification, which can diminish test power and bias odds ratio estimates. We propose to instead directly model the algorithm-derived probability of being a case. We demonstrate how our approach improves test power and effect estimation in simulation studies, and we describe its performance in a study of rheumatoid arthritis. Our work provides an easily implemented solution to a major practical challenge that arises in the use of EMR data, which can facilitate the use of EMR infrastructure for more powerful, cost-effective, and diverse genetic studies.

Genetic basis of autoantibody positive and negative rheumatoid arthritis risk in a multi-ethnic cohort derived from electronic health records.

Discovering and following up on genetic associations with complex phenotypes require large patient cohorts. This is particularly true for patient cohorts of diverse ancestry and clinically relevant subsets of disease. The ability to mine the electronic health records (EHRs) of patients followed as part of routine clinical care provides a potential opportunity to efficiently identify affected cases and unaffected controls for appropriate-sized genetic studies. Here, we demonstrate proof-of-concept that it is possible to use EHR data linked with biospecimens to establish a multi-ethnic case-control cohort for genetic research of a complex disease, rheumatoid arthritis (RA). In 1,515 EHR-derived RA cases and 1,480 controls matched for both genetic ancestry and disease-specific autoantibodies (anti-citrullinated protein antibodies [ACPA]), we demonstrate that the odds ratios and aggregate genetic risk score (GRS) of known RA risk alleles measured in individuals of European ancestry within our EHR cohort are nearly identical to those derived from a genome-wide association study (GWAS) of 5,539 autoantibody-positive RA cases and 20,169 controls. We extend this approach to other ethnic groups and identify a large overlap in the GRS among individuals of European, African, East Asian, and Hispanic ancestry. We also demonstrate that the distribution of a GRS based on 28 non-HLA risk alleles in ACPA+ cases partially overlaps with ACPA- subgroup of RA cases. Our study demonstrates that the genetic basis of rheumatoid arthritis risk is similar among cases of diverse ancestry divided into subsets based on ACPA status and emphasizes the utility of linking EHR clinical data with biospecimens for genetic studies.

Thromboprophylaxis is associated with reduced post-hospitalization venous thromboembolic events in patients with inflammatory bowel diseases.

BACKGROUND & AIMS: Patients with inflammatory bowel diseases (IBDs) have increased risk for venous thromboembolism (VTE); those who require hospitalization have particularly high risk. Few hospitalized patients with IBD receive thromboprophylaxis. We analyzed the frequency of VTE after IBD-related hospitalization, risk factors for post-hospitalization VTE, and the efficacy of prophylaxis in preventing post-hospitalization VTE. METHODS: In a retrospective study, we analyzed data from a multi-institutional cohort of patients with Crohn's disease or ulcerative colitis and at least 1 IBD-related hospitalization. Our primary outcome was a VTE event. All patients contributed person-time from the date of the index hospitalization to development of VTE, subsequent hospitalization, or end of follow-up. Our main predictor variable was pharmacologic thromboprophylaxis. Cox proportional hazard models adjusting for potential confounders were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs). RESULTS: From a cohort of 2788 patients with at least 1 IBD-related hospitalization, 62 patients developed VTE after discharge (2%). Incidences of VTE at 30, 60, 90, and 180 days after the index hospitalization were 3.7/1000, 4.1/1000, 5.4/1000, and 9.4/1000 person-days, respectively. Pharmacologic thromboprophylaxis during the index hospital stay was associated with a significantly lower risk of post-hospitalization VTE (HR, 0.46; 95% CI, 0.22-0.97). Increased numbers of comorbidities (HR, 1.30; 95% CI, 1.16-1.47) and need for corticosteroids before hospitalization (HR, 1.71; 95% CI, 1.02-2.87) were also independently associated with risk of VTE. Length of hospitalization or surgery during index hospitalization was not associated with post-hospitalization VTE. CONCLUSIONS: Pharmacologic thromboprophylaxis during IBD-related hospitalization is associated with reduced risk of post-hospitalization VTE.

Serum inflammatory markers and risk of colorectal cancer in patients with inflammatory bowel diseases.

BACKGROUND & AIMS: Patients with inflammatory bowel diseases (IBDs) (Crohn's disease, ulcerative colitis) are at increased risk of colorectal cancer (CRC). Persistent inflammation is hypothesized to increase risk of CRC in patients with IBD; however, the few studies in this area have been restricted to cross-sectional assessments of histologic severity. No prior studies have examined association between C-reactive protein (CRP) or erythrocyte sedimentation rate (ESR) elevation and risk of CRC in an IBD cohort. METHODS: From a multi-institutional validated IBD cohort, we identified all patients with at least one measured CRP or ESR value. Patients were stratified into quartiles of severity of inflammation on the basis of their median CRP or ESR value, and subsequent diagnosis of CRC was ascertained. Logistic regression adjusting for potential confounders was used to identify the independent association between CRP or ESR elevation and risk of CRC. RESULTS: Our study included 3145 patients with at least 1 CRP value (CRP cohort) and 4008 with at least 1 ESR value (ESR cohort). Thirty-three patients in the CRP cohort and 102 patients in the ESR cohort developed CRC during a median follow-up of 5 years at a median age of 55 years. On multivariate analysis, there was a significant increase in risk of CRC across quartiles of CRP elevation (P(trend) = .017; odds ratio for quartile 4 vs quartile 1, 2.72; 95% confidence interval, 0.95-7.76). Similarly higher median ESR was also independently associated with risk of CRC across the quartiles (odds ratio, 2.06; 95% confidence interval, 1.14-3.74) (P(trend) = .007). CONCLUSIONS: An elevated CRP or ESR is associated with increased risk of CRC in patients with IBD.

Association between reduced plasma 25-hydroxy vitamin D and increased risk of cancer in patients with inflammatory bowel diseases.

BACKGROUND & AIMS: Vitamin D deficiency is common among patients with inflammatory bowel diseases (IBD) (Crohn's disease or ulcerative colitis). The effects of low plasma 25-hydroxy vitamin D (25[OH]D) on outcomes other than bone health are understudied in patients with IBD. We examined the association between plasma level of 25(OH)D and risk of cancers in patients with IBD. METHODS: From a multi-institutional cohort of patients with IBD, we identified those with at least 1 measurement of plasma 25(OH)D. The primary outcome was development of any cancer. We examined the association between plasma 25(OH)D and risk of specific subtypes of cancer, adjusting for potential confounders in a multivariate regression model. RESULTS: We analyzed data from 2809 patients with IBD and a median plasma level of 25(OH)D of 26 ng/mL. Nearly one-third had deficient levels of vitamin D (<20 ng/mL). During a median follow-up period of 11 years, 196 patients (7%) developed cancer, excluding nonmelanoma skin cancer (41 cases of colorectal cancer). Patients with vitamin D deficiency had an increased risk of cancer (adjusted odds ratio, 1.82; 95% confidence interval, 1.25-2.65) compared with those with sufficient levels. Each 1-ng/mL increase in plasma 25(OH)D was associated with an 8% reduction in risk of colorectal cancer (odds ratio, 0.92; 95% confidence interval, 0.88-0.96). A weaker inverse association was also identified for lung cancer. CONCLUSIONS: In a large multi-institutional IBD cohort, a low plasma level of 25(OH)D was associated with an increased risk of cancer, especially colorectal cancer.

Association between low density lipoprotein and rheumatoid arthritis genetic factors with low density lipoprotein levels in rheumatoid arthritis and non-rheumatoid arthritis controls.

OBJECTIVES: While genetic determinants of low density lipoprotein (LDL) cholesterol levels are well characterised in the general population, they are understudied in rheumatoid arthritis (RA). Our objective was to determine the association of established LDL and RA genetic alleles with LDL levels in RA cases compared with non-RA controls. METHODS: Using data from electronic medical records, we linked validated RA cases and non-RA controls to discarded blood samples. For each individual, we extracted data on: first LDL measurement, age, gender and year of LDL measurement. We genotyped subjects for 11 LDL and 44 non-HLA RA alleles, and calculated RA and LDL genetic risk scores (GRS). We tested the association between each GRS and LDL level using multivariate linear regression models adjusted for age, gender, year of LDL measurement and RA status. RESULTS: Among 567 RA cases and 979 controls, 80% were female and mean age at the first LDL measurement was 55 years. RA cases had significantly lower mean LDL levels than controls (117.2 vs 125.6 mg/dl, respectively, p<0.0001). Each unit increase in LDL GRS was associated with 0.8 mg/dl higher LDL levels in both RA cases and controls (p=3.0×10(-7)). Each unit increase in RA GRS was associated with 4.3 mg/dl lower LDL levels in both groups (p=0.01). CONCLUSIONS: LDL alleles were associated with higher LDL levels in RA. RA alleles were associated with lower LDL levels in both RA cases and controls. As RA cases carry more RA alleles, these findings suggest a genetic basis for epidemiological observations of lower LDL levels in RA.