Emotional problems across development: examining measurement invariance across childhood, adolescence and early adulthood.

Emotional problems (anxiety, depression) are prevalent in children, adolescents and young adults with varying ages at onset. Studying developmental changes in emotional problems requires repeated assessments using the same or equivalent measures. The parent-rated Strengths and Difficulties Questionnaire is commonly used to assess emotional problems in childhood and adolescence, but there is limited research about whether it captures a similar construct across these developmental periods. Our study addressed this by investigating measurement invariance in the scales' emotional problems subscale (SDQ-EP) across childhood, adolescence and early adulthood. Data from two UK population cohorts were utilised: the Millennium Cohort Study (ages 3-17 years) and the Avon Longitudinal Study of Parents and Children (4-25 years). In both samples we observed weak (metric) measurement invariance by age, suggesting that the parent-rated SDQ-EP items contribute to the underlying construct of emotional problems similarly across age. This supports the validity of using the subscale to rank participants on their levels of emotional problems in childhood, adolescence and early adulthood. However strong (scalar) measurement invariance was not observed, suggesting that the same score may correspond to different levels of emotional problems across developmental periods. Comparisons of mean parent-rated SDQ-EP scores across age may therefore not be valid.

Counterfactual Mediation Analysis with a Latent Class Exposure.

Latent classes are a useful tool in developmental research, however there are challenges associated with embedding them within a counterfactual mediation model. We develop and test a new method "updated pseudo class draws (uPCD)" to examine the association between a latent class exposure and distal outcome that could easily be extended to allow the use of any counterfactual mediation method. UPCD extends an existing group of methods (based on pseudo class draws) that assume that the true values of the latent class variable are missing, and need to be multiply imputed using class membership probabilities. We simulate data based on the Avon Longitudinal Study of Parents and Children, examine performance for existing techniques to relate a latent class exposure to a distal outcome ("one-step," "bias-adjusted three-step," "modal class assignment," "non-inclusive pseudo class draws," and "inclusive pseudo class draws") and compare bias in parameter estimates and their precision to uPCD when estimating counterfactual mediation effects. We found that uPCD shows minimal bias when estimating counterfactual mediation effects across all levels of entropy. UPCD performs similarly to recommended methods (one-step and bias-adjusted three-step), but provides greater flexibility and scope for incorporating the latent grouping within any commonly-used counterfactual mediation approach.

Depressive symptoms in adolescence and adult educational and employment outcomes: a structured life course analysis.

BACKGROUND: Depression is a common mental health disorder that often starts during adolescence, with potentially important future consequences including 'Not in Education, Employment or Training' (NEET) status. METHODS: We took a structured life course modeling approach to examine how depressive symptoms during adolescence might be associated with later NEET status, using a high-quality longitudinal data resource. We considered four plausible life course models: (1) an early adolescent sensitive period model where depressive symptoms in early adolescence are more associated with later NEET status relative to exposure at other stages; (2) a mid adolescent sensitive period model where depressive symptoms during the transition from compulsory education to adult life might be more deleterious regarding NEET status; (3) a late adolescent sensitive period model, meaning that depressive symptoms around the time when most adults have completed their education and started their careers are the most strongly associated with NEET status; and (4) an accumulation of risk model which highlights the importance of chronicity of symptoms. RESULTS: Our analysis sample included participants with full information on NEET status (N = 3951), and the results supported the accumulation of risk model, showing that the odds of NEET increase by 1.015 (95% CI 1.012-1.019) for an increase of 1 unit in depression at any age between 11 and 24 years. CONCLUSIONS: Given the adverse implications of NEET status, our results emphasize the importance of supporting mental health during adolescence and early adulthood, as well as considering specific needs of young people with re-occurring depressed mood.

Non-linear Mendelian randomization: detection of biases using negative controls with a focus on BMI, Vitamin D and LDL cholesterol.

Mendelian randomisation (MR) is an established technique in epidemiological investigation, using the principle of random allocation of genetic variants at conception to estimate the causal linear effect of an exposure on an outcome. Extensions to this technique include non-linear approaches that allow for differential effects of the exposure on the outcome depending on the level of the exposure. A widely used non-linear method is the residual approach, which estimates the causal effect within different strata of the non-genetically predicted exposure (i.e. the "residual" exposure). These "local" causal estimates are then used to make inferences about non-linear effects. Recent work has identified that this method can lead to estimates that are seriously biased, and a new method-the doubly-ranked method-has been introduced as a possibly more robust approach. In this paper, we perform negative control outcome analyses in the MR context. These are analyses with outcomes onto which the exposure should have no predicted causal effect. Using both methods we find clearly biased estimates in certain situations. We additionally examined a situation for which there are robust randomised controlled trial estimates of effects-that of low-density lipoprotein cholesterol (LDL-C) reduction onto myocardial infarction, where randomised trials have provided strong evidence of the shape of the relationship. The doubly-ranked method did not identify the same shape as the trial data, and for LDL-C and other lipids they generated some highly implausible findings. Therefore, we suggest there should be extensive simulation and empirical methodological examination of performance of both methods for NLMR under different conditions before further use of these methods. In the interim, use of NLMR methods needs justification, and a number of sanity checks (such as analysis of negative and positive control outcomes, sensitivity analyses excluding removal of strata at the extremes of the distribution, examination of biological plausibility and triangulation of results) should be performed.

Regression discontinuity design for the study of health effects of exposures acting early in life.

Regression discontinuity design (RDD) is a quasi-experimental approach to study the causal effect of an exposure on later outcomes by exploiting the discontinuity in the exposure probability at an assignment variable cut-off. With the intent of facilitating the use of RDD in the Developmental Origins of Health and Disease (DOHaD) research, we describe the main aspects of the study design and review the studies, assignment variables and exposures that have been investigated to identify short- and long-term health effects of early life exposures. We also provide a brief overview of some of the methodological considerations for the RDD identification using an example of a DOHaD study. An increasing number of studies investigating the effects of early life environmental stressors on health outcomes use RDD, mostly in the context of education, social and welfare policies, healthcare organization and insurance, and clinical management. Age and calendar time are the mostly used assignment variables to study the effects of various early life policies and programs, shock events and guidelines. Maternal and newborn characteristics, such as age, birth weight and gestational age are frequently used assignment variables to study the effects of the type of neonatal care, health insurance, and newborn benefits, while socioeconomic measures have been used to study the effects of social and welfare programs. RDD has advantages, including intuitive interpretation, and transparent and simple graphical representation. It provides valid causal estimates if the assumptions, relatively weak compared to other non-experimental study designs, are met. Its use to study health effects of exposures acting early in life has been limited to studies based on registries and administrative databases, while birth cohort data has not been exploited so far using this design. Local causal effect around the cut-off, difficulty in reaching high statistical power compared to other study designs, and the rarity of settings outside of policy and program evaluations hamper the widespread use of RDD in the DOHaD research. Still, the assignment variables' cut-offs for exposures applied in previous studies can be used, if appropriate, in other settings and with additional outcomes to address different research questions.

A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts.

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

A tool to assess risk of bias in non-randomized follow-up studies of exposure effects (ROBINS-E).

BACKGROUND: Observational epidemiologic studies provide critical data for the evaluation of the potential effects of environmental, occupational and behavioural exposures on human health. Systematic reviews of these studies play a key role in informing policy and practice. Systematic reviews should incorporate assessments of the risk of bias in results of the included studies. OBJECTIVE: To develop a new tool, Risk Of Bias In Non-randomized Studies - of Exposures (ROBINS-E) to assess risk of bias in estimates from cohort studies of the causal effect of an exposure on an outcome. METHODS AND RESULTS: ROBINS-E was developed by a large group of researchers from diverse research and public health disciplines through a series of working groups, in-person meetings and pilot testing phases. The tool aims to assess the risk of bias in a specific result (exposure effect estimate) from an individual observational study that examines the effect of an exposure on an outcome. A series of preliminary considerations informs the core ROBINS-E assessment, including details of the result being assessed and the causal effect being estimated. The assessment addresses bias within seven domains, through a series of 'signalling questions'. Domain-level judgements about risk of bias are derived from the answers to these questions, then combined to produce an overall risk of bias judgement for the result, together with judgements about the direction of bias. CONCLUSION: ROBINS-E provides a standardized framework for examining potential biases in results from cohort studies. Future work will produce variants of the tool for other epidemiologic study designs (e.g. case-control studies). We believe that ROBINS-E represents an important development in the integration of exposure assessment, evidence synthesis and causal inference.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

Background: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. Methods: We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR). Results: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. Conclusions: We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Using the global randomization test as a Mendelian randomization falsification test for the exclusion restriction assumption.

Mendelian randomization may give biased causal estimates if the instrument affects the outcome not solely via the exposure of interest (violating the exclusion restriction assumption). We demonstrate use of a global randomization test as a falsification test for the exclusion restriction assumption. Using simulations, we explored the statistical power of the randomization test to detect an association between a genetic instrument and a covariate set due to (a) selection bias or (b) horizontal pleiotropy, compared to three approaches examining associations with individual covariates: (i) Bonferroni correction for the number of covariates, (ii) correction for the effective number of independent covariates, and (iii) an r2 permutation-based approach. We conducted proof-of-principle analyses in UK Biobank, using CRP as the exposure and coronary heart disease (CHD) as the outcome. In simulations, power of the randomization test was higher than the other approaches for detecting selection bias when the correlation between the covariates was low (r2 < 0.1), and at least as powerful as the other approaches across all simulated horizontal pleiotropy scenarios. In our applied example, we found strong evidence of selection bias using all approaches (e.g., global randomization test p < 0.002). We identified 51 of the 58 CRP genetic variants as horizontally pleiotropic, and estimated effects of CRP on CHD attenuated somewhat to the null when excluding these from the genetic risk score (OR = 0.96 [95% CI: 0.92, 1.00] versus 0.97 [95% CI: 0.90, 1.05] per 1-unit higher log CRP levels). The global randomization test can be a useful addition to the MR researcher's toolkit.

Meta-regression of genome-wide association studies to estimate age-varying genetic effects.

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Investigation of genetic determinants of cognitive change in later life.

Cognitive decline is a major health concern and identification of genes that may serve as drug targets to slow decline is important to adequately support an aging population. Whilst genetic studies of cross-sectional cognition have been carried out, cognitive change is less well-understood. Here, using data from the TOMMORROW trial, we investigate genetic associations with cognitive change in a cognitively normal older cohort. We conducted a genome-wide association study of trajectories of repeated cognitive measures (using generalised estimating equation (GEE) modelling) and tested associations with polygenic risk scores (PRS) of potential risk factors. We identified two genetic variants associated with change in attention domain scores, rs534221751 (p = 1 × 10-8 with slope 1) and rs34743896 (p = 5 × 10-10 with slope 2), implicating NCAM2 and CRIPT/ATP6V1E2 genes, respectively. We also found evidence for the association between an education PRS and baseline cognition (at >65 years of age), particularly in the language domain. We demonstrate the feasibility of conducting GWAS of cognitive change using GEE modelling and our results suggest that there may be novel genetic associations for cognitive change that have not previously been associated with cross-sectional cognition. We also show the importance of the education PRS on cognition much later in life. These findings warrant further investigation and demonstrate the potential value of using trial data and trajectory modelling to identify genetic variants associated with cognitive change.

Two sample Mendelian Randomisation using an outcome from a multilevel model of disease progression.

Identifying factors that are causes of disease progression, especially in neurodegenerative diseases, is of considerable interest. Disease progression can be described as a trajectory of outcome over time-for example, a linear trajectory having both an intercept (severity at time zero) and a slope (rate of change). A technique for identifying causal relationships between one exposure and one outcome in observational data whilst avoiding bias due to confounding is two sample Mendelian Randomisation (2SMR). We consider a multivariate approach to 2SMR using a multilevel model for disease progression to estimate the causal effect an exposure has on the intercept and slope. We carry out a simulation study comparing a naïve univariate 2SMR approach to a multivariate 2SMR approach with one exposure that effects both the intercept and slope of an outcome that changes linearly with time since diagnosis. The simulation study results, across six different scenarios, for both approaches were similar with no evidence against a non-zero bias and appropriate coverage of the 95% confidence intervals (for intercept 93.4-96.2% and the slope 94.5-96.0%). The multivariate approach gives a better joint coverage of both the intercept and slope effects. We also apply our method to two Parkinson's cohorts to examine the effect body mass index has on disease progression. There was no strong evidence that BMI affects disease progression, however the confidence intervals for both intercept and slope were wide.

Multiple imputation strategies for missing event times in a multi-state model analysis.

In clinical studies, multi-state model (MSM) analysis is often used to describe the sequence of events that patients experience, enabling better understanding of disease progression. A complicating factor in many MSM studies is that the exact event times may not be known. Motivated by a real dataset of patients who received stem cell transplants, we considered the setting in which some event times were exactly observed and some were missing. In our setting, there was little information about the time intervals in which the missing event times occurred and missingness depended on the event type, given the analysis model covariates. These additional challenges limited the usefulness of some missing data methods (maximum likelihood, complete case analysis, and inverse probability weighting). We show that multiple imputation (MI) of event times can perform well in this setting. MI is a flexible method that can be used with any complete data analysis model. Through an extensive simulation study, we show that MI by predictive mean matching (PMM), in which sampling is from a set of observed times without reliance on a specific parametric distribution, has little bias when event times are missing at random, conditional on the observed data. Applying PMM separately for each sub-group of patients with a different pathway through the MSM tends to further reduce bias and improve precision. We recommend MI using PMM methods when performing MSM analysis with Markov models and partially observed event times.

Modeling Risk Factors for Intraindividual Variability: A Mixed-Effects Beta-Binomial Model Applied to Cognitive Function in Older People in the English Longitudinal Study of Ageing.

Cognitive functioning in older age profoundly impacts quality of life and health. While most research on cognition in older age has focused on mean levels, intraindividual variability (IIV) around this may have risk factors and outcomes independent of the mean value. Investigating risk factors associated with IIV has typically involved deriving a summary statistic for each person from residual error around a fitted mean. However, this ignores uncertainty in the estimates, prohibits exploring associations with time-varying factors, and is biased by floor/ceiling effects. To address this, we propose a mixed-effects location scale beta-binomial model for estimating average probability and IIV in a word recall test in the English Longitudinal Study of Ageing. After adjusting for mean performance, an analysis of 9,873 individuals across 7 (mean = 3.4) waves (2002-2015) found IIV to be greater at older ages, with lower education, in females, with more difficulties in activities of daily living, in later birth cohorts, and when interviewers recorded issues potentially affecting test performance. Our study introduces a novel method for identifying groups with greater IIV in bounded discrete outcomes. Our findings have implications for daily functioning and care, and further work is needed to identify the impact for future health outcomes.

Young Adult ADHD Symptoms in the General Population and Neurocognitive Impairment.

OBJECTIVE: Neurocognitive impairments are associated with child and adult ADHD in clinical settings. However, it is unknown whether adult ADHD symptoms in the general population are associated with the same pattern of cognitive impairment. We examined this using a prospective, population-based cohort spanning birth to age 25 years. METHODS: We examined associations between self-reported adult ADHD symptoms and cognitive task performance (attention and response inhibition) in adulthood and childhood. RESULTS: Self-rated ADHD symptoms at age 25 were associated with poorer performance in age 25 cognitive tasks capturing ADHD-related functioning (attention B = -0.03, 95% CI [0.05, -0.01], p = .005; response inhibition B = -0.03, 95% CI [-0.05, -0.01], p = .002). CONCLUSIONS: Neurocognitive impairments linked to adult ADHD symptoms in the general population, are similar to those found in people with childhood ADHD symptoms and are consistent with findings in adult ADHD clinical samples.

An epigenome-wide analysis of DNA methylation, racialized and economic inequities, and air pollution.

Importance: DNA methylation (DNAm) provides a plausible mechanism by which adverse exposures become embodied and contribute to health inequities, due to its role in genome regulation and responsiveness to social and biophysical exposures tied to societal context. However, scant epigenome-wide association studies (EWAS) have included structural and lifecourse measures of exposure, especially in relation to structural discrimination. Objective: Our study tests the hypothesis that DNAm is a mechanism by which racial discrimination, economic adversity, and air pollution become biologically embodied. Design: A series of cross-sectional EWAS, conducted in My Body My Story (MBMS, biological specimens collected 2008-2010, DNAm assayed in 2021); and the Multi Ethnic Study of Atherosclerosis (MESA; biological specimens collected 2010-2012, DNAm assayed in 2012-2013); using new georeferenced social exposure data for both studies (generated in 2022). Setting: MBMS was recruited from four community health centers in Boston; MESA was recruited from four field sites in: Baltimore, MD; Forsyth County, NC; New York City, NY; and St. Paul, MN. Participants: Two population-based samples of US-born Black non-Hispanic (Black NH), white non-Hispanic (white NH), and Hispanic individuals (MBMS; n=224 Black NH and 69 white NH) and (MESA; n=229 Black NH, n=555 white NH and n=191 Hispanic). Exposures: Eight social exposures encompassing racial discrimination, economic adversity, and air pollution. Main outcome: Genome-wide changes in DNAm, as measured using the Illumina EPIC BeadChip (MBMS; using frozen blood spots) and Illumina 450k BeadChip (MESA; using purified monocytes). Our hypothesis was formulated after data collection. Results: We observed the strongest associations with traffic-related air pollution (measured via black carbon and nitrogen oxides exposure), with evidence from both studies suggesting that air pollution exposure may induce epigenetic changes related to inflammatory processes. We also found suggestive associations of DNAm variation with measures of structural racial discrimination (e.g., for Black NH participants, born in a Jim Crow state; adult exposure to racialized economic residential segregation) situated in genes with plausible links to effects on health. Conclusions and Relevance: Overall, this work suggests that DNAm is a biological mechanism through which structural racism and air pollution become embodied and may lead to health inequities.

Childhood socioeconomic position and sex-specific trajectories of metabolic traits across early life: prospective cohort study.

BACKGROUND: Socioeconomic inequalities in cardiovascular disease risk begin early in life and are more pronounced in females than males later in life. Causal atherogenic traits explaining this are not well understood. We explored sex-specific associations between childhood socioeconomic position (SEP) and molecular measures of systemic metabolism across early life. METHODS: Data were from the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based birth cohort in southwest England. Pregnant women with an expected delivery date between 1991 and 1992 were invited to participate. Maternal education was the primary indicator of SEP. Concentrations of 148 metabolic traits from targeted metabolomics (nuclear magnetic resonance spectroscopy) from research clinics at ages 7, 15, 18 and 25 years were analysed. The sex-specific slope index of inequality (SII) in trajectories of metabolic traits was estimated using multilevel models. FINDINGS: Total number of participants included was 6537 (12,543 repeated measures). Lower maternal education was associated with more adverse levels of several atherogenic lipids and key metabolic traits among females at age 7 years, but not males. For instance, SII for very small very-low-density lipoprotein (VLDL) concentrations was 0.16SD (95% CI: 0.01, 0.30) among females and -0.02SD (95% CI: -0.16, 0.13) among males. Between 7 and 25 years, inequalities widened among females and emerged among males particularly for VLDL particle concentrations, apolipoprotein-B concentrations, and inflammatory glycoprotein acetyls. For instance, at 25 years, SII for very small VLDL concentrations was 0.36SD (95% CI: 0.20, 0.52) and 0.22SD (95% CI: 0.04, 0.40) among females and males respectively. INTERPRETATION: Prevention of socioeconomic inequalities in cardiovascular disease risk requires a life course approach beginning at the earliest opportunity, especially among females. FUNDING: The UK Medical Research Council and Wellcome (grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). KON is supported by a Health Research Board (HRB) of Ireland Investigator Led Award (ILP-PHR-2022-008). JB, GDS and KT work in a unit funded by the UK MRC (MC_UU_00011/1 and MC UU 00011/3) and the University of Bristol. OR is supported by a UKRI Future Leaders Fellowship (MR/S03532X/1). These funding sources had no role in the design and conduct of this study. This publication is the work of the authors and KON will serve as guarantor for the contents of this paper.

Methodological approaches, challenges, and opportunities in the application of Mendelian randomisation to lifecourse epidemiology: A systematic literature review.

Diseases diagnosed in adulthood may have antecedents throughout (including prenatal) life. Gaining a better understanding of how exposures at different stages in the lifecourse influence health outcomes is key to elucidating the potential benefits of disease prevention strategies. Mendelian randomisation (MR) is increasingly used to estimate causal effects of exposures across the lifecourse on later life outcomes. This systematic literature review explores MR methods used to perform lifecourse investigations and reviews previous work that has utilised MR to elucidate the effects of factors acting at different stages of the lifecourse. We conducted searches in PubMed, Embase, Medline and MedRXiv databases. Thirteen methodological studies were identified. Four studies focused on the impact of time-varying exposures in the interpretation of "standard" MR techniques, five presented methods for repeat measures of the same exposure, and four described methodological approaches to handling multigenerational exposures. A further 127 studies presented the results of an applied research question. Over half of these estimated effects in a single generation and were largely confined to the exploration of questions regarding body composition. The remaining mostly estimated maternal effects. There is a growing body of research focused on the development and application of MR methods to address lifecourse research questions. The underlying assumptions require careful consideration and the interpretation of results rely on select conditions. Whilst we do not advocate for a particular strategy, we encourage practitioners to make informed decisions on how to approach a research question in this field with a solid understanding of the limitations present and how these may be affected by the research question, modelling approach, instrument selection, and data availability.

Multiple imputation of missing data under missing at random: including a collider as an auxiliary variable in the imputation model can induce bias.

Epidemiological studies often have missing data, which are commonly handled by multiple imputation (MI). In MI, in addition to those required for the substantive analysis, imputation models often include other variables ("auxiliary variables"). Auxiliary variables that predict the partially observed variables can reduce the standard error (SE) of the MI estimator and, if they also predict the probability that data are missing, reduce bias due to data being missing not at random. However, guidance for choosing auxiliary variables is lacking. We examine the consequences of a poorly chosen auxiliary variable: if it shares a common cause with the partially observed variable and the probability that it is missing (i.e., it is a "collider"), its inclusion can induce bias in the MI estimator and may increase the SE. We quantify, both algebraically and by simulation, the magnitude of bias and SE when either the exposure or outcome is incomplete. When the substantive analysis outcome is partially observed, the bias can be substantial, relative to the magnitude of the exposure coefficient. In settings in which a complete records analysis is valid, the bias is smaller when the exposure is partially observed. However, bias can be larger if the outcome also causes missingness in the exposure. When using MI, it is important to examine, through a combination of data exploration and considering plausible casual diagrams and missingness mechanisms, whether potential auxiliary variables are colliders.

Association between perinatal mortality and morbidity and customised and non-customised birthweight centiles in Denmark, Finland, Norway, Wales, and England: comparative, population based, record linkage study.

Objectives: To compare the risk of adverse perinatal outcomes according to infants who are born small for gestational age (SGA; <10th centile) or large for gestational age (LGA; >90th centile), as defined by birthweight centiles that are non-customised (ie, standardised by sex and gestational age only) and customised (by sex, gestational age, maternal weight, height, parity, and ethnic group). Design: Comparative, population based, record linkage study with meta-analysis of results. Setting: Denmark, Finland, Norway, Wales, and England (city of Bradford), 1986-2019. Participants: 2 129 782 infants born at term in birth registries. Main outcome measures: Stillbirth, neonatal death, infant death, admission to neonatal intensive care unit, and low Apgar score (<7) at 5 minutes. Results: Relative to those infants born average for gestational age (AGA), both SGA and LGA births were at increased risk of all five outcomes, but observed relative risks were similar irrespective of whether non-customised or customised charts were used. For example, for SGA versus AGA births, when non-customised and customised charts were used, relative risks pooled over countries were 3.60 (95% confidence interval 3.29 to 3.93) versus 3.58 (3.02 to 4.24) for stillbirth, 2.83 (2.18 to 3.67) versus 3.32 (2.05 to 5.36) for neonatal death, 2.82 (2.07 to 3.83) versus 3.17 (2.20 to 4.56) for infant death, 1.66 (1.49 to 1.86) versus 1.54 (1.30 to 1.81) for low Apgar score at 5 minutes, and (based on Bradford data only) 1.97 (1.74 to 2.22) versus 1.94 (1.70 to 2.21) for admission to the neonatal intensive care unit. The estimated sensitivity of combined SGA or LGA births to identify the three mortality outcomes ranged from 31% to 34% for non-customised charts and from 34% to 38% for customised charts, with a specificity of 82% and 80% with non-customised and customised charts, respectively. Conclusions: These results suggest an increased risk of adverse perinatal outcomes of a similar magnitude among SGA or LGA term infants when customised and non-customised centiles are used. Use of customised charts for SGA/LGA births-over and above use of non-customised charts for SGA/LGA births-is unlikely to provide benefits in terms of identifying term births at risk of these outcomes.