Search

Malaria-elimination interventions aim to extinguish hotspots and prevent transmission to nearby areas. Here, we re-analyzed a cluster-randomized trial of reactive, focal interventions (chemoprevention using artemether-lumefantrine and/or indoor residual spraying with pirimiphos-methyl) delivered within 500 m of confirmed malaria index cases in Namibia to measure direct effects (among intervention recipients within 500 m) and spillover effects (among non-intervention recipients within 3 km) on incidence, prevalence and seroprevalence. There was no or weak evidence of direct effects, but the sample size of intervention recipients was small, limiting statistical power. There was the strongest evidence of spillover effects of combined chemoprevention and indoor residual spraying. Among non-recipients within 1 km of index cases, the combined intervention reduced malaria incidence by 43% (95% confidence interval, 20-59%). In analyses among non-recipients within 3 km of interventions, the combined intervention reduced infection prevalence by 79% (6-95%) and seroprevalence, which captures recent infections and has higher statistical power, by 34% (20-45%). Accounting for spillover effects increased the cost-effectiveness of the combined intervention by 42%. Targeting hotspots with combined chemoprevention and vector-control interventions can indirectly benefit non-recipients up to 3 km away.

BACKGROUND: While many patients seem to recover from SARS-CoV-2 infections, many patients report experiencing SARS-CoV-2 symptoms for weeks or months after their acute COVID-19 ends, even developing new symptoms weeks after infection. These long-term effects are called post-acute sequelae of SARS-CoV-2 (PASC) or, more commonly, Long COVID. The overall prevalence of Long COVID is currently unknown, and tools are needed to help identify patients at risk for developing long COVID. METHODS: A working group of the Rapid Acceleration of Diagnostics-radical (RADx-rad) program, comprised of individuals from various NIH institutes and centers, in collaboration with REsearching COVID to Enhance Recovery (RECOVER) developed and organized the Long COVID Computational Challenge (L3C), a community challenge aimed at incentivizing the broader scientific community to develop interpretable and accurate methods for identifying patients at risk of developing Long COVID. From August 2022 to December 2022, participants developed Long COVID risk prediction algorithms using the National COVID Cohort Collaborative (N3C) data enclave, a harmonized data repository from over 75 healthcare institutions from across the United States (U.S.). FINDINGS: Over the course of the challenge, 74 teams designed and built 35 Long COVID prediction models using the N3C data enclave. The top 10 teams all scored above a 0.80 Area Under the Receiver Operator Curve (AUROC) with the highest scoring model achieving a mean AUROC of 0.895. Included in the top submission was a visualization dashboard that built timelines for each patient, updating the risk of a patient developing Long COVID in response to clinical events. INTERPRETATION: As a result of L3C, federal reviewers identified multiple machine learning models that can be used to identify patients at risk for developing Long COVID. Many of the teams used approaches in their submissions which can be applied to future clinical prediction questions. FUNDING: Research reported in this RADx® Rad publication was supported by the National Institutes of Health. Timothy Bergquist, Johanna Loomba, and Emily Pfaff were supported by Axle Subcontract: NCATS-STSS-P00438.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P

BackgroundEach year 25 000-32 000 children develop rifampicin- or multidrug-resistant tuberculosis (RR/MDR-TB), and many more require preventive treatment. Levofloxacin is a key component of RR/MDR-TB treatment and prevention, but the existing pharmacokinetic data in children have not yet been comprehensively summarized. We aimed to characterize levofloxacin pharmacokinetics through an individual patient data meta-analysis of available studies and to determine optimal dosing in children.MethodsLevofloxacin concentration and demographic data were pooled from 5 studies and analyzed using nonlinear mixed effects modeling. Simulations were performed using current World Health Organization (WHO)-recommended and model-informed optimized doses. Optimal levofloxacin doses were identified to target median adult area under the time-concentration curve (AUC)24 of 101 mg·h/L given current standard adult doses.ResultsData from 242 children (2.8 years [0.2-16.8] was used). Apparent clearance was 3.16 L/h for a 13-kg child. Age affected clearance, reaching 50% maturation at birth and 90% maturation at 8 months. Nondispersible tablets had 29% lower apparent oral bioavailability compared to dispersible tablets. Median exposures at current WHO-recommended doses were below the AUC target for children weighing 20 mg/kg are required to ensure adequate exposure. Further studies are needed to determine safety and tolerability of these higher doses.

Clofazimine is recommended for the treatment of rifampicin-resistant tuberculosis (RR-TB), but there is currently no verified dosing guideline for its use in children. There is only limited safety and no pharmacokinetic (PK) data available for children. We aimed to characterize clofazimine PK and its relationship with QT-interval prolongation in children. An observational cohort study of South African children

BackgroundLong COVID, also known as post-acute sequelae of COVID-19 (PASC), is a poorly understood condition with symptoms across a range of biological domains that often have debilitating consequences. Some have recently suggested that lingering SARS-CoV-2 virus particles in the gut may impede serotonin production and that low serotonin may drive many Long COVID symptoms across a range of biological systems. Therefore, selective serotonin reuptake inhibitors (SSRIs), which increase synaptic serotonin availability, may be used to prevent or treat Long COVID. SSRIs are commonly prescribed for depression, therefore restricting a study sample to only include patients with depression can reduce the concern of confounding by indication.MethodsIn an observational sample of electronic health records from patients in the National COVID Cohort Collaborative (N3C) with a COVID-19 diagnosis between September 1, 2021, and December 1, 2022, and a comorbid depressive disorder, the leading indication for SSRI use, we evaluated the relationship between SSRI use during acute COVID-19 and subsequent 12-month risk of Long COVID (defined by ICD-10 code U09.9). We defined SSRI use as a prescription for SSRI medication beginning at least 30 days before acute COVID-19 and not ending before SARS-CoV-2 infection. To minimize bias, we estimated relationships using nonparametric targeted maximum likelihood estimation to aggressively adjust for high-dimensional covariates.ResultsWe analyzed a sample (n = 302,626) of patients with a diagnosis of a depressive condition before COVID-19 diagnosis, where 100,803 (33%) were using an SSRI. We found that SSRI users had a significantly lower risk of Long COVID compared to nonusers (adjusted causal relative risk 0.92, 95% CI (0.86, 0.99)) and we found a similar relationship comparing new SSRI users (first SSRI prescription 1 to 4 months before acute COVID-19 with no prior history of SSRI use) to nonusers (adjusted causal relative risk 0.89, 95% CI (0.80, 0.98)).ConclusionsThese findings suggest that SSRI use during acute COVID-19 may be protective against Long COVID, supporting the hypothesis that serotonin may be a key mechanistic biomarker of Long COVID.

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

In biomedical science, analyzing treatment effect heterogeneity plays an essential role in assisting personalized medicine. The main goals of analyzing treatment effect heterogeneity include estimating treatment effects in clinically relevant subgroups and predicting whether a patient subpopulation might benefit from a particular treatment. Conventional approaches often evaluate the subgroup treatment effects via parametric modeling and can thus be susceptible to model mis-specifications. In this paper, we take a model-free semiparametric perspective and aim to efficiently evaluate the heterogeneous treatment effects of multiple subgroups simultaneously under the one-step targeted maximum-likelihood estimation (TMLE) framework. When the number of subgroups is large, we further expand this path of research by looking at a variation of the one-step TMLE that is robust to the presence of small estimated propensity scores in finite samples. From our simulations, our method demonstrates substantial finite sample improvements compared to conventional methods. In a case study, our method unveils the potential treatment effect heterogeneity of rs12916-T allele (a proxy for statin usage) in decreasing Alzheimers disease risk.

Across research disciplines, cluster randomized trials (CRTs) are commonly implemented to evaluate interventions delivered to groups of participants, such as communities and clinics. Despite advances in the design and analysis of CRTs, several challenges remain. First, there are many possible ways to specify the causal effect of interest (eg, at the individual-level or at the cluster-level). Second, the theoretical and practical performance of common methods for CRT analysis remain poorly understood. Here, we present a general framework to formally define an array of causal effects in terms of summary measures of counterfactual outcomes. Next, we provide a comprehensive overview of CRT estimators, including the t-test, generalized estimating equations (GEE), augmented-GEE, and targeted maximum likelihood estimation (TMLE). Using finite sample simulations, we illustrate the practical performance of these estimators for different causal effects and when, as commonly occurs, there are limited numbers of clusters of different sizes. Finally, our application to data from the Preterm Birth Initiative (PTBi) study demonstrates the real-world impact of varying cluster sizes and targeting effects at the cluster-level or at the individual-level. Specifically, the relative effect of the PTBi intervention was 0.81 at the cluster-level, corresponding to a 19% reduction in outcome incidence, and was 0.66 at the individual-level, corresponding to a 34% reduction in outcome risk. Given its flexibility to estimate a variety of user-specified effects and ability to adaptively adjust for covariates for precision gains while maintaining Type-I error control, we conclude TMLE is a promising tool for CRT analysis.

BACKGROUND: The Targeted Learning roadmap provides a systematic guide for generating and evaluating real-world evidence (RWE). From a regulatory perspective, RWE arises from diverse sources such as randomized controlled trials that make use of real-world data, observational studies, and other study designs. This paper illustrates a principled approach to assessing the validity and interpretability of RWE. METHODS: We applied the roadmap to a published observational study of the dose-response association between ritodrine hydrochloride and pulmonary edema among women pregnant with twins in Japan. The goal was to identify barriers to causal effect estimation beyond unmeasured confounding reported by the studys authors, and to explore potential options for overcoming the barriers that robustify results. RESULTS: Following the roadmap raised issues that led us to formulate alternative causal questions that produced more reliable, interpretable RWE. The process revealed a lack of information in the available data to identify a causal dose-response curve. However, under explicit assumptions the effect of treatment with any amount of ritodrine versus none, albeit a less ambitious parameter, can be estimated from data. CONCLUSIONS: Before RWE can be used in support of clinical and regulatory decision-making, its quality and reliability must be systematically evaluated. The TL roadmap prescribes how to carry out a thorough, transparent, and realistic assessment of RWE. We recommend this approach be a routine part of any decision-making process.

Distinct morphologies of the mitochondrial network support divergent metabolic and regulatory processes that determine cell function and fate1-3. The mechanochemical GTPase optic atrophy 1 (OPA1) influences the architecture of cristae and catalyses the fusion of the mitochondrial inner membrane4,5. Despite its fundamental importance, the molecular mechanisms by which OPA1 modulates mitochondrial morphology are unclear. Here, using a combination of cellular and structural analyses, we illuminate the molecular mechanisms that are key to OPA1-dependent membrane remodelling and fusion. Human OPA1 embeds itself into cardiolipin-containing membranes through a lipid-binding paddle domain. A conserved loop within the paddle domain inserts deeply into the bilayer, further stabilizing the interactions with cardiolipin-enriched membranes. OPA1 dimerization through the paddle domain promotes the helical assembly of a flexible OPA1 lattice on the membrane, which drives mitochondrial fusion in cells. Moreover, the membrane-bending OPA1 oligomer undergoes conformational changes that pull the membrane-inserting loop out of the outer leaflet and contribute to the mechanics of membrane remodelling. Our findings provide a structural framework for understanding how human OPA1 shapes mitochondrial morphology and show us how human disease mutations compromise OPA1 functions.

Questions remain regarding the effect of baseline host and exposure factors on vaccine efficacy (VE) across pathogens and vaccine platforms. We report placebo-controlled data from four Phase 3 COVID-19 trials during the early period of the pandemic. This was a cross-protocol analysis of four randomized, placebo-controlled efficacy trials (Moderna/mRNA1273, AstraZeneca/AZD1222, Janssen/Ad26.COV2.S, and Novavax/NVX-CoV2373) using a harmonized design. Trials were conducted in the United States and international sites in adults ≥ 18 years of age. VE was assessed for symptomatic and severe COVID-19. We analyzed 114,480 participants from both placebo and vaccine arms, enrolled July 2020 to February 2021, with follow up through July 2021. VE against symptomatic COVID-19 showed little heterogeneity across baseline socio-demographic, clinical or exposure characteristics, in either univariate or multivariate analysis, regardless of vaccine platform. Similarly, VE against severe COVID-19 in the single trial (Janssen) with sufficient endpoints for analysis showed little evidence of heterogeneity. COVID-19 VE is not influenced by baseline host or exposure characteristics across efficacy trials of different vaccine platforms and countries when well matched to circulating virus strains. This supports use of these vaccines, regardless of platform type, as effective tools in the near term for reducing symptomatic and severe COVID-19, particularly for older individuals and those with common co-morbidities during major variant shifts. Clinical trial registration numbers: NCT04470427, NCT04516746, NCT04505722, and NCT04611802.

The optimal dynamic treatment rule (ODTR) framework offers an approach for understanding which kinds of patients respond best to specific treatments - in other words, treatment effect heterogeneity. Recently, there has been a proliferation of methods for estimating the ODTR. One such method is an extension of the SuperLearner algorithm - an ensemble method to optimally combine candidate algorithms extensively used in prediction problems - to ODTRs. Following the ``causal roadmap," we causally and statistically define the ODTR and provide an introduction to estimating it using the ODTR SuperLearner. Additionally, we highlight practical choices when implementing the algorithm, including choice of candidate algorithms, metalearners to combine the candidates, and risk functions to select the best combination of algorithms. Using simulations, we illustrate how estimating the ODTR using this SuperLearner approach can uncover treatment effect heterogeneity more effectively than traditional approaches based on fitting a parametric regression of the outcome on the treatment, covariates and treatment-covariate interactions. We investigate the implications of choices in implementing an ODTR SuperLearner at various sample sizes. Our results show the advantages of: (1) including a combination of both flexible machine learning algorithms and simple parametric estimators in the library of candidate algorithms; (2) using an ensemble metalearner to combine candidates rather than selecting only the best-performing candidate; (3) using the mean outcome under the rule as a risk function. Finally, we apply the ODTR SuperLearner to the ``Interventions" study, an ongoing randomized controlled trial, to identify which justice-involved adults with mental illness benefit most from cognitive behavioral therapy to reduce criminal re-offending.

Cluster randomized trials (CRTs) randomly assign an intervention to groups of individuals (e.g., clinics or communities) and measure outcomes on individuals in those groups. While offering many advantages, this experimental design introduces challenges that are only partially addressed by existing analytic approaches. First, outcomes are often missing for some individuals within clusters. Failing to appropriately adjust for differential outcome measurement can result in biased estimates and inference. Second, CRTs often randomize limited numbers of clusters, resulting in chance imbalances on baseline outcome predictors between arms. Failing to adaptively adjust for these imbalances and other predictive covariates can result in efficiency losses. To address these methodological gaps, we propose and evaluate a novel two-stage targeted minimum loss-based estimator to adjust for baseline covariates in a manner that optimizes precision, after controlling for baseline and postbaseline causes of missing outcomes. Finite sample simulations illustrate that our approach can nearly eliminate bias due to differential outcome measurement, while existing CRT estimators yield misleading results and inferences. Application to real data from the SEARCH community randomized trial demonstrates the gains in efficiency afforded through adaptive adjustment for baseline covariates, after controlling for missingness on individual-level outcomes.

BackgroundAdverse childhood experiences (ACEs) increase the risk of poor health outcomes later in life. Psychosocial stressors may also have intergenerational health effects by which parental ACEs are associated with mental and physical health of children. Epigenetic programming may be one mechanism linking parental ACEs to child health. This study aimed to investigate epigenome-wide associations of maternal preconception ACEs with DNA methylation patterns of children. In the Center for the Health Assessment of Mothers and Children of Salinas study, cord blood DNA methylation was measured using the Illumina HumanMethylation450 BeadChip. Preconception ACEs, which occurred during the mothers' childhoods, were collected using a standard ACE questionnaire including 10 ACE indicators. Maternal ACE exposures were defined in this study as (1) the total number of ACEs; (2) the total number of ACEs categorized as 0, 1-3, and > 4; and (3) individual ACEs. Associations of ACE exposures with differential methylated positions, regions, and CpG modules determined using weighted gene co-expression network analysis were evaluated adjusting for covariates.ResultsData on maternal ACEs and cord blood DNA methylation were available for 196 mother/newborn pairs. One differential methylated position was associated with maternal experience of emotional abuse (cg05486260/FAM135B gene; q value

Gestational age (GA) is an important determinant of child health and disease risk. Two epigenetic GA clocks have been developed using DNA methylation (DNAm) patterns in cord blood. We investigate the accuracy of GA clocks and determinants of epigenetic GA acceleration (GAA), a biomarker of biological ageing. We hypothesize that prenatal and birth characteristics are associated with altered GAA, thereby disrupting foetal biological ageing. We examined 372 mother-child pairs from the Center for the Health Assessment of Mothers and Children of Salinas study of primarily Latino farmworkers in California. Chronological GA was robustly correlated with epigenetic GA (DNAm GA) estimated by the Knight (r = 0.48, p

Measuring immune correlates of disease acquisition and protection in the context of a clinical trial is a prerequisite for improved vaccine design. We analysed binding and neutralizing antibody measurements 4 weeks post vaccination as correlates of risk of moderate to severe-critical COVID-19 through 83 d post vaccination in the phase 3, double-blind placebo-controlled phase of ENSEMBLE, an international randomized efficacy trial of a single dose of Ad26.COV2.S. We also evaluated correlates of protection in the trial cohort. Of the three antibody immune markers we measured, we found most support for 50% inhibitory dilution (ID50) neutralizing antibody titre as a correlate of risk and of protection. The outcome hazard ratio was 0.49 (95% confidence interval 0.29, 0.81; P = 0.006) per 10-fold increase in ID50; vaccine efficacy was 60% (43%, 72%) at non-quantifiable ID50 (

BackgroundPharmacokinetic data for bedaquiline in children are limited. We described the pharmacokinetics and safety of bedaquiline in South African children and adolescents receiving treatment for multidrug/rifampicin-resistant tuberculosis (MDR/RR-TB) in routine care.MethodsIn this observational cohort study, children aged 6-17 years receiving bedaquiline at recommended doses as part of MDR/RR-TB treatment underwent semi-intensive pharmacokinetic sampling. Bedaquiline and the M2 metabolite plasma concentrations were quantified, and nonlinear mixed-effects modeling performed. Pediatric data were described using a pre-established model of bedaquiline pharmacokinetics in adults. The exposure reference was 187 µg ⋅ h/mL, the median weekly area under the curve (AUC) of adults at week 24 of treatment with bedaquiline. Safety was assessed through monthly clinical, blood and electrocardiogram monitoring, and treatment outcomes described.ResultsFifteen children (3 human immunodeficiency virus [HIV]-positive) with median age 13.3 years (range 6.5-16.3) were included. A bedaquiline pharmacokinetic model was adapted to be allometrically scaled in clearance and volume, centered in the median child population weight. Bedaquiline bioavailability was 57% of that in adults. Overall bedaquiline exposures were below target, and AUC reference attainment was achieved in only 3 (20%) children. Ten children experienced 27 adverse events at least possibly related to bedaquiline; no adverse events led to bedaquiline withdrawal. Two adverse events (arthritis and arthralgia) were considered severe, and 2 children had mild QT interval corrected for heart rate using Fridericia's formula (QT) prolongation.ConclusionsThe evaluated doses of bedaquiline in children ≥ 6 years of age were safe but achieved slightly lower plasma concentrations compared to adults receiving the recommended dose, possibly due to delayed food intake relative to bedaquiline administration.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value

Several recently developed methods have the potential to harness machine learning in the pursuit of target quantities inspired by causal inference, including inverse weighting, doubly robust estimating equations and substitution estimators like targeted maximum likelihood estimation. There are even more recent augmentations of these procedures that can increase robustness, by adding a layer of cross-validation (cross-validated targeted maximum likelihood estimation and double machine learning, as applied to substitution and estimating equation approaches, respectively). While these methods have been evaluated individually on simulated and experimental data sets, a comprehensive analysis of their performance across real data based simulations have yet to be conducted. In this work, we benchmark multiple widely used methods for estimation of the average treatment effect using ten different nutrition intervention studies data. A nonparametric regression method, undersmoothed highly adaptive lasso, is used to generate the simulated distribution which preserves important features from the observed data and reproduces a set of true target parameters. For each simulated data, we apply the methods above to estimate the average treatment effects as well as their standard errors and resulting confidence intervals. Based on the analytic results, a general recommendation is put forth for use of the cross-validated variants of both substitution and estimating equation estimators. We conclude that the additional layer of cross-validation helps in avoiding unintentional over-fitting of nuisance parameter functionals and leads to more robust inferences.

BackgroundMoxifloxacin is a recommended drug for rifampin-resistant tuberculosis (RR-TB) treatment, but there is limited pediatric pharmacokinetic and safety data, especially in young children. We characterize moxifloxacin population pharmacokinetics and QT interval prolongation and evaluate optimal dosing in children with RR-TB.MethodsPharmacokinetic data were pooled from 2 observational studies in South African children with RR-TB routinely treated with oral moxifloxacin once daily. The population pharmacokinetics and Fridericia-corrected QT (QTcF)-interval prolongation were characterized in NONMEM. Pharmacokinetic simulations were performed to predict expected exposure and optimal weight-banded dosing.ResultsEighty-five children contributed pharmacokinetic data (median [range] age of 4.6 [0.8-15] years); 16 (19%) were aged

Epigenetic aging biomarkers are associated with increased morbidity and mortality. We evaluated if occupational exposure to three established chemical carcinogens is associated with acceleration of epigenetic aging. We studied workers in China occupationally exposed to benzene, trichloroethylene (TCE) or formaldehyde by measuring personal air exposures prior to blood collection. Unexposed controls matched by age and sex were selected from nearby factories. We measured leukocyte DNA methylation (DNAm) in peripheral white blood cells using the Infinium HumanMethylation450 BeadChip to calculate five epigenetic aging clocks and DNAmTL, a biomarker associated with leukocyte telomere length and cell replication. We tested associations between exposure intensity and epigenetic age acceleration (EAA), defined as the residuals of regressing the DNAm aging biomarker on chronological age, matching factors and potential confounders. Median differences in EAA between exposure groups were tested using a permutation test with exact p-values. Epigenetic clocks were strongly correlated with age (Spearman r > 0.8) in all three occupational studies. There was a positive exposure-response relationship between benzene and the Skin-Blood Clock EAA biomarker: median EAA was -0.91 years in controls (n = 44), 0.78 years in workers exposed to

Emerging research suggests associations of physical and psychosocial stressors with epigenetic aging. Although this work has included early-life exposures, data on maternal exposures and epigenetic aging of their children remain sparse. Using longitudinally collected data from the California, Salinas Valley CHAMACOS study, we examined relationships between maternal Adverse Childhood Experiences (ACEs) reported up to 18 years of life, prior to pregnancy, with eight measures (Horvath, Hannum, SkinBloodClock, Intrinsic, Extrinsic, PhenoAge, GrimAge, and DNAm telomere length) of blood leukocyte epigenetic age acceleration (EAA) in their children at ages 7, 9, and 14 years (N = 238 participants with 483 observations). After adjusting for maternal chronological age at delivery, pregnancy smoking/alcohol use, parity, child gestational age, and estimated leukocyte proportions, higher maternal ACEs were significantly associated with at least a 0.76-year increase in child Horvath and Intrinsic EAA. Higher maternal ACEs were also associated with a 0.04 kb greater DNAm estimate of telomere length of children. Overall, our data suggests that maternal preconception ACEs are associated with biological aging in their offspring in childhood and that preconception ACEs have differential relationships with EAA measures, suggesting different physiologic utilities of EEA measures. Studies are necessary to confirm these findings and to elucidate potential pathways to explain these relationships, which may include intergenerational epigenetic inheritance and persistent physical and social exposomes.

Metformin and weight loss relationships with epigenetic age measures-biological aging biomarkers-remain understudied. We performed a post-hoc analysis of a randomized controlled trial among overweight/obese breast cancer survivors (N = 192) assigned to metformin, placebo, weight loss with metformin, or weight loss with placebo interventions for 6 months. Epigenetic age was correlated with chronological age (r = 0.20-0.86; P

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

BackgroundIn generalized epidemic settings, strategies are needed to prioritize individuals at higher risk of human immunodeficiency virus (HIV) acquisition for prevention services. We used population-level HIV testing data from rural Kenya and Uganda to construct HIV risk scores and assessed their ability to identify seroconversions.MethodsDuring 2013-2017, >75% of residents in 16 communities in the SEARCH study were tested annually for HIV. In this population, we evaluated 3 strategies for using demographic factors to predict the 1-year risk of HIV seroconversion: membership in ≥1 known "risk group" (eg, having a spouse living with HIV), a "model-based" risk score constructed with logistic regression, and a "machine learning" risk score constructed with the Super Learner algorithm. We hypothesized machine learning would identify high-risk individuals more efficiently (fewer persons targeted for a fixed sensitivity) and with higher sensitivity (for a fixed number targeted) than either other approach.ResultsA total of 75 558 persons contributed 166 723 person-years of follow-up; 519 seroconverted. Machine learning improved efficiency. To achieve a fixed sensitivity of 50%, the risk-group strategy targeted 42% of the population, the model-based strategy targeted 27%, and machine learning targeted 18%. Machine learning also improved sensitivity. With an upper limit of 45% targeted, the risk-group strategy correctly classified 58% of seroconversions, the model-based strategy 68%, and machine learning 78%.ConclusionsMachine learning improved classification of individuals at risk of HIV acquisition compared with a model-based approach or reliance on known risk groups and could inform targeting of prevention strategies in generalized epidemic settings.Clinical trials registrationNCT01864603.

Astronauts undertaking long-duration space missions may be vulnerable to unique stressors that can impact human aging. Nevertheless, few studies have examined the relationship of mission duration with DNA-methylation-based biomarkers of aging in astronauts. Using data from the six participants of the Mars-500 mission, a high-fidelity 520-day ground simulation experiment, we tested relationships of mission duration with five longitudinally measured blood DNA-methylation-based metrics: DNAmGrimAge, DNAmPhenoAge, DNA-methylation-based estimator of telomere length (DNAmTL), mitotic divisions (epigenetic mitotic clock [epiTOC2]), and pace of aging (PoA). We provide evidence that, relative to baseline, mission duration was associated with significant decreases in epigenetic aging. However, only decreases in DNAmPhenoAge remained significant 7 days post-mission. We also observed significant changes in estimated proportions of plasmablasts, CD4T, CD8 naive, and natural killer (NK) cells. Only decreases in NK cells remained significant post-mission. If confirmed more broadly, these findings contribute insights to improve the understanding of the biological aging implications for individuals experiencing long-duration space travel.

BackgroundPopulation-level estimates of disease prevalence and control are needed to assess prevention and treatment strategies. However, available data often suffer from differential missingness. For example, population-level HIV viral suppression is the proportion of all HIV-positive persons with suppressed viral replication. Individuals with measured HIV status, and among HIV-positive individuals those with measured viral suppression, likely differ from those without such measurements.MethodsWe discuss three sets of assumptions to identify population-level suppression in the intervention arm of the SEARCH Study (NCT01864603), a community randomized trial in rural Kenya and Uganda (2013-2017). Using data on nearly 100,000 participants, we compare estimates from (1) an unadjusted approach assuming data are missing-completely-at-random (MCAR); (2) stratification on age group, sex, and community; and (3) targeted maximum likelihood estimation to adjust for a larger set of baseline and time-updated variables.ResultsDespite high measurement coverage, estimates of population-level viral suppression varied by identification assumption. Unadjusted estimates were most optimistic: 50% (95% confidence interval [CI] = 46%, 54%) of HIV-positive persons suppressed at baseline, 80% (95% CI = 78%, 82%) at year 1, 85% (95% CI = 83%, 86%) at year 2, and 85% (95% CI = 83%, 87%) at year 3. Stratifying on baseline predictors yielded slightly lower estimates, and full adjustment reduced estimates meaningfully: 42% (95% CI = 37%, 46%) of HIV-positive persons suppressed at baseline, 71% (95% CI = 69%, 73%) at year 1, 76% (95% CI = 74%, 78%) at year 2, and 79% (95% CI = 77%, 81%) at year 3.ConclusionsEstimation of population-level disease burden and control requires appropriate adjustment for missing data. Even in large studies with limited missingness, estimates relying on the MCAR assumption or baseline stratification should be interpreted cautiously.

We propose a method for summarizing the strength of association between a set of variables and a multivariate outcome. Classical summary measures are appropriate when linear relationships exist between covariates and outcomes, while our approach provides an alternative that is useful in situations where complex relationships may be present. We utilize machine learning to detect nonlinear relationships and covariate interactions and propose a measure of association that captures these relationships. A hypothesis test about the proposed associative measure can be used to test the strong null hypothesis of no association between a set of variables and a multivariate outcome. Simulations demonstrate that this hypothesis test has greater power than existing methods against alternatives where covariates have nonlinear relationships with outcomes. We additionally propose measures of variable importance for groups of variables, which summarize each groups' association with the outcome. We demonstrate our methodology using data from a birth cohort study on childhood health and nutrition in the Philippines.

A number of sophisticated estimators of longitudinal effects have been proposed for estimating the intervention-specific mean outcome. However, there is a relative paucity of research comparing these methods directly to one another. In this study, we compare various approaches to estimating a causal effect in a longitudinal treatment setting using both simulated data and data measured from a human immunodeficiency virus cohort. Six distinct estimators are considered: (i) an iterated conditional expectation representation, (ii) an inverse propensity weighted method, (iii) an augmented inverse propensity weighted method, (iv) a double robust iterated conditional expectation estimator, (v) a modified version of the double robust iterated conditional expectation estimator, and (vi) a targeted minimum loss-based estimator. The details of each estimator and its implementation are presented along with nuisance parameter estimation details, which include potentially pooling the observed data across all subjects regardless of treatment history and using data adaptive machine learning algorithms. Simulations are constructed over six time points, with each time point steadily increasing in positivity violations. Estimation is carried out for both the simulations and applied example using each of the six estimators under both stratified and pooled approaches of nuisance parameter estimation. Simulation results show that double robust estimators remained without meaningful bias as long as at least one of the two nuisance parameters were estimated with a correctly specified model. Under full misspecification, the bias of the double robust estimators remained better than that of the inverse propensity estimator under misspecification, but worse than the iterated conditional expectation estimator. Weighted estimators tended to show better performance than the covariate estimators. As positivity violations increased, the mean squared error and bias of all estimators considered became worse, with covariate-based double robust estimators especially susceptible. Applied analyses showed similar estimates at most time points, with the important exception of the inverse propensity estimator which deviated markedly as positivity violations increased. Given its efficiency, ability to respect the parameter space, and observed performance, we recommend the pooled and weighted targeted minimum loss-based estimator.

The assumption that no subject's exposure affects another subject's outcome, known as the no-interference assumption, has long held a foundational position in the study of causal inference. However, this assumption may be violated in many settings, and in recent years has been relaxed considerably. Often this has been achieved with either the aid of a known underlying network, or the assumption that the population can be partitioned into separate groups, between which there is no interference, and within which each subject's outcome may be affected by all the other subjects in the group via the proportion exposed (the stratified interference assumption). In this article, we instead consider a complete interference setting, in which each subject affects every other subject's outcome. In particular, we make the stratified interference assumption for a single group consisting of the entire sample. We show that a targeted maximum likelihood estimator for the i.i.d. setting can be used to estimate a class of causal parameters that includes direct effects and overall effects under certain interventions. This estimator remains doubly-robust, semiparametric efficient, and continues to allow for incorporation of machine learning under our model. We conduct a simulation study, and present results from a data application where we study the effect of a nurse-based triage system on the outcomes of patients receiving HIV care in Kenyan health clinics.

BackgroundUniversal antiretroviral therapy (ART) with annual population testing and a multidisease, patient-centered strategy could reduce new human immunodeficiency virus (HIV) infections and improve community health.MethodsWe randomly assigned 32 rural communities in Uganda and Kenya to baseline HIV and multidisease testing and national guideline-restricted ART (control group) or to baseline testing plus annual testing, eligibility for universal ART, and patient-centered care (intervention group). The primary end point was the cumulative incidence of HIV infection at 3 years. Secondary end points included viral suppression, death, tuberculosis, hypertension control, and the change in the annual incidence of HIV infection (which was evaluated in the intervention group only).ResultsA total of 150,395 persons were included in the analyses. Population-level viral suppression among 15,399 HIV-infected persons was 42% at baseline and was higher in the intervention group than in the control group at 3 years (79% vs. 68%; relative prevalence, 1.15; 95% confidence interval [CI], 1.11 to 1.20). The annual incidence of HIV infection in the intervention group decreased by 32% over 3 years (from 0.43 to 0.31 cases per 100 person-years; relative rate, 0.68; 95% CI, 0.56 to 0.84). However, the 3-year cumulative incidence (704 incident HIV infections) did not differ significantly between the intervention group and the control group (0.77% and 0.81%, respectively; relative risk, 0.95; 95% CI, 0.77 to 1.17). Among HIV-infected persons, the risk of death by year 3 was 3% in the intervention group and 4% in the control group (0.99 vs. 1.29 deaths per 100 person-years; relative risk, 0.77; 95% CI, 0.64 to 0.93). The risk of HIV-associated tuberculosis or death by year 3 among HIV-infected persons was 4% in the intervention group and 5% in the control group (1.19 vs. 1.50 events per 100 person-years; relative risk, 0.79; 95% CI, 0.67 to 0.94). At 3 years, 47% of adults with hypertension in the intervention group and 37% in the control group had hypertension control (relative prevalence, 1.26; 95% CI, 1.15 to 1.39).ConclusionsUniversal HIV treatment did not result in a significantly lower incidence of HIV infection than standard care, probably owing to the availability of comprehensive baseline HIV testing and the rapid expansion of ART eligibility in the control group. (Funded by the National Institutes of Health and others; SEARCH ClinicalTrials.gov number, NCT01864603.).

We often seek to estimate the impact of an exposure naturally occurring or randomly assigned at the cluster-level. For example, the literature on neighborhood determinants of health continues to grow. Likewise, community randomized trials are applied to learn about real-world implementation, sustainability, and population effects of interventions with proven individual-level efficacy. In these settings, individual-level outcomes are correlated due to shared cluster-level factors, including the exposure, as well as social or biological interactions between individuals. To flexibly and efficiently estimate the effect of a cluster-level exposure, we present two targeted maximum likelihood estimators (TMLEs). The first TMLE is developed under a non-parametric causal model, which allows for arbitrary interactions between individuals within a cluster. These interactions include direct transmission of the outcome (i.e. contagion) and influence of one individual's covariates on another's outcome (i.e. covariate interference). The second TMLE is developed under a causal sub-model assuming the cluster-level and individual-specific covariates are sufficient to control for confounding. Simulations compare the alternative estimators and illustrate the potential gains from pairing individual-level risk factors and outcomes during estimation, while avoiding unwarranted assumptions. Our results suggest that estimation under the sub-model can result in bias and misleading inference in an observational setting. Incorporating working assumptions during estimation is more robust than assuming they hold in the underlying causal model. We illustrate our approach with an application to HIV prevention and treatment.

ObjectivePrior studies around the relationship between smoking and rheumatoid arthritis (RA) disease activity have reported inconsistent findings, which may be ascribed to heterogeneous study designs or biases in statistical analyses. We examined the association between smoking and RA outcomes using statistical methods that account for time-varying confounding and loss to followup.MethodsWe included 282 individuals with an RA diagnosis using electronic health record data collected at a public hospital between 2013 and 2017. Current smoking status and disease activity were assessed at each visit; covariates included sex, race/ethnicity, age, obesity, and medication use. We used longitudinal targeted maximum likelihood estimation to estimate the causal effect of smoking on disease activity measures at 27 months, and compared results to conventional longitudinal methods.ResultsSmoking was associated with an increase of 0.64 units in the patient global score compared to nonsmoking (p = 0.01), and with 2.58 more swollen joints (p < 0.001). While smoking was associated with a higher clinical disease activity score (2.11), the difference was not statistically significant (p = 0.22). We found no association between smoking and physician global score, or C-reactive protein levels, and an inverse association between smoking and tender joint count (p = 0.05). Analyses using conventional methods showed a null relationship for all outcomes.ConclusionSmoking is associated with higher levels of disease activity in RA. Causal methods may be useful for investigations of additional exposures on longitudinal outcome measures in rheumatologic disease.

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

Binary classification problems are ubiquitous in health and social sciences. In many cases, one wishes to balance two competing optimality considerations for a binary classifier. For instance, in resource-limited settings, an human immunodeficiency virus prevention program based on offering pre-exposure prophylaxis (PrEP) to select high-risk individuals must balance the sensitivity of the binary classifier in detecting future seroconverters (and hence offering them PrEP regimens) with the total number of PrEP regimens that is financially and logistically feasible for the program. In this article, we consider a general class of constrained binary classification problems wherein the objective function and the constraint are both monotonic with respect to a threshold. These include the minimization of the rate of positive predictions subject to a minimum sensitivity, the maximization of sensitivity subject to a maximum rate of positive predictions, and the Neyman-Pearson paradigm, which minimizes the type II error subject to an upper bound on the type I error. We propose an ensemble approach to these binary classification problems based on the Super Learner methodology. This approach linearly combines a user-supplied library of scoring algorithms, with combination weights and a discriminating threshold chosen to minimize the constrained optimality criterion. We then illustrate the application of the proposed classifier to develop an individualized PrEP targeting strategy in a resource-limited setting, with the goal of minimizing the number of PrEP offerings while achieving a minimum required sensitivity. This proof of concept data analysis uses baseline data from the ongoing Sustainable East Africa Research in Community Health study. Copyright © 2017 John Wiley & Sons, Ltd.

Standard statistical practice used for determining the relative importance of competing causes of disease typically relies on ad hoc methods, often byproducts of machine learning procedures (stepwise regression, random forest, etc.). Causal inference framework and data-adaptive methods may help to tailor parameters to match the clinical question and free one from arbitrary modeling assumptions. Our focus is on implementations of such semiparametric methods for a variable importance measure (VIM). We propose a fully automated procedure for VIM based on collaborative targeted maximum likelihood estimation (cTMLE), a method that optimizes the estimate of an association in the presence of potentially numerous competing causes. We applied the approach to data collected from traumatic brain injury patients, specifically a prospective, observational study including three US Level-1 trauma centers. The primary outcome was a disability score (Glasgow Outcome Scale - Extended (GOSE)) collected three months post-injury. We identified clinically important predictors among a set of risk factors using a variable importance analysis based on targeted maximum likelihood estimators (TMLE) and on cTMLE. Via a parametric bootstrap, we demonstrate that the latter procedure has the potential for robust automated estimation of variable importance measures based upon machine-learning algorithms. The cTMLE estimator was associated with substantially less positivity bias as compared to TMLE and larger coverage of the 95% CI. This study confirms the power of an automated cTMLE procedure that can target model selection via machine learning to estimate VIMs in complicated, high-dimensional data.

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74g

The increasing burden and the continued suboptimal outcomes for patients with heart failure underlines the importance of continued research to develop novel therapeutics for this disorder. This can only be accomplished with successful translation of basic science discoveries into direct human application through effective clinical trial design and execution that results in a substantially improved clinical course and outcomes. In this respect, phase II clinical trials play a pivotal role in determining which of the multitude of potential basic science discoveries should move to the large and expansive registration trials in humans. A critical examination of the phase II trials in heart failure reveals multiple shortcomings in their concept, design, execution, and interpretation. To further a dialogue on the challenges and potential for improvement and the role of phase II trials in patients with heart failure, the Food and Drug Administration facilitated a meeting on October 17, 2016, represented by clinicians, researchers, industry members, and regulators. This document summarizes the discussion from this meeting and provides key recommendations for future directions.

In randomized trials, adjustment for measured covariates during the analysis can reduce variance and increase power. To avoid misleading inference, the analysis plan must be pre-specified. However, it is often unclear a priori which baseline covariates (if any) should be adjusted for in the analysis. Consider, for example, the Sustainable East Africa Research in Community Health (SEARCH) trial for HIV prevention and treatment. There are 16 matched pairs of communities and many potential adjustment variables, including region, HIV prevalence, male circumcision coverage, and measures of community-level viral load. In this paper, we propose a rigorous procedure to data-adaptively select the adjustment set, which maximizes the efficiency of the analysis. Specifically, we use cross-validation to select from a pre-specified library the candidate targeted maximum likelihood estimator (TMLE) that minimizes the estimated variance. For further gains in precision, we also propose a collaborative procedure for estimating the known exposure mechanism. Our small sample simulations demonstrate the promise of the methodology to maximize study power, while maintaining nominal confidence interval coverage. We show how our procedure can be tailored to the scientific question (intervention effect for the study sample vs. for the target population) and study design (pair-matched or not). Copyright © 2016 John Wiley & Sons, Ltd.

In cluster randomized trials, the study units usually are not a simple random sample from some clearly defined target population. Instead, the target population tends to be hypothetical or ill-defined, and the selection of study units tends to be systematic, driven by logistical and practical considerations. As a result, the population average treatment effect (PATE) may be neither well defined nor easily interpretable. In contrast, the sample average treatment effect (SATE) is the mean difference in the counterfactual outcomes for the study units. The sample parameter is easily interpretable and arguably the most relevant when the study units are not sampled from some specific super-population of interest. Furthermore, in most settings, the sample parameter will be estimated more efficiently than the population parameter. To the best of our knowledge, this is the first paper to propose using targeted maximum likelihood estimation (TMLE) for estimation and inference of the sample effect in trials with and without pair-matching. We study the asymptotic and finite sample properties of the TMLE for the sample effect and provide a conservative variance estimator. Finite sample simulations illustrate the potential gains in precision and power from selecting the sample effect as the target of inference. This work is motivated by the Sustainable East Africa Research in Community Health (SEARCH) study, a pair-matched, community randomized trial to estimate the effect of population-based HIV testing and streamlined ART on the 5-year cumulative HIV incidence (NCT01864603). The proposed methodology will be used in the primary analysis for the SEARCH trial. Copyright © 2016 John Wiley & Sons, Ltd.

BackgroundSo far, more than 170 loci have been associated with circulating lipid levels through genome-wide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels.MethodsWe used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ∼60 000 individuals in the discovery stage and ∼90 000 samples in the replication stage.ResultsOur study resulted in the identification of five new associations with circulating lipid levels at four loci. All four loci are within genes that can be linked biologically to lipid metabolism. One of the variants, rs116843064, is a damaging missense variant within the ANGPTL4 gene.ConclusionsThis study illustrates that GWAS with high-scale imputation may still help us unravel the biological mechanism behind circulating lipid levels.

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.

In social and health sciences, many research questions involve understanding the causal effect of a longitudinal treatment on mortality (or time-to-event outcomes in general). Often, treatment status may change in response to past covariates that are risk factors for mortality, and in turn, treatment status may also affect such subsequent covariates. In these situations, Marginal Structural Models (MSMs), introduced by Robins (1997. Marginal structural models Proceedings of the American Statistical Association. Section on Bayesian Statistical Science, 1-10), are well-established and widely used tools to account for time-varying confounding. In particular, a MSM can be used to specify the intervention-specific counterfactual hazard function, i. e. the hazard for the outcome of a subject in an ideal experiment where he/she was assigned to follow a given intervention on their treatment variables. The parameters of this hazard MSM are traditionally estimated using the Inverse Probability Weighted estimation Robins (1999. Marginal structural models versus structural nested models as tools for causal inference. In: Statistical models in epidemiology: the environment and clinical trials. Springer-Verlag, 1999:95-134), Robins et al. (2000), (IPTW, van der Laan and Petersen (2007. Causal effect models for realistic individualized treatment and intention to treat rules. Int J Biostat 2007;3:Article 3), Robins et al. (2008. Estimaton and extrapolation of optimal treatment and testing strategies. Statistics in Medicine 2008;27(23):4678-721)). This estimator is easy to implement and admits Wald-type confidence intervals. However, its consistency hinges on the correct specification of the treatment allocation probabilities, and the estimates are generally sensitive to large treatment weights (especially in the presence of strong confounding), which are difficult to stabilize for dynamic treatment regimes. In this paper, we present a pooled targeted maximum likelihood estimator (TMLE, van der Laan and Rubin (2006. Targeted maximum likelihood learning. The International Journal of Biostatistics 2006;2:1-40)) for MSM for the hazard function under longitudinal dynamic treatment regimes. The proposed estimator is semiparametric efficient and doubly robust, offering bias reduction over the incumbent IPTW estimator when treatment probabilities may be misspecified. Moreover, the substitution principle rooted in the TMLE potentially mitigates the sensitivity to large treatment weights in IPTW. We compare the performance of the proposed estimator with the IPTW and a on-targeted substitution estimator in a simulation study.

ObjectiveConsistent estimation of causal effects with inverse probability weighting estimators is known to rely on consistent estimation of propensity scores. To alleviate the bias expected from incorrect model specification for these nuisance parameters in observational studies, data-adaptive estimation and in particular an ensemble learning approach known as Super Learning has been proposed as an alternative to the common practice of estimation based on arbitrary model specification. While the theoretical arguments against the use of the latter haphazard estimation strategy are evident, the extent to which data-adaptive estimation can improve inferences in practice is not. Some practitioners may view bias concerns over arbitrary parametric assumptions as academic considerations that are inconsequential in practice. They may also be wary of data-adaptive estimation of the propensity scores for fear of greatly increasing estimation variability due to extreme weight values. With this report, we aim to contribute to the understanding of the potential practical consequences of the choice of estimation strategy for the propensity scores in real-world comparative effectiveness research.MethodWe implement secondary analyses of Electronic Health Record data from a large cohort of type 2 diabetes patients to evaluate the effects of four adaptive treatment intensification strategies for glucose control (dynamic treatment regimens) on subsequent development or progression of urinary albumin excretion. Three Inverse Probability Weighting estimators are implemented using both model-based and data-adaptive estimation strategies for the propensity scores. Their practical performances for proper confounding and selection bias adjustment are compared and evaluated against results from previous randomized experiments.ConclusionResults suggest both potential reduction in bias and increase in efficiency at the cost of an increase in computing time when using Super Learning to implement Inverse Probability Weighting estimators to draw causal inferences.

Comparing the relative fit of competing models can be used to address many different scientific questions. In classical statistics one can, if appropriate, use likelihood ratio tests and information based criterion, whereas clinical medicine has tended to rely on comparisons of fit metrics like C-statistics. However, for many data adaptive modelling procedures such approaches are not suitable. In these cases, statisticians have used cross-validation, which can make inference challenging. In this paper we propose a general approach that focuses on the "conditional" risk difference (conditional on the model fits being fixed) for the improvement in prediction risk. Specifically, we derive a Wald-type test statistic and associated confidence intervals for cross-validated test sets utilizing the independent validation within cross-validation in conjunction with a test for multiple comparisons. We show that this test maintains proper Type I Error under the null fit, and can be used as a general test of relative fit for any semi-parametric model alternative. We apply the test to a candidate gene study to test for the association of a set of genes in a genetic pathway.

An individualized treatment rule (ITR) is a treatment rule which assigns treatments to individuals based on (a subset of) their measured covariates. An optimal ITR is the ITR which maximizes the population mean outcome. Previous works in this area have assumed that treatment is an unlimited resource so that the entire population can be treated if this strategy maximizes the population mean outcome. We consider optimal ITRs in settings where the treatment resource is limited so that there is a maximum proportion of the population which can be treated. We give a general closed-form expression for an optimal stochastic ITR in this resource-limited setting, and a closed-form expression for the optimal deterministic ITR under an additional assumption. We also present an estimator of the mean outcome under the optimal stochastic ITR in a large semiparametric model that at most places restrictions on the probability of treatment assignment given covariates. We give conditions under which our estimator is efficient among all regular and asymptotically linear estimators. All of our results are supported by simulations.

BackgroundWe investigated the incidence of ischemic heart disease (IHD) in relation to accumulated exposure to particulate matter (PM) in a cohort of aluminum workers. We adjusted for time varying confounding characteristic of the healthy worker survivor effect, using a recently introduced method for the estimation of causal target parameters.MethodsApplying longitudinal targeted minimum loss-based estimation, we estimated the difference in marginal cumulative risk of IHD in the cohort comparing counterfactual outcomes if always exposed above to always exposed below a PM2.5 exposure cut-off. Analyses were stratified by sub-cohort employed in either smelters or fabrication facilities. We selected two exposure cut-offs a priori, at the median and 10th percentile in each sub-cohort.ResultsIn smelters, the estimated IHD risk difference after 15 years of accumulating PM2.5 exposure during follow-up was 2.9% (0.6%, 5.1%) using the 10th percentile cut-off of 0.10 mg/m. For fabrication workers, the difference was 2.5% (0.8%, 4.1%) at the 10th percentile of 0.06 mg/m. Using the median exposure cut-off, results were similar in direction but smaller in size. We present marginal incidence curves describing the cumulative risk of IHD over the course of follow-up for each sub-cohort under each intervention regimen.ConclusionsThe accumulation of exposure to PM2.5 appears to result in higher risks of IHD in both aluminum smelter and fabrication workers. This represents the first longitudinal application of targeted minimum loss-based estimation, a method for generating doubly robust semi-parametric efficient substitution estimators of causal parameters, in the fields of occupational and environmental epidemiology.

Each year, 500,000 patients are treated with radiotherapy for head and neck cancer, resulting in relatively high survival rates. However, in 40% of patients, quality of life is severely compromised because of radiation-induced impairment of salivary gland function and consequent xerostomia (dry mouth). New radiation treatment technologies enable sparing of parts of the salivary glands. We have determined the parts of the major salivary gland, the parotid gland, that need to be spared to ensure that the gland continues to produce saliva after irradiation treatment. In mice, rats, and humans, we showed that stem and progenitor cells reside in the region of the parotid gland containing the major ducts. We demonstrated in rats that inclusion of the ducts in the radiation field led to loss of regenerative capacity, resulting in long-term gland dysfunction with reduced saliva production. Then we showed in a cohort of patients with head and neck cancer that the radiation dose to the region of the salivary gland containing the stem/progenitor cells predicted the function of the salivary glands one year after radiotherapy. Finally, we showed that this region of the salivary gland could be spared during radiotherapy, thus reducing the risk of post-radiotherapy xerostomia.