Estimating the effect of realistic improvements of metformin adherence on COVID-19 mortality using targeted machine learning.

Background: Type 2 diabetes elevates the risk of severe outcomes in COVID-19 patients, with multiple studies reporting higher case fatality rates. Metformin is a widely used medication for glycemic management. We hypothesize that improved adherence to metformin may lower COVID-19 post-infection mortality risk in this group. Utilizing data from the Mexican Social Security Institute (IMSS), we investigate the relationship between metformin adherence and mortality following COVID-19 infection in patients with chronic metformin prescriptions. Methods: This is a retrospective cohort study consisting of 61,180 IMSS beneficiaries who received a positive polymerase chain reaction (PCR) or rapid test for SARS-CoV-2 and had at least two consecutive months of metformin prescriptions prior to the positive test. The hypothetical intervention is improved adherence to metformin, measured by proportion of days covered (PDC), with the comparison being the observed metformin adherence values. The primary outcome is all-cause mortality following COVID-19 infection. We defined the causal parameter using shift intervention, an example of modified treatment policies. We used the targeted learning framework for estimation of the target estimand. Findings: Among COVID-19 positive patients with chronic metformin prescriptions, we found that a 5% and 10% absolute increase in metformin adherence is associated with a respective 0.26% (95% CI: -0.28%, 0.79%) and 1.26% (95% CI: 0.72%, 1.80%) absolute decrease in mortality risk. Interpretation: Subject to the limitations of a real-world data study, our results indicate a causal association between improved metformin adherence and reduced COVID-19 post-infection mortality risk.

Treatment Heterogeneity of Water, Sanitation, Hygiene, and Nutrition Interventions on Child Growth by Environmental Enteric Dysfunction and Pathogen Status for Young Children in Bangladesh.

Background: Water, sanitation, hygiene (WSH), nutrition (N), and combined (N+WSH) interventions are often implemented by global health organizations, but WSH interventions may insufficiently reduce pathogen exposure, and nutrition interventions may be modified by environmental enteric dysfunction (EED), a condition of increased intestinal permeability and inflammation. This study investigated the heterogeneity of these treatments' effects based on individual pathogen and EED biomarker status with respect to child linear growth. Methods: We applied cross-validated targeted maximum likelihood estimation and super learner ensemble machine learning to assess the conditional treatment effects in subgroups defined by biomarker and pathogen status. We analyzed treatment (N+WSH, WSH, N, or control) randomly assigned in-utero, child pathogen and EED data at 14 months of age, and child LAZ at 28 months of age. We estimated the difference in mean child length for age Z-score (LAZ) under the treatment rule and the difference in stratified treatment effect (treatment effect difference) comparing children with high versus low pathogen/biomarker status while controlling for baseline covariates. Results: We analyzed data from 1,522 children, who had median LAZ of -1.56. We found that myeloperoxidase (N+WSH treatment effect difference 0.0007 LAZ, WSH treatment effect difference 0.1032 LAZ, N treatment effect difference 0.0037 LAZ) and Campylobacter infection (N+WSH treatment effect difference 0.0011 LAZ, WSH difference 0.0119 LAZ, N difference 0.0255 LAZ) were associated with greater effect of all interventions on growth. In other words, children with high myeloperoxidase or Campylobacter infection experienced a greater impact of the interventions on growth. We found that a treatment rule that assigned the N+WSH (LAZ difference 0.23, 95% CI (0.05, 0.41)) and WSH (LAZ difference 0.17, 95% CI (0.04, 0.30)) interventions based on EED biomarkers and pathogens increased predicted child growth compared to the randomly allocated intervention. Conclusions: These findings indicate that EED biomarker and pathogen status, particularly Campylobacter and myeloperoxidase (a measure of gut inflammation), may be related to impact of N+WSH, WSH, and N interventions on child linear growth.

A cluster-randomized trial of water, sanitation, handwashing and nutritional interventions on stress and epigenetic programming.

A regulated stress response is essential for healthy child growth and development trajectories. We conducted a cluster-randomized trial in rural Bangladesh (funded by the Bill & Melinda Gates Foundation, ClinicalTrials.gov NCT01590095) to assess the effects of an integrated nutritional, water, sanitation, and handwashing intervention on child health. We previously reported on the primary outcomes of the trial, linear growth and caregiver-reported diarrhea. Here, we assessed additional prespecified outcomes: physiological stress response, oxidative stress, and DNA methylation (N = 759, ages 1-2 years). Eight neighboring pregnant women were grouped into a study cluster. Eight geographically adjacent clusters were block-randomized into the control or the combined nutrition, water, sanitation, and handwashing (N + WSH) intervention group (receiving nutritional counseling and lipid-based nutrient supplements, chlorinated drinking water, upgraded sanitation, and handwashing with soap). Participants and data collectors were not masked, but analyses were masked. There were 358 children (68 clusters) in the control group and 401 children (63 clusters) in the intervention group. We measured four F2-isoprostanes isomers (iPF(2α)-III; 2,3-dinor-iPF(2α)-III; iPF(2α)-VI; 8,12-iso-iPF(2α)-VI), salivary alpha-amylase and cortisol, and methylation of the glucocorticoid receptor (NR3C1) exon 1F promoter including the NGFI-A binding site. Compared with control, the N + WSH group had lower concentrations of F2-isoprostanes isomers (differences ranging from -0.16 to -0.19 log ng/mg of creatinine, P < 0.01), elevated post-stressor cortisol (0.24 log µg/dl; P < 0.01), higher cortisol residualized gain scores (0.06 µg/dl; P = 0.023), and decreased methylation of the NGFI-A binding site (-0.04; P = 0.037). The N + WSH intervention enhanced adaptive responses of the physiological stress system in early childhood.

Stress biomarkers and child development in young children in Bangladesh.

BACKGROUND: Hundreds of millions of children in low- and middle-income countries are exposed to chronic stressors, such as poverty, poor sanitation and hygiene, and sub-optimal nutrition. These stressors can have physiological consequences for children and may ultimately have detrimental effects on child development. This study explores associations between biological measures of chronic stress in early life and developmental outcomes in a large cohort of young children living in rural Bangladesh. METHODS: We assessed physiologic measures of stress in the first two years of life using measures of the hypothalamic-pituitary-adrenal (HPA) axis (salivary cortisol and glucocorticoid receptor gene methylation), the sympathetic-adrenal-medullary (SAM) system (salivary alpha-amylase, heart rate, and blood pressure), and oxidative status (F2-isoprostanes). We assessed child development in the first two years of life with the MacArthur-Bates Communicative Development Inventories (CDI), the WHO gross motor milestones, and the Extended Ages and Stages Questionnaire (EASQ). We compared development outcomes of children at the 75th and 25th percentiles of stress biomarker distributions while adjusting for potential confounders using generalized additive models, which are statistical models where the outcome is predicted by a potentially non-linear function of predictor variables. RESULTS: We analyzed data from 684 children (49% female) at both 14 and 28 months of age; we included an additional 765 children at 28 months of age. We detected a significant relationship between HPA axis activity and child development, where increased HPA axis activity was associated with poor development outcomes. Specifically, we found that cortisol reactivity (coefficient -0.15, 95% CI (-0.29, -0.01)) and post-stressor levels (coefficient -0.12, 95% CI (-0.24, -0.01)) were associated with CDI comprehension score, post-stressor cortisol was associated with combined EASQ score (coefficient -0.22, 95% CI (-0.41, -0.04), and overall glucocorticoid receptor methylation was associated with CDI expression score (coefficient -0.09, 95% CI (-0.17, -0.01)). We did not detect a significant relationship between SAM activity or oxidative status and child development. CONCLUSIONS: Our observations reveal associations between the physiological evidence of stress in the HPA axis with developmental status in early childhood. These findings add to the existing evidence exploring the developmental consequences of early life stress.

Improved Child Feces Management Mediates Reductions in Childhood Diarrhea from an On-Site Sanitation Intervention: Causal Mediation Analysis of a Cluster-Randomized Trial in Rural Bangladesh.

BACKGROUND: The WASH benefits Bangladesh trial multi-component sanitation intervention reduced diarrheal disease among children < 5 years. Intervention components included latrine upgrades, child feces management tools, and behavioral promotion. It remains unclear which components most impacted diarrhea. METHODS: We conducted mediation analysis within a subset of households (n = 720) from the sanitation and control arms. Potential mediators were categorized into indicators of latrine quality, latrine use practices, and feces management practices. We estimated average causal mediation effects (ACME) as prevalence differences (PD), defined as the intervention's effect on diarrhea through its effect on the mediator. RESULTS: The intervention improved all indicators compared to controls. We found significant mediation through multiple latrine use and feces management practice indicators. The strongest mediators during monsoon seasons were reduced open defecation among children aged < 3 and 3-8 years, and increased disposal of child feces into latrines. The strongest mediators during dry seasons were access to a flush/pour-flush latrine, reduced open defecation among children aged 3-8 years, and increased disposal of child feces into latrines. Individual mediation effects were small (PD = 0.5-2 percentage points) compared to the overall intervention effect but collectively describe significant mediation pathways. DISCUSSION: The effect of the WASH Benefits Bangladesh sanitation intervention on diarrheal disease was mediated through improved child feces management and reduced child open defecation. Although the intervention significantly improved latrine quality, relatively high latrine quality at baseline may have limited benefits from additional improvements. Targeting safe child feces management may increase the health benefits of rural sanitation interventions.

Geographic pair matching in large-scale cluster randomized trials.

Cluster randomized trials are often used to study large-scale public health interventions. In large trials, even small improvements in statistical efficiency can have profound impacts on the required sample size and cost. Location integrates many socio-demographic and environmental characteristics into a single, readily available feature. Here we show that pair matching by geographic location leads to substantial gains in statistical efficiency for 14 child health outcomes that span growth, development, and infectious disease through a re-analysis of two large-scale trials of nutritional and environmental interventions in Bangladesh and Kenya. Relative efficiencies from pair matching are ≥1.1 for all outcomes and regularly exceed 2.0, meaning an unmatched trial would need to enroll at least twice as many clusters to achieve the same level of precision as the geographically pair matched design. We also show that geographically pair matched designs enable estimation of fine-scale, spatially varying effect heterogeneity under minimal assumptions. Our results demonstrate broad, substantial benefits of geographic pair matching in large-scale, cluster randomized trials.

Sex differences in ascending aortic size reporting and growth on chest computed tomography and magnetic resonance imaging.

PURPOSE: Diameter-based guidelines for prophylactic repair of ascending aortic aneurysms have led to routine aortic evaluation in chest imaging. Despite sex differences in aneurysm outcomes, there is little understanding of sex-specific aortic growth rates. Our objective was to evaluate sex-specific temporal changes in radiologist-reported aortic size as well as sex differences in aortic reporting. METHOD: In this cohort study, we queried radiology reports of chest computed tomography or magnetic resonance imaging at an academic medical center from 1994 to 2022, excluding type A dissection. Aortic diameter was extracted using a custom text-processing algorithm. Growth rates were estimated using mixed-effects modeling with fixed terms for sex, age, and imaging modality, and patient-level random intercepts. Sex, age, and modality were evaluated as predictors of aortic reporting by logistic regression. RESULTS: This study included 89,863 scans among 46,622 patients (median [interquartile range] age, 64 [52-73]; 22,437 women [48%]). Aortic diameter was recorded in 14% (12,722/89,863 reports). Temporal trends were analyzed in 7194 scans among 1998 patients (age, 68 [60-75]; 677 women [34%]) with ≥2 scans. Aortic growth rate was significantly higher in women (0.22 mm/year [95% confidence interval 0.17-0.28] vs. 0.09 mm/year [0.06-0.13], respectively). Aortic reporting was significantly less common in women (odds ratio, 0.54; 95% CI, 0.52-0.56; p < 0.001). CONCLUSIONS: While aortic growth rates were small overall, women had over twice the growth rate of men. Aortic dimensions were much less frequently reported in women than men. Sex-specific standardized assessment of aortic measurements may be needed to address sex differences in aneurysm outcomes.

Who is most at risk of dying if infected with SARS-CoV-2? A mortality risk factor analysis using machine learning of patients with COVID-19 over time: a large population-based cohort study in Mexico.

OBJECTIVE: COVID-19 would kill fewer people if health programmes can predict who is at higher risk of mortality because resources can be targeted to protect those people from infection. We predict mortality in a very large population in Mexico with machine learning using demographic variables and pre-existing conditions. DESIGN: Cohort study. SETTING: March 2020 to November 2021 in Mexico, nationally represented. PARTICIPANTS: 1.4 million laboratory-confirmed patients with COVID-19 in Mexico at or over 20 years of age. PRIMARY AND SECONDARY OUTCOME MEASURES: Analysis is performed on data from March 2020 to November 2021 and over three phases: (1) from March to October in 2020, (2) from November 2020 to March 2021 and (3) from April to November 2021. We predict mortality using an ensemble machine learning method, super learner, and independently estimate the adjusted mortality relative risk of each pre-existing condition using targeted maximum likelihood estimation. RESULTS: Super learner fit has a high predictive performance (C-statistic: 0.907), where age is the most predictive factor for mortality. After adjusting for demographic factors, renal disease, hypertension, diabetes and obesity are the most impactful pre-existing conditions. Phase analysis shows that the adjusted mortality risk decreased over time while relative risk increased for each pre-existing condition. CONCLUSIONS: While age is the most important predictor of mortality, younger individuals with hypertension, diabetes and obesity are at comparable mortality risk as individuals who are 20 years older without any of the three conditions. Our model can be continuously updated to identify individuals who should most be protected against infection as the pandemic evolves.

Child wasting and concurrent stunting in low- and middle-income countries.

Sustainable Development Goal 2.2-to end malnutrition by 2030-includes the elimination of child wasting, defined as a weight-for-length z-score that is more than two standard deviations below the median of the World Health Organization standards for child growth1. Prevailing methods to measure wasting rely on cross-sectional surveys that cannot measure onset, recovery and persistence-key features that inform preventive interventions and estimates of disease burden. Here we analyse 21 longitudinal cohorts and show that wasting is a highly dynamic process of onset and recovery, with incidence peaking between birth and 3 months. Many more children experience an episode of wasting at some point during their first 24 months than prevalent cases at a single point in time suggest. For example, at the age of 24 months, 5.6% of children were wasted, but by the same age (24 months), 29.2% of children had experienced at least one wasting episode and 10.0% had experienced two or more episodes. Children who were wasted before the age of 6 months had a faster recovery and shorter episodes than did children who were wasted at older ages; however, early wasting increased the risk of later growth faltering, including concurrent wasting and stunting (low length-for-age z-score), and thus increased the risk of mortality. In diverse populations with high seasonal rainfall, the population average weight-for-length z-score varied substantially (more than 0.5 z in some cohorts), with the lowest mean z-scores occurring during the rainiest months; this indicates that seasonally targeted interventions could be considered. Our results show the importance of establishing interventions to prevent wasting from birth to the age of 6 months, probably through improved maternal nutrition, to complement current programmes that focus on children aged 6-59 months.

Early-childhood linear growth faltering in low- and middle-income countries.

Globally, 149 million children under 5 years of age are estimated to be stunted (length more than 2 standard deviations below international growth standards)1,2. Stunting, a form of linear growth faltering, increases the risk of illness, impaired cognitive development and mortality. Global stunting estimates rely on cross-sectional surveys, which cannot provide direct information about the timing of onset or persistence of growth faltering-a key consideration for defining critical windows to deliver preventive interventions. Here we completed a pooled analysis of longitudinal studies in low- and middle-income countries (n = 32 cohorts, 52,640 children, ages 0-24 months), allowing us to identify the typical age of onset of linear growth faltering and to investigate recurrent faltering in early life. The highest incidence of stunting onset occurred from birth to the age of 3 months, with substantially higher stunting at birth in South Asia. From 0 to 15 months, stunting reversal was rare; children who reversed their stunting status frequently relapsed, and relapse rates were substantially higher among children born stunted. Early onset and low reversal rates suggest that improving children's linear growth will require life course interventions for women of childbearing age and a greater emphasis on interventions for children under 6 months of age.

Causes and consequences of child growth faltering in low-resource settings.

Growth faltering in children (low length for age or low weight for length) during the first 1,000 days of life (from conception to 2 years of age) influences short-term and long-term health and survival1,2. Interventions such as nutritional supplementation during pregnancy and the postnatal period could help prevent growth faltering, but programmatic action has been insufficient to eliminate the high burden of stunting and wasting in low- and middle-income countries. Identification of age windows and population subgroups on which to focus will benefit future preventive efforts. Here we use a population intervention effects analysis of 33 longitudinal cohorts (83,671 children, 662,763 measurements) and 30 separate exposures to show that improving maternal anthropometry and child condition at birth accounted for population increases in length-for-age z-scores of up to 0.40 and weight-for-length z-scores of up to 0.15 by 24 months of age. Boys had consistently higher risk of all forms of growth faltering than girls. Early postnatal growth faltering predisposed children to subsequent and persistent growth faltering. Children with multiple growth deficits exhibited higher mortality rates from birth to 2 years of age than children without growth deficits (hazard ratios 1.9 to 8.7). The importance of prenatal causes and severe consequences for children who experienced early growth faltering support a focus on pre-conception and pregnancy as a key opportunity for new preventive interventions.

Development of a Machine Learning Model of Postoperative Acute Kidney Injury Using Non-Invasive Time-Sensitive Intraoperative Predictors.

Acute kidney injury (AKI) is a major postoperative complication that lacks established intraoperative predictors. Our objective was to develop a prediction model using preoperative and high-frequency intraoperative data for postoperative AKI. In this retrospective cohort study, we evaluated 77,428 operative cases at a single academic center between 2016 and 2022. A total of 11,212 cases with serum creatinine (sCr) data were included in the analysis. Then, 8519 cases were randomly assigned to the training set and the remainder to the validation set. Fourteen preoperative and twenty intraoperative variables were evaluated using elastic net followed by hierarchical group least absolute shrinkage and selection operator (LASSO) regression. The training set was 56% male and had a median [IQR] age of 62 (51-72) and a 6% AKI rate. Retained model variables were preoperative sCr values, the number of minutes meeting cutoffs for urine output, heart rate, perfusion index intraoperatively, and the total estimated blood loss. The area under the receiver operator characteristic curve was 0.81 (95% CI, 0.77-0.85). At a score threshold of 0.767, specificity was 77% and sensitivity was 74%. A web application that calculates the model score is available online. Our findings demonstrate the utility of intraoperative time series data for prediction problems, including a new potential use of the perfusion index. Further research is needed to evaluate the model in clinical settings.

What Happens on Call Doesn't Stay on Call. The Effects of In-house Call on Acute Care Surgeons' Sleep and Burnout: Results of the Surgeon Performance (SuPer) Trial.

OBJECTIVE: We sought to quantify the effects of in-house call(IHC) on sleep patterns and burnout among acute care surgeons (ACS). BACKGROUND: Many ACS take INC, which leads to disrupted sleep and high levels of stress and burnout. METHODS: Physiological and survey data of 224 ACS with IHC were collected over 6 months. Participants continuously wore a physiological tracking device and responded to daily electronic surveys. Daily surveys captured work and life events as well as feelings of restfulness and burnout. The Maslach Burnout Inventory (MBI) was administered at the beginning and end of the study period. RESULTS: Physiological data were recorded for 34,135 days, which includes 4389 nights of IHC. Feelings of moderate, very, or extreme burnout occurred 25.7% of days and feelings of being moderately, slightly, or not at all rested occurred 75.91% of days. Decreased amount of time since the last IHC, reduced sleep duration, being on call, and having a bad outcome all contribute to greater feelings of daily burnout ( P <0.001). Decreased time since last call also exacerbates the negative effect of IHC on burnout ( P <0.01). CONCLUSIONS: ACS exhibit lower quality and reduced amount of sleep compared with an age-matched population. Furthermore, reduced sleep and decreased time since the last call led to increased feelings of daily burnout, accumulating in emotional exhaustion as measured on the MBI. A reevaluation of IHC requirements and patterns as well as identification of countermeasures to restore homeostatic wellness in ACS is essential to protect and optimize our workforce.

Geographic pair-matching in large-scale cluster randomized trials.

Custer randomized trials are often used to study large-scale public health interventions. In large trials, even small improvements in statistical efficiency can have profound impacts on the required sample size and cost. Pair matched randomization is one strategy with potential to increase trial efficiency, but to our knowledge there have been no empirical evaluations of pair-matching in large-scale, epidemiologic field trials. Location integrates many socio-demographic and environmental characteristics into a single feature. Here, we show that geographic pair-matching leads to substantial gains in statistical efficiency for 14 child health outcomes that span growth, development, and infectious disease through a re-analysis of two large-scale trials of nutritional and environmental interventions in Bangladesh and Kenya. We estimate relative efficiencies ≥1.1 for all outcomes assessed and relative efficiencies regularly exceed 2.0, meaning an unmatched trial would have needed to enroll at least twice as many clusters to achieve the same level of precision as the geographically pair-matched design. We also show that geographically pair-matched designs enable estimation of fine-scale, spatially varying effect heterogeneity under minimal assumptions. Our results demonstrate broad, substantial benefits of geographic pair-matching in large-scale, cluster randomized trials.

Personalized online ensemble machine learning with applications for dynamic data streams.

In this work we introduce the personalized online super learner (POSL), an online personalizable ensemble machine learning algorithm for streaming data. POSL optimizes predictions with respect to baseline covariates, so personalization can vary from completely individualized, that is, optimization with respect to subject ID, to many individuals, that is, optimization with respect to common baseline covariates. As an online algorithm, POSL learns in real time. As a super learner, POSL is grounded in statistical optimality theory and can leverage a diversity of candidate algorithms, including online algorithms with different training and update times, fixed/offline algorithms that are not updated during POSL's fitting procedure, pooled algorithms that learn from many individuals' time series, and individualized algorithms that learn from within a single time series. POSL's ensembling of the candidates can depend on the amount of data collected, the stationarity of the time series, and the mutual characteristics of a group of time series. Depending on the underlying data-generating process and the information available in the data, POSL is able to adapt to learning across samples, through time, or both. For a range of simulations that reflect realistic forecasting scenarios and in a medical application, we examine the performance of POSL relative to other current ensembling and online learning methods. We show that POSL is able to provide reliable predictions for both short and long time series, and it's able to adjust to changing data-generating environments. We further cultivate POSL's practicality by extending it to settings where time series dynamically enter and exit.

A generalization of moderated statistics to data adaptive semiparametric estimation in high-dimensional biology.

The widespread availability of high-dimensional biological data has made the simultaneous screening of many biological characteristics a central problem in computational and high-dimensional biology. As the dimensionality of datasets continues to grow, so too does the complexity of identifying biomarkers linked to exposure patterns. The statistical analysis of such data often relies upon parametric modeling assumptions motivated by convenience, inviting opportunities for model misspecification. While estimation frameworks incorporating flexible, data adaptive regression strategies can mitigate this, their standard variance estimators are often unstable in high-dimensional settings, resulting in inflated Type-I error even after standard multiple testing corrections. We adapt a shrinkage approach compatible with parametric modeling strategies to semiparametric variance estimators of a family of efficient, asymptotically linear estimators of causal effects, defined by counterfactual exposure contrasts. Augmenting the inferential stability of these estimators in high-dimensional settings yields a data adaptive approach for robustly uncovering stable causal associations, even when sample sizes are limited. Our generalized variance estimator is evaluated against appropriate alternatives in numerical experiments, and an open source R/Bioconductor package, biotmle, is introduced. The proposal is demonstrated in an analysis of high-dimensional DNA methylation data from an observational study on the epigenetic effects of tobacco smoking.

A comparison of COVID-19 outbreaks across US Combined Statistical Areas using new methods for estimating R0 and social distancing behaviour.

We investigated the initial outbreak rates and subsequent social distancing behaviour over the initial phase of the COVID-19 pandemic across 29 Combined Statistical Areas (CSAs) of the United States. We used the Numerus Model Builder Data and Simulation Analysis (NMB-DASA) web application to fit the exponential phase of a SCLAIV+D (Susceptible, Contact, Latent, Asymptomatic infectious, symptomatic Infectious, Vaccinated, Dead) disease classes model to outbreaks, thereby allowing us to obtain an estimate of the basic reproductive number R0 for each CSA. Values of R0 ranged from 1.9 to 9.4, with a mean and standard deviation of 4.5±1.8. Fixing the parameters from the exponential fit, we again used NMB-DASA to estimate a set of social distancing behaviour parameters to compute an epidemic flattening index cflatten. Finally, we applied hierarchical clustering methods using this index to divide CSA outbreaks into two clusters: those presenting a social distancing response that was either weaker or stronger. We found cflatten to be more influential in the clustering process than R0. Thus, our results suggest that the behavioural response after a short initial exponential growth phase is likely to be more determinative of the rise of an epidemic than R0 itself.

Influence of community-level sanitation coverage and population density on environmental fecal contamination and child health in a longitudinal cohort in rural Bangladesh.

BACKGROUND: Household-level sanitation interventions have had limited effects on child health or environmental contamination, potentially due to low community coverage. Higher community-level coverage with safely managed sanitation can reduce opportunities for disease transmission. METHODS: We estimated associations between community sanitation coverage, environmental fecal contamination, and child health among 360 compounds in the control arm of the WASH Benefits trial in rural Bangladesh (NCT01590095). In each compound, we enumerated E. coli in environmental samples and recorded the 7-day prevalence of caregiver-reported diarrheal disease and acute respiratory infections (ARI) in children under five. We observed indicators of latrine access and quality among all neighboring compounds within 100 m of study compounds. We defined community coverage as the proportion of neighboring compounds with (1) at least one latrine, and (2) exclusively hygienic latrines (improved facility observed to safely contain feces), within both 50 m and 100 m of study compounds. We assessed effect modification by population density and season. RESULTS: Adjusted for confounders, study compounds surrounded by 100% coverage of at least one latrine per compound within 50 m had slightly lower log10E. coli counts in stored water (Δlog = -0.13, 95% CI -0.26, -0.01), child hand rinses (Δlog = -0.13, 95% CI -0.24, -0.02), and caregiver hand rinses (Δlog = -0.16, 95% CI -0.29, -0.03) and marginally lower prevalence of diarrheal disease (prevalence ratio [PR] = 0.82, 95% CI 0.64, 1.04) and ARI (PR = 0.84, 95% CI 0.69, 1.03) compared to compounds surrounded by <100% coverage. Effects were similar but less pronounced at 100 m. At higher population densities, community latrine coverage was associated with larger reductions in E. coli on child and caregiver hands and prevalence of diarrheal disease. Coverage with exclusively hygienic latrines was not associated with any outcome. CONCLUSION: Higher community sanitation coverage was associated with reduced fecal contamination and improved child health, with stronger effects at highly local scales (50m) and at high population densities. Our findings indicate that the relationship between community sanitation coverage, environmental contamination, and child health varies by definition of coverage, distance, and population density. This work highlights significant uncertainty around how to best measure sanitation coverage and the expected health effects of increasing sanitation coverage using a specific metric. Better understanding of community-level sanitation access is needed to inform policy for implementing sanitation systems that effectively protect community health.