Confounding and exposure measurement error can introduce bias when drawing inference about the marginal effect of an exposure on an outcome of interest. While there are broad methodologies for addressing each source of bias individually, confounding and exposure measurement error frequently co-occur, and there is a need for methods that address them simultaneously. In this paper, corrected score methods are derived under classical additive measurement error to draw inference about marginal exposure effects using only measured variables. Three estimators are proposed based on g-formula, inverse probability weighting, and doubly-robust estimation techniques. The estimators are shown to be consistent and asymptotically normal, and the doubly-robust estimator is shown to exhibit its namesake property. The methods, which are implemented in the R package mismex, perform well in finite samples under both confounding and measurement error as demonstrated by simulation studies. The proposed doubly-robust estimator is applied to study the effects of two biomarkers on HIV-1 infection using data from the HVTN 505 preventative vaccine trial.

Receiver operating characteristic (ROC) analysis is a tool to evaluate the capacity of an ordinal or numeric measure to distinguish between states, often employed in the evaluation of diagnostic tests. Assumptions regarding the nature of the data or the structure of the ROC curve can improve efficiency of estimation for the ROC curve itself or its associated area under the curve (a measure of stochastic ordering closely related to the target parameter of a Mann-Whitney U test). Despite recent interest in methods to conduct comparisons by way of stochastic ordering, ROC modeling strategies are not widely known but may be appreciated by practitioners of these methods. The major goals of this manuscript are to (1) provide an overview of non-parametric, parametric, and semiparametric frameworks for ROC curve estimation, (2) guide the practice of ROC analysis by offering considerations regarding methodological trade-offs, and (3) supply sample R code to further guide implementation. We utilize simulations to demonstrate various features of the methods we discuss. As an illustrative example, we analyze data from a recent cohort study in order to compare responses to SARS-CoV-2 vaccination between solid organ transplant recipients and healthy individuals.

Understanding whether and how treatment effects vary across subgroups is crucial to inform clinical practice and recommendations. Accordingly, the assessment of heterogeneous treatment effects based on pre-specified potential effect modifiers has become a common goal in modern randomized trials. However, when one or more potential effect modifiers are missing, complete-case analysis may lead to bias and under-coverage. While statistical methods for handling missing data have been proposed and compared for individually randomized trials with missing effect modifier data, few guidelines exist for the cluster-randomized setting, where intracluster correlations in the effect modifiers, outcomes, or even missingness mechanisms may introduce further threats to accurate assessment of heterogeneous treatment effect. In this article, the performance of several missing data methods are compared through a simulation study of cluster-randomized trials with continuous outcome and missing binary effect modifier data, and further illustrated using real data from the Work, Family, and Health Study. Our results suggest that multilevel multiple imputation and Bayesian multilevel multiple imputation have better performance than other available methods, and that Bayesian multilevel multiple imputation has lower bias and closer to nominal coverage than standard multilevel multiple imputation when there are model specification or compatibility issues.

Estimands can help clarify the interpretation of treatment effects and ensure that estimators are aligned with the study's objectives. Cluster-randomised trials require additional attributes to be defined within the estimand compared to individually randomised trials, including whether treatment effects are marginal or cluster-specific, and whether they are participant- or cluster-average. In this paper, we provide formal definitions of estimands encompassing both these attributes using potential outcomes notation and describe differences between them. We then provide an overview of estimators for each estimand, describe their assumptions, and show consistency (i.e. asymptotically unbiased estimation) for a series of analyses based on cluster-level summaries. Then, through a re-analysis of a published cluster-randomised trial, we demonstrate that the choice of both estimand and estimator can affect interpretation. For instance, the estimated odds ratio ranged from 1.38 ( p = 0.17) to 1.83 ( p = 0.03) depending on the target estimand, and for some estimands, the choice of estimator affected the conclusions by leading to smaller treatment effect estimates. We conclude that careful specification of the estimand, along with an appropriate choice of estimator, is essential to ensuring that cluster-randomised trials address the right question.

BACKGROUND/OBJECTIVE:Osteosarcoma treatment incorporates chemotherapy and surgery. Resection of the primary tumor usually occurs after induction chemotherapy. Occasionally, scheduling challenges and medical complications result in delay. The goal of this study is to determine if an increased interval between completion of neoadjuvant therapy and surgical resection correlates with decreased tumor necrosis and inferior outcomes in children and young adults with osteosarcoma. DESIGN/METHOD:We conducted a retrospective chart review of 121 patients age less than 40 years diagnosed with osteosarcoma treated at a single tertiary medical center between 2000 and 2022. Inclusion criteria included receipt of two cycles of neoadjuvant methotrexate, cisplatin, and doxorubicin. Association of the interval from completion of induction chemotherapy to resection with tumor necrosis (Spearman's correlation) and outcomes (multivariable Cox hazard regression) were analyzed. RESULTS:There was no significant correlation between interval length and tumor necrosis. However, patients with an interval greater than 16 days had lower 5-year event-free survival (p = 0.019). Multivariable adjusted analysis of patients with initially localized disease revealed that each day increase in interval length corresponds with a 1.1 times greater hazard of having an event (95% CI: 1.0-1.2; p = 0.016). CONCLUSION:Delays in local control were not associated with tumor necrosis. This is consistent with the hypothesis that tumor necrosis is a biologic marker of a tumor's sensitivity to chemotherapy and may not be affected by minor regimen aberrations. However, surgical delay from completion of induction chemotherapy may confer worse outcomes. Longer intervals generally confer worse outcomes in patients with initially localized disease.

Postpartum depression (PPD) affects nearly 20% of postpartum women in Sub-Saharan Africa (SSA), where HIV prevalence is high. Depression is associated with worse HIV outcomes in non-pregnant adults and mental health disorders may worsen HIV outcomes for postpartum women and their infants. PPD is effectively treated with psychosocial or pharmacologic interventions; however, few studies have evaluated the acceptability of treatment modalities in SSA. We analyzed interviews with 23 postpartum women with HIV to assess the acceptability of two depression treatments provided in the context of a randomized trial. Most participants expressed acceptability of treatment randomization and study visit procedures. Participants shared perceptions of high treatment efficacy of their assigned intervention. They reported ongoing HIV and mental health stigma in their communities and emphasized the importance of social support from clinic staff. Our findings suggest a full-scale trial of PPD treatment will be acceptable among women with HIV in Zambia.

BACKGROUND:Data on long-term neurodevelopmental outcomes of normocephalic children (born with normal head circumference) exposed to Zika virus in utero are scarce. We aimed to compare neurodevelopmental outcomes in normocephalic children up to age 48 months with and without Zika virus exposure in utero. METHODS:In this prospective cohort study, we included infants from two cohorts of normocephalic children born in León and Managua, Nicaragua during the 2016 Zika epidemic. In León, all women pregnant during the two enrolment periods were eligible. In Managua, mother-child pairs were included from three districts in the municipality of Managua: all women who became pregnant before June 15, 2016, and had a due date of Sept 15, 2016 or later were eligible. Infants were serologically classified as Zika virus-exposed or Zika virus-unexposed in utero and were followed up prospectively until age 48 months. At 36 months and 48 months of age, the Mullen Scales of Early Learning (MSEL) assessment was administered. Primary outcomes were MSEL early learning composite (ELC) scores at 30-48 months in León and 36-48 months in Managua. We used an inverse probability weighting generalised estimating equations model to assess the effect of Zika virus exposure on individual MSEL cognitive domain scores and ELC scores, adjusted for maternal education and age, poverty status, and infant sex. FINDINGS:The initial enrolment period for the León cohort was between Jan 31 and April 5, 2017 and the second was between Aug 30, 2017, and Feb 22, 2018. The enrolment period for the Managua cohort was between Oct 24, 2019, and May 5, 2020. 478 mothers (482 infants) from the León cohort and 615 mothers (609 infants) from the Managua cohort were enrolled, of whom 622 children (303 from the León cohort; 319 from the Managua cohort) were included in the final analysis; four children had microcephaly at birth and thus were excluded from analyses, two from each cohort. 33 (11%) of 303 children enrolled in León and 219 (69%) of 319 children enrolled in Managua were exposed to Zika virus in utero. In both cohorts, no significant differences were identified in adjusted mean ELC scores between Zika virus-exposed and unexposed infants at 36 months (between-group difference 1·2 points [95% CI -4·2 to 6·5] in the León cohort; 2·8 [-2·4 to 8·1] in the Managua cohort) or at 48 months (-0·9 [-10·8 to 8·8] in the León cohort; 0·1 [-5·1 to 5·2] in the Managua cohort). No differences in ELC scores between Zika virus-exposed and unexposed infants exceeded 6 points at any time between 30 months and 48 months in León or between 36 months and 48 months in Managua, which was considered clinically significant in other settings. INTERPRETATION:We found no significant differences in neurodevelopmental scores between normocephalic children with in-utero Zika virus exposure and Zika virus-unexposed children at age 36 months or 48 months. These findings are promising, supporting typical neurodevelopment in Zika virus-exposed normocephalic children, although additional follow-up and research is warranted. FUNDING:National Institute of Child Health and Development, National Institute of Allergy and Infectious Diseases, and Fogarty International Center. TRANSLATION:For the Spanish translation of the abstract see Supplementary Materials section.

Key Points Question Are patients’ individual health expectations associated with how they subsequently report health-related quality of life? Findings In this cohort study of 207 outpatients with severe chronic obstructive pulmonary disease, most patients were overoptimistic about their future symptom burdens. Overoptimistic expectations of dyspnea burden and negative emotions were associated with worse health-related quality of life over a 24-month study period. Meaning These results suggest that strategies to improve patients’ understanding of their illness and likely future symptom burden may improve decision-making and quality of life in serious illness care.

Background Recent work has shown that cluster-randomised trials can estimate two distinct estimands: the participant-average and cluster-average treatment effects. These can differ when participant outcomes or the treatment effect depends on the cluster size (termed informative cluster size). In this case, estimators that target one estimand (such as the analysis of unweighted cluster-level summaries, which targets the cluster-average effect) may be biased for the other. Furthermore, commonly used estimators such as mixed-effects models or generalised estimating equations with an exchangeable correlation structure can be biased for both estimands. However, there has been little empirical research into whether informative cluster size is likely to occur in practice. Method We re-analysed a cluster-randomised trial comparing two different thresholds for red blood cell transfusion in patients with acute upper gastrointestinal bleeding to explore whether estimates for the participant- and cluster-average effects differed, to provide empirical evidence for whether informative cluster size may be present. For each outcome, we first estimated a participant-average effect using independence estimating equations, which are unbiased under informative cluster size. We then compared this to two further methods: (1) a cluster-average effect estimated using either weighted independence estimating equations or unweighted cluster-level summaries, and (2) estimates from a mixed-effects model or generalised estimating equations with an exchangeable correlation structure. We then performed a small simulation study to evaluate whether observed differences between cluster- and participant-average estimates were likely to occur even if no informative cluster size was present. Results For most outcomes, treatment effect estimates from different methods were similar. However, differences of >10% occurred between participant- and cluster-average estimates for 5 of 17 outcomes (29%). We also observed several notable differences between estimates from mixed-effects models or generalised estimating equations with an exchangeable correlation structure and those based on independence estimating equations. For example, for the EQ-5D VAS score, the independence estimating equation estimate of the participant-average difference was 4.15 (95% confidence interval: −3.37 to 11.66), compared with 2.84 (95% confidence interval: −7.37 to 13.04) for the cluster-average independence estimating equation estimate, and 3.23 (95% confidence interval: −6.70 to 13.16) from a mixed-effects model. Similarly, for thromboembolic/ischaemic events, the independence estimating equation estimate for the participant-average odds ratio was 0.43 (95% confidence interval: 0.07 to 2.48), compared with 0.33 (95% confidence interval: 0.06 to 1.77) from the cluster-average estimator. Conclusion In this re-analysis, we found that estimates from the various approaches could differ, which may be due to the presence of informative cluster size. Careful consideration of the estimand and the plausibility of assumptions underpinning each estimator can help ensure an appropriate analysis methods are used. Independence estimating equations and the analysis of cluster-level summaries (with appropriate weighting for each to correspond to either the participant-average or cluster-average treatment effect) are a desirable choice when informative cluster size is deemed possible, due to their unbiasedness in this setting.

Background Targeting short-term improvements in multicomponent risk scores for mortality in patients with pulmonary arterial hypertension (PAH) could result in improved long-term outcomes. We aimed to determine whether PAH risk scores were adequate surrogates for clinical worsening or mortality outcomes in PAH randomised clinical trials (RCTs).Methods We performed an individual participant data meta-analysis of RCTs selected from PAH trials provided by the US Food and Drug Administration (FDA). We calculated predicted risk using the COMPERA, COMPERA 2.0, non-invasive FPHR, REVEAL 2.0, and REVEAL Lite 2 risk scores. The primary outcome of interest was time to clinical worsening, a composite endpoint composed of any of the following events: all-cause death, hospitalisation for worsening PAH, lung transplantation, atrial septostomy, discontinuation of study treatment (or study withdrawal) for worsening PAH, initiation of parenteral prostacyclin analogue therapy, or decrease of at least 15% in 6-min walk distance from baseline, combined with either worsening of WHO functional class from baseline or the addition of an approved PAH treatment. The secondary outcome of interest was time to all-cause mortality. We assessed the surrogacy of these risk scores, parameterised as attainment of low-risk status by 16 weeks, for improvement in long-term clinical worsening and survival using mediation and meta-analysis frameworks.Findings Of 28 trials received from the FDA, three RCTs (AMBITION, GRIPHON, and SERAPHIN; n=2508) had the data necessary to assess long-term surrogacy. The mean age was 49 years (SD 16), 1956 (78%) participants were women, 1704 (68%) were classified as White, and 280 (11%) were Hispanic or Latino. 1388 (55%) of 2503 participants with available data had idiopathic PAH and 776 (31%) of 2503 had PAH associated with connective tissue disease. In a mediation analysis, the proportions of treatment effects explained by attainment of low-risk status ranged only from 7% to 13%. In a meta-analysis of trial-regions, the treatment effects on low-risk status were not predictive of the treatment effects on time to clinical worsening (R2 values 0 center dot 01-0 center dot 19) nor the treatment effects on time to all-cause mortality (R2 values 0-0 center dot 2). A leave-one-out analysis suggested that the use of these risk scores as surrogates might lead to biased inferences regarding the effect of therapies on clinical outcomes in PAH RCTs. Results were similar when using absolute risk scores at 16 weeks as the potential surrogates.Interpretation Multicomponent risk scores have utility for the prediction of outcomes in patients with PAH. Clinical surrogacy for long-term outcomes cannot be inferred from observational studies of outcomes. Our analyses of three PAH trials with long-term follow-up suggest that further study is necessary before using these or other scores as surrogate outcomes in PAH RCTs or clinical care. Copyright (c) 2023 Elsevier Ltd. All rights reserved.

Rationale: Small trials and professional recommendations support mobilization interventions to improve recovery among critically ill patients, but their real-world effectiveness is unknown. Objective: To evaluate a low-cost, multifaceted mobilization intervention. Methods: We conducted a stepped-wedge cluster-randomized trial across 12 ICUs with diverse case mixes. The primary and secondary samples included patients mechanically ventilated for >= 48 hours who were ambulatory before admission, and all patients with ICU stays >= 48 hours, respectively. The mobilization intervention included 1) designation and posting of daily mobilization goals; 2) interprofessional closed-loop communication coordinated by each ICU's facilitator; and 3) performance feedback. Measurements and Main Results: From March 4, 2019 through March 15, 2020, 848 and 1,069 patients were enrolled in the usual care and intervention phases in the primary sample, respectively. The intervention did not increase the primary outcome, patient's maximal Intensive Care Mobility Scale (range, 0-10) score within 48 hours before ICU discharge (estimated mean difference, 0.16; 95% confidence interval, 20.31 to 0.63; P = 0.51). More patients in the intervention (37.2%) than usual care (30.7%) groups achieved the prespecified secondary outcome of ability to stand before ICU discharge (odds ratio, 1.48; 95% confidence interval, 1.02 to 2.15; P = 0.04). Similar results were observed among the 7,115 patients in the secondary sample. The percentage of days on which patients received physical therapy mediated 90.1% of the intervention effect on standing. ICU mortality (31.5% vs. 29.0%), falls (0.7% vs. 0.4%), and unplanned extubations (2.0% vs. 1.8%) were similar between groups (all P. 0.3). Conclusions: A low-cost, multifaceted mobilization intervention did not improve overall mobility but improved patients' odds of standing and was safe.

The currently recommended dose of dexamethasone for patients with severe or critical COVID-19 is 6 mg per day (mg/d) regardless of patient features and variation. However, patients with severe or critical COVID-19 are heterogenous in many ways (e.g., age, weight, comorbidities, disease severity, and immune features). Thus, it is conceivable that a standardized dosing protocol may not be optimal. We assessed treatment effect heterogeneity in the COVID STEROID 2 trial, which compared 6 mg/d to 12 mg/d, using a causal inference framework with Bayesian Additive Regression Trees, a flexible modeling method that detects interactive effects and nonlinear relationships among multiple patient characteristics simultaneously. We found that 12 mg/d of dexamethasone, relative to 6 mg/d, was probably associated with better long-term outcomes (days alive without life support and mortality after 90 days) among the entire trial population (i.e., no signals of harm), and probably more beneficial among those without diabetes mellitus, that were older, were not using IL-6 inhibitors at baseline, weighed less, or had higher level respiratory support at baseline. This adds more evidence supporting the use of 12 mg/d in practice for most patients not receiving other immunosuppressants and that additional study of dosing could potentially optimize clinical outcomes.

The COVE trial randomized participants to receive two doses of mRNA-1273 vaccine or placebo on Days 1 and 29 (D1, D29). Anti-SARS-CoV-2 Spike IgG binding antibodies (bAbs), anti-receptor binding domain IgG bAbs, 50% inhibitory dilution neutralizing antibody (nAb) titers, and 80% inhibitory dilution nAb titers were measured at D29 and D57. We assessed these markers as correlates of protection (CoPs) against COVID-19 using stochastic interventional vaccine efficacy (SVE) analysis and principal surrogate (PS) analysis, frameworks not used in our previous COVE immune correlates analyses. By SVE analysis, hypothetical shifts of the D57 Spike IgG distribution from a geometric mean concentration (GMC) of 2737 binding antibody units (BAU)/mL (estimated vaccine efficacy (VE): 92.9% (95% CI: 91.7%, 93.9%)) to 274 BAU/mL or to 27,368 BAU/mL resulted in an overall estimated VE of 84.2% (79.0%, 88.1%) and 97.6% (97.4%, 97.7%), respectively. By binary marker PS analysis of Low and High subgroups (cut-point: 2094 BAU/mL), the ignorance interval (IGI) and estimated uncertainty interval (EUI) for VE were [85%, 90%] and (78%, 93%) for Low compared to [95%, 96%] and (92%, 97%) for High. By continuous marker PS analysis, the IGI and 95% EUI for VE at the 2.5th percentile (519.4 BAU/mL) vs. at the 97.5th percentile (9262.9 BAU/mL) of D57 Spike IgG concentration were [92.6%, 93.4%] and (89.2%, 95.7%) vs. [94.3%, 94.6%] and (89.7%, 97.0%). Results were similar for other D29 and D57 markers. Thus, the SVE and PS analyses additionally support all four markers at both time points as CoPs.

Background: Norovirus and sapovirus cause a large burden of acute gastroenteritis (AGE) in young children. We assessed protection conferred by norovirus and sapovirus AGE episodes against future episodes. Methods: Between June 2017 and July 2018, we recruited 444 newborns in León, Nicaragua. Weekly household surveys identified AGE episodes over 36 months, and AGE stools were tested by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) for norovirus genogroup (G)I/GII and sapovirus. We used recurrent-event Cox models and negative control methods to estimate protection conferred by first episodes, controlling for observed and unobserved risk factors, respectively. Results: Sapovirus episodes conferred a 69% reduced hazard of subsequent episodes using the negative control method. Norovirus GI (hazard ratio [HR] = 0.67; 95% confidence interval [CI] = 0.31, 1.3) and GII (HR = 0.20; 95% CI = 0.04, 0.44) episodes also appeared highly protective. Protection against norovirus GII was enhanced following two episodes. Conclusions: Evidence of natural immunity in early childhood provides optimism for the future success of pediatric norovirus and sapovirus vaccines.

BACKGROUND: Pediatric out-of-hospital cardiac arrest results in high morbidity and mortality. Currently, there are no recommended therapies beyond supportive care. The THAPCA-OH (Therapeutic Hypothermia after Pediatric Cardiac Arrest Out-of-Hospital) trial compared hypothermia (33.0°C) with normothermia (36.8°C) in 295 children. Good neurobehavioral outcome and survival at 1 year were higher in the hypothermia group (20 vs. 12% and 38 vs. 29%, respectively). These differences did not meet the planned statistical threshold of P75% for all informative prior integrations with the THAPCA-OH results, except those with the most pessimistic priors. CONCLUSIONS: There is a high probability that hypothermia provides a modest benefit in neurobehavioral outcome and survival at 1 year. (ClinicalTrials.gov number, NCT00878644.)

Introduction Postpartum depression (PPD) is a prevalent and debilitating disease that may affect medication adherence and thus maternal health and vertical transmission among women with HIV. We assessed the feasibility of a trial of interpersonal psychotherapy (IPT) versus antidepressant medication (ADM) to treat PPD and/or anxiety among postpartum women with HIV in Lusaka, Zambia. Methods Between 29 October 2019 and 8 September 2020, we pre-screened women 6-8 weeks after delivery with the Edinburgh Postnatal Depression Scale (EPDS) and diagnosed PPD or anxiety with the Mini International Neuropsychiatric Interview. Consenting participants were randomized 1:1 to up to 11 sessions of IPT or daily self-administered sertraline and followed for 24 weeks. We assessed EPDS score, Clinical Global Impression-Severity of Illness (CGI-S) and medication side effects at each visit and measured maternal HIV viral load at baseline and final study visit. Retention, visit adherence, change in EPDS, CGI-S and log viral load were compared between groups with t-tests and Wilcoxon signed rank tests; we report mean differences, relative risks and 95% confidence intervals. A participant satisfaction survey assessed trial acceptability. Results 78/80 (98%) participants were retained at the final study visit. In the context of the COVID-19 pandemic, visit adherence was greater among women allocated to ADM (9.9 visits, SD 2.2) versus IPT (8.9 visits, SD 2.4; p = 0.06). EPDS scores decreased from baseline to final visit overall, though mean change was greater in the IPT group (-13.8 points, SD 4.7) compared to the ADM group (-11.4 points, SD 5.5; p = 0.04). Both groups showed similar changes in mean log viral load from baseline to final study visit (mean difference -0.43, 95% CI -0.32, 1.18; p = 0.48). In the IPT group, viral load decreased significantly from baseline (0.9 log copies/ml, SD 1.7) to final visit (0.2 log copies/ml, SD 0.9; p = 0.01). Conclusions This pilot study demonstrates that a trial of two forms of PPD treatment is feasible and acceptable among women with HIV in Zambia. IPT and ADM both improved measures of depression severity; however, a full-scale trial is required to determine whether treatment of PPD and anxiety improves maternal-infant HIV outcomes.

The American Zika virus (ZIKV) epidemic has highlighted the need to gain a better understanding of this emerging virus. The goal of this study was to describe the clinical symptoms, laboratory findings, and risk factors for symptomatic ZIKV infection in an area with ongoing transmission of other arboviral infections. We recruited patients at least 2 years of age seeking care at public health centers in León, Nicaragua, between January 2016 and August 2017, for fever, maculopapular rash, and/or nonsuppurative conjunctivitis with a duration of less than 1 week. A laboratory diagnosis of ZIKV was established using a combination of molecular and serological tests. Clinical and laboratory findings and potential risk factors were compared between participants with and without acute ZIKV infection. Fifty-eight (26%) of the 225 participants included in the analysis were found to have acute ZIKV infection. Pregnancy and reports of previous arboviral infection were associated with a higher risk of ZIKV infection. Rash, conjunctivitis, sore throat, and lower absolute neutrophil counts were associated with acute ZIKV infection. The clinical characteristics and risk factors identified were consistent with those identified by previous studies; however, we found sore throat to be a feature of ZIKV infection. We also found that neutrophil counts were lower in ZIKV-infected subjects. These clinical symptoms and laboratory data may help clinicians suspect ZIKV infection during future outbreaks.

BACKGROUND:Sapovirus is increasingly recognized as an important cause of acute gastroenteritis (AGE) in children. We identified risk factors and characterized the clinical profile of sapovirus AGE in a birth cohort in León, Nicaragua.METHODS:We conducted a case-control study nested within a birth cohort (n = 444). Fieldworkers conducted weekly household AGE surveillance. AGE stools were tested for sapovirus by reverse transcriptase quantitative polymerase chain reaction. For each first sapovirus episode, we selected 2 healthy age-matched controls and estimated independent risk factors of sapovirus AGE using conditional logistic regression. We compared clinical characteristics of sapovirus AGE episodes with episodes associated with other etiologies and identified co-infections with other enteric pathogens.RESULTS:From June 2017 to July 2019, we identified 63 first sapovirus AGE episodes and selected 126 controls. Having contact with an individual with AGE symptoms and vaginal delivery were independent risk factors for sapovirus AGE. All cases experienced diarrhea, lasting a median 6 days; 23% experienced vomiting. Compared with children with AGE due to another etiology, sapovirus AGE was similar in severity, with less reported fever. Most cases experienced co-infections and were more likely than controls to be infected with diarrheagenic Escherichia coli or astrovirus.CONCLUSIONS:Sapovirus was a commonly identified AGE etiology in this Central American setting, and symptoms were similar to AGE associated with other etiologies. The association between vaginal delivery and sapovirus is a novel finding. Gut microbiome composition might mediate this relationship, or vaginal delivery might be a proxy for other risk factors. Further investigation into more specific biological mechanisms is warranted.