IncidencePrevalence: An R package to calculate population-level incidence rates and prevalence using the OMOP common data model.

PURPOSE: Real-world data (RWD) offers a valuable resource for generating population-level disease epidemiology metrics. We aimed to develop a well-tested and user-friendly R package to compute incidence rates and prevalence in data mapped to the observational medical outcomes partnership (OMOP) common data model (CDM). MATERIALS AND METHODS: We created IncidencePrevalence, an R package to support the analysis of population-level incidence rates and point- and period-prevalence in OMOP-formatted data. On top of unit testing, we assessed the face validity of the package. To do so, we calculated incidence rates of COVID-19 using RWD from Spain (SIDIAP) and the United Kingdom (CPRD Aurum), and replicated two previously published studies using data from the Netherlands (IPCI) and the United Kingdom (CPRD Gold). We compared the obtained results to those previously published, and measured execution times by running a benchmark analysis across databases. RESULTS: IncidencePrevalence achieved high agreement to previously published data in CPRD Gold and IPCI, and showed good performance across databases. For COVID-19, incidence calculated by the package was similar to public data after the first-wave of the pandemic. CONCLUSION: For data mapped to the OMOP CDM, the IncidencePrevalence R package can support descriptive epidemiological research. It enables reliable estimation of incidence and prevalence from large real-world data sets. It represents a simple, but extendable, analytical framework to generate estimates in a reproducible and timely manner.

Trends of use and characterisation of anti-dementia drugs users: a large multinational-network population-based study.

BACKGROUND: An updated time-trend analysis of anti-dementia drugs (ADDs) is lacking. The aim of this study is to assess the incident rate (IR) of ADD in individuals with dementia using real-world data. SETTING: Primary care data (country/database) from the UK/CPRD-GOLD (2007-20), Spain/SIDIAP (2010-20) and the Netherlands/IPCI (2008-20), standardised to a common data model. METHODS: Cohort study. Participants: dementia patients ≥40 years old with ≥1 year of previous data. Follow-up: until the end of the study period, transfer out of the catchment area, death or incident prescription of rivastigmine, galantamine, donepezil or memantine. Other variables: age/sex, type of dementia, comorbidities. Statistics: overall and yearly age/sex IR, with 95% confidence interval, per 100,000 person-years (IR per 105 PY (95%CI)). RESULTS: We identified a total of (incident anti-dementia users/dementia patients) 41,024/110,642 in UK/CPRD-GOLD, 51,667/134,927 in Spain/SIDIAP and 2,088/17,559 in the Netherlands/IPCI.In the UK, IR (per 105 PY (95%CI)) of ADD decreased from 2007 (30,829 (28,891-32,862)) to 2010 (17,793 (17,083-18,524)), then increased up to 2019 (31,601 (30,483 to 32,749)) and decrease in 2020 (24,067 (23,021-25,148)). In Spain, IR (per 105 PY (95%CI)) of ADD decreased by 72% from 2010 (51,003 (49,199-52,855)) to 2020 (14,571 (14,109-15,043)). In the Netherlands, IR (per 105 PY (95%CI)) of ADD decreased by 77% from 2009 (21,151 (14,967-29,031)) to 2020 (4763 (4176-5409)). Subjects aged ≥65-79 years and men (in the UK and the Netherlands) initiated more frequently an ADD. CONCLUSIONS: Treatment of dementia remains highly heterogeneous. Further consensus in the pharmacological management of patients living with dementia is urgently needed.

Impact of the COVID-19 pandemic on eating disorders diagnoses among adolescents and young adults in Catalonia: A population-based cohort study.

OBJECTIVE: To describe the impact of the COVID-19 pandemic on trends in incidence rates (IR) of diagnoses of eating disorders (ED) among adolescents and young adults. METHODS: Population-based cohort study using primary care records of people aged 10-24 years between January, 2016 and December, 2021 in Catalonia, Spain. IRs were calculated monthly and grouped by the different stages of the COVID-19 pandemic in Catalonia: (1) the pre-lockdown (January, 2016-February, 2020), (2) lockdown (March-June, 2020) and, (3) post-lockdown (July, 2020-December, 2021) periods. Incidence rate ratios (IRR) relative to the corresponding periods in 2018-2019 were calculated. RESULTS: A total of 1,179,009 individuals were included. The IR was 9.2 per 100,000 person-months (95% confidence intervals [CI]: 8.9-9.5) during the pre-lockdown period. It decreased during the lockdown period (6.3 per 100,000 person-months [5.5-7.3]), but substantially increased during the following period (19.4. per 100,000 person-months [18.7-20.1]). While large reductions in IRs were observed for both sexes during the lockdown period (IRR 95% CI: 0.65 [0.54-0.78] in females and 0.46 [0.29-0.71] in males), substantial increases during the post-lockdown period were limited to females, and were particularly pronounced among those aged 10-14 and 15-19 years (2.50 [2.23-2.80] and 2.29 [2.07-2.54], respectively). DISCUSSION: The COVID-19 pandemic has resulted in a substantial increase in ED diagnoses, primarily driven by higher rates among adolescent females. PUBLIC SIGNIFICANCE: This population-based cohort study demonstrated a substantial increase in incidence rates of eating disorders in primary care following the end of lockdown in Catalonia, Spain, with adolescent girls seen to be most affected.

Socioeconomic Inequalities in COVID-19 Vaccination and Infection in Adults, Catalonia, Spain.

Evidence on the impact of the COVID-19 vaccine rollout on socioeconomic COVID-19-related inequalities is scarce. We analyzed associations between socioeconomic deprivation index (SDI) and COVID-19 vaccination, infection, and hospitalization before and after vaccine rollout in Catalonia, Spain. We conducted a population-based cohort study during September 2020-June 2021 that comprised 2,297,146 adults >40 years of age. We estimated odds ratio of nonvaccination and hazard ratios (HRs) of infection and hospitalization by SDI quintile relative to the least deprived quintile, Q1. Six months after rollout, vaccination coverage differed by SDI quintile in working-age (40-64 years) persons: 81% for Q1, 71% for Q5. Before rollout, we found a pattern of increased HR of infection and hospitalization with deprivation among working-age and retirement-age (>65 years) persons. After rollout, infection inequalities decreased in both age groups, whereas hospitalization inequalities decreased among retirement-age persons. Our findings suggest that mass vaccination reduced socioeconomic COVID-19-related inequalities.

Cancer and the risk of coronavirus disease 2019 diagnosis, hospitalisation and death: A population-based multistate cohort study including 4 618 377 adults in Catalonia, Spain.

The relationship between cancer and coronavirus disease 2019 (COVID-19) infection and severity remains poorly understood. We conducted a population-based cohort study between 1 March and 6 May 2020 describing the associations between cancer and risk of COVID-19 diagnosis, hospitalisation and COVID-19-related death. Data were obtained from the Information System for Research in Primary Care (SIDIAP) database, including primary care electronic health records from ~80% of the population in Catalonia, Spain. Cancer was defined as any primary invasive malignancy excluding non-melanoma skin cancer. We estimated adjusted hazard ratios (aHRs) for the risk of COVID-19 (outpatient) clinical diagnosis, hospitalisation (with or without a prior COVID-19 diagnosis) and COVID-19-related death using Cox proportional hazard regressions. Models were estimated for the overall cancer population and by years since cancer diagnosis (<1 year, 1-5 years and ≥5 years), sex, age and cancer type; and adjusted for age, sex, smoking status, deprivation and comorbidities. We included 4 618 377 adults, of which 260 667 (5.6%) had a history of cancer. A total of 98 951 individuals (5.5% with cancer) were diagnosed, and 6355 (16.4% with cancer) were directly hospitalised with COVID-19. Of those diagnosed, 6851 were subsequently hospitalised (10.7% with cancer), and 3227 died without being hospitalised (18.5% with cancer). Among those hospitalised, 1963 (22.5% with cancer) died. Cancer was associated with an increased risk of COVID-19 diagnosis (aHR: 1.08; 95% confidence interval [1.05-1.11]), direct COVID-19 hospitalisation (1.33 [1.24-1.43]) and death following hospitalisation (1.12 [1.01-1.25]). These associations were stronger for patients recently diagnosed with cancer, aged <70 years, and with haematological cancers. These patients should be prioritised in COVID-19 vaccination campaigns and continued non-pharmaceutical interventions.

Adherence of long-term breast cancer survivors to follow-up care guidelines: a study based on real-world data from the SURBCAN cohort.

PURPOSE: To identify adherence to follow-up recommendations in long-term breast cancer survivors (LTBCS) of the SURBCAN cohort and to identify its determinants, using real-world data. METHODS: We conducted a retrospective study using electronic health records from 2012 to 2016 of women diagnosed with incident breast cancer in Spain between 2000 and 2006 and surviving at least 5 years. Adherence to basic follow-up recommendations, adherence according to risk of recurrence, and overall adherence were calculated based on attendance at medical appointments and imaging surveillance, by year of survivorship. Logistic regression models were fitted to depict the association between adherence and its determinants. RESULTS: A total of 2079 LTBCS were followed up for a median of 4.97 years. Of them, 23.6% had survived ≥ 10 years at baseline. We estimated that 79.5% of LTBCS were overall adherent to at least one visit and one imaging test. Adherence to recommendations decreased over time and no differences were found according to recurrence risk. Determinants of better overall adherence were diagnosis in middle age (50-69 years old), living in a more-deprived area, having fewer years of survival, receiving primary treatment, and being alive at the end of follow-up. CONCLUSION: We identified women apparently not complying with surveillance visits and tests. Special attention should be paid to the youngest and eldest women at diagnosis and to those with longer survival.

Knee and hip osteoarthritis and risk of nine cancers in a large real-world matched cohort study.

OBJECTIVES: Joint replacement due to end-stage OA has been linked to incidence of several cancers. We aimed to estimate the association between newly diagnosed knee and hip OA and incidence of nine common cancer types. METHODS: We identified persons with incident knee or hip OA, aged ≥40 years, between 2009 and 2015 in the SIDIAP database in Catalonia, Spain. We matched up to three OA-free controls on age, sex and general practitioner. We followed participants from 1 year after OA diagnosis until migration, death, end of study at 31 December 2017 or incident cancer of: stomach, colorectal, liver, pancreas, lung, skin, breast, prostate and bladder. We used flexible parametric survival models, adjusted for confounders. Estimates were corrected for misclassification using probabilistic bias analysis. RESULTS: We included 117 750 persons with knee OA and matched 309 913 persons without, with mean (s.d.) age of 67.5 (11.1) years and 63% women. The hip cohort consisted of 39 133 persons with hip OA and 116 713 controls. For most of the included cancers, the hazard ratios (HRs) were close to 1. The HR of lung cancer for knee OA exposure was 0.80 (95% CI: 0.71, 0.89) and attenuated to 0.98 (0.76, 1.27) in non-smokers. The hazard of colorectal cancer was lower in persons with both knee and hip OA by 10-20%. CONCLUSIONS: Knee and hip OA are not associated with studied incident cancers, apart from lower risk of colorectal cancer. The often-reported protective association of knee OA with lung cancer is explained by residual confounding.

Seek COVER: using a disease proxy to rapidly develop and validate a personalized risk calculator for COVID-19 outcomes in an international network.

BACKGROUND: We investigated whether we could use influenza data to develop prediction models for COVID-19 to increase the speed at which prediction models can reliably be developed and validated early in a pandemic. We developed COVID-19 Estimated Risk (COVER) scores that quantify a patient's risk of hospital admission with pneumonia (COVER-H), hospitalization with pneumonia requiring intensive services or death (COVER-I), or fatality (COVER-F) in the 30-days following COVID-19 diagnosis using historical data from patients with influenza or flu-like symptoms and tested this in COVID-19 patients. METHODS: We analyzed a federated network of electronic medical records and administrative claims data from 14 data sources and 6 countries containing data collected on or before 4/27/2020. We used a 2-step process to develop 3 scores using historical data from patients with influenza or flu-like symptoms any time prior to 2020. The first step was to create a data-driven model using LASSO regularized logistic regression, the covariates of which were used to develop aggregate covariates for the second step where the COVER scores were developed using a smaller set of features. These 3 COVER scores were then externally validated on patients with 1) influenza or flu-like symptoms and 2) confirmed or suspected COVID-19 diagnosis across 5 databases from South Korea, Spain, and the United States. Outcomes included i) hospitalization with pneumonia, ii) hospitalization with pneumonia requiring intensive services or death, and iii) death in the 30 days after index date. RESULTS: Overall, 44,507 COVID-19 patients were included for model validation. We identified 7 predictors (history of cancer, chronic obstructive pulmonary disease, diabetes, heart disease, hypertension, hyperlipidemia, kidney disease) which combined with age and sex discriminated which patients would experience any of our three outcomes. The models achieved good performance in influenza and COVID-19 cohorts. For COVID-19 the AUC ranges were, COVER-H: 0.69-0.81, COVER-I: 0.73-0.91, and COVER-F: 0.72-0.90. Calibration varied across the validations with some of the COVID-19 validations being less well calibrated than the influenza validations. CONCLUSIONS: This research demonstrated the utility of using a proxy disease to develop a prediction model. The 3 COVER models with 9-predictors that were developed using influenza data perform well for COVID-19 patients for predicting hospitalization, intensive services, and fatality. The scores showed good discriminatory performance which transferred well to the COVID-19 population. There was some miscalibration in the COVID-19 validations, which is potentially due to the difference in symptom severity between the two diseases. A possible solution for this is to recalibrate the models in each location before use.

Air pollution and green spaces in relation to breast cancer risk among pre and postmenopausal women: A mega cohort from Catalonia.

BACKGROUND: The association between air pollution and green spaces with breast cancer risk stratified by menopausal status has not been frequently investigated despite its importance given the different impact of risk factors on breast cancer risk depending on menopausal status. OBJECTIVES: To study the association between air pollution, green spaces and pre and postmenopausal breast cancer risk. METHODS: We conducted a population-based cohort study using electronic primary care records in Catalonia. We included women aged 17-85 years free of cancer at study entry between 2009 and 2017. Our exposures were particulate matter <2.5 µm (PM2.5) & <10 µm (PM10), nitrogen dioxide (NO2), normalized difference vegetation index (NDVI), and percentage of green spaces estimated at the census tract level. Breast cancer was identified with ICD-10 code C50. We estimated cause-specific hazard ratios (HR) for the relationship between each individual exposure and pre and postmenopausal breast cancer risk, using linear and non-linear models. RESULTS: Of the 1,054,180 pre and 744,658 postmenopausal women followed for a median of 10 years, 6,126 and 17,858 developed breast cancer, respectively. Among premenopausal women, only very high levels of PM10 (≥46 µg/m3) were associated with increased cancer risk (compared to lower levels) in non-linear models. Among postmenopausal women, an interquartile range increase in PM2.5 (HR:1.03; 95%CI:1.01-1.04), PM10 (1.03; 1.01-1.05), and NO2 (1.05; 1.02-1.08) were associated with higher cancer risk. NDVI was negatively associated with decreased cancer risk only among postmenopausal women who did not change residence during follow-up (0.84; 0.71-0.99) or who were followed for at least three years (0.82; 0.69-0.98). DISCUSSION: Living in areas with high concentrations of PM2.5, PM10, and NO2 increases breast cancer risk in postmenopausal women while long-term exposure to green spaces may decrease this risk. Only very high concentrations of PM10 increase breast cancer risk in premenopausal women.

Multimorbidity clusters among long-term breast cancer survivors in Spain: Results of the SURBCAN study.

The disease management of long-term breast cancer survivors (BCS) is hampered by the scarce knowledge of multimorbidity patterns. The aim of our study was to identify multimorbidity clusters among long-term BCS and assess their impact on mortality and health services use. We conducted a retrospective study using electronic health records of 6512 BCS from Spain surviving at least 5 years. Hierarchical cluster analysis was used to identify groups of similar patients based on their chronic diagnoses, which were assessed using the Clinical Classifications Software. As a result, multimorbidity clusters were obtained, clinically defined and named according to the comorbidities with higher observed/expected prevalence ratios. Multivariable Cox and negative binomial regression models were fitted to estimate overall mortality risk and probability of contacting health services according to the clusters identified. 83.7% of BCS presented multimorbidity, essential hypertension (34.5%) and obesity and other metabolic disorders (27.4%) being the most prevalent chronic diseases at the beginning of follow-up. Five multimorbidity clusters were identified: C1-unspecific (29.9%), C2-metabolic and neurodegenerative (28.3%), C3-anxiety and fractures (9.7%), C4-musculoskeletal and cardiovascular (9.6%) and C5-thyroid disorders (5.3%). All clusters except C5-thyroid disorders were associated with higher mortality compared to BCS without comorbidities. The risk of mortality in C4 was increased by 64% (adjusted hazard ratio 1.64, 95% confidence interval 1.52-2.07). Stratified analysis showed an increased risk of death among BCS with 5 to 10 years of survival in all clusters. These results help to identify subgroups of long-term BCS with specific needs and mortality risks and to guide BCS clinical practice regarding multimorbidity.

Metabolic perturbations prior to hepatocellular carcinoma diagnosis: Findings from a prospective observational cohort study.

Hepatocellular carcinoma (HCC) development entails changes in liver metabolism. Current knowledge on metabolic perturbations in HCC is derived mostly from case-control designs, with sparse information from prospective cohorts. Our objective was to apply comprehensive metabolite profiling to detect metabolites whose serum concentrations are associated with HCC development, using biological samples from within the prospective European Prospective Investigation into Cancer and Nutrition (EPIC) cohort (>520 000 participants), where we identified 129 HCC cases matched 1:1 to controls. We conducted high-resolution untargeted liquid chromatography-mass spectrometry-based metabolomics on serum samples collected at recruitment prior to cancer diagnosis. Multivariable conditional logistic regression was applied controlling for dietary habits, alcohol consumption, smoking, body size, hepatitis infection and liver dysfunction. Corrections for multiple comparisons were applied. Of 9206 molecular features detected, 220 discriminated HCC cases from controls. Detailed feature annotation revealed 92 metabolites associated with HCC risk, of which 14 were unambiguously identified using pure reference standards. Positive HCC-risk associations were observed for N1-acetylspermidine, isatin, p-hydroxyphenyllactic acid, tyrosine, sphingosine, l,l-cyclo(leucylprolyl), glycochenodeoxycholic acid, glycocholic acid and 7-methylguanine. Inverse risk associations were observed for retinol, dehydroepiandrosterone sulfate, glycerophosphocholine, γ-carboxyethyl hydroxychroman and creatine. Discernible differences for these metabolites were observed between cases and controls up to 10 years prior to diagnosis. Our observations highlight the diversity of metabolic perturbations involved in HCC development and replicate previous observations (metabolism of bile acids, amino acids and phospholipids) made in Asian and Scandinavian populations. These findings emphasize the role of metabolic pathways associated with steroid metabolism and immunity and specific dietary and environmental exposures in HCC development.

Body mass index and waist circumference in relation to the risk of 26 types of cancer: a prospective cohort study of 3.5 million adults in Spain.

BACKGROUND: A high body mass index (BMI) has been associated with increased risk of several cancers; however, whether BMI is related to a larger number of cancers than currently recognized is unclear. Moreover, whether waist circumference (WC) is more strongly associated with specific cancers than BMI is not well established. We aimed to investigate the associations between BMI and 26 cancers accounting for non-linearity and residual confounding by smoking status as well as to compare cancer risk estimates between BMI and WC. METHODS: Prospective cohort study with population-based electronic health records from Catalonia, Spain. We included 3,658,417 adults aged ≥ 18 years and free of cancer at baseline between 2006 and 2017. Our main outcome measures were cause-specific hazard ratios (HRs) with 99% confidence intervals (CIs) for incident cancer at 26 anatomical sites. RESULTS: After a median follow-up time of 8.3 years, 202,837 participants were diagnosed with cancer. A higher BMI was positively associated with risk of nine cancers (corpus uteri, kidney, gallbladder, thyroid, colorectal, breast post-menopausal, multiple myeloma, leukemia, non-Hodgkin lymphoma) and was positively associated with three additional cancers among never smokers (head and neck, brain and central nervous system, Hodgkin lymphoma). The respective HRs (per 5 kg/m2 increment) ranged from 1.04 (99%CI 1.01 to 1.08) for non-Hodgkin lymphoma to 1.49 (1.45 to 1.53) for corpus uteri cancer. While BMI was negatively associated to five cancer types in the linear analyses of the overall population, accounting for non-linearity revealed that BMI was associated to prostate cancer in a U-shaped manner and to head and neck, esophagus, larynx, and trachea, bronchus and lung cancers in an L-shaped fashion, suggesting that low BMIs are an approximation of heavy smoking. Of the 291,305 participants with a WC measurement, 27,837 were diagnosed with cancer. The 99%CIs of the BMI and WC point estimates (per 1 standard deviation increment) overlapped for all cancers. CONCLUSIONS: In this large Southern European study, a higher BMI was associated with increased risk of twelve cancers, including four hematological and head and neck (only among never smokers) cancers. Furthermore, BMI and WC showed comparable estimates of cancer risk associated with adiposity.

Ambient air pollution and the development of overweight and obesity in children: a large longitudinal study.

BACKGROUND: Ambient air pollution may play a role in childhood obesity development, but evidence is scarce, and the modifying role of socioeconomic status (SES) is unclear. We aimed to examine the association between exposure to air pollution during early childhood and subsequent risk of developing overweight and obesity, and to evaluate whether SES is a modifier of this association. METHODS: This longitudinal study included 416,955 children identified as normal weight between 2-5 years old and registered in an electronic primary healthcare record between 2006 and 2016 in Catalonia (Spain). Children were followed-up until they developed overweight or obesity, reached 15 years of age, died, transferred out, or end of study period (31/12/2018). Overweight and obesity were defined following the WHO reference obtained from height and weight measures. We estimated annual residential census levels of nitrogen dioxide (NO2) and particulate matter <10 µm (PM10), <2.5 µm (PM2.5), and 2.5-10 µm (PMcoarse) at study entry. We estimated the risk of developing overweight and obesity per interquartile range increase in air pollution exposure with Cox proportional hazard models. RESULTS: A total of 142,590 (34.2%) children developed overweight or obesity. Increased exposure to NO2, PM10, and PMcoarse was associated with a 2-3% increased risk of developing overweight and obesity (hazard ratio [HR] per 21.8 µg/m3 NO2 = 1.03 [95% CI: 1.02-1.04]; HR per 6.4 µg/m3 PM10 = 1.02 [95% CI: 1.02-1.03]; HR per 4.6 µg/m3 PMcoarse = 1.02, [95% CI: 1.01-1.02]). For all air pollutants, associations were stronger among children living in most compared to least deprived areas. CONCLUSIONS: This study suggests that early life exposure to air pollution may be associated with a small increase in the risk of developing overweight and obesity in childhood, and that this association may be exacerbated in the most deprived areas. Even these small associations are of potential global health importance because air pollution exposure is widespread and the long-term health consequences of childhood obesity are clear.

Characteristics and outcomes of 627 044 COVID-19 patients living with and without obesity in the United States, Spain, and the United Kingdom.

BACKGROUND: A detailed characterization of patients with COVID-19 living with obesity has not yet been undertaken. We aimed to describe and compare the demographics, medical conditions, and outcomes of COVID-19 patients living with obesity (PLWO) to those of patients living without obesity. METHODS: We conducted a cohort study based on outpatient/inpatient care and claims data from January to June 2020 from Spain, the UK, and the US. We used six databases standardized to the OMOP common data model. We defined two non-mutually exclusive cohorts of patients diagnosed and/or hospitalized with COVID-19; patients were followed from index date to 30 days or death. We report the frequency of demographics, prior medical conditions, and 30-days outcomes (hospitalization, events, and death) by obesity status. RESULTS: We included 627 044 (Spain: 122 058, UK: 2336, and US: 502 650) diagnosed and 160 013 (Spain: 18 197, US: 141 816) hospitalized patients with COVID-19. The prevalence of obesity was higher among patients hospitalized (39.9%, 95%CI: 39.8-40.0) than among those diagnosed with COVID-19 (33.1%; 95%CI: 33.0-33.2). In both cohorts, PLWO were more often female. Hospitalized PLWO were younger than patients without obesity. Overall, COVID-19 PLWO were more likely to have prior medical conditions, present with cardiovascular and respiratory events during hospitalization, or require intensive services compared to COVID-19 patients without obesity. CONCLUSION: We show that PLWO differ from patients without obesity in a wide range of medical conditions and present with more severe forms of COVID-19, with higher hospitalization rates and intensive services requirements. These findings can help guiding preventive strategies of COVID-19 infection and complications and generating hypotheses for causal inference studies.

Risk of depression, suicide and psychosis with hydroxychloroquine treatment for rheumatoid arthritis: a multinational network cohort study.

OBJECTIVES: Concern has been raised in the rheumatology community regarding recent regulatory warnings that HCQ used in the coronavirus disease 2019 pandemic could cause acute psychiatric events. We aimed to study whether there is risk of incident depression, suicidal ideation or psychosis associated with HCQ as used for RA. METHODS: We performed a new-user cohort study using claims and electronic medical records from 10 sources and 3 countries (Germany, UK and USA). RA patients ≥18 years of age and initiating HCQ were compared with those initiating SSZ (active comparator) and followed up in the short (30 days) and long term (on treatment). Study outcomes included depression, suicide/suicidal ideation and hospitalization for psychosis. Propensity score stratification and calibration using negative control outcomes were used to address confounding. Cox models were fitted to estimate database-specific calibrated hazard ratios (HRs), with estimates pooled where I2 <40%. RESULTS: A total of 918 144 and 290 383 users of HCQ and SSZ, respectively, were included. No consistent risk of psychiatric events was observed with short-term HCQ (compared with SSZ) use, with meta-analytic HRs of 0.96 (95% CI 0.79, 1.16) for depression, 0.94 (95% CI 0.49, 1.77) for suicide/suicidal ideation and 1.03 (95% CI 0.66, 1.60) for psychosis. No consistent long-term risk was seen, with meta-analytic HRs of 0.94 (95% CI 0.71, 1.26) for depression, 0.77 (95% CI 0.56, 1.07) for suicide/suicidal ideation and 0.99 (95% CI 0.72, 1.35) for psychosis. CONCLUSION: HCQ as used to treat RA does not appear to increase the risk of depression, suicide/suicidal ideation or psychosis compared with SSZ. No effects were seen in the short or long term. Use at a higher dose or for different indications needs further investigation. TRIAL REGISTRATION: Registered with EU PAS (reference no. EUPAS34497; http://www.encepp.eu/encepp/viewResource.htm? id=34498). The full study protocol and analysis source code can be found at https://github.com/ohdsi-studies/Covid19EstimationHydroxychloroquine2.

COVID-19 in patients with autoimmune diseases: characteristics and outcomes in a multinational network of cohorts across three countries.

OBJECTIVE: Patients with autoimmune diseases were advised to shield to avoid coronavirus disease 2019 (COVID-19), but information on their prognosis is lacking. We characterized 30-day outcomes and mortality after hospitalization with COVID-19 among patients with prevalent autoimmune diseases, and compared outcomes after hospital admissions among similar patients with seasonal influenza. METHODS: A multinational network cohort study was conducted using electronic health records data from Columbia University Irving Medical Center [USA, Optum (USA), Department of Veterans Affairs (USA), Information System for Research in Primary Care-Hospitalization Linked Data (Spain) and claims data from IQVIA Open Claims (USA) and Health Insurance and Review Assessment (South Korea). All patients with prevalent autoimmune diseases, diagnosed and/or hospitalized between January and June 2020 with COVID-19, and similar patients hospitalized with influenza in 2017-18 were included. Outcomes were death and complications within 30 days of hospitalization. RESULTS: We studied 133 589 patients diagnosed and 48 418 hospitalized with COVID-19 with prevalent autoimmune diseases. Most patients were female, aged ≥50 years with previous comorbidities. The prevalence of hypertension (45.5-93.2%), chronic kidney disease (14.0-52.7%) and heart disease (29.0-83.8%) was higher in hospitalized vs diagnosed patients with COVID-19. Compared with 70 660 hospitalized with influenza, those admitted with COVID-19 had more respiratory complications including pneumonia and acute respiratory distress syndrome, and higher 30-day mortality (2.2-4.3% vs 6.32-24.6%). CONCLUSION: Compared with influenza, COVID-19 is a more severe disease, leading to more complications and higher mortality.

Use of hydrochlorothiazide and risk of skin cancer in a large nested case-control study in Spain.

PURPOSE: Hydrochlorothiazide (HCTZ) use has been linked to skin cancer in northern European countries. We assessed the association between HCTZ exposure and risk of malignant melanoma (MM) and keratinocyte carcinoma (KC) in a European Mediterranean population. METHODS: Two parallel nested case-control studies were conducted in Spain using two electronic primary healthcare databases, each one providing data on both exposure and outcomes: SIDIAP and BIFAP. Cancer cases were matched to 10 controls by age and gender through risk-set sampling. The ORs and 95% CI for MM and KC associated with previous HCTZ use were estimated using conditional logistic regression. In BIFAP, KC cases were further identified as basal cell carcinoma (BCC) or squamous cell carcinoma (SCC). RESULTS: In adjusted analyses, both ever and cumulative high (≥50,000 mg) use of HCTZ were associated with an increased risk of KC. The risk estimates for high use were 1.30 (1.26-1.34) in SIDIAP and 1.20 (1.12-1.30) in BIFAP, with a lower risk for BCC (1.11 [1.02-1.21]) than for SCC (1.71 [1.45-2.02]). A dose-response relationship was observed between cumulative doses of HCTZ and KC risk. Inconsistent results were found for high use of HCTZ and risk of MM: 1.25 (1.09-1.43) in SIDIAP and 0.85 (0.64-1.13) in BIFAP. CONCLUSIONS: In this European Mediterranean population, a high cumulative use of HCTZ was related to an increased risk of KC with a clear dose-response pattern.

The relationship between residential proximity to outdoor play spaces and children's mental and behavioral health: The importance of neighborhood socio-economic characteristics.

Urban outdoor play spaces are reported to improve children's health. However, there is little empirical evidence on the impact of outdoor play spaces on childhood mental and behavioral health. To fill this gap, we investigated the associations between residential proximity to outdoor play spaces and the prevalence of diagnosed mental and behavioral disorders. We explored whether these associations differ by individual and area-level socio-economic status (SES). This cross-sectional study included 151 110 children who were 0-12 years old in 2014 and were visited in public primary health care centers in Barcelona (Spain). Each child's demographic and mental and behavioral disorders information was extracted for 2005-2014, including diagnoses on disorders of psychological development together with other four types of mental and behavioral disorders. The pediatrician diagnosed mental and behavioral disorders we explored in this study were: mood/affective; neurotic, stress-related and somatoform; psychological development; behavioral and emotional; and overall mental and behavioral disorders. We assessed 300 m network buffer residential proximity to overall outdoor play spaces (i.e., the overall sum of play spaces of any type), outdoor green play spaces, and to a diversity of outdoor play spaces. We used robust Poisson regression models to investigate the association between proximity to outdoor play spaces indicators and each health outcome. We tested interaction terms for indicators of proximity to outdoor play spaces and individual and area SES. For measures with significant interaction terms, we conducted stratified models. We found residential proximity to outdoor play spaces to be protective of disorders of psychological development. Proximity to overall outdoor play spaces, proximity to outdoor green play spaces and proximity to a greater diversity of outdoor play spaces were associated with a 4% (95% CI: 1,7), 4% (95% CI: 1,7) and 5% (95% CI: 2,9) lower prevalence rates of disorders of psychological development respectively. Most of the associations were found to be in the same direction-although more pronounced-in low SES areas, but in the opposite direction for children living in high SES areas. No differences in these associations were found by individual SES. Residential proximity to outdoor play spaces is protective of children's mental and behavioral health living in low SES areas.

Impact of a single-age cohort human papillomavirus vaccination strategy in Catalonia, Spain: Population-based analysis of anogenital warts in men and women.

Extensive multiple-age cohort human papillomavirus (HPV) vaccination has proved to be highly effective. We aimed to determine the 8-year population impact of a female single-age cohort HPV vaccination programme on the incidence of anogenital warts (AGW). In 2008, Catalonia initiated a school-based quadrivalent HPV vaccination programme targeting 11-year-old girls, achieving coverage over 80%. Data on diagnoses of AGW and genital herpes were obtained from a population-based database of electronic health records covering 74% of the population. The annual incidence rates from 2009 to 2016 were calculated, stratified by age and sex using Joinpoint regression to estimate trends and annual percentage changes (APC). Among women aged 16-19 years, the AGW incidence decreased by 61% from 2012 to 2016 (APC -19.4%; 95% CI: -30.0 to -7.3). In contrast, the incidence of genital herpes in same-aged women increased throughout the study period (APC 11.1%; 95% CI: 7.2-15.2). Among men aged 20-22 years, the increasing incidence of AGW shifted to a downward trend in 2013 (APC 2009-2013: 17.0%; 95% CI: 8.2-26.5; and APC 2013-2016: -4.5%; 95% CI: -14.6 to 6.9). A similar pattern was observed among men aged 23-25 years (APC 2009-2014: 16.0%; 95% CI: 12.0-20.2; and APC 2014-2016: -6.0%; 95% CI: -18.4 to 8.3). In contrast to AGW, among men aged 20-25 years, the incidence of genital herpes increased over this period. Our study strongly suggests that a single-cohort HPV vaccination strategy with high vaccine uptake not only provides direct benefit in the vaccinated cohorts but also extends protection through a herd effect to unvaccinated men.

Risks and clinical predictors of cirrhosis and hepatocellular carcinoma diagnoses in adults with diagnosed NAFLD: real-world study of 18 million patients in four European cohorts.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is a common condition that progresses in some patients to steatohepatitis (NASH), cirrhosis and hepatocellular carcinoma (HCC). Here we used healthcare records of 18 million adults to estimate risk of acquiring advanced liver disease diagnoses in patients with NAFLD or NASH compared to individually matched controls. METHODS: Data were extracted from four European primary care databases representing the UK, Netherlands, Italy and Spain. Patients with a recorded diagnosis of NAFLD or NASH (NAFLD/NASH) were followed up for incident cirrhosis and HCC diagnoses. Each coded NAFLD/NASH patient was matched to up to 100 "non-NAFLD" patients by practice site, gender, age ± 5 years and visit recorded within ± 6 months. Hazard ratios (HR) were estimated using Cox models adjusted for age and smoking status and pooled across databases by random effects meta-analyses. RESULTS: Out of 18,782,281 adults, we identified 136,703 patients with coded NAFLD/NASH. Coded NAFLD/NASH patients were more likely to have diabetes, hypertension and obesity than matched controls. HR for cirrhosis in patients compared to controls was 4.73 (95% CI 2.43-9.19) and for HCC, 3.51 (95% CI 1.72-7.16). HR for either outcome was higher in patients with NASH and those with high-risk Fib-4 scores. The strongest independent predictor of a diagnosis of HCC or cirrhosis was baseline diagnosis of diabetes. CONCLUSIONS: Real-world population data show that recorded diagnosis of NAFLD/NASH increases risk of life-threatening liver outcomes. Diabetes is an independent predictor of advanced liver disease diagnosis, emphasising the need to identify specific groups of patients at highest risk.