https://orcid.org/0000-0001-8778-4132
Figures
Abstract
Background
Insomnia is common and associated with adverse pregnancy and perinatal outcomes in observational studies. However, those associations could be vulnerable to residual confounding or reverse causality. Our aim was to estimate the association of insomnia with stillbirth, miscarriage, gestational diabetes (GD), hypertensive disorders of pregnancy (HDP), perinatal depression, preterm birth (PTB), and low/high offspring birthweight (LBW/HBW).
Methods and findings
We used 2-sample mendelian randomization (MR) with 81 single-nucleotide polymorphisms (SNPs) instrumenting for a lifelong predisposition to insomnia. Our outcomes included ever experiencing stillbirth, ever experiencing miscarriage, GD, HDP, perinatal depression, PTB (gestational age <37 completed weeks), LBW (<2,500 grams), and HBW (>4,500 grams). We used data from women of European descent (N = 356,069, mean ages at delivery 25.5 to 30.0 years) from UK Biobank (UKB), FinnGen, Avon Longitudinal Study of Parents and Children (ALSPAC), Born in Bradford (BiB), and the Norwegian Mother, Father and Child Cohort (MoBa). Main MR analyses used inverse variance weighting (IVW), with weighted median and MR-Egger as sensitivity analyses. We compared MR estimates with multivariable regression of insomnia in pregnancy on outcomes in ALSPAC (N = 11,745). IVW showed evidence of an association of genetic susceptibility to insomnia with miscarriage (odds ratio (OR): 1.60, 95% confidence interval (CI): 1.18, 2.17, p = 0.002), perinatal depression (OR 3.56, 95% CI: 1.49, 8.54, p = 0.004), and LBW (OR 3.17, 95% CI: 1.69, 5.96, p < 0.001). IVW results did not support associations of insomnia with stillbirth, GD, HDP, PTB, and HBW, with wide CIs including the null. Associations of genetic susceptibility to insomnia with miscarriage, perinatal depression, and LBW were not observed in weighted median or MR-Egger analyses. Results from these sensitivity analyses were directionally consistent with IVW results for all outcomes, with the exception of GD, perinatal depression, and PTB in MR-Egger. Multivariable regression showed associations of insomnia at 18 weeks of gestation with perinatal depression (OR 2.96, 95% CI: 2.42, 3.63, p < 0.001), but not with LBW (OR 0.92, 95% CI: 0.69, 1.24, p = 0.60). Multivariable regression with miscarriage and stillbirth was not possible due to small numbers in index pregnancies. Key limitations are potential horizontal pleiotropy (particularly for perinatal depression) and low statistical power in MR, and residual confounding in multivariable regression.
Conclusions
In this study, we observed some evidence in support of a possible causal relationship between genetically predicted insomnia and miscarriage, perinatal depression, and LBW. Our study also found observational evidence in support of an association between insomnia in pregnancy and perinatal depression, with no clear multivariable evidence of an association with LBW. Our findings highlight the importance of healthy sleep in women of reproductive age, though replication in larger studies, including with genetic instruments specific to insomnia in pregnancy are important.
Author summary
Why was this study done?
- Insomnia in pregnancy was associated with higher risks of adverse pregnancy and perinatal outcomes in observational studies.
- It is currently not clear whether insomnia causes adverse pregnancy and perinatal outcomes or whether the unfavourable associations are explained by confounding.
- To the best of our knowledge, mendelian randomization (MR) has not been used to explore whether there is evidence to support a causal association between insomnia and adverse pregnancy and perinatal outcomes.
What did the researchers do and find?
- We used data on up to 356,069 women from UK Biobank (UKB), FinnGen, and 3 birth cohorts and assessed whether genetic susceptibility to insomnia was associated with stillbirth, miscarriage, gestational diabetes (GD), hypertensive disorders of pregnancy (HDP), perinatal depression, preterm birth (PTB), low offspring birthweight (LBW), and high offspring birthweight (HBW) in 2-sample MR.
- To triangulate with our MR estimates, we conducted multivariable regression in 11,745 women from the Avon Longitudinal Study of Parents and Children (ALSPAC), where insomnia was measured in pregnancy for all outcomes except miscarriage and stillbirth for which there were too few cases in the index pregnancy.
- We found evidence from MR and multivariable regression that insomnia was associated with a higher risk of perinatal depression, and MR analyses also suggested evidence for an association between genetically predicted insomnia and risks of miscarriage and LBW.
What do these findings mean?
- These findings raise the possibility that insomnia maybe related to adverse pregnancy outcomes, implying that interventions to improve healthy sleep may be beneficial to a healthy pregnancy.
- Key limitations of our study are potential horizontal pleiotropy (particularly for perinatal depression) and low statistical power in MR and residual confounding in multivariable regression. Replication in larger MR studies would be valuable.
Citation: Yang Q, Borges MC, Sanderson E, Magnus MC, Kilpi F, Collings PJ, et al. (2022) Associations between insomnia and pregnancy and perinatal outcomes: Evidence from mendelian randomization and multivariable regression analyses. PLoS Med 19(9): e1004090. https://doi.org/10.1371/journal.pmed.1004090
Academic Editor: Sarah J. Stock, The University of Edinburgh Usher Institute of Population Health Sciences and Informatics, UNITED KINGDOM
Received: October 7, 2021; Accepted: August 15, 2022; Published: September 6, 2022
Copyright: © 2022 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: We used both individual participant cohort data and publicly available summary statistics. We present summary statistics that we generated from those individual participant cohort data in S4 and S5 Tables. Full information on how to access UKB data can be found at its website (https://www.ukbiobank.ac.uk/researchers/). All ALSPAC data are available to scientists on request to the ALSPAC Executive via this website (http://www.bristol.ac.uk/alspac/researchers/), which also provides full details and distributions of the ALSPAC study variables. Similarly, data from BiB are available on request to the BiB Executive (https://borninbradford.nhs.uk/research/how-to-access-data/). Data from MoBa are available from the Norwegian Institute of Public Health after application to the MoBa Scientific Management Group (see its website https://www.fhi.no/en/op/data-access-from-health-registries-health-studies-and-biobanks/data-access/applying-for-access-to-data/ for details). Summary statistics from FinnGen are publicly available on its website (https://finngen.gitbook.io/documentation/data-download).
Funding: QY, MCB, ES, FK, ALS, KT and DAL work in a unit that is supported by the University of Bristol and UK Medical Research Council (MRC, MM_UU_00011/1, MM_UU_00011/3 to KT and MM_UU_00011/6 to DAL). This work was supported by China Scholarship Council PhD Scholarship (CSC201807060273 to QY), UK MRC Skilled Development Fellowship (MR/P014054/1 to MCB), Research Council of Norway through its Centres of Excellence funding scheme (262700 to MCM and SHE), European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (947684 to MCM), and British Heart Foundation (BHF) Immediate Postdoctoral Basic Science Research Fellowship (FS/17/37/32937 to PJC). DAL is a BHF Chair (CH/F/20/9003) and National Institute of Health Research (NIHR) Senior Investigator (NF-0616-10102). Core funding for ALSPAC is provided by UK MRC and University of Bristol (217065/Z/19/Z). Genotyping of maternal samples was funded by the Wellcome Trust (WT088806), and offspring samples were genotyped by Sample Logistics and Genotyping Facilities at the Wellcome Trust Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). BiB is supported by a Wellcome programme grant (WT223601/Z/21/Z: Age of Wonder), an infrastructure grant (WT101597MA), a UK MRC and Economic and Social Science Research Council programme grant (MR/N024397/1), a BHF Clinical Study grant (CS/16/4/32482), and NIHR under its Applied Health Research Collaboration Yorkshire and Humber (NIHR200166) and the NIHR Clinical Research Network. Further supports for genome-wide and multiple omics measurements in BiB are from UK MRC (G0600705), NIHR (NF-SI-0611010196), US National Institute of Health (R01DK10324), and ERC via Advanced Grant (669545) and under the European Union’s Seventh Framework Programme (FP7/2007-2013). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: KT has acted as a consultant for CHDI Foundation, and Expert Witness to the High Court in England, called by the UK Medicines and Healthcare products Regulatory Agency, defendants in a case on hormonal pregnancy tests and congenital anomalies 2021/22. DAL has received support from Medtronic LTD and Roche Diagnostics for biomarker research that is not related to the study presented in this paper. The other authors report no conflicts.
Abbreviations: ALSPAC, Avon Longitudinal Study of Parents and Children; BiB, Born in Bradford; CI, confidence interval; GD, gestational diabetes; GWAS, genome-wide association study; HBW, high offspring birthweight; HDP, hypertensive disorders of pregnancy; IV, instrumental variable; IVW, inverse variance weighting; LBW, low offspring birthweight; MoBa, the Norwegian Mother, Father and Child Cohort; MR, mendelian randomization; OR, odds ratio; PC, principal component; PTB, preterm birth; SNP, single-nucleotide polymorphism; STROBE, Strengthening the Reporting of Observational Studies in Epidemiology; UKB, UK Biobank
Introduction
Insomnia, which affects approximately 10% to 20% of the adult population, is usually defined as a difficulty in getting to sleep or remaining asleep, or having a nonrestorative sleep, and such sleep impairment can be associated with daytime sleepiness [1,2]. Physical and hormonal changes during pregnancy increase susceptibility to insomnia [3,4].
Most evidence on the relationship between insomnia during pregnancy and adverse pregnancy and perinatal outcomes has come from observational studies. The most recently updated systematic reviews of observational studies suggest that pregnancy-related insomnia and poor sleep quality are associated with higher risks of gestational diabetes (GD) [5,6], hypertensive disorders of pregnancy (HDP) [6], perinatal depression [7], and preterm birth (PTB) [6]. Other observational studies have shown that specific conditions that relate to insomnia are also associated with adverse pregnancy and perinatal outcomes. Sleep-disordered breathing, obstructive sleep apnoea, and restless legs syndrome have also been shown to associate with higher risks of GD, HDP, large-for-gestational age, and low offspring birthweight (LBW) [6]. However, it remains unclear whether insomnia causes adverse pregnancy outcomes or whether these associations are explained by confounding, e.g., due to socioeconomic status and lifestyle factors. It is also possible that some of these studies reflect reverse causation. For example, all 4 studies included in the systematic review for perinatal depression were cross-sectional [7], in which disturbed sleep could be either a symptom of or a risk factor for depression. Furthermore, most individual studies focus on just 1 or 2 outcomes. Examining potential effects on a range of adverse pregnancy and perinatal outcomes is important to understand the overall health impact of insomnia during pregnancy.
Three randomized control trials assessing the effects of interventions to prevent insomnia on adverse pregnancy and perinatal outcomes have been published [8–10]. All 3 of these used cognitive behavioural interventions targeted at reducing insomnia, with the primary outcome being Edinburgh Postnatal Depression Scale scores. The small number of randomized control trials, their small sample sizes, and directional inconsistency, but overlapping 95% confidence intervals (CIs), make it difficult to draw conclusions, and none of them explored other adverse pregnancy or perinatal outcomes.
Mendelian randomization (MR) provides an alternative way to assess the impact of insomnia on adverse pregnancy and perinatal outcomes by using genetic variants (mostly single-nucleotide polymorphisms [SNPs]) as instrumental variables (IVs) for insomnia [11,12]. MR is less prone to confounding than observational studies, as genetic variants are randomly allocated at meiosis and cannot be influenced by the wide range of sociodemographic or behavioural factors which conventionally confound observational studies nor can they be influenced by health status [11,12]. Under key assumptions (discussed in Methods), MR can be used to estimate a causal association from the SNPs-exposure and SNPs-outcome associations, if the underlying assumptions (in Discussion) are true. In 2-sample MR, the SNP-exposure and SNP-outcome associations are estimated using different (ideally independent) studies from the same underlying population [13]. This approach has previously been used to evaluate causal associations of insomnia with type 2 diabetes [14,15], hypertension [16], and cardiovascular disease [15,17,18] in non-pregnant populations, but to the best of our knowledge not pregnancy and perinatal outcomes.
The aims of this study are to (I) explore the causal associations of maternal genetic susceptibility to insomnia with stillbirth, miscarriage, GD, HDP, perinatal depression, PTB, LBW, and high offspring birthweight (HBW), using 2-sample MR; and (II) compare MR findings with conventional multivariable regression analyses of self-reported insomnia during pregnancy with these outcomes, where possible.
Methods
Study populations
This study was undertaken using data from the MR-PREG collaboration, which aims to explore causes and consequences of different pregnancy and perinatal outcomes [19]. We used individual-level data from UK Biobank (UKB) women (N = 208,140, recruited between 2006 to 2010) and mother-offspring pairs from Avon Longitudinal Study of Parents and Children (ALSPAC, N = 6,826, recruited between 1991 to 1992), Born in Bradford (BiB, N = 2,940, recruited between 2007 to 2010), and the Norwegian Mother, Father and Child Cohort (MoBa, N = 14,584, recruited between 1999 to 2009). To be comparable across all cohorts, only genetically unrelated women of European descent with qualified genotype data (and with singleton offspring in birth cohorts) were eligible for inclusion in our analyses (S1 Fig). We also used summary-level genetic association data from FinnGen—the national wide network of Finnish biobanks (N = up to 123,579 women) [20]. All studies had ethical approval from relevant national or local bodies and participants provided written informed consent. Details of the recruitment, information on genetic data, and measurements of baseline characteristics of each cohort are described in S1 Text. This study was initiated using UKB in January 2018, with extra exploration of insomnia IVs and MR sensitivity analyses completed in February 2020 [21]. We searched for additional cohorts till July 2021, and harmonization across the cohorts had to be made continuously. Therefore, we did not have a prespecified analysis plan.
Outcomes measures
We explored potential effects of insomnia on 8 binary outcomes: ever experiencing stillbirth, ever experiencing miscarriage, GD, HDP, perinatal depression, PTB (gestational age <37 completed weeks), LBW (<2,500 grams), and HBW (>4,500 grams). Full details about how these outcomes were measured and derived in each participating study and how we harmonised them across studies can be found in S1 Table. We were not able to measure pre-eclampsia and gestational hypertension separately, because of the small number of definite cases of pre-eclampsia, and because of differences between studies in data collection and definitions.
In UKB, gestational age was only available for a small subset of women (N = 7,280) who delivered a child during or after 1989, the earliest date for which linked hospital labour and perinatal data are available [22]. As a result, numbers with data on PTB are smaller than for any other outcome, and we a priori decided to examine associations with LBW and HBW rather than small-for-gestational age and large-for-gestational age. For most outcomes in UKB, women reported their experience retrospectively in a questionnaire completed at recruitment when they were aged 40 to 60 years.
In the 3 birth cohorts, most outcomes were prospectively obtained (from self-report or clinical records) during an index pregnancy and the perinatal period. The 2 exceptions were history of stillbirth and miscarriage, which were retrospectively reported at the time of the index pregnancy when women were asked if they had ever experienced a (previous) stillbirth or miscarriage. We explored the possibility of examining associations with miscarriage and stillbirth in the index pregnancy. However, numbers were too small for reliable results, and for miscarriage, we were concerned about misclassification or selection bias due to women who had experienced a miscarriage prior to recruitment. Therefore, we used MR to explore the association of susceptibility to insomnia on a history of miscarriage and stillbirth and did not undertake any multivariable regression analyses for these 2 outcomes as suggested during peer review. There were a small proportion of women who contributed more than 1 pregnancy (<5% of total N for each outcome). Given that choosing the first pregnancy could introduce selection towards younger age, lower parity, and higher morbidity of HDP [23], we followed EGG consortium convention [24] to randomly select 1 pregnancy per woman [25].
Data from FinnGen were available for 4 of our outcomes: ever experiencing miscarriage, GD, HDP, and PTB, which were defined based on International Classification of Diseases codes.
Insomnia measures
Self-reported information on insomnia was obtained from 2 of the studies. In UKB, information on lifetime insomnia was used to generate SNP-insomnia associations in women for use in MR analyses in UKB and the birth cohorts. ALSPAC collected data on insomnia during pregnancy, and this was used for conventional confounder-adjusted multivariable regression.
In UKB, insomnia was self-reported at recruitment via the question “Do you have trouble falling asleep at night or do you wake up in the middle of the night?” with responses “never/rarely,” “sometimes,” “usually,” and “prefer not to answer.” For our analyses, we collapsed these categories to generate a binary variable of usually experiencing insomnia (i.e., “usually” [cases] versus “sometimes” + “never/rarely” [controls]) as this was how the responses were categorised in the published genome-wide association study (GWAS) that we have used to select genetic IVs [15].
In ALSPAC, insomnia in pregnancy was self-reported, at 18 and 32 weeks of gestation, using the question “Can you get off to sleep alright?” with options “Very often,” “Often,” “Not very often,” and “Never.” At each time point, we compared “Not very often” + “Never” [cases] versus “Very often” + “Often” [controls]. We acknowledge that the 2 studies are using different questions and that definitions of insomnia vary across published literature [2]. For ease of reading throughout the paper, we refer to results reflecting genetic susceptibility to insomnia (MR) and reporting insomnia in pregnancy (multivariable regression).
SNP selection and SNP-insomnia associations
To identify genetic IVs for insomnia, we searched the GWAS published between January 2017 and February 2021 on PubMed and Neale Lab website [26]. We found 7 insomnia GWAS reporting genome-wide significant SNPs (details in S2 Table). Of these, we selected SNPs from the largest GWAS (total N = 709,986 women, 29% from UKB, and 71% from 23andMe), which provided female-specific results [15]. This GWAS identified 83 loci containing 87 lead SNPs that were robustly associated with insomnia (P-value < 5 × 10−8) after pooling UKB and 23andMe women together. We removed 6 SNPs that were correlated to other SNPs (linkage disequilibrium) at an R2 threshold of 0.01 or higher, based on all European samples from the 1,000 genome project [27]. Associations (reported in log odds ratios [ORs]) of the remaining 81 lead SNPs from the women only GWAS were extracted and listed in S3 Table.
We followed the standard IV approach [28], first using linear regression with individual-level data from 208,140 UKB women to obtain SNP-insomnia association summary data for 2-sample MR analyses. This provides estimates on the risk difference scale, which is more interpretable and comparable to our multivariable regression results [29]. We adjusted the linear models for genotyping batch, top 40 principal components (PCs) and women’s age. During peer review, we were asked to regenerate SNP-insomnia associations using logistic regression to repeat MR analyses. Therefore, we reconducted: (I) split-sample analyses in UKB by generating SNP-insomnia and SNP-outcome associations in logistic regression; (II) 2-sample MR using SNP-insomnia associations generated in logistic regression by the GWAS where we selected our IVs [15], and the pooled SNP-outcome associations combining ALSPAC, BiB, MoBa, and FinnGen; and (III) a meta-analysis of MR estimates from (I) and (II) using fixed-effects (with inverse variance weights) for each insomnia-outcome pair. Consistent with a previous MR study of binary exposures [30], our MR estimates were reported as odds ratios per 1 unit higher log-odds of insomnia.
SNP-outcome associations
We estimated the associations between maternal SNPs and outcomes (log OR and standard errors) for each of the 81 insomnia-related SNPs. In UKB, we randomly separating women in half (giving 2 datasets, A and B) for our split cross-over 2-sample MR [31], given UKB was also included in the GWAS of insomnia. We then estimated SNP-outcome associations in each split sample using logistic regression, adjusting for genotyping batch, top 40 PCs, and women’s age. In the birth cohorts, we estimated the SNP-outcome associations using logistic regression, adjusting for (I) top 20 PCs and women’s age in ALSPAC; (II) top 10 PCs and women’s age in BiB; and (III) genotyping batch, top 10 PCs, and women’s age in MoBa. We extracted associations of the 81 SNPs with the following from FinnGen (words in brackets are the outcome labels from FinnGen): miscarriage (O15_ABORT_SPONTAN), GD (GEST_DIABETES), HDP (O15_GESTAT_HYPERT), and PTB (O15_PRETERM). These summary data were generated by FinnGen using the R-package called SAIGE that fits mixed-effects logistic regression [32], adjusting for genotyping batch, top 10 PCs, and women’s age [20]. Then, we meta-analysed associations from ALSPAC, BiB, MoBa, and FinnGen using fixed-effects with inverse variance weights. Two SNPs (i.e., rs10947428 and rs117037340) were excluded from BiB analyses due to their minor allele frequency lower than 1%.
Assessment of confounders in ALSPAC for multivariable regression
We considered maternal age at time of delivery, education, body mass index at 12 weeks of gestation, smoking status in pregnancy, alcohol intake in the first 3 months of pregnancy, and household occupational social class as potential confounders based on their known or plausible associations with maternal insomnia and pregnancy and perinatal outcomes. Details of confounders were based on maternal self-report and are fully described in S1 Text.
Statistical analyses
Two-sample MR.
As shown in Fig 1, we conducted 2-sample MR analyses of maternal insomnia on pregnancy and perinatal outcomes. In UKB, we conducted a split cross-over 2-sample MR [31]. Specifically, we used SNP-insomnia associations from dataset A and SNP-outcomes associations from dataset B (A on B) and vice-versa (B on A), and then meta-analysed the MR estimates from the 2 together for each insomnia-outcome pair using fixed-effects (with inverse variance weights). For the 2-sample MR using the rest of the cohorts, we used SNP-insomnia associations from UKB women and the pooled SNP-outcome associations combining ALSPAC, BiB, MoBa, and FinnGen. For each outcome, we pooled MR estimates from all cohorts using fixed-effects (with inverse variance weights) and used leave-one (study)-out analysis (initially across all cohorts and then among non-UKB cohorts during peer review) to assess the degree of heterogeneity between cohorts.
(a) Two-sample MR methods include: IVW, MR-Egger, weighted median, and leave-one-out analysis. (b) Multivariable regression analysis adjusted for maternal age at time of delivery, social class, education, body mass index at 12 weeks of gestation, smoking status in pregnancy, and alcohol intake in the first 3 months of ALSPAC pregnancy. ALSPAC, Avon Longitudinal Study of Parents and Children; BiB, Born in Bradford; GWAS, genome-wide association study; IVW, inverse variance weighted; MoBa, Norwegian Mother, Father and Child Cohort Study; MR, mendelian randomization; SNP, single-nucleotide polymorphism; UKB, UK Biobank.
In the main analyses, we used the MR inverse variance weighting (IVW) method, which is a regression of the estimates for SNP-outcomes associations on SNP-insomnia associations weighted by the inverse of the SNP-outcome associations variances, with the intercept of the regression line forced through zero [33]. The IVW estimates should provide an unbiased estimate of a causal effect in the absence of unbalanced horizontal pleiotropy [33]. To explore potential unbalanced horizontal pleiotropy, our sensitivity analyses included (I) estimating between-SNP heterogeneity (which if present may be due to one or more SNPs having horizontal pleiotropic effects on the outcome) using Cochran’s Q-statistic and leave-one (SNP)-out analysis; and (II) undertaking analyses with weighted median [34] and MR-Egger [35], which are more likely to be robust in the presence of invalid IVs. The weighted median method is unbiased so long as less than 50% of the weight is from invalid instruments (i.e., if 1 SNP contributing more than 50% of the weight across the SNP-insomnia associations or several SNPs that contribute more than 50% introduce horizontal pleiotropy the effect estimate is likely to be biased) [34]. MR-Egger is similar to IVW except it does not constrain the regression line to go through zero; if the MR-Egger intercept is not null, it suggests the presence of unbalanced horizontal pleiotropy, and the MR-Egger slope provides an effect estimate corrected for unbalanced horizontal pleiotropy [35]. However, MR-Egger has considerably less statistical power than IVW. Further details of these MR methods are provided in our previous study [21]. When using MR to assess the effect of maternal exposures in pregnancy on offspring outcomes, results might be biased via a path from maternal genotypes to maternal/offspring outcomes due to fetal genotype [36]. To explore this, we compared SNP-outcome associations with versus without adjustments for fetal genotypes in the pooled birth cohort analyses.
We evaluated the strength of IVs using both proportion of variances of maternal insomnia explained by the 81 SNPs (R2) and F-statistic [37]. We selected SNPs robustly related to insomnia in the general female population rather than in pregnant women. Therefore, we explored associations of the 81 SNPs with woman’s insomnia measured at 18 and 32 weeks of gestation in ALSPAC using logistic regressions to determine whether those SNPs related similarly to insomnia in pregnancy. We adjusted for the top 20 PCs and women’s age. As suggested during peer review, we used Steiger filtering to identify SNPs explaining substantially more of the variance in an outcome than in insomnia (i.e., P-value < 0.05) [38] and reconducted MR IVW after removing those SNPs (listed in S3 Table).
Multivariable regression in ALSPAC.
In ALSPAC, we explored the observational associations of insomnia at 18 weeks of gestation with binary outcomes using logistic regression, with adjustment for measured confounders. During peer review, insomnia at 32 weeks of gestation was not considered in the analysis due to potential reverse causality for some outcomes.
All analyses were performed using R 3.5.1 (R Foundation for Statistical Computing, Vienna, Austria). Two-sample MR analyses were conducted using the “TwoSampleMR” R package [27]. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline, specific for MR (S1 STROBE Checklist) [39].
Results
Table 1 summarizes the characteristics of included women from UKB, ALSPAC, BiB, MoBa, and FinnGen. The SNP-insomnia associations in UKB and ALSPAC are listed in S4 Table. The 81 SNPs explained approximately 0.42% of the variance of insomnia among the 208,140 UKB women included in this study (S4 Table), and the mean F-statistic of the 81 SNPs was 11. The pooled 81 SNP-insomnia associations at 18 (OR 1.02 per effect allele, 95% CI: 1.01, 1.03, p = 0.004) and 32 (OR 1.02 per effect allele, 95% CI: 1.01, 1.03, p < 0.001) weeks of gestation in ALSPAC were in the same direction as (but weaker than) the pooled association in the original GWAS of UKB plus 23andMe women (OR 1.05 per effect allele, 95% CI: 1.05, 1.06, p < 0.001). The SNP-outcome associations in UKB, ALSPAC, BiB, and MoBa are listed in S5 Table.
Two-sample MR
In MR IVW combining all cohorts, point estimates for associations between lifetime susceptibility to insomnia (versus no insomnia) and outcomes ranging from ORs of 1.20 (95% CI: 0.52, 2.77, p = 0.67) for GD, to 3.56 (95% CI: 1.49, 8.54, p = 0.004) for perinatal depression (Fig 2). Despite combining data from the largest genetic studies available estimates were imprecise, with 95% CIs for all but 3 outcomes including the null. The 3 that did not include the null were miscarriage, perinatal depression, and LBW (Fig 2). S2 Fig shows IVW results for leave-one (study)-out analysis. Results were broadly consistent but dominated by large cohorts (e.g., UKB and FinnGen), with the point estimates inflated and very wide CIs in small birth cohorts. We further removed 26, 1, and 7 SNPs from analyses for stillbirth, perinatal depression, and LBW, respectively (S3 Table), because Steiger filtering suggested these SNPs potentially more associated with the respective outcome than with susceptibility to insomnia (see Methods). MR IVW estimates after Steiger filtering were consistent for perinatal depression, slightly attenuated for LBW, and in the opposite direction for stillbirth with overlapped CIs both including the null (Fig 2).
(a) For those outcomes that FinnGen did not contribute to. (b) Steiger filtering did not suggest a removal of any SNPs from MR analyses (details shown in S3 Table). (c) p for Cochran’s Q-statistic <0.05 suggests between-study heterogeneity in the meta-analysis. CI, confidence interval; IVW, inverse variance weighted; MR, mendelian randomization; OR, odds ratio; SNP, single-nucleotide polymorphism; UKB, UK Biobank.
Sensitivity analyses using weighted median and MR-Egger for all outcomes were directionally consistent with IVW but attenuated to the null for all outcomes (except stillbirth), and MR-Egger results for GD, perinatal depression and PTB were attenuated to the null (Fig 2). Between-SNP heterogeneity for MR analyses was observed with LBW and HDP (S6A Table), but leave-one (SNP)-out analyses were consistent with the main IVW estimates including all SNPs for all outcomes (S3–S5 Figs). The MR-Egger intercept p-value indicated unbalanced horizontal pleiotropy only for perinatal depression in UKB (S6A Table). Adjusting for fetal genotype (only possible in the birth cohorts) did not alter the SNP-outcome associations with stillbirth, miscarriage, LBW, or HBW; SNP-outcome associations with GD, HDP, and perinatal depression were slightly attenuated; SNP-PTB associations moved slightly away from the null (S6 Fig).
After combining all cohorts, most MR estimates based on SNP-insomnia associations from linear (Fig 2) versus logistic (S6B Table) regression were in the same directions (S7 Table). An association of lifetime susceptibility to insomnia with HBW was observed using IVW (S6B Table), which previously had a wide 95% CI including the null (Fig 2).
Multivariable regression in ALSPAC
Tables 2 and S7 summarize the characteristics of women from ALSPAC. After adjusting for maternal age, education, BMI, smoking, alcohol intake, and household occupational social class, there was an association of insomnia (versus no insomnia) at 18 weeks of gestation with perinatal depression (OR 2.96, 95% CI: 2.42, 3.63, p < 0.001, Fig 3). Associations with other outcomes had imprecise 95% CIs including the null, although their point estimates were in similar magnitudes to those seen in MR (Fig 3).
(a) We adjusted for maternal age at time of delivery, education, body mass index at 12 weeks of gestation, smoking status in pregnancy and alcohol intake in the first 3 month of ALSPAC pregnancy, and household occupational social class. (b) The numbers of women in adjusted models are slightly smaller than those in crude models due to missingness (<8%) in these covariates. ALSPAC, Avon Longitudinal Study of Parents and Children; CI, confidence interval; OR, odds ratio.
Discussion
To the best of our knowledge, this is the first MR study to explore the relationship of insomnia with pregnancy and perinatal outcomes. We interpreted the MR results as reflecting a lifetime susceptibility to insomnia on the basis that SNPs are determined at conception, and evidence suggested that with similar analyses of other exposures (e.g., blood pressure and C-reactive protein) this is the case [40,41]. We interpreted the multivariable regression results as reflecting associations of insomnia during pregnancy, though we could not distinguish this from preexisting insomnia as we did not have information on sleep traits before conception. The associations of the insomnia genetic IVs with reported insomnia during pregnancy in ALSPAC provided some support that the exposures in our MR and multivariable regression analyses had some consistency with each other. Overall, our MR results provide some evidence that a lifetime susceptibility to insomnia might be associated with higher risks of miscarriage, perinatal depression, and LBW. We did not observe evidence to support associations between genetically predicted insomnia and stillbirth, GD, HDP, PTB, and HBW. In multivariable regression, we were unable to assess associations with miscarriage in the index pregnancy. Result for perinatal depression were consistent with the MR results, but this was not the case for LBW, for which no significant association with insomnia reported at 18 weeks gestation was observed.
Our findings in both MR and multivariable regression of an association of insomnia with perinatal depression are consistent with the systematic review and meta-analysis of observational studies [7], and with randomized control trials suggesting that pregnancy intervention with cognitive behavioural therapy to reduce insomnia decreases perinatal depression [8,9]. Recent systematic reviews have only identified 1 cross-sectional study of the association of insomnia with stillbirth [6,42]. This cross-sectional study compared outcomes between 190 women reporting poor sleep quality and 30 women who did not and found no association with stillbirth, although this was not the main focus of the paper [43]. We did not identify any previous studies of insomnia associations with miscarriage. Thus, our novel finding of a possible association of insomnia with miscarriage in MR warrants replication, and larger studies that support analyses with both miscarriage and stillbirth would be valuable. Previous systematic reviews of observational associations of insomnia with GD (OR 1.37, 95% CI: 1.12, 1.69), HDP (OR 1.72, 95% CI: 1.16, 2.56), and PTB (OR 1.49, 95% CI: 1.17, 1.90) are directionally consistent but with stronger ORs than our main MR results [6]. These stronger associations could be due to insufficient adjustment of potential confounders and reverse causality, as many cross-sectional studies and unadjusted associations were included in the meta-analyses.
Several mechanisms have been suggested for why insomnia might influence pregnancy and perinatal outcomes, including insomnia resulting in increased risks of adiposity and insulin resistance that could then influence related pregnancy outcomes (GD, HDP, and HBW). Insomnia has also been suggested to influence maternal blood pressure and placental function which in turn would increase risks of HDP, miscarriage, stillbirth, and PTB. MR analyses support causal associations of insomnia with coronary heart disease, higher glycated haemoglobin, and higher glycoprotein acetyls (an inflammatory marker) in general populations of women and men [17,44,45]. Thus, an increase in cardio-metabolic risk and inflammation may mediate effects of insomnia on miscarriage and LBW, and outcomes for which our MR analyses are currently imprecise. Similarly, MR analyses have found a causal association of insomnia with depressive symptoms [17], which is coherent with our findings in relation to perinatal depression.
Key strengths of our study are that (I) to the best of our knowledge, it is the first study to use MR to explore associations of insomnia with pregnancy and perinatal outcomes; (II) we conducted confounder-adjusted multivariable regression of insomnia in pregnancy in ALSPAC—a larger sample than most previous studies; and (III) we explored a range of pregnancy and perinatal outcomes in 1 paper.
Our MR analyses may be biased by horizontal pleiotropy, particularly given our previous research showing that SNPs for insomnia are also associated with several factors that could influence pregnancy and perinatal outcomes, including education, age at first live birth, and smoking [21]. We explored this potential with a range of sensitivity analyses, including exploring between-SNP heterogeneity and using weighted median and MR-Egger methods that are more robust to such bias than IVW [33]. Results from these sensitivity analyses were broadly consistent with point estimates from IVW; however, the associations between insomnia and miscarriage, perinatal depression, and LBW no longer reached statistical significance. The wider 95% CIs observed could be attribute to the fact that those sensitivity analyses are known to have less statistical power [46]. Those attenuations towards the null could be due to weak IVs, and MR-Egger point estimates are known to be attenuated more severely than weighted median ones [13,46]. Further MR studies in larger samples with more cases would be needed for all outcomes. Specially, our results for perinatal depression require further validation using multivariable MR to account for unbalanced horizontal pleiotropy. Adjusting for fetal genotype did not alter results suggesting that bias due to fetal genotypic effects is unlikely. We did not further adjust for paternal genotype because of limited data with paternal, maternal, and offspring genotype. Furthermore, the most plausible mechanism for paternal genotype to affect pregnancy outcomes is via fetal genotype, which we have adjusted for. Interpretation of our MR estimates requires a further assumption of monotonicity in the SNP-insomnia associations. This requires that all of the women with genetic IVs related to higher susceptibility to insomnia symptoms should report more symptoms (compared to those with fewer alleles related to insomnia)—i.e., that they are “compliers” [47]. The monotonicity assumption cannot be tested. A previous study indicated potential bias when the standard IV approach (see [28]) was used for a nonlinear model [48]. In our study, using linear versus logistic regression to obtain SNP-insomnia associations showed consistent directions between MR estimates. However, magnitudes of MR estimates cannot be compared directly due to their different scales. Further MR studies of binary exposures could apply both approaches to explore an association.
Both our MR and multivariable regression estimates could be vulnerable to selection bias, which has been extensively discussed in previous papers [25,49,50]. UKB is a selective sample (5.5% response to invitation) of adults who are healthier and better educated than the general UK adult population of the same age [51]. Information on perinatal depression and PTB was only available in a subsample of UKB women and such missingness might not be at random [52,53]. By definition our study only includes women who have experienced at least 1 pregnancy, and if insomnia influences fertility then our results might be biased [54]. However, we are not aware of robust evidence of insomnia (or SNPs related to insomnia) influencing infertility or number of children [55,56], suggesting any selection bias through only including pregnant women is unlikely to have a meaningful impact on our MR estimates [54,57].
Insomnia was measured via one self-administrated question in both UKB and ALSPAC, which could mean the binary exposure is misclassified. Non-differential misclassification of insomnia would be expected to bias MR results away from the null (given the attenuated genetic IVs-insomnia associations is the denominator), but multivariable regression results towards the null [58,59]. Similarly, there may be misclassification in some of our outcomes because of the absence of universal testing (e.g., GD in ALSPAC [60]), assessment via self-report questionnaires (e.g., birthweight in UKB), or differences between studies in definitions (e.g., in older women in UKB the gestational age thresholds for defining stillbirth and miscarriage would have differed from those used in the more contemporary birth cohorts). Non-differential misclassification of our binary outcomes would be expected to bias both MR and multivariable regression results towards the null [58,59]. Moreover, the first live-born babies of UKB women are known to be lighter than babies with various birth orders from the more contemporary birth cohorts [61,62].
Although we examined the possibility of reverse causality for individual SNPs using Steiger filtering in MR, this test could be influenced by measurement errors in insomnia and our outcomes and by confounding with opposite directions for insomnia and the outcomes [38]. Multivariable regression results for maternal outcomes could also be vulnerable to reverse causality, as tendency towards the outcomes might have influenced insomnia reported at 18 weeks of gestation. As the sources of bias in our 2 methods (MR and multivariable regression) differ, consistent results between them could strength confidence in the findings even considering different timings of an exposure [63–67]. Our previous study discussed how timings affected the interpretation of MR estimates for insomnia [21]. The similarity of the multivariable regression and MR results for GD, HDP, and perinatal depression suggests it is unlikely that residual confounding has biased regression results, horizontal pleiotropy has substantially affect MR results, or different sources of selection bias in the 2 have a strong impact, for these outcomes. The associations for PTB, LBW, and HBW were attenuated to the null compared to MR results. These suggest possible masking confounding, other biases specific to the multivariable regression, or that the MR is estimating an accumulative effect of insomnia across the life course [68], whereas the observational analyses reflect exposure only from 18 weeks of gestation to occurrence of outcome. MR analyses are statistically inefficient and despite combining relevant studies in order to increase sample size, several of our MR and multivariable regression estimates are imprecise due to small numbers of cases. Our study is limited to women of European ancestry, and we cannot assume that our results generalize to other populations.
Our findings provide some evidence for associations between insomnia and adverse pregnancy outcomes, raising the possibility that interventions to improve healthy sleep (e.g., cognitive behavioural therapy) in women of reproductive age might be beneficial to a healthy pregnancy. However, we acknowledge the need for further MR studies based on larger GWAS of pregnancy and perinatal outcomes, larger observational studies, and studies in women from ethnic backgrounds other than white European. Further studies on the association of insomnia with recurrent miscarriage would help policy makers decide whether to allocate sleep interventions to women with a history of miscarriage when they prepare to be pregnant again.
In conclusion, our study raises the possibility of associations between insomnia and miscarriage, perinatal depression, and LBW.
Supporting information
S1 STROBE Checklist. STROBE-MR checklist of recommended items to address in reports of mendelian randomization studies.
https://doi.org/10.1371/journal.pmed.1004090.s001
(DOCX)
S2 Fig. Leave-one (study)-out analyses of MR IVW estimates.
https://doi.org/10.1371/journal.pmed.1004090.s004
(DOCX)
S3 Fig. Leave-one-out sensitivity analysis in UK Biobank (dataset A on dataset B).
https://doi.org/10.1371/journal.pmed.1004090.s005
(DOCX)
S4 Fig. Leave-one-out sensitivity analysis in UK Biobank (dataset B on dataset A).
https://doi.org/10.1371/journal.pmed.1004090.s006
(DOCX)
S5 Fig. Leave-one-out sensitivity analysis in ALSPAC, BiB, MoBa, and FinnGen.
https://doi.org/10.1371/journal.pmed.1004090.s007
(DOCX)
S6 Fig. Maternal SNP-outcome associations comparing unadjusted to adjusted for fetal genotypes.
https://doi.org/10.1371/journal.pmed.1004090.s008
(DOCX)
S1 Table. Definitions of pregnancy and perinatal outcomes.
https://doi.org/10.1371/journal.pmed.1004090.s009
(XLSX)
S2 Table. Key characteristics of recently conducted genome-wide association studies of insomnia.
https://doi.org/10.1371/journal.pmed.1004090.s010
(XLSX)
S3 Table. Characteristics of genome-wide significant genetic variants for insomnia in women.
https://doi.org/10.1371/journal.pmed.1004090.s011
(XLSX)
S4 Table. Associations of the 81 SNPs with insomnia in UKB and ALSPAC.
https://doi.org/10.1371/journal.pmed.1004090.s012
(XLSX)
S5 Table. Associations of the 81 SNPs with outcomes in UKB, ALSPAC, BiB, and MoBa.
https://doi.org/10.1371/journal.pmed.1004090.s013
(XLSX)
S6 Table. Two-sample MR estimates in UKB and in the birth cohorts and FinnGen.
https://doi.org/10.1371/journal.pmed.1004090.s014
(XLSX)
S7 Table. Characteristics of women in ALSPAC by insomnia at 32 weeks of gestation.
https://doi.org/10.1371/journal.pmed.1004090.s015
(XLSX)
Acknowledgments
This research has been conducted using the UKB Resources under application number 23938. The authors would like to thank the participants and researchers from UKB who contributed or collected data. We are extremely grateful to all the families who took part in ALSPAC, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, and nurses. BiB is only possible because of the enthusiasm and commitment of the Children and Parents in BiB. We are grateful to all the participants, teachers, school staff, health professionals, and researchers who have made BiB happen. This research has been conducted using MoBa data using application number 2552. MoBa is supported by the Norwegian Ministry of Health and Care services and the Ministry of Education and Research. We are grateful to all the participating families in Norway who take part in this ongoing cohort study. We thank the Norwegian Institute of Public Health (NIPH) for generating high-quality genomic data. This research is part of the HARVEST collaboration, supported by the Research Council of Norway (#229624). We also thank the NORMENT Centre for providing genotype data, funded by the Research Council of Norway (#223273), South East Norway Health Authority, and KG Jebsen Stiftelsen. We further thank the Center for Diabetes Research, the University of Bergen for providing genotype data and performing quality control and imputation of the data funded by the ERC AdG project SELECTionPREDISPOSED, Stiftelsen Kristian Gerhard Jebsen, Trond Mohn Foundation, the Research Council of Norway, the Novo Nordisk Foundation, the University of Bergen, and the Western Norway health Authorities (Helse Vest). The authors thank FinnGen investigators for sharing their summary-level data.
References
- 1. Morin CM, Benca R. Chronic insomnia. Lancet. 2012;379:1129–1141. pmid:22265700
- 2. Roth T. Insomnia: definition, prevalence, etiology, and consequences. J Clin Sleep Med. 2007;3:S7–10. pmid:17824495
- 3. Pengo MF, Won CH, Bourjeily G. Sleep in Women Across the Life Span. Chest. 2018;154:196–206. pmid:29679598
- 4. Sedov ID, Cameron EE, Madigan S, Tomfohr-Madsen LM. Sleep quality during pregnancy: A meta-analysis. Sleep Med Rev. 2018;38:168–176. pmid:28866020
- 5. Zhang X, Zhang R, Cheng L, Wang Y, Ding X, Fu J, et al. The effect of sleep impairment on gestational diabetes mellitus: a systematic review and meta-analysis of cohort studies. Sleep Med. 2020;74:267–277. pmid:32862011
- 6. Lu Q, Zhang X, Wang Y, Li J, Xu Y, Song X, et al. Sleep disturbances during pregnancy and adverse maternal and fetal outcomes: a systematic review and meta-analysis. Sleep Med Rev. 2021;58:101436. pmid:33571887
- 7. Yang Z, Zhu Z, Wang C, Zhang F, Zeng H. Association between adverse perinatal outcomes and sleep disturbances during pregnancy: a systematic review and meta-analysis. J Matern Fetal Neonatal Med. 2020:1–9. pmid:31902261
- 8. Felder JN, Epel ES, Neuhaus J, Krystal AD, Prather AA. Efficacy of Digital Cognitive Behavioral Therapy for the Treatment of Insomnia Symptoms Among Pregnant Women: A Randomized Clinical Trial. JAMA Psychiat. 2020;77:484–492. pmid:31968068
- 9. Manber R, Bei B, Simpson N, Asarnow L, Rangel E, Sit A, et al. Cognitive Behavioral Therapy for Prenatal Insomnia: A Randomized Controlled Trial. Obstet Gynecol. 2019;133:911–919. pmid:30969203
- 10. Kalmbach DA, Cheng P, O’Brien LM, Swanson LM, Sangha R, Sen S, et al. A randomized controlled trial of digital cognitive behavioral therapy for insomnia in pregnant women. Sleep Med. 2020;72:82–92. pmid:32559716
- 11. Davey Smith G, Lawlor DA, Harbord R, Timpson N, Day I, Ebrahim S. Clustered environments and randomized genes: a fundamental distinction between conventional and genetic epidemiology. PLoS Med. 2007;4:e352. pmid:18076282
- 12. Davey Smith G, Ebrahim S. ’Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22. pmid:12689998
- 13. Lawlor DA. Commentary: Two-sample Mendelian randomization: opportunities and challenges. Int J Epidemiol. 2016;45:908–915. pmid:27427429
- 14. Yuan S, Larsson SC. An atlas on risk factors for type 2 diabetes: a wide-angled Mendelian randomisation study. Diabetologia. 2020;63:2359–2371. pmid:32895727
- 15. Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, et al. Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways. Nat Genet. 2019;51:394–403. pmid:30804565
- 16. van Oort S, Beulens JWJ, van Ballegooijen AJ, Grobbee DE, Larsson SC. Association of Cardiovascular Risk Factors and Lifestyle Behaviors With Hypertension: A Mendelian Randomization Study. Hypertension. 2020;76:1971–1979. pmid:33131310
- 17. Lane JM, Jones SE, Dashti HS, Wood AR, Aragam KG, van Hees VT, et al. Biological and clinical insights from genetics of insomnia symptoms. Nat Genet. 2019;51:387–393. pmid:30804566
- 18. Larsson SC, Markus HS. Genetic Liability to Insomnia and Cardiovascular Disease Risk. Circulation. 2019;140:796–798. pmid:31422675
- 19. Brand JS, Gaillard R, West J, McEachan RRC, Wright J, Voerman E, et al. Associations of maternal quitting, reducing, and continuing smoking during pregnancy with longitudinal fetal growth: Findings from Mendelian randomization and parental negative control studies. PLoS Med. 2019;16:e1002972. pmid:31721775
- 20. FinnGen. FinnGen Documentation of R5 release. 2021 May 11 [cited 2021 Jul 27]. Available from: https://finngen.gitbook.io/documentation/v/r5/.
- 21. Yang Q, Sanderson E, Tilling K, Borges MC, Lawlor DA. Exploring and mitigating potential bias when genetic instrumental variables are associated with multiple non-exposure traits in Mendelian randomization. Eur J Epidemiol. 2022;37:683–700. pmid:35622304
- 22.
UK Biobank Hospital inpatient data Version 3.0. 2020 Aug [cited 2021 May 11]. Available from: https://biobank.ctsu.ox.ac.uk/crystal/crystal/docs/HospitalEpisodeStatistics.pdf.
- 23. Hernández-Díaz S, Toh S, Cnattingius S. Risk of pre-eclampsia in first and subsequent pregnancies: prospective cohort study. BMJ. 2009;338:b2255. pmid:19541696
- 24. Beaumont RN, Warrington NM, Cavadino A, Tyrrell J, Nodzenski M, Horikoshi M, et al. Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics. Hum Mol Genet. 2018;27:742–756. pmid:29309628
- 25. Hughes RA, Davies NM, Davey Smith G, Tilling K. Selection Bias When Estimating Average Treatment Effects Using One-sample Instrumental Variable Analysis. Epidemiology. 2019;30:350–357. pmid:30896457
- 26. NealLab. We’re thrilled to announce an updated GWAS analysis of the UK Biobank. 2017 Sep 20 [cited 2021 Feb 28]. Available from: http://www.nealelab.is/uk-biobank.
- 27. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7. pmid:29846171
- 28. Greenland S. An introduction to instrumental variables for epidemiologists. Int J Epidemiol. 2000;29:722–729. pmid:10922351
- 29. Burgess S, Labrecque JA. Mendelian randomization with a binary exposure variable: interpretation and presentation of causal estimates. Eur J Epidemiol. 2018;33:947–952. pmid:30039250
- 30. Østergaard SD, Mukherjee S, Sharp SJ, Proitsi P, Lotta LA, Day F, et al. Associations between Potentially Modifiable Risk Factors and Alzheimer Disease: A Mendelian Randomization Study. PLoS Med. 2015;12:e1001841. pmid:26079503
- 31. Henry A, Katsoulis M, Masi S, Fatemifar G, Denaxas S, Acosta D, et al. The relationship between sleep duration, cognition and dementia: a Mendelian randomization study. Int J Epidemiol. 2019;48:849–860. pmid:31062029
- 32. Zhou W, Nielsen JB, Fritsche LG, Dey R, Gabrielsen ME, Wolford BN, et al. Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies. Nat Genet. 2018;50:1335–1341. pmid:30104761
- 33. Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum Mol Genet. 2018;27:R195–208. pmid:29771313
- 34. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40:304–314. pmid:27061298
- 35. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–525. pmid:26050253
- 36. Lawlor D, Richmond R, Warrington N, McMahon G, Davey Smith G, Bowden J, et al. Using Mendelian randomization to determine causal effects of maternal pregnancy (intrauterine) exposures on offspring outcomes: Sources of bias and methods for assessing them. Wellcome Open Res. 2017;2:11. pmid:28405635
- 37. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21:223–242. pmid:21216802
- 38. Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13:e1007081. pmid:29149188
- 39. Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ. 2021;375:n2233. pmid:34702754
- 40. Ference BA, Bhatt DL, Catapano AL, Packard CJ, Graham I, Kaptoge S, et al. Association of Genetic Variants Related to Combined Exposure to Lower Low-Density Lipoproteins and Lower Systolic Blood Pressure With Lifetime Risk of Cardiovascular Disease. JAMA. 2019;322:1381–1391. pmid:31475726
- 41. Kivimäki M, Lawlor DA, Smith GD, Eklund C, Hurme M, Lehtimäki T, et al. Variants in the CRP gene as a measure of lifelong differences in average C-reactive protein levels: the Cardiovascular Risk in Young Finns Study, 1980–2001. Am J Epidemiol. 2007;166:760–764. pmid:17641153
- 42. Warland J, Dorrian J, Morrison JL, O’Brien LM. Maternal sleep during pregnancy and poor fetal outcomes: A scoping review of the literature with meta-analysis. Sleep Med Rev. 2018;41:197–219. pmid:29910107
- 43. Owusu JT, Anderson FJ, Coleman J, Oppong S, Seffah JD, Aikins A, et al. Association of maternal sleep practices with pre-eclampsia, low birth weight, and stillbirth among Ghanaian women. Int J Gynaecol Obstet. 2013;121:261–265. pmid:23507553
- 44. Bos MM, Goulding NJ, Lee MA, Hofman A, Bot M, Pool R, et al. Investigating the relationships between unfavourable habitual sleep and metabolomic traits: evidence from multi-cohort multivariable regression and Mendelian randomization analyses. BMC Med. 2021;19:69. pmid:33731105
- 45. Liu J, Richmond RC, Bowden J, Barry C, Dashti HS, Daghlas I, et al. Assessing the Causal Role of Sleep Traits on Glycated Hemoglobin: A Mendelian Randomization Study. Diabetes Care. 2022;45:772–781. pmid:35349659
- 46. Burgess S, Foley CN, Allara E, Staley JR, Howson JMM. A robust and efficient method for Mendelian randomization with hundreds of genetic variants. Nat Commun. 2020;11:376. pmid:31953392
- 47. Labrecque J, Swanson SA. Understanding the Assumptions Underlying Instrumental Variable Analyses: a Brief Review of Falsification Strategies and Related Tools. Curr Epidemiol Rep. 2018;5:214–220. pmid:30148040
- 48. Terza JV, Bradford WD, Dismuke CE. The use of linear instrumental variables methods in health services research and health economics: a cautionary note. Health Serv Res. 2008;43:1102–1120. pmid:18546544
- 49. Munafò MR, Tilling K, Taylor AE, Evans DM, Davey Smith G. Collider scope: when selection bias can substantially influence observed associations. Int J Epidemiol. 2018;47:226–235. pmid:29040562
- 50. Paternoster L, Tilling K, Davey Smith G. Genetic epidemiology and Mendelian randomization for informing disease therapeutics: Conceptual and methodological challenges. PLoS Genet. 2017;13:e1006944. pmid:28981501
- 51. Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, et al. Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population. Am J Epidemiol. 2017;186:1026–1034. pmid:28641372
- 52. Adams MJ, Hill WD, Howard DM, Dashti HS, Davis KAS, Campbell A, et al. Factors associated with sharing e-mail information and mental health survey participation in large population cohorts. Int J Epidemiol. 2020;49:410–421. pmid:31263887
- 53. Davis KAS, Coleman JRI, Adams M, Allen N, Breen G, Cullen B, et al. Mental health in UK Biobank—development, implementation and results from an online questionnaire completed by 157 366 participants: a reanalysis. BJPsych Open. 2020;6:e18. pmid:32026800
- 54. Mahmoud O, Dudbridge F, Davey Smith G, Munafo M, Tilling K. A robust method for collider bias correction in conditional genome-wide association studies. Nat Commun. 2022;13:619. pmid:35110547
- 55. Caetano G, Bozinovic I, Dupont C, Léger D, Lévy R, Sermondade N. Impact of sleep on female and male reproductive functions: a systematic review. Fertil Steril. 2021;115:715–731. pmid:33054981
- 56. Barban N, Jansen R, de Vlaming R, Vaez A, Mandemakers JJ, Tropf FC, et al. Genome-wide analysis identifies 12 loci influencing human reproductive behavior. Nat Genet. 2016;48:1462–1472. pmid:27798627
- 57. Gkatzionis A, Burgess S. Contextualizing selection bias in Mendelian randomization: how bad is it likely to be? Int J Epidemiol. 2019;48:691–701. pmid:30325422
- 58. van Smeden M, Lash TL, Groenwold RHH. Reflection on modern methods: five myths about measurement error in epidemiological research. Int J Epidemiol. 2020;49:338–347. pmid:31821469
- 59. Pierce BL, VanderWeele TJ. The effect of non-differential measurement error on bias, precision and power in Mendelian randomization studies. Int J Epidemiol. 2012;41:1383–1393. pmid:23045203
- 60. Lawlor DA, Fraser A, Lindsay RS, Ness A, Dabelea D, Catalano P, et al. Association of existing diabetes, gestational diabetes and glycosuria in pregnancy with macrosomia and offspring body mass index, waist and fat mass in later childhood: findings from a prospective pregnancy cohort. Diabetologia. 2010;53:89–97. pmid:19841891
- 61. Hyppönen E, Power C. An intergenerational study of birthweight: investigating the birth order effect. BJOG. 2004;111:377–379. pmid:15008776
- 62. Ghosh RE, Berild JD, Sterrantino AF, Toledano MB, Hansell AL. Birth weight trends in England and Wales (1986–2012): babies are getting heavier. Arch Dis Child Fetal Neonatal Ed. 2018;103:F264–f270. pmid:28780501
- 63. Lawlor DA, Tilling K, Davey Smith G. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016;45:1866–1886. pmid:28108528
- 64. Richmond RC, Timpson NJ, Felix JF, Palmer T, Gaillard R, McMahon G, et al. Using Genetic Variation to Explore the Causal Effect of Maternal Pregnancy Adiposity on Future Offspring Adiposity: A Mendelian Randomisation Study. PLoS Med. 2017;14:e1002221. pmid:28118352
- 65. Bond TA, Richmond RC, Karhunen V, Cuellar-Partida G, Borges MC, Zuber V, et al. Exploring the causal effect of maternal pregnancy adiposity on offspring adiposity: Mendelian randomisation using polygenic risk scores. BMC Med. 2022;20:34. pmid:35101027
- 66. Hernáez Á, Rogne T, Skåra KH, Håberg SE, Page CM, Fraser A, et al. Body mass index and subfertility: multivariable regression and Mendelian randomization analyses in the Norwegian Mother, Father and Child Cohort Study. Hum Reprod. 2021;36:3141–3151. pmid:34668019
- 67. Magnus MC, Miliku K, Bauer A, Engel SM, Felix JF, Jaddoe VWV, et al. Vitamin D and risk of pregnancy related hypertensive disorders: mendelian randomisation study. BMJ. 2018;361:k2167. pmid:29925546
- 68. Sanderson E, Glymour MM, Holmes MV, Kang H, Morrison J, Munafò MR, et al. Mendelian randomization. Nat Rev Methods Primers. 2022;2:6.