Community profiles in northeastern and central Pennsylvania characterized by distinct social, natural, food, and physical activity environments and their relation to type 2 diabetes.

Background: Understanding geographic disparities in type 2 diabetes (T2D) requires approaches that account for communities' multidimensional nature. Methods: In an electronic health record nested case-control study, we identified 15,884 cases of new-onset T2D from 2008 to 2016, defined using encounter diagnoses, medication orders, and laboratory test results, and frequency-matched controls without T2D (79,400; 65,069 unique persons). We used finite mixture models to construct community profiles from social, natural, physical activity, and food environment measures. We estimated T2D odds ratios (OR) with 95% confidence intervals (CI) using logistic generalized estimating equation models, adjusted for sociodemographic variables. We examined associations with the profiles alone and combined them with either community type based on administrative boundaries or Census-based urban/rural status. Results: We identified four profiles in 1069 communities in central and northeastern Pennsylvania along a rural-urban gradient: "sparse rural," "developed rural," "inner suburb," and "deprived urban core." Urban areas were densely populated with high physical activity resources and food outlets; however, they also had high socioeconomic deprivation and low greenness. Compared with "developed rural," T2D onset odds were higher in "deprived urban core" (1.24, CI = 1.16-1.33) and "inner suburb" (1.10, CI = 1.04-1.17). These associations with model-based community profiles were weaker than when combined with administrative boundaries or urban/rural status. Conclusions: Our findings suggest that in urban areas, diabetogenic features overwhelm T2D-protective features. The community profiles support the construct validity of administrative-community type and urban/rural status, previously reported, to evaluate geographic disparities in T2D onset in this geography.

Associations between PM2.5 and O3 exposures and new onset type 2 diabetes in regional and national samples in the United States.

BACKGROUND: Exposure to particulate matter ≤2.5 µm in diameter (PM2.5) and ozone (O3) has been linked to numerous harmful health outcomes. While epidemiologic evidence has suggested a positive association with type 2 diabetes (T2D), there is heterogeneity in findings. We evaluated exposures to PM2.5 and O3 across three large samples in the US using a harmonized approach for exposure assignment and covariate adjustment. METHODS: Data were obtained from the Veterans Administration Diabetes Risk (VADR) cohort (electronic health records [EHRs]), the Reasons for Geographic and Racial Disparities in Stroke (REGARDS) cohort (primary data collection), and the Geisinger health system (EHRs), and reflect the years 2003-2016 (REGARDS) and 2008-2016 (VADR and Geisinger). New onset T2D was ascertained using EHR information on medication orders, laboratory results, and T2D diagnoses (VADR and Geisinger) or report of T2D medication or diagnosis and/or elevated blood glucose levels (REGARDS). Exposure was assigned using pollutant annual averages from the Downscaler model. Models stratified by community type (higher density urban, lower density urban, suburban/small town, or rural census tracts) evaluated likelihood of new onset T2D in each study sample in single- and two-pollutant models of PM2.5 and O3. RESULTS: In two pollutant models, associations of PM2.5, and new onset T2D were null in the REGARDS cohort except for in suburban/small town community types in models that also adjusted for NSEE, with an odds ratio (95% CI) of 1.51 (1.01, 2.25) per 5 µg/m3 of PM2.5. Results in the Geisinger sample were null. VADR sample results evidenced nonlinear associations for both pollutants; the shape of the association was dependent on community type. CONCLUSIONS: Associations between PM2.5, O3 and new onset T2D differed across three large study samples in the US. None of the results from any of the three study populations found strong and clear positive associations.

Sinus inflammation and chronic rhinosinusitis are associated with a diagnosis of new onset asthma in the following year.

BACKGROUND: Chronic rhinosinusitis (CRS) and asthma commonly co-occur. No studies have leveraged large samples needed to formally address whether preexisting CRS is associated with new onset asthma over time. METHODS: We evaluated whether prevalent CRS [identified in two ways: validated text algorithm applied to sinus computerized tomography (CT) scan or two diagnoses] was associated with new onset adult asthma in the following year. We used electronic health record data from Geisinger from 2008 to 2019. For each year we removed persons with any evidence of asthma through the end of the year, then identified those with new diagnosis of asthma in the following year. Complementary log-log regression was used to adjust for confounding variables (e.g., sociodemographic, contact with the health system, comorbidities), and hazard ratios (HRs) and 95% confidence intervals (CI) were calculated. RESULTS: A total of 35,441 persons were diagnosed with new onset asthma and were compared to 890,956 persons who did not develop asthma. Persons with new onset asthma tended to be female (69.6%) and younger (mean [SD] age 45.9 [17.0] years). Both CRS definitions were associated (HR, 95% CI) with new onset asthma, with 2.21 (1.93, 2.54) and 1.48 (1.38, 1.59) for CRS based on sinus CT scan and two diagnoses, respectively. New onset asthma was uncommonly observed in persons with a history of sinus surgery. CONCLUSION: Prevalent CRS identified with two complementary approaches was associated with a diagnosis of new onset asthma in the following year. The findings may have clinical implications for the prevention of asthma.

Associations of four indexes of social determinants of health and two community typologies with new onset type 2 diabetes across a diverse geography in Pennsylvania.

Evaluation of geographic disparities in type 2 diabetes (T2D) onset requires multidimensional approaches at a relevant spatial scale to characterize community types and features that could influence this health outcome. Using Geisinger electronic health records (2008-2016), we conducted a nested case-control study of new onset T2D in a 37-county area of Pennsylvania. The study included 15,888 incident T2D cases and 79,435 controls without diabetes, frequency-matched 1:5 on age, sex, and year of diagnosis or encounter. We characterized patients' residential census tracts by four dimensions of social determinants of health (SDOH) and into a 7-category SDOH census tract typology previously generated for the entire United States by dimension reduction techniques. Finally, because the SDOH census tract typology classified 83% of the study region's census tracts into two heterogeneous categories, termed rural affordable-like and suburban affluent-like, to further delineate geographies relevant to T2D, we subdivided these two typology categories by administrative community types (U.S. Census Bureau minor civil divisions of township, borough, city). We used generalized estimating equations to examine associations of 1) four SDOH indexes, 2) SDOH census tract typology, and 3) modified typology, with odds of new onset T2D, controlling for individual-level confounding variables. Two SDOH dimensions, higher socioeconomic advantage and higher mobility (tracts with fewer seniors and disabled adults) were independently associated with lower odds of T2D. Compared to rural affordable-like as the reference group, residence in tracts categorized as extreme poverty (odds ratio [95% confidence interval] = 1.11 [1.02, 1.21]) or multilingual working (1.07 [1.03, 1.23]) were associated with higher odds of new onset T2D. Suburban affluent-like was associated with lower odds of T2D (0.92 [0.87, 0.97]). With the modified typology, the strongest association (1.37 [1.15, 1.63]) was observed in cities in the suburban affluent-like category (vs. rural affordable-like-township), followed by cities in the rural affordable-like category (1.20 [1.05, 1.36]). We conclude that in evaluating geographic disparities in T2D onset, it is beneficial to conduct simultaneous evaluation of SDOH in multiple dimensions. Associations with the modified typology showed the importance of incorporating governmentally, behaviorally, and experientially relevant community definitions when evaluating geographic health disparities.

Urban and rural differences in new onset type 2 diabetes: Comparisons across national and regional samples in the diabetes LEAD network.

Introduction: Geographic disparities in diabetes burden exist throughout the United States (US), with many risk factors for diabetes clustering at a community or neighborhood level. We hypothesized that the likelihood of new onset type 2 diabetes (T2D) would differ by community type in three large study samples covering the US. Research design and methods: We evaluated the likelihood of new onset T2D by a census tract-level measure of community type, a modification of RUCA designations (higher density urban, lower density urban, suburban/small town, and rural) in three longitudinal US study samples (REGARDS [REasons for Geographic and Racial Differences in Stroke] cohort, VADR [Veterans Affairs Diabetes Risk] cohort, Geisinger electronic health records) representing the CDC Diabetes LEAD (Location, Environmental Attributes, and Disparities) Network. Results: In the REGARDS sample, residing in higher density urban community types was associated with the lowest odds of new onset T2D (OR [95% CI]: 0.80 [0.66, 0.97]) compared to rural community types; in the Geisinger sample, residing in higher density urban community types was associated with the highest odds of new onset T2D (OR [95% CI]: 1.20 [1.06, 1.35]) compared to rural community types. In the VADR sample, suburban/small town community types had the lowest hazard ratios of new onset T2D (HR [95% CI]: 0.99 [0.98, 1.00]). However, in a regional stratified analysis of the VADR sample, the likelihood of new onset T2D was consistent with findings in the REGARDS and Geisinger samples, with highest likelihood of T2D in the rural South and in the higher density urban communities of the Northeast and West regions; likelihood of T2D did not differ by community type in the Midwest. Conclusions: The likelihood of new onset T2D by community type varied by region of the US. In the South, the likelihood of new onset T2D was higher among those residing in rural communities.

Urban-Rural Differences in Health Care Utilization and COVID-19 Outcomes in Patients With Type 2 Diabetes.

INTRODUCTION: Two studies in Pennsylvania aimed to determine whether community type and community socioeconomic deprivation (CSD) 1) modified associations between type 2 diabetes (hereinafter, diabetes) and COVID-19 hospitalization outcomes, and 2) influenced health care utilization among individuals with diabetes during the COVID-19 pandemic. METHODS: The hospitalization study evaluated a retrospective cohort of patients hospitalized with COVID-19 through 2020 for COVID-19 outcomes: death, intensive care unit (ICU) admission, mechanical ventilation, elevated D-dimer, and elevated troponin level. We used adjusted logistic regression models, adding interaction terms to evaluate effect modification by community type (township, borough, or city census tract) and CSD. The utilization study included patients with diabetes and a clinical encounter between 2017 and 2020. Autoregressive integrated moving average time-series models evaluated changes in weekly rates of emergency department and outpatient visits, hemoglobin A1c (HbA1c) laboratory tests, and antihyperglycemic medication orders from 2018 to 2020. RESULTS: In the hospitalization study, of 2,751 patients hospitalized for COVID-19, 1,020 had diabetes, which was associated with ICU admission and elevated troponin. Associations did not differ by community type or CSD. In the utilization study, among 93,401 patients with diabetes, utilization measures decreased in March 2020. Utilization increased in July, and then began to stabilize or decline through the end of 2020. Changes in HbA1c tests and medication order trends during the pandemic differed by community type and CSD. CONCLUSION: Diabetes was associated with selected outcomes among individuals hospitalized for COVID-19, but these did not differ by community features. Utilization trajectories among individuals with diabetes during the pandemic were influenced by community type and CSD and could be used to identify individuals at risk of gaps in diabetes care.

Strong and consistent associations of precedent chronic rhinosinusitis with risk of non-cystic fibrosis bronchiectasis.

BACKGROUND: Chronic rhinosinusitis (CRS) and bronchiectasis commonly co-occur, but most prior studies were not designed to evaluate temporality and causality. OBJECTIVES: In a sample representing the general population in 37 counties in Pennsylvania, and thus the full spectrum of sinonasal and relevant lung diseases, we aimed to evaluate the temporality and strength of associations of CRS with non-cystic fibrosis bronchiectasis. METHODS: We completed case-control analyses for each of 3 primary bronchiectasis case finding methods. We used electronic health records to identify CRS and bronchiectasis with diagnoses, procedure orders, and/or specific text in sinus or chest computerized tomography scan radiology reports. The controls never had any indication of bronchiectasis and were frequency-matched to the 3 bronchiectasis groups on the basis of age, sex, and encounter year. There were 5,329 unique persons with bronchiectasis and 33,363 without bronchiectasis in the 3 analyses. Important co-occurring conditions were identified with diagnoses, medication orders, and encounter types. Logistic regression was used to evaluate associations (odds ratios [ORs] and 95% CIs) of CRS with bronchiectasis while adjusting for confounding variables. RESULTS: In adjusted analyses, CRS was consistently and strongly associated with all 3 bronchiectasis definitions. The strongest associations for CRS (ORs and 95% CIs) were those that were based on the text of sinus computerized tomography scan reports; the associations were generally stronger for CRS without nasal polyps (eg, OR = 4.46 [95% CI = 2.09-9.51] for diagnosis-based bronchiectasis). On average, CRS was identified more than 6 years before bronchiectasis. CONCLUSION: Precedent CRS was strongly and consistently associated with increased risk of bronchiectasis.

Risk factors for Lyme disease stage and manifestation using electronic health records.

BACKGROUND: Little is known about risk factors for early (e.g., erythema migrans) and disseminated Lyme disease manifestations, such as arthritis, neurological complications, and carditis. No study has used both diagnoses and free text to classify Lyme disease by disease stage and manifestation. METHODS: We identified Lyme disease cases in 2012-2016 in the electronic health record (EHR) of a large, integrated health system in Pennsylvania. We developed a rule-based text-matching algorithm using regular expressions to extract clinical data from free text. Lyme disease cases were then classified by stage and manifestation using data from both diagnoses and free text. Among cases classified by stage, we evaluated individual, community, and health care variables as predictors of disseminated stage (vs. early) disease using Poisson regression models with robust errors. Final models adjusted for sociodemographic factors, receipt of Medical Assistance (i.e., Medicaid, a proxy for low socioeconomic status), primary care contact, setting of diagnosis, season of diagnosis, and urban/rural status. RESULTS: Among 7310 cases of Lyme disease, we classified 62% by stage. Overall, 23% were classified using both diagnoses and text, 26% were classified using diagnoses only, and 13% were classified using text only. Among the staged diagnoses (n = 4530), 30% were disseminated stage (762 arthritis, 426 neurological manifestations, 76 carditis, 95 secondary erythema migrans, and 76 other manifestations). In adjusted models, we found that persons on Medical Assistance at least 50% of time under observation, compared to never users, had a higher risk (risk ratio [95% confidence interval]) of disseminated Lyme disease (1.20 [1.05, 1.37]). Primary care contact (0.59 [0.54, 0.64]) and diagnosis in the urgent care (0.22 [0.17, 0.29]), compared to the outpatient setting, were associated with lower risk of disseminated Lyme disease. CONCLUSIONS: The associations between insurance payor, primary care status, and diagnostic setting with disseminated Lyme disease suggest that lower socioeconomic status and less health care access could be linked with disseminated stage Lyme disease. Intervening on these factors could reduce the individual and health care burden of disseminated Lyme disease. Our findings demonstrate the value of both diagnostic and narrative text data to identify Lyme disease manifestations in the EHR.

Proximity to freshwater blue space and type 2 diabetes onset: the importance of historical and economic context.

Salutogenic effects of living near aquatic areas (blue space) remain underexplored, particularly in non-coastal and non-urban areas. We evaluated associations of residential proximity to inland freshwater blue space with new onset type 2 diabetes (T2D) in central and northeast Pennsylvania, USA, using medical records to conduct a nested case-control study. T2D cases (n=15,888) were identified from diabetes diagnoses, medication orders, and laboratory test results and frequency-matched on age, sex, and encounter year to diabetes-free controls (n=79,435). We calculated distance from individual residences to the nearest lake, river, tributary, or large stream, and residence within the 100-year floodplain. Logistic regression models adjusted for community socioeconomic deprivation and other confounding variables and stratified by community type (townships [rural/suburban], boroughs [small towns], city census tracts). Compared to individuals living ≥1.25 miles from blue space, those within 0.25 miles had 8% and 17% higher odds of T2D onset in townships and boroughs, respectively. Among city residents, T2D odds were 38-39% higher for those living 0.25 to <0.75 miles from blue space. Residing within the floodplain was associated with 16% and 14% higher T2D odds in townships and boroughs. A post-hoc analysis demonstrated patterns of lower residential property values with nearer distance to the region's predominant waterbody, suggesting unmeasured confounding by socioeconomic disadvantage. This may explain our unexpected findings of higher T2D odds with closer proximity to blue space. Our findings highlight the importance of historic and economic context and interrelated factors such as flood risk and lack of waterfront development in blue space research.

Relations of peri-residential temperature and humidity in tick-life-cycle-relevant time periods with human Lyme disease risk in Pennsylvania, USA.

How weather affects tick development and behavior and human Lyme disease remains poorly understood. We evaluated relations of temperature and humidity during critical periods for the tick lifecycle with human Lyme disease. We used electronic health records from 479,344 primary care patients in 38 Pennsylvania counties in 2006-2014. Lyme disease cases (n = 9657) were frequency-matched (5:1) by year, age, and sex. Using daily weather data at ~4 km2 resolution, we created cumulative metrics hypothesized to promote (warm and humid) or inhibit (hot and dry) tick development or host-seeking during nymph development (March 1-May 31), nymph activity (May 1-July 30), and prior year larva activity (Aug 1-Sept 30). We estimated odds ratios (ORs) of Lyme disease by quartiles of each weather variable, adjusting for demographic, clinical, and other weather variables. Exposure-response patterns were observed for higher cumulative same-year temperature, humidity, and hot and dry days (nymph-relevant), and prior year hot and dry days (larva-relevant), with same-year hot and dry days showing the strongest association (4th vs. 1st quartile OR = 0.40; 95% confidence interval [CI] = 0.36, 0.43). Changing temperature and humidity could increase or decrease human Lyme disease risk.

Association of Greenness with Blood Pressure among Individuals with Type 2 Diabetes across Rural to Urban Community Types in Pennsylvania, USA.

Greenness may impact blood pressure (BP), though evidence is limited among individuals with type 2 diabetes (T2D), for whom BP management is critical. We evaluated associations of residential greenness with BP among individuals with T2D in geographically diverse communities in Pennsylvania. To address variation in greenness type, we evaluated modification of associations by percent forest. We obtained systolic (SBP) and diastolic (DBP) BP measurements from medical records of 9593 individuals following diabetes diagnosis. Proximate greenness was estimated within 1250-m buffers surrounding individuals' residences using the normalized difference vegetation index (NDVI) prior to blood pressure measurement. Percent forest was calculated using the U.S. National Land Cover Database. Linear mixed models with robust standard errors accounted for spatial clustering; models were stratified by community type (townships/boroughs/cities). In townships, the greenest communities, an interquartile range increase in NDVI was associated with reductions in SBP of 0.87 mmHg (95% CI: -1.43, -0.30) and in DBP of 0.41 mmHg (95% CI: -0.78, -0.05). No significant associations were observed in boroughs or cities. Evidence for modification by percent forest was weak. Findings suggest a threshold effect whereby high greenness may be necessary to influence BP in this population and support a slight beneficial impact of greenness on cardiovascular disease risk.

Association of community types and features in a case-control analysis of new onset type 2 diabetes across a diverse geography in Pennsylvania.

OBJECTIVES: To evaluate associations of community types and features with new onset type 2 diabetes in diverse communities. Understanding the location and scale of geographic disparities can lead to community-level interventions. DESIGN: Nested case-control study within the open dynamic cohort of health system patients. SETTING: Large, integrated health system in 37 counties in central and northeastern Pennsylvania, USA. PARTICIPANTS AND ANALYSIS: We used electronic health records to identify persons with new-onset type 2 diabetes from 2008 to 2016 (n=15 888). Persons with diabetes were age, sex and year matched (1:5) to persons without diabetes (n=79 435). We used generalised estimating equations to control for individual-level confounding variables, accounting for clustering of persons within communities. Communities were defined as (1) townships, boroughs and city census tracts; (2) urbanised area (large metro), urban cluster (small cities and towns) and rural; (3) combination of the first two; and (4) county. Community socioeconomic deprivation and greenness were evaluated alone and in models stratified by community types. RESULTS: Borough and city census tract residence (vs townships) were associated (OR (95% CI)) with higher odds of type 2 diabetes (1.10 (1.04 to 1.16) and 1.34 (1.25 to 1.44), respectively). Urbanised areas (vs rural) also had increased odds of type 2 diabetes (1.14 (1.08 to 1.21)). In the combined definition, the strongest associations (vs townships in rural areas) were city census tracts in urban clusters (1.41 (1.22 to 1.62)) and city census tracts in urbanised areas (1.33 (1.22 to 1.45)). Higher community socioeconomic deprivation and lower greenness were each associated with increased odds. CONCLUSIONS: Urban residence was associated with higher odds of type 2 diabetes than for other areas. Higher community socioeconomic deprivation in city census tracts and lower greenness in all community types were also associated with type 2 diabetes.

Unconventional Natural Gas Development and Hospitalization for Heart Failure in Pennsylvania.

BACKGROUND: Growing literature linking unconventional natural gas development (UNGD) to adverse health has implicated air pollution and stress pathways. Persons with heart failure (HF) are susceptible to these stressors. OBJECTIVES: This study sought to evaluate associations between UNGD activity and hospitalization among HF patients, stratified by both ejection fraction (EF) status (reduced [HFrEF], preserved [HFpEF], not classifiable) and HF severity. METHODS: We evaluated the odds of hospitalization among patients with HF seen at Geisinger from 2008 to 2015 using electronic health records. We assigned metrics of UNGD activity by phase (pad preparation, drilling, stimulation, and production) 30 days before hospitalization or a frequency-matched control selection date. We assigned phenotype status using a validated algorithm. RESULTS: We identified 9,054 patients with HF with 5,839 hospitalizations (mean age 71.1 ± 12.7 years; 47.7% female). Comparing 4th to 1st quartiles, adjusted odds ratios (95% confidence interval) for hospitalization were 1.70 (1.35 to 2.13), 0.97 (0.75 to 1.27), 1.80 (1.35 to 2.40), and 1.62 (1.07 to 2.45) for pad preparation, drilling, stimulation, and production metrics, respectively. We did not find effect modification by HFrEF or HFpEF status. Associations of most UNGD metrics with hospitalization were stronger among those with more severe HF at baseline. CONCLUSIONS: Three of 4 phases of UNGD activity were associated with hospitalization for HF in a large sample of patients with HF in an area of active UNGD, with similar findings by HFrEF versus HFpEF status. Older patients with HF seem particularly vulnerable to adverse health impacts from UNGD activity.

Peridomestic and community-wide landscape risk factors for Lyme disease across a range of community contexts in Pennsylvania.

Land use and forest fragmentation are thought to be major drivers of Lyme disease incidence and its geographic distribution. We examined the association between landscape composition and configuration and Lyme disease in a population-based case control study in the Geisinger health system in Pennsylvania. Lyme disease cases (n = 9657) were identified using a combination of diagnosis codes, laboratory codes, and antibiotic orders from electronic health records (EHRs). Controls (5:1) were randomly selected and frequency matched on year, age, and sex. We measured six landscape variables based on prior literature, derived from the National Land Cover Database and MODIS satellite imagery: greenness (normalized difference vegetation index), percent forest, percent herbaceous, forest edge density, percent forest-herbaceous edge, and mean forest patch size. We assigned landscape variables within two spatial contexts (community and ½-mile [805 m] Euclidian residential buffer). In models stratified by community type, landscape variables were modeled as tertiles and flexible splines and associations were adjusted for demographic and clinical covariates. In general, we observed positive associations between landscape metrics and Lyme disease, except for percent herbaceous, where associations differed by community type. For example, compared to the lowest tertile, individuals with highest tertile of greenness in residential buffers had higher odds of Lyme disease (odds ratio: 95% confidence interval [CI]) in townships (1.73: 1.55, 1.93), boroughs (1.70: 1.40, 2.07), and cities (3.71: 1.74, 7.92). Similarly, corresponding odds ratios (95% CI) for forest edge density were 1.34 (1.22, 1.47), 1.56 (1.33, 1.82), and 1.90 (1.13, 3.18). Associations were generally higher in residential buffers, compared to community, and in cities, compared to boroughs or townships. Our results reinforce the importance of peridomestic landscape in Lyme disease risk, particularly measures that reflect human interaction with tick habitat. Linkage of EHR data to public data on residential and community context may lead to new health system-based approaches for improving Lyme disease diagnosis, treatment, and prevention.

Associations of multidimensional socioeconomic and built environment factors with body mass index trajectories among youth in geographically heterogeneous communities.

Understanding contextual influences on obesity requires comparison of heterogeneous communities and concurrent assessment of multiple contextual domains. We used a theoretically-based measurement model to assess multidimensional socioeconomic and built environment factors theorized to influence childhood obesity across a diverse geography ranging from rural to urban. Confirmatory factor analysis specified four factors-community socioeconomic deprivation (CSED), food outlet abundance (FOOD), fitness and recreational assets (FIT), and utilitarian physical activity favorability (UTIL)-which were assigned to communities (townships, boroughs, city census tracts) in 37 Pennsylvania counties. Using electronic health records from 2001 to 2012 from 163,820 youth aged 3-18 years from 1288 communities, we conducted multilevel linear regression analyses with factor quartiles and their cross products with age, age2, and age3 to test whether community factors impacted body mass index (BMI) growth trajectories. Models controlled for sex, age, race/ethnicity, and Medical Assistance. Factor scores were lowest in townships, indicating less deprivation, fewer food and physical activity outlets, and lower utilitarian physical activity favorability. BMI at average age was lower in townships versus boroughs (beta [SE]) (0.217 [0.027], P < 0.001) and cities (0.378 [0.036], P < 0.001), as was BMI growth over time. Factor distributions across community types lacked overlap, requiring stratified analyses to avoid extrapolation. In townships, FOOD, UTIL, and FIT were inversely associated with BMI trajectories. Across community types, youth in the lowest (versus higher) CSED quartiles had lower BMI at average age and slower BMI growth, signifying the importance of community deprivation to the obesogenicity of environments.

Household Food Insecurity and Home Food Availability in Relation to Youth Diet, Body Mass Index, and Adiposity.

BACKGROUND: Food security status is related to food types available in the home, which may shape youth dietary patterns, with implications for obesity. OBJECTIVE: Investigate whether household food insecurity and home food availability (HFA) are associated with youth fruit and vegetable (F/V) consumption and anthropometric outcomes. DESIGN: Cross-sectional study. Youth and parents completed questionnaires during in-home visits (2013-2014). Research staff obtained anthropometric measures. PARTICIPANTS/SETTING: Medical record data for 10- to 15-year-old Pennsylvania youths were used to identify 434 parent-youth dyads, with 408 evaluated after excluding missing data. MAIN OUTCOME MEASURES: Parent-reported household food security was assessed with the six-item US Department of Agriculture Food Security Scale (dichotomized as high vs low). Healthy and obesogenic HFA scales assessed parent report of how frequently particular foods were present in the home. Youth self-reported daily average F/V consumption. Anthropometric outcomes included age- and sex-standardized z scores for body mass index (BMIz), waist circumference (WCz), and percent body fat (PBFz). STATISTICAL ANALYSES: Associations were evaluated with multivariable linear regression adjusted for youth age, sex, and race or ethnicity, and parent age and income. RESULTS: Compared with food secure counterparts, youth from food insecure households had higher mean (beta [standard error]) BMIz (.30 [.15]), WCz (.27 [.12]), and PBFz (.43 [.16]). Food insecure households had lower mean healthy HFA scores (-1.23 [.54]); there was no evidence obesogenic HFA differed between food secure and insecure households. Youth from lower healthy HFA or higher obesogenic HFA households reported fewer mean daily F/V servings (healthy HFA: .08 [.02]; obesogenic HFA: -.06 [.02]). Food security status was not associated with F/V consumption, nor was there evidence HFA modified associations between food insecurity and anthropometric outcomes. CONCLUSIONS: Despite an observed association between healthy HFA and youth F/V consumption, this study did not provide evidence that HFA explained associations between food insecurity and youth anthropometric outcomes.

Epidemiology of Lyme disease in Pennsylvania 2006-2014 using electronic health records.

Lyme disease is the most common vector-borne disease in the United States. Electronic health record (EHR)-based research on Lyme disease is limited. We used Geisinger EHR data from 479,344 primary care patients in 38 Pennsylvania counties in 2006-2014 to compare EHR-based Lyme disease incidence rates to surveillance incidence rates, evaluate individual and community risk factors for incident Lyme disease, and to characterize the proportion of cases with diagnoses consistent with post-treatment Lyme disease syndrome in the EHR (PTLDSEHR). We primarily identified Lyme disease cases using diagnosis codes, serologic testing order codes, and medication orders but also completed subgroup analyses among those with positive serology and those with both diagnosis code and antibiotic treatment. We compared annual incidence rates from the EHR to surveillance by age, sex, and county. In case-control analyses, we compared cases to randomly selected controls (5:1) frequency-matched on year, age, and sex. We identified 9657 cases of Lyme disease, including 1791 cases with positive serology and 4992 cases with both diagnosis code and antibiotic treatment. Annual incidence rates in the EHR were 4.25-7.43 times higher than surveillance. In adjusted analyses, white non-Hispanic race/ethnicity (vs. black, Hispanic, or other) was associated with higher odds of Lyme disease (odds ratio [OR]: 2.06, 95% confidence interval [CI]: 1.73-2.44). Medical Assistance insurance use (always vs. never; OR: 0.77, 95% CI: 0.68-0.88), and higher community-level socioeconomic deprivation (quartile 4 vs. 1 OR: 0.50 (95% CI: 0.42-0.59) were associated with lower odds of Lyme disease. Within 4-52 weeks after Lyme disease diagnosis, 20.8% (n = 735) of cases with a diagnosis code and treatment had a diagnosis of malaise or fatigue, pain, or cognitive difficulties not present in the past 26 weeks. These results highlight the utility of EHR data for epidemiologic research on Lyme disease for case-finding, surveillance, risk factor evaluation, and characterization of PTLDS using EHR data.

Associations of unconventional natural gas development with depression symptoms and disordered sleep in Pennsylvania.

Environmental and community factors may influence the development or course of depression and sleep problems. We evaluated the association of unconventional natural gas development (UNGD) with depression symptoms and disordered sleep diagnoses using the Patient Health Questionnaire-8 and electronic health record data among Geisinger adult primary care patients in Pennsylvania. Participants received a retrospective metric for UNGD at their residence (very low, low, medium, and high) that incorporated dates and durations of well development, distance from patient homes to wells, and well characteristics. Analyses included 4,762 participants with no (62%), mild (23%), moderate (10%), and moderately severe or severe (5%) depression symptoms in 2014-2015 and 3,868 disordered sleep diagnoses between 2009-2015. We observed associations between living closer to more and bigger wells and depression symptoms, but not disordered sleep diagnoses in models weighted to account for sampling design and participation. High UNGD (vs. very low) was associated with depression symptoms in an adjusted negative binomial model (exponentiated coefficient = 1.18, 95% confidence interval [CI]: 1.04-1.34). High and low UNGD (vs. very low) were associated with depression symptoms (vs. none) in an adjusted multinomial logistic model. Our findings suggest that UNGD may be associated with adverse mental health in Pennsylvania.

Measures of SES for Electronic Health Record-based Research.

INTRODUCTION: Although infrequently recorded in electronic health records (EHRs), measures of SES are essential to describe health inequalities and account for confounding in epidemiologic research. Medical Assistance (i.e., Medicaid) is often used as a surrogate for SES, but correspondence between conventional SES and Medical Assistance has been insufficiently studied. METHODS: Geisinger Clinic EHR data from 2001 to 2014 and a 2014 questionnaire were used to create six SES measures: EHR-derived Medical Assistance and proportion of time under observation on Medical Assistance; educational attainment, income, and marital status; and area-level poverty. Analyzed in 2016-2017, associations of SES measures with obesity, hypertension, type 2 diabetes, chronic rhinosinusitis, fatigue, and migraine headache were assessed using weighted age- and sex-adjusted logistic regression. RESULTS: Among 5,550 participants (interquartile range, 39.6-57.5 years, 65.9% female), 83% never used Medical Assistance. All SES measures were correlated (Spearman's p≤0.4). Medical Assistance was significantly associated with all six health outcomes in adjusted models. For example, the OR for prevalent type 2 diabetes associated with Medical Assistance was 1.7 (95% CI=1.3, 2.2); the OR for high school versus college graduates was 1.7 (95% CI=1.2, 2.5). Medical Assistance was an imperfect proxy for SES: associations between conventional SES measures and health were attenuated <20% after adjustment for Medical Assistance. CONCLUSIONS: Because systematically collected SES measures are rarely available in EHRs and are unlikely to appear soon, researchers can use EHR-based Medical Assistance to describe inequalities. As SES has many domains, researchers who use Medical Assistance to evaluate the association of SES with health should expect substantial unmeasured confounding.

Comparing objective measures of the built environment in their associations with youth physical activity and sedentary behavior across heterogeneous geographies.

We compared two strategies for measuring built environment features in their associations with youth physical activity and sedentary behavior across heterogeneous geographies of Pennsylvania. Physical activity environments of communities representing a rural-to-urban gradient were characterized through direct observation and spatially referenced archival data subjected to confirmatory factor analysis. Stratified regression analyses assessed associations between environmental measures and behavioral outcomes by community type. Neither strategy was consistently associated with behavior across communities. Findings highlight the importance of differentiating community type in evaluating associations of the built environment, and the challenge of measuring meaningful differences that influence youth behavior across heterogeneous geographies.