Investigating the relationships between unfavorable sleep and metabolomic traits: evidence from multi-cohort multivariable regression and Mendelian randomization analyses
==========================================================================================================================================================================

* Maxime M Bos
* Neil J Goulding
* Matthew A Lee
* Amy Hofman
* Mariska Bot
* René Pool
* Lisanne S Vijfhuizen
* Xiang Zhang
* Chihua Li
* Rima Mustafa
* Matt J Neville
* Ruifang Li-Gao
* Stella Trompet
* Marian Beekman
* Nienke R Biermasz
* Dorret I Boomsma
* Irene de Boer
* Constantinos Christodoulides
* Abbas Dehghan
* Ko Willems van Dijk
* Ian Ford
* He Gao
* Mohsen Ghanbari
* Bastiaan T Heijmans
* M Arfan Ikram
* J Wouter Jukema
* Dennis O Mook-Kanamori
* Fredrik Karpe
* Annemarie I Luik
* L.H. Lumey
* Arn M.J.M. van den Maagdenberg
* Simon P Mooijaart
* Renée de Mutsert
* Brenda W.J.H. Penninx
* Patrick CN Rensen
* Rebecca C Richmond
* Frits R Rosendaal
* Naveed Sattar
* Robert A Schoevers
* P Eline Slagboom
* Gisela M Terwindt
* Carisha S Thesing
* Kaitlin H Wade
* Carolien A Wijsman
* Gonneke Willemsen
* Aeilko H. Zwinderman
* Diana van Heemst
* Raymond Noordam
* Deborah A Lawlor

## Abstract

**Background** Sleep traits are associated with cardiometabolic disease risk, with evidence from Mendelian randomization (MR) suggesting that insomnia symptoms and shorter sleep duration increase coronary artery disease risk. We combined adjusted multivariable regression (AMV) and MR analyses of phenotypes of unfavourable sleep on 113 metabolomic traits to investigate possible biochemical mechanisms linking sleep to cardiovascular disease.

**Methods** We used AMV (N=17,370) combined with two-sample MR (N=38,618) to examine effects of self-reported insomnia symptoms, total habitual sleep duration, and chronotype on 113 metabolomic traits. The AMV analyses were conducted on data from 10 cohorts of mostly Europeans, adjusted for age, sex and body mass index. For the MR analyses, we used summary results from published European-ancestry genome-wide association studies of self-reported sleep traits and of nuclear magnetic resonance (NMR) serum metabolites. We used the inverse-variance weighted (IVW) method and complemented this with sensitivity analyses to assess MR assumptions.

**Results** We found consistent evidence from AMV and MR analyses for associations of usual vs. sometimes/rare/never insomnia symptoms with lower citrate (−0.08 standard deviation (SD)[95% confidence interval (CI): −0.12, −0.03] in AMV and −0.03SD [−0.07, −0.003] in MR), higher glycoprotein acetyls (0.08SD [95%CI: 0.03, 0.12] in AMV and 0.06SD [0.03, 0.10) in MR]), lower total very large HDL particles (−0.04SD [−0.08, 0.00] in AMV and - 0.05SD [−0.09, −0.02] in MR) and lower phospholipids in very large HDL particles (−0.04SD [−0.08, 0.002] in AMV and −0.05SD [−0.08, −0.02] in MR). Longer total sleep duration associated with higher creatinine concentrations using both methods (0.02SD per 1-hour [0.01, 0.03] in AMV and 0.15SD [0.02, 0.29] in MR) and with isoleucine in MR analyses (0.22SD [0.08, 0.35]). No consistent evidence was observed for effects of chronotype on metabolomic measures.

**Conclusions** Whilst our results suggested that unfavourable sleep traits may not cause widespread metabolic disruption, some notable effects were observed. The evidence for possible effects of insomnia symptoms on glycoprotein acetyls and citrate and longer total sleep duration on creatinine and isoleucine might explain some of the effects, found in MR analyses of these sleep traits on coronary heart disease, which warrant further investigation.

## Introduction

Several systematic reviews and large biobank studies have reported associations of self-reported insomnia symptoms, short and long sleep duration, and chronotype (i.e., having an evening rather than morning preference) with increased risk of cardiovascular disease, type 2 diabetes and risk factors for these1-9. The mechanisms underlying these associations are unclear, and it is plausible that specific sleep traits may contribute to the misalignment of various behavioral and internal physiological processes, including aspects of metabolism that causes adverse cardiometabolic health.

There is some evidence of poor sleep quality, shorter sleep duration and having an evening chronotype being associated with higher triglyceride, total cholesterol and low- density lipoprotein cholesterol (LDL-C) levels and lower high-density lipoprotein cholesterol (HDL-C) concentrations10-12. However, the extent to which these associations are explained by confounding factors, such as body mass index11, is unclear. Beyond conventional multivariable-adjusted regression analyses, we have previously demonstrated that sleep duration modifies the associations of genetic variation with triglycerides, LDL-C and HDL-C in a large sleep-gene interaction analysis, suggesting that possible different biological mechanisms underlie the associations of short and long sleep duration with these lipid traits13. However, these genetic interaction analyses do not assess causality and, like previous multivariable-adjusted regression analyses, have focused on a limited number of lipid traits.

Mendelian randomization (MR) uses genetic variants that are robustly associate with an exposure as an instrumental variable to obtain unconfounded effects of that exposure on an outcome of interest14-16. Recent MR analyses have suggested a causal effect of insomnia symptoms on coronary heart disease17 and of short (<6 hours) sleep duration on myocardial infarction risk18.

The aim of this study was to determine the possible causal effect of sleep traits on metabolomic traits. We compared findings from adjusted multivariable regression (AMV) and MR analysis, to determine the relationships between self-reported insomnia symptoms (usually vs sometimes/rare/never), total habitual sleep duration (per 1 hour longer) and chronotype (evening vs morning preference) and 113 nuclear magnetic resonance (NMR) metabolomic traits. Cross-sectional AMV was performed with adjustment for age, sex and BMI in 17,370 individuals from 10 cohorts of mostly Europeans. Two-sample MR used summary results from genome-wide association studies (GWAS) of different sleep traits in 1,331,010 (insomnia)19, 446,118 (sleep duration)20, and 651,295 (chronotype)21 European adults and summary results from four GWAS of 113 circulating metabolomic measures from NMR in 38,618 European adults. In secondary analyses, we explored effects of short (<7 vs 7-<9 hours) and long (≥9 vs 7-<9 hours) sleep duration on the metabolomic traits. We highlight results that were consistent across both methods, as the different key sources of bias of the two methods (e.g., residual confounding in AMV and unbalanced horizontal pleiotropy in MR, respectively) mean that, where there is consistency, this is more likely to reflect a causal effect22.

## Methods

### Studies used for AMV

Cross-sectional AMV analyses were performed using data from 10 cohorts: the Active and Healthy Ageing (AGO) study23, the Dutch Hunger Winter Families Study (DWFS)24, the Healthy Life in an Urban Setting (HELIUS) Study25, the Leiden University Migraine Neuro-Analysis (LUMINA)26, Netherlands Study of Depression and Anxiety (NESDA)27, the Netherlands Twin Register (NTR)28, the Netherlands Epidemiology of Obesity (NEO)

Study29, and the Rotterdam Study cohorts 1, 2 and 3 (RS1, RS2 and RS3)30. Details, and distributions for characteristics, of each study are given in the **Supplemental Material** and **Supplemental Table S1**. Each participating study obtained written informed consent from all participants and received approval from the appropriate local institutional review boards. Before the analyses, we excluded all participants with diabetes (defined as self-report/hospital record, fasting plasma glucose >7 mmol/L and/or use of hypoglycaemic medication) given the known disturbances on many metabolomic traits.

### Studies used for MR analyses

We performed two-sample MR analyses using publicly available summary-level data14 from the following GWAS:

### Sleep trait GWAS

We selected genome-wide significant (p-value<5e-8) variants as instrumental variables from the following GWAS:

*   **Insomnia**: A GWAS that pooled data from two large biobanks (UK Biobank and 23andMe) and included 1,331,010 unrelated European-ancestry adults. This GWAS identified 248 variants for experience of insomnia symptoms (usually vs sometimes/rare/never)19.

*   **Sleep duration**: A GWAS undertaken in UK Biobank of 446,118 unrelated European-ancestry adults20. This GWAS identified 78 variants for total sleep duration (mean 7.2h; SD1.1h). In addition, this GWAS identified 27 variants for short sleep duration (<7 h vs 7 to <9 hours; N = 106,192 cases) and 8 for long sleep duration (≥9 h vs 7 to <9 hours; N = 34,184 cases).

*   **Chronotype**: A GWAS that pooled data from two large biobanks (UK Biobank and 23andMe) and included 697,828 unrelated European-ancestry adults (651,295 of whom were in the combined (both biobanks) GWAS of morning versus evening preference that we have used in this two-sample MR study. This GWAS identified 351 variants for chronotype21. Because previous observational studies have found increased risk of cardiometabolic diseases and risk factors in those with an evening preference, we transformed the GWAS results to reflect alleles associated with evening preference.

### NMR Metabolite GWAS

*   MAGNETIC consortium (N = 24,925)31 with summary-level GWAS data downloaded from [http://www.computationalmedicine.fi/data#NMRGWAS](http://www.computationalmedicine.fi/data#NMRGWAS)

*   In addition, to increase statistical power in the MR analyses, we generated new summary-level GWAS data from three cohorts using similar analyses procedures to the MAGNETIC consortium: 
    
    *   Oxford Biobank (N = 6,616)32;
    
    *   NEO (N = 4,734)29;
    
    *   Pravastatin in Elderly Individuals at Risk of Vascular Disease (PROSPER) (N = 2,343; placebo arm only)33.

All GWASs were undertaken in participants of European ancestry. Further details of the MAGNETIC consortium and the three additional cohorts are provided in **Supplemental Methods**. There was no overlap between the cohorts included in the sleep trait GWAS and those included in the NMR GWAS.

### Sleep traits

In both AMV and MR analyses, sleep traits were self-reported and analysed in the same units/categories. Insomnia symptoms were assessed with a question similar to “Do you have trouble falling asleep at night or do you wake up in the middle of the night?” with the following answers possible: “never/rarely”, “sometimes”, “usually”, or “prefer not to answer”. In the AMV and GWAS analyses participants who answered “usually” were defined as having insomnia symptoms and were compared to those answering “never/rarely” or “sometimes”. Habitual sleep duration was assessed using a question similar to “On an average day, how many hours of sleep do you get?”. For our main analyses, we examined effects of total sleep duration (per 1 hour longer) on metabolomic measures. In secondary analyses, we explored associations of short (<7 vs 7-<9) and long (≥9 vs 7-<9 hours) habitual sleep. These latter two analyses were considered exploratory because of lower statistical power and possible weak instrument bias in the MR analyses. For chronotype, a question similar to “Are you naturally a night person or a morning person?” with the possible responses “Night owl/night person”, “Early bird/morning person”, “Neither/not sure” was used in most studies. A variation on the question in UK Biobank included more responses: “Definitely a morning person”, “More a morning than evening person”, “More an evening than a morning person”, “Definitely an evening person”, “Do not know”. Participants were classified as having a ‘morning preference’ (“Early bird/morning person”, “Definitely a morning person” or “More a morning than evening person”), the reference group, or an ‘evening preference’ (“Night owl/night person”, “More an evening than a morning person” or “Definitely an evening person”). For all traits those responding “do not know”, “unsure” or “prefer not to answer” were excluded.

### NMR-based metabolomic profiling

In both the metabolite GWAS and studies included in the AMV meta-analysis, metabolites were quantified using a high-throughput proton (1H) NMR metabolomics platform34 ([https://nightingalehealth.com/](https://nightingalehealth.com/)) to quantify a maximum of 148 (excluding ratios) lipid and lipoprotein and metabolite concentrations in fasting serum or plasma samples. The quantitative NMR measures include numerous lipid species and fatty acids, as well as some amino acids, markers of glucose homeostasis, fluid balance and an inflammatory marker. This platform has been used widely in population-based studies of cardiometabolic diseases, and has been described in detail elsewhere34-36. There were 113 metabolomic trait measurements that were available for both AMV and MR analyses.

### Statistical Analyses

In both AMV and MR analyses, we estimated the same effect: the difference in mean NMR metabolites (SD units of the natural log-transformed metabolomic traits; as dependent variables) comparing (i) usually experiencing insomnia symptoms to sometimes, rarely or never, (ii) per 1 hour longer habitual sleep duration and (iii) an evening to a morning preference. All analyses were performed in R (v3.6.1)37.

### Multivariable-adjusted regression meta-analysis

Cross-sectional AMV was performed by each of the individual cohorts according to a pre-specified analysis plan and standardized analysis script. Results were collected centrally for quality control subsequent fixed-effect meta-analyses using the R “rmeta” package, using similar procedures as described previously38. AMV analyses adjusted for age, sex and BMI.

### MR analyses

We excluded all palindromic single nucleotide polymorphisms (SNPs) and those in linkage disequilibrium at R2>0.001 (based on the 1000genomes (phase 1) panel). After these exclusions, we searched for all remaining independent sleep-associated variants (149 for insomnia, 57 for total sleep duration [and an additional 25 and 71 variants for short and long habitual sleep duration, respectively] and 208 for chronotype) in the GWAS of NMR metabolomic measures, and the directions of the summary data were harmonised (i.e. making sure that each effect estimate was coded in the same direction with respect to the effect allele as SNP associations from the summary sleep trait data) with those of the sleep trait summary data.

The MRCIEU/TwoSampleMR package was used for harmonization of the exposure and outcome SNPs and to perform the MR analyses16. For our main analyses, we used the multiplicative random effects inverse variance-weighted (IVW) approach39. This method generates a causal estimate of the sleep traits on metabolomic traits by regressing the SNP- sleep trait association on the SNP-metabolomic measure association, weighted by the inverse of the SNP-metabolomic measure association, and constraining the intercept of this regression to zero. Standard errors are corrected to take into account any between SNP heterogeneity and assumes that there is no directional horizontal pleiotropy. To explore this assumption further, we performed sensitivity analyses using MR-Egger40 and weighted- median estimator41 methods. MR-Egger is similar to the IVW method but does not force the regression line (i.e., of the SNP-sleep trait association on the SNP-metabolomic measure association) through an intercept of zero. It is statistically less efficient (providing wider confidence intervals) but provides a causal estimate (i.e., the regression slope) that is corrected for directional horizontal pleiotropy, and a non-zero intercept is an indication of the existence of directional pleiotropy. The weighted-median estimator is valid if more than 50% of the weight of the genetic instrument is from valid variants (i.e., if one single SNP or several SNPs jointly contributing 50% or more of the weight in the MR analysis exhibit horizontal pleiotropy the calculated effect estimate may be biased).

We performed MR analyses for all sleep traits with each of the 4 metabolomic GWAS data sources (MAGNETIC, Oxford Biobank, NEO and PROSPER) and the results were subsequently meta-analysed using fixed-effect meta-analyses as implemented in the R package rmeta.

### Comparing multivariable regression and MR analysis results

Circos plots were used to summarise and visually compare the AMV and the IVW MR results. These were generated using R package EpiViz (version 0.1.0) (further details in the **Supplemental material**). We also generated scatter plots of the AMV vs MR results for each metabolite and compared the linear fit across all metabolites to a slope of perfect concordance and used R2 as a measure of goodness of fit (agreement) between the two methods across all 113 metabolomic traits.

Having compared results for the AMV and IVW MR methods across all metabolites, we then selected all sleep trait-metabolite associations that reached a pre-defined p-value threshold in AMV or IVW MR. We then compared results across AMV, IVW MR, MR- Egger and weighted median MR for those selected associations. Whilst we focus on results reaching a pre-defined p-value threshold in either AMV or IVW MR in the main paper and our conclusion, a full set of all results (AMV, unadjusted MV, IVW MR and all MR sensitivity analyses) are presented in Supplemental **Tables S2 to S11** in the **Supplemental Tables File**. We applied the same Bonferroni multiple testing corrected p-value threshold separately to the AMV and MR analyses. The threshold was determined taking into account the correlation structure of the metabolomic measures by using information from previous studies that have identified 17 principal components, which explain 95% of the metabolomic traits data variance42. Therefore, the two-sided threshold of *P <* 0.05 adjusted for multiple testing becomes *P <* 0.0029 (0.05/17). For any association that passed this threshold with either AMV or IVW MR we considered the result from the second method to be consistent if the point estimate had a similar direction of effect and the p-value for the second association was <0.05. This was justified on the basis that once one method passed the Bonferroni threshold, we were treating that result as a hypothesised effect and seeking replication and triangulation in the second method.

## Results

Full results of all AMV and MR analysis, including MR sensitivity analysis results, are presented in Supplemental **Tables S2 to S11** in the **Supplemental Tables File**.

### Insomnia symptoms

Visually inspecting the circos plot shows there was directional consistency between the AMV and IVW MR results for most of the metabolomic traits (**Figure 1**). With both methods, insomnia symptoms were associated with higher concentrations of small and medium very large density lipoprotein (VLDL) particles, small HDL particles and glycoprotein acetyls, and with lower concentrations of large HDL particles. Across all 113 metabolomic traits, there was good concordance of effect size and direction (**Figure 2**; R2 = 0.57).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F1)

Figure 1: IVW Mendelian randomization estimates and age-, sex- and BMI-adjusted multivariable regression estimates for the associations between insomnia symptoms and 113 NMR-derived metabolomic measures.
Results are expressed as the difference in mean metabolite concentrations (in standard deviation units) between those reporting usually versus sometimes/rarely/never experiencing insomnia symptoms. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; IDL, intermediate density lipoprotein; IVW MR, Inverse variance weighted Mendelian randomization; LDL, low density lipoprotein; NMR, nuclear magnetic resonance; VLDL, very large density lipoprotein.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F2)

Figure 2: Comparison of the point estimates of the IVW Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses for the associations between insomnia symptoms and 113 NMR-derived metabolomic measures.
Each green dot in the scatter plot represents a metabolic trait and the positions of the dots are determined by the differences in mean metabolite concentrations (in standard deviation units) between those reporting usually versus sometimes/rarely/never experiencing insomnia symptoms. These are estimated by Inverse variance weighted (IVW) Mendelian randomization (vertical axes) and age, sex and BMI adjusted multivariable regression (horizontal axes). The vertical grey lines for each dot indicate the 95% confidence intervals (CI) for the Mendelian randomization estimates and the horizontal grey lines for each dot indicate the 95% CI for the adjusted multivariable regression estimates. A linear fit (red dashed line) summarizes the similarity between the two estimates. A slope of 1 with an intercept of 0 (dashed grey line), with all green dots sitting on that line (R2 = 1), would indicate identical magnitude and direction between the two methods. R indicates goodness of linear fit and is a measure of the consistency between the two estimates. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; CI, confidence interval; DHA, 22:6, docosahexaenoic acid; IVW MR, Inverse variance weighted Mendelian randomization, SD, standard deviation.

Associations of insomnia symptoms with the 113 metabolomic traits passed the multiple testing threshold *(P* < 0.0029) for 13 in the AMV analyses, and for 3 in the MR analyses (glycoprotein acetyls passed the threshold in both). Based on our pre-specified definition of consistency (i.e. same direction and p-value <0.05 in MR for any AMV results reaching the corrected p-value, and vice versa) we found consistent evidence from AMV and MR analyses for 4 associations. Specifically, usual vs. sometimes/rare/never insomnia symptoms lowered citrate (−0.08SD [95%CI: −0.12, −0.03] in AMV and −0.03SD [−0.07, - 0.003] in MR), increased glycoprotein acetyls (0.08SD [95%CI: 0.03, 0.12] in AMV and 0.06 [0.03, 0.10] in MR), and lowered total very large HDL particles (−0.04SD [−0.08, 0.00] in AMV and −0.05SD [−0.09, −0.02] in MR) and phospholipids in very large HDL particles (− 0.04SD [−0.08, 0.002] in AMV and −0.05SD [−0.08, −0.02] in MR) (**Figure 3**). MR sensitivity analyses were generally consistent with the main IVW analyses though confidence intervals were wide for the MR-Egger results. (**Figure 3**).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F3)

Figure 3: Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses results for select associations of insomnia symptoms with NMR- derived metabolomic measures.
Figure shows inverse variance weighted (IVW) Mendelian randomization, Mendelian randomization sensitivity (weighted median (WM) and MR-Egger) and adjusted multivariable (AMV) regression analysis results. Results presented were selected on the basis of passing multiple testing threshold for either IVW or AMV (p-values < 0.0029) Theestimates are the difference in mean metabolite (in standard deviation units) between those reporting usually versus sometimes/rarely/never experiencing insomnia symptoms. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; IVW MR, Inverse variance weighted Mendelian randomization; NMR, nuclear magnetic resonance; SD, Standard Error; VLDL, very low density lipoprotein; WM, Weighted Median.

### Sleep duration

Several associations of total sleep duration with metabolomics traits were directionally consistent in the the AMV and MR analyses (**Figure 4**). Where association directions were consistent, the MR results often had a stronger magnitude of association than the AMV results. Consistency of magnitude (as well as direction) was poor to moderate between the two methods (**Figure 5**, R2 = 0.37).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F4)

Figure 4: IVW Mendelian randomization estimates and age-, sex- and BMI-adjusted multivariable regression estimates for the associations between total sleep duration and 113 NMR-derived metabolomic measures.
Results are expressed as the difference in mean metabolite concentrations (in standard deviation units) for each 1 hour greater reported total sleep duration. For visualization purposes the axes have unequal scaling. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; IDL, intermediate density lipoprotein; IVW MR, Inverse variance weighted Mendelian randomization; LDL, low density lipoprotein; NMR, nuclear magnetic resonance; VLDL, very large density lipoprotein.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F5)

Figure 5: Comparison of the point estimates of the IVW Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses for the associations between total sleep duration and 113 NMR-derived metabolomic measures
Each green dot in the scatter plot represents a metabolic trait and the positions of the dots are determined by the differences in mean metabolite concentrations (in standard deviation units) for each 1 hour greater reported total sleep duration. These are estimated by Inverse variance weighted (IVW) Mendelian randomization (vertical axes) and age, sex and BMI adjusted multivariable regression (horizontal axes). The vertical grey lines for each dot indicate the 95% confidence intervals (CI) for the Mendelian randomization estimates and the horizontal grey lines for each dot indicate the 95% CI for the adjusted multivariable regression estimates. A linear fit (red dashed line) summarizes the similarity between the two estimates. A slope of1 with an intercept of 0 (dashed grey line), with all green dots sitting on that line (R2 = 1), would indicate identical magnitude and direction between the two methods. R indicates goodness of linear fit and is a measure of the consistency between the two estimates. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; CI, confidence interval; IVW MR., Inverse variance weighted Mendelian randomization, SD, standard deviation.

Associations for total sleep duration passed the multiple testing threshold for 8 of the 113 metabolomic trait associations in AMV analyses and one in MR. Only one of the 8 AMV associations replicated in the MR analyses (difference in mean creatinine for a 1 hour longer sleep was 0.02SD [0.01, 0.03] in AMV and 0.15 [0.02, 0.29] in MR). Isoleucine was the one metabolite to pass the multiple testing threshold in IVW MR analyses but it did not replicate in AMV analyses (0.01SD [−0.001, 0.02] in AMV and 0.22 [0.08, 0.36] in MR analyses]) (**Figure 6**). For the associations with creatinine, the weighted median MR result was consistent with that of the main (IVW) results but MR-Egger was in the opposite direction (though with very wide confidence intervals). For isoleucine, both MR sensitivity analyses had point estimates that were directionally, and in magnitude, similar to the main IVW MR results (**Figure 6**).

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F6)

Figure 6: Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses results for selected associations of total sleep duration with NMR- derived metabolomic measures.
Figure shows inverse variance weighted (IVW) Mendelian randomization, Mendelian randomization sensitivity (weighted median (WM) and MR-Egger) and adjusted multivariable (AMV) regression analysis results. Results presented were selected on the basis of passing multiple testing threshold for either IVW or AMV (p-values < 0.0029). The estimates are the difference in mean metabolite (in standard deviation units) per 1 hour greater total sleep duration. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; DHA, 22:6, docosahexaenoic acid; HDL, high density lipoprotein; IVW MR, Inverse variance weighted Mendelian randomization; NMR, nuclear magnetic resonance; SD, Standard Error; WM, Weighted Median.

In exploratory analyses, most associations of short sleep duration (<7 hours) were close to the null in both AMV and MR analyses, with very little overall agreement between the two methods (R2=0.09, **Supplemental Figures S1** and **S2**). Two associations of short sleep passed the multiple testing corrected p-value in AMV analyses (22:6 docosahexonic acid (DHA) and omega-3 fatty acids), with short sleep duration associated with lower levels for both of these; none passed the multiple testing threshold in the MR analyses. For docosahexonic acid (DHA) and omega-3 fatty acids there was an inverse association in IVW MR analyses that had a larger effect estimate than in the AMV analyses but with wide confidence intervals that included the null (**Supplemental Figure S3**).

A total of 31 of the 113 metabolites passed the multiple testing threshold in the AMV analyses of long sleep duration (≥9 hours), including higher concentrations of most extremely large, large and medium VLDL, triglycerides, and concentrations of glycoprotein acetyls and isoleucine (**Supplemental Figures S4**). MR analyses did not support a causal effect for any of these, with MR analysis point estimate close to the null or in the opposite direction (**Figure S6**). We did not identify any metabolic traits passing the multiple testing threshold in IVW MR.

### Chronotype

There was very little consistency in direction and magnitude of association between AMV and MR analyses of chronotype with the metabolomic traits (**Figures 7 and 8**, R2 = 0.17). Chronotype was associated with isoleucine after multiple testing correction in the AMV analyses (difference in mean comparing evening to morning preference (0.13SD [0.04, 0.21]), but this was not supported in MR analyses (−0.02 [−0.05, 0.02])) (**Figure 9**). No associations of chronotype with the metabolomics traits passed the multiple testing threshold in the IVW MR analyses.

![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F7.medium.gif)

[Figure 7:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F7)

Figure 7: IVW Mendelian randomization estimates and age-, sex- and BMI-adjusted multivariable regression estimates for the associations between chronotype and 113 NMR-derived metabolomic measures.
Results are the difference in mean metabolite concentrations (in standard deviation units) between those reporting an evening versus morning preference. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; IDL, intermediate Results

![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F8.medium.gif)

[Figure 8:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F8)

Figure 8: Comparison of the point estimates of the IVW Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses for the associations
Each green dot in the scatter plot represents a metabolic trait and the positions of the dots are determined by the differences in mean metabolite concentrations (in standard deviation units) comparing those reporting an evening preference versus morning preference. These are estimated by Inverse variance weighted (IVW) Mendelian randomization (vertical axes) and age, sex and BMI adjusted multivariable regression (horizontal axes). The vertical grey lines for each dot indicate the 95% confidence intervals (CI) for the Mendelian randomization estimates and the horizontal grey lines for each dot indicate the 95% CI for the adjusted multivariable regression estimates. A linear fit (red dashed line) summarizes the similarity between the two estimates. A slope of 1 with an intercept of 0 (dashedgrey line), with all green dots sitting on that line (R = 1), would indicate identical magnitude and direction between the two methods. R indicates goodness of linear fit and is a measure of the consistency between the two estimates. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; CI, confidence interval; IVW MR, Inverse variance weighted Mendelian randomization, SD, standard deviation.

![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/01/2020.08.27.20173518/F9.medium.gif)

[Figure 9:](http://medrxiv.org/content/early/2020/09/01/2020.08.27.20173518/F9)

Figure 9: Mendelian randomization and age-, sex- and BMI-adjusted multivariable regression analyses results for selected associations of chronotype with NMR-derived metabolomic measures
Figure shows inverse variance weighted (IVW) Mendelian randomization, Mendelian randomization sensitivity (weighted median (WM) and MR-Egger) and adjusted multivariable (AMV) regression analysis results. Results presented were selected on the basis of passing multiple testing threshold for either IVW or AMV(p-values < 0.0029). The estimates are the difference in mean metabolite (in standard deviation units) comparing report of having an evening versus morning preference. Abbreviations: AMV, adjusted (age, sex, BMI) multivariable regression; BMI, body mass index; IVW MR, Inverse variance weighted Mendelian randomization; NMR., nuclear magnetic resonance; SD, standard deviation; WM, Weighted Median.

## Discussion

With the present multi-cohort effort, we intended to identify the potential biochemical mechanisms linking sleep to cardiometabolic disease risk. We found consistent evidence with both AMV and MR that usually (vs. sometimes, rarely or never) experiencing insomnia symptoms causes lower concentrations of citrate, total very large HDL particles and phospholipids in very large HDL particles and higher concentrations of glycoprotein acetyls. There was little consistency between AMV and MR results for total habitual sleep duration across all metabolomic traits, though a longer total sleep duration was associated with higher concentrations of creatinine in both methods. For chronotype, whilst having an evening preference was associated with higher isoleucine concentrations at our multiple-testing threshold in the AMV analyses, MR analyses did not support causality. Chronotype did not pass multiple testing with any other metabolites. Therefore, our findings do not support the notion that sleep traits have widespread effects on the investigated metabolomic traits. Nevertheless, they suggest that insomnia symptoms may influence cardiometabolic disease (as previously shown in MR17) through increased inflammation, and also result in lower citrate levels.

The lack of a more widespread impact of sleep traits on multiple metabolomic traits is in contrast with some experimental sleep studies, although direct comparisons are not possible. For example, targeted and untargeted mass spectrometry measurements performed in frequently sampled blood (every 2 hours) from 12 healthy men revealed that 109 out of 171 metabolites exhibited a circadian rhythm43. Furthermore, in controlled experimental conditions this circadian variation was maintained for 78 out of these 109 metabolites over a 24-hour period of total sleep deprivation. For 27 metabolites, including some lipids (13 glycerophospholipids and 3 sphingolipids), as well as tryptophan, serotonin, taurine and 8 acylcarnitines, marked acute increases in concentrations were observed during 24 hours of sleep deprivation compared with the 24 hours of habitual sleep43. Importantly, the MR analyses assessed long-term (lifelong), rather than acute, effects of a predisposition for unfavourable quality or quantity of sleep on metabolic disturbances, which could explain the generally stronger effects in the total sleep duration MR analyses.

Glycoprotein acetyls, which we identified as a novel trait potentially influenced by insomnia symptoms, are elevated in response to infection and inflammation. C-reactive protein (CRP) is the most widely recognized marker of acute and chronic inflammation in epidemiological studies. Whilst observational studies have shown that higher circulating CRP is associated with increased cardiovascular disease risk, MR studies suggest this is not a causal relationship44, 45. Glycoprotein acetyls have emerged as a potentially better measure of cumulative inflammation than CRP, since glycoprotein acetyls increase late in the inflammatory process and levels are relatively stable within individuals over many years46, 47. In AMV analyses in prospective cohorts, glycoprotein acetyls were positively associated with cardiovascular diseases and type 2 diabetes, independently of established risk factors and CRP47. If these associations are shown to be causal, then it is possible that cumulative chronic inflammation, as measured by glycoprotein acetyls, mediates the effect of insomnia on coronary heart disease identified in MR analyses17. However, we acknowledge that our results for the effect of insomnia on glycoprotein acetyls require replication in independent and larger studies and testing in ancestries other than Europeans.

The inverse association of insomnia symptoms with citrate in both AMV and MR analyses is novel. A recent narrative review highlighted the physiological control of plasma citrate concentrations in health and disease48. One possible mechanism through which insomnia might influence citrate is via the relationship of insomnia with night-time eating (which is also accompanied with higher night physical activity)49, which would result in higher TCA cycle activity and consequently lower plasma citrate concentrations. However despite a plausible role, there is a paucity of clinical and epidemiological studies of the effect of citrate levels on disease outcomes48. Therefore, the meaning of a possible effect of insomnia on citrate levels is hard to discern.

We found evidence for associations of experiencing insomnia symptoms with higher concentrations of very large total HDL particles and phospholipids in very large HDL particles. MR and randomized controlled trials suggest that circulating HDL cholesterol is not causally related to cardiovascular disease50-52. The amount of cholesterol carried in HDL particles increases with increasing particle size and emerging evidence highlights the importance of considering size, structure and composition of lipoprotein particles when exploring their effects on cardiovascular disease 53. In AMV analyses, inverse associations of very large, large, medium and small HDL particles with cardiovascular disease have been observed, but these attenuated to the null with adjustment for lipids used by clinicians36. Thus, the relevance of possible insomnia effects on very large HDL particle concentrations, and specifically phospholipids in these particles, is unclear and require additional studies.

We found evidence in both AMV and MR analyses of a possible association of longer total sleep duration with higher creatine concentrations, a biomarker used to estimate kidney function. Established cardiovascular risk factors, such as high blood pressure and type 2 diabetes, are associated with higher creatinine concentrations54. Findings from multivariable regression suggest that the association of kidney function with cardiovascular disease largely reflects confounding and/or reverse causality55. Thus, our observations possibly suggest that longer sleep duration is an additional risk factor for chronic kidney disease rather than cardiovascular diseases, though we acknowledge MR sensitivity analyses did not support a causal effect. We also found a novel association of longer total sleep duration with the branched-chain amino acid isoleucine in MR analyses. Higher concentrations of branched- chain amino acids, including isoleucine, are associated with increased risk of cardiovascular disease36, though this has not been explored in MR studies. MR analyses supports a causal effect of the branched-chain amino acids on type 2 diabetes56 and our results suggest that longer total sleep duration may mediate some of this effect.

Key strengths of our study are its novelty and the comparison of results from the largest AMV study of sleep traits with multiple circulating metabolomic measures22 with equivalent results from MR. We harmonised questionnaire-based sleep data across all contributing studies and the NMR metabolomic platform was consistent across studies in both the AMV and MR analyses. We were able to increase the power of our two-sample MR analyses by combining unpublished summary-level GWAS results from three cohorts (total N=13,693) with those of the largest published GWAS of the same NMR platform (N=24,925) to date31. Two-sample MR assumes that the two samples are from the same underlying population and independent of each other. Given all GWAS were undertaken in adults of European ancestry and the lack of overlap in studies contributing to the metabolite GWAS with any of the sleep trait GWAS, we are confident this assumption is largely met. Most observed differences in mean metabolomic concentrations were close to the null, and in general (true) null results are less subject to bias than non-null results57.

Important limitations include the lack of statistical power, particularly to explore possible non-linear associations for sleep duration. The platform misses a high proportion of currently quantifiable metabolites in human serum/plasma, including markers of energy balance, microbiota metabolism, vitamins, co-factors and xenobiotics, that may be influenced by sleep traits43. Still the NMR platform used in the analyses covers considerably more of the lipidome than conventional clinical chemistry measures (total cholesterol, LDL-C, HDL-C and triglycerides) that have previously been explored, and in addition includes amino acids, glycolysis metabolites, ketone bodies and an inflammatory marker. Whilst we adjusted for age, sex and BMI, the results obtained in multivariable-adjusted regression may be exaggerated by residual confounding from unobserved confounders such as socioeconomic position, smoking, and physical activity. As the AMV results were cross-sectional it is also possible that variation in metabolomic traits influences sleep patterns, and some of the multivariable regression results not verified in MR are due to reverse causality. In addition, we restricted the analyses to cohorts containing mostly European participants (one cohort contributing to AMV meta-analysis, HELIUS, included non-European participants, whereas all MR analyses were restricted to Europeans). This reduces the potential for population stratification to bias our MR analyses, but hampers generalization of our findings to other ancestry groups. The MR results may have been influenced by weak instrument bias, which, if present, would be expected to bias results towards the null. Sensitivity analyses exploring possible bias due to directional horizontal pleiotropy were mostly consistent with the main IVW findings, though MR-Egger estimates, were imprecise as expected with this method which is statistically less efficient than the main IVW method.

Taken together, our findings do not suggest widespread metabolic disruption caused by sleep traits. However, the evidence for possible effects of insomnia symptoms on glycoprotein acetyls and citrate, and longer total sleep duration on creatinine and isoleucine might explain some of the effects, found in MR analyses, of these sleep traits on cardiometabolic diseases. These warrant further investigation.

## Data Availability

RN and DAL act as guarantors for the paper.

## Sources of Funding

This research was funded by the British Heart Foundation (AA/18/7/34219), Diabetes UK (17/0005700) and the European Research Council (DevelopObese; 669545), which funds NJG’s salary. KHW is supported by the Elizabeth Blackwell Institute for Health Research, University of Bristol and the Wellcome Trust Institutional Strategic Support Fund (204813/Z/16/Z). NJG, MAL, RCR, KHW and DAL work in a unit that receives support from the University of Bristol and UK Medical Research Council (MC_UU_00011/6). AD is funded by the Wellcome Trust (206046/Z/17/Z). RM is supported by the President’s PhD scholarship from Imperial College, London. MAL is funded by a UK Medical Research Council PhD studentship (MR/R502340/1). DAL is a UK National Institute for Health Research Senior Investigator (NF-0616-10102). Sources of funding for individual studies used in this paper are listed in **Supplemental Material**.

None of the funders influenced data collection, analyses or interpretation, or the decision to publish this research. The views expressed are those of the authors and not necessarily those of the sponsors.

## Disclosures

DOMK is a part-time research consultant at Metabolon, Inc. DAL has received research support from national and international government and charitable funders, as well as Roche Diagnostics and Medtronic for research unrelated to that presented here. All others declare no conflicts of interest.

## Acknowledgements

We are most grateful to the participants of the different participating studies. In addition, we are grateful to the management team, research nurses, interviewers, research assistants and other staff who have taken part in gathering the data of the different studies.

*   Received August 27, 2020.
*   Revision received August 27, 2020.
*   Accepted September 1, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## Reference List

1.  1.Wu Y, Zhai L and Zhang D. Sleep duration and obesity among adults: a meta-analysis of prospective studies. Sleep Med. 2014;15:1456–62.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.sleep.2014.07.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25450058&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

2.  2.Xi B, He D, Zhang M, Xue J and Zhou D. Short sleep duration predicts risk of metabolic syndrome: a systematic review and meta-analysis. Sleep Med Rev. 2014;18:293–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.smrv.2013.06.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23890470&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

3.  3.Itani O, Jike M, Watanabe N and Kaneita Y. Short sleep duration and health outcomes: a systematic review, meta-analysis, and meta-regression. Sleep Med. 2017;32:246–256.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.sleep.2016.08.006&link_type=DOI) 

4.  4.Jike M, Itani O, Watanabe N, Buysse DJ and Kaneita Y. Long sleep duration and health outcomes: A systematic review, meta-analysis and meta-regression. Sleep Med Rev. 2018;39:25–36.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

5.  5.Cappuccio FP, Cooper D, D'Elia L, Strazzullo P and Miller MA. Sleep duration predicts cardiovascular outcomes: a systematic review and meta-analysis of prospective studies. Eur Heart J. 2011;32:1484–92.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehr007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21300732&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295679100012&link_type=ISI) 

6.  6.Ge L, Guyatt G, Tian J, Pan B, Chang Y, Chen Y, Li H, Zhang J, Li Y, Ling J and Yang K. Insomnia and risk of mortality from all-cause, cardiovascular disease, and cancer: Systematic review and meta-analysis of prospective cohort studies. Sleep Med Rev. 2019;48:101215.
    
    

7.  7.Sofi F, Cesari F, Casini A, Macchi C, Abbate R and Gensini GF. Insomnia and risk of cardiovascular disease: a meta-analysis. Eur J Prev Cardiol. 2014;21:57–64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/2047487312460020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22942213&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

8.  8.Zheng B, Yu C, Lv J, Guo Y, Bian Z, Zhou M, Yang L, Chen Y, Li X, Zou J, Ning F, Chen J, Chen Z, Li L and China Kadoorie Biobank Collaborative G. Insomnia symptoms and risk of cardiovascular diseases among 0.5 million adults: A 10-year cohort. Neurology. 2019;93:e2110–e2120.
    
    

9.  9.Knutson KL and von Schantz M. Associations between chronotype, morbidity and mortality in the UK Biobank cohort. Chronobiol Int. 2018;35:1045–1053.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/07420528.2018.1454458&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29642757&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

10. 10.Petrov ME, Kim Y, Lauderdale D, Lewis CE, Reis JP, Carnethon MR, Knutson K and Glasser SJ. Longitudinal associations between objective sleep and lipids: the CARDIA study. Sleep. 2013;36:1587–95.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24179290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

11. 11.Bos MM, Noordam R, van den Berg R, de Mutsert R, Rosendaal FR, Blauw GJ, Rensen PCN, Biermasz NR and van Heemst D. Associations of sleep duration and quality with serum and hepatic lipids: The Netherlands Epidemiology of Obesity Study. J Sleep Res. 2019;28:e12776.
    
    

12. 12.Kaneita Y, Uchiyama M, Yoshiike N and Ohida T. Associations of usual sleep duration with serum lipid and lipoprotein levels. Sleep. 2008;31:645–52.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18517035&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255451600009&link_type=ISI) 

13. 13.Noordam R, Bos MM, Wang H, Winkler TW, Bentley AR, Kilpelainen TO, de Vries PS, Sung YJ, Schwander K, Cade BE, Manning A, Aschard H, Brown MR, Chen H, Franceschini N, Musani SK, Richard M, Vojinovic D, Aslibekyan S, Bartz TM, de Las Fuentes L, Feitosa M, Horimoto AR, Ilkov M, Kho M, Kraja A, Li C, Lim E, Liu Y, Mook- Kanamori DO, Rankinen T, Tajuddin SM, van der Spek A, Wang Z, Marten J, Laville V, Alver M, Evangelou E, Graff ME, He M, Kuhnel B, Lyytikainen LP, Marques-Vidal P, Nolte IM, Palmer ND, Rauramaa R, Shu XO, Snieder H, Weiss S, Wen W, Yanek LR, Adolfo C, Ballantyne C, Bielak L, Biermasz NR, Boerwinkle E, Dimou N, Eiriksdottir G, Gao C, Gharib SA, Gottlieb DJ, Haba-Rubio J, Harris TB, Heikkinen S, Heinzer R, Hixson JE, Homuth G, Ikram MA, Komulainen P, Krieger JE, Lee J, Liu J, Lohman KK, Luik AI, Magi R, Martin LW, Meitinger T, Metspalu A, Milaneschi Y, Nalls MA, O'Connell J, Peters A, Peyser P, Raitakari OT, Reiner AP, Rensen PCN, Rice TK, Rich SS, Roenneberg T, Rotter JI, Schreiner PJ, Shikany J, Sidney SS, Sims M, Sitlani CM, Sofer T, Strauch K, Swertz MA, Taylor KD, Uitterlinden AG, van Duijn CM, Volzke H, Waldenberger M, Wallance RB, van Dijk KW, Yu C, Zonderman AB, Becker DM, Elliott P, Esko T, Gieger C, Grabe HJ, Lakka TA, Lehtimaki T, North KE, Penninx B, Vollenweider P, Wagenknecht LE, Wu T, Xiang YB, Zheng W, Arnett DK, Bouchard C, Evans MK, Gudnason V, Kardia S, Kelly TN, Kritchevsky SB, Loos RJF, Pereira AC, Province M, Psaty BM, Rotimi C, Zhu X, Amin N, Cupples LA, Fornage M, Fox EF, Guo X, Gauderman WJ, Rice K, Kooperberg C, Munroe PB, Liu CT, Morrison AC, Rao DC, van Heemst D and Redline S. Multi-ancestry sleep-by- SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration. Nat Commun. 2019;10:5121.
    
    

14. 14.Lawlor DA. Commentary: Two-sample Mendelian randomization: opportunities and challenges. Int J Epidemiol. 2016;45:908–15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw127&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27427429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

15. 15.Lawlor DA, Harbord RM, Sterne JA, Timpson N and Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27:1133–63.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.3034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17886233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

16. 16.Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J, Langdon R, Tan VY, Yarmolinsky J, Shihab HA, Timpson NJ, Evans DM, Relton C, Martin RM, Davey Smith G, Gaunt TR and Haycock PC. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7.
    
    

17. 17.Lane JM, Jones SE, Dashti HS, Wood AR, Aragam KG, van Hees VT, Strand LB, Winsvold BS, Wang H, Bowden J, Song Y, Patel K, Anderson SG, Beaumont RN, Bechtold DA, Cade BE, Haas M, Kathiresan S, Little MA, Luik AI, Loudon AS, Purcell S, Richmond RC, Scheer F, Schormair B, Tyrrell J, Winkelman JW, Winkelmann J, Sleep HAI, Hveem K, Zhao C, Nielsen JB, Willer CJ, Redline S, Spiegelhalder K, Kyle SD, Ray DW, Zwart JA, Brumpton B, Frayling TM, Lawlor DA, Rutter MK, Weedon MN and Saxena R. Biological and clinical insights from genetics of insomnia symptoms. Nat Genet. 2019;51:387–393.
    
    

18. 18.Daghlas I, Dashti HS, Lane J, Aragam KG, Rutter MK, Saxena R and Vetter C. Sleep Duration and Myocardial Infarction. J Am Coll Cardiol. 2019;74:1304–1314.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjEwOiI3NC8xMC8xMzA0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

19. 19.Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, de Leeuw CA, Benjamins JS, Munoz-Manchado AB, Nagel M, Savage JE, Tiemeier H, White T, and Me Research T, Tung JY, Hinds DA, Vacic V, Wang X, Sullivan PF, van der Sluis S, Polderman TJC, Smit AB, Hjerling-Leffler J, Van Someren EJW and Posthuma D. Genome- wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways. Nat Genet. 2019;51:394–403.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0333-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30804565&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

20. 20.Dashti HS, Jones SE, Wood AR, Lane JM, van Hees VT, Wang H, Rhodes JA, Song Y, Patel K, Anderson SG, Beaumont RN, Bechtold DA, Bowden J, Cade BE, Garaulet M, Kyle SD, Little MA, Loudon AS, Luik AI, Scheer F, Spiegelhalder K, Tyrrell J, Gottlieb DJ, Tiemeier H, Ray DW, Purcell SM, Frayling TM, Redline S, Lawlor DA, Rutter MK, Weedon MN and Saxena R. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat Commun. 2019;10:1100.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-08917-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30846698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

21. 21.Jones SE, Lane JM, Wood AR, van Hees VT, Tyrrell J, Beaumont RN, Jeffries AR, Dashti HS, Hillsdon M, Ruth KS, Tuke MA, Yaghootkar H, Sharp SA, Jie Y, Thompson WD, Harrison JW, Dawes A, Byrne EM, Tiemeier H, Allebrandt KV, Bowden J, Ray DW, Freathy RM, Murray A, Mazzotti DR, Gehrman PR, Lawlor DA, Frayling TM, Rutter MK, Hinds DA, Saxena R and Weedon MN. Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat Commun. 2019;10:343.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-08259-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30696823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

22. 22.Lawlor DA, Tilling K and Davey Smith G. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016;45:1866–1886.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28108528&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

23. 23.Wijsman CA, Westendorp RG, Verhagen EA, Catt M, Slagboom PE, de Craen AJ, Broekhuizen K, van Mechelen W, van Heemst D, van der Ouderaa F and Mooijaart SP. Effects of a web-based intervention on physical activity and metabolism in older adults: randomized controlled trial. J Med Internet Res. 2013;15:e233.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2196/jmir.2843&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24195965&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

24. 24.Lumey LH, Stein AD, Kahn HS, van der Pal-de Bruin KM, Blauw GJ, Zybert PA and Susser ES. Cohort profile: the Dutch Hunger Winter families study. Int J Epidemiol. 2007;36:1196–204.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dym126&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17591638&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252247400008&link_type=ISI) 

25. 25.Snijder MB, Galenkamp H, Prins M, Derks EM, Peters RJG, Zwinderman AH and Stronks K. Cohort profile: the Healthy Life in an Urban Setting (HELIUS) study in Amsterdam, The Netherlands. BMJ Open. 2017;7:e017873.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiNy8xMi9lMDE3ODczIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

26. 26.van Oosterhout WP, Weller CM, Stam AH, Bakels F, Stijnen T, Ferrari MD and Terwindt GM. Validation of the web-based LUMINA questionnaire for recruiting large cohorts of migraineurs. Cephalalgia. 2011;31:1359–67.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0333102411418846&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21914734&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295887500005&link_type=ISI) 

27. 27.Penninx BW, Beekman AT, Smit JH, Zitman FG, Nolen WA, Spinhoven P, Cuijpers P, De Jong PJ, Van Marwijk HW, Assendelft WJ, Van Der Meer K, Verhaak P, Wensing M, De Graaf R, Hoogendijk WJ, Ormel J, Van Dyck R and Consortium NR. The Netherlands Study of Depression and Anxiety (NESDA): rationale, objectives and methods. Int J Methods Psychiatr Res. 2008;17:121–40.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mpr.256&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18763692&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259296400001&link_type=ISI) 

28. 28.Willemsen G, de Geus EJ, Bartels M, van Beijsterveldt CE, Brooks AI, Estourgie-van Burk GF, Fugman DA, Hoekstra C, Hottenga JJ, Kluft K, Meijer P, Montgomery GW, Rizzu P, Sondervan D, Smit AB, Spijker S, Suchiman HE, Tischfield JA, Lehner T, Slagboom PE and Boomsma DI. The Netherlands Twin Register biobank: a resource for genetic epidemiological studies. Twin Res Hum Genet. 2010;13:231–45.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1375/twin.13.3.231&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20477721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

29. 29.de Mutsert R, den Heijer M, Rabelink TJ, Smit JW, Romijn JA, Jukema JW, de Roos A, Cobbaert CM, Kloppenburg M, le Cessie S, Middeldorp S and Rosendaal FR. The Netherlands Epidemiology of Obesity (NEO) study: study design and data collection. Eur J Epidemiol. 2013;28:513–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10654-013-9801-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23576214&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

30. 30.Ikram MA, Brusselle GGO, Murad SD, van Duijn CM, Franco OH, Goedegebure A, Klaver CCW, Nijsten TEC, Peeters RP, Stricker BH, Tiemeier H, Uitterlinden AG, Vernooij MW and Hofman A. The Rotterdam Study: 2018 update on objectives, design and main results. Eur J Epidemiol. 2017;32:807–850.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10654-017-0321-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

31. 31.Kettunen J, Demirkan A, Wurtz P, Draisma HH, Haller T, Rawal R, Vaarhorst A, Kangas AJ, Lyytikainen LP, Pirinen M, Pool R, Sarin AP, Soininen P, Tukiainen T, Wang Q, Tiainen M, Tynkkynen T, Amin N, Zeller T, Beekman M, Deelen J, van Dijk KW, Esko T, Hottenga JJ, van Leeuwen EM, Lehtimaki T, Mihailov E, Rose RJ, de Craen AJ, Gieger C, Kahonen M, Perola M, Blankenberg S, Savolainen MJ, Verhoeven A, Viikari J, Willemsen G, Boomsma DI, van Duijn CM, Eriksson J, Jula A, Jarvelin MR, Kaprio J, Metspalu A, Raitakari O, Salomaa V, Slagboom PE, Waldenberger M, Ripatti S and Ala-Korpela M. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms11122&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27005778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

32. 32.Karpe F, Vasan SK, Humphreys SM, Miller J, Cheeseman J, Dennis AL and Neville MJ. Cohort Profile: The Oxford Biobank. Int J Epidemiol. 2018;47:21–21g.
    
    

33. 33.Shepherd J, Blauw GJ, Murphy MB, Cobbe SM, Bollen EL, Buckley BM, Ford I, Jukema JW, Hyland M, Gaw A, Lagaay AM, Perry IJ, Macfarlane PW, Meinders AE, Sweeney BJ, Packard CJ, Westendorp RG, Twomey C and Stott DJ. The design of a prospective study of Pravastatin in the Elderly at Risk (PROSPER). PROSPER Study Group. PROspective Study of Pravastatin in the Elderly at Risk. Am J Cardiol. 1999;84:1192–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0002-9149(99)00533-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10569329&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000083534700009&link_type=ISI) 

34. 34.Soininen P, Kangas AJ, Würtz P, Suna T and Ala-Korpela M. Quantitative Serum Nuclear Magnetic Resonance Metabolomics in Cardiovascular Epidemiology and Genetics. Circulation: Cardiovascular Genetics. 2015;8:192–206.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiY2lyY2N2ZyI7czo1OiJyZXNpZCI7czo3OiI4LzEvMTkyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35.Sliz E, Kettunen J, Holmes MV, Williams CO, Boachie C, Wang Q, Mannikko M, Sebert S, Walters R, Lin K, Millwood LY, Clarke R, Li L, Rankin N, Welsh P, Delles C, Jukema JW, Trompet S, Ford I, Perola M, Salomaa V, Jarvelin MR, Chen Z, Lawlor DA, Ala-Korpela M, Danesh J, Smith GD, Sattar N, Butterworth A and Wurtz P. Metabolomic Consequences of Genetic Inhibition of PCSK9 Compared With Statin Treatment. Circulation. 2018;138:2499–2512.
    
    

36. 36.Wurtz P, Havulinna AS, Soininen P, Tynkkynen T, Prieto-Merino D, Tillin T, Ghorbani A, Artati A, Wang Q, Tiainen M, Kangas AJ, Kettunen J, Kaikkonen J, Mikkila V, Jula A, Kahonen M, Lehtimaki T, Lawlor DA, Gaunt TR, Hughes AD, Sattar N, Illig T, Adamski J, Wang TJ, Perola M, Ripatti S, Vasan RS, Raitakari OT, Gerszten RE, Casas JP, Chaturvedi N, Ala-Korpela M and Salomaa V. Metabolite profiling and cardiovascular event risk: a prospective study of 3 population-based cohorts. Circulation. 2015;131:774–85.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjEzMS85Lzc3NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzAxLzIwMjAuMDguMjcuMjAxNzM1MTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

37. 37.Team RDC. R: A Language and Environment for Statistical Computing. 2008.
    
    

38. 38.Noordam R, Young WJ, Salman R, Kanters JK, van den Berg ME, van Heemst D, Lin HJ, Barreto SM, Biggs ML, Biino G, Catamo E, Concas MP, Ding J, Evans DS, Foco L, Grarup N, Lyytikainen LP, Mangino M, Mei H, van der Most PJ, Muller-Nurasyid M, Nelson CP, Qian Y, Repetto L, Said MA, Shah N, Schramm K, Vidigal PG, Weiss S, Yao J, Zilhao NR, Brody JA, Braund PS, Brumat M, Campana E, Christofidou P, Caulfield MJ, De Grandi A, Dominiczak AF, Doney ASF, Eiriksdottir G, Ellervik C, Giatti L, Gogele M, Graff C, Guo X, van der Harst P, Joshi PK, Kahonen M, Kestenbaum B, Lima-Costa MF, Linneberg A, Maan AC, Meitinger T, Padmanabhan S, Pattaro C, Peters A, Petersmann A, Sever P, Sinner MF, Shen X, Stanton A, Strauch K, Soliman EZ, Tarasov KV, Taylor KD, Thio CHL, Uitterlinden AG, Vaccargiu S, Waldenberger M, Robino A, Correa A, Cucca F, Cummings SR, Dorr M, Girotto G, Gudnason V, Hansen T, Heckbert SR, Juhl CR, Kaab S, Lehtimaki T, Liu Y, Lotufo PA, Palmer CNA, Pirastu M, Pramstaller PP, Ribeiro ALP, Rotter JI, Samani NJ, Snieder H, Spector TD, Stricker BH, Verweij N, Wilson JF, Wilson JG, Jukema JW, Tinker A, Newton-Cheh CH, Sotoodehnia N, Mook-Kanamori DO, Munroe PB and Warren HR. Effects of Calcium, Magnesium, and Potassium Concentrations on Ventricular Repolarization in Unselected Individuals. J Am Coll Cardiol. 2019;73:3118–3131.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjEwOiI3My8yNC8zMTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

39. 39.Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N and Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36:1783–1802.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.7221&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28114746&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

40. 40.Bowden J, Davey Smith G and Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyv080&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26050253&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

41. 41.Bowden J, Davey Smith G, Haycock PC and Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40:304–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21965&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27061298&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

42. 42.Santos Ferreira DL, Williams DM, Kangas AJ, Soininen P, Ala-Korpela M, Smith GD, Jarvelin MR and Lawlor DA. Association of pre-pregnancy body mass index with offspring metabolic profile: Analyses of 3 European prospective birth cohorts. PLoS Med. 2017;14:e1002376.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002376&link_type=DOI) 

43. 43.Davies SK, Ang JE, Revell VL, Holmes B, Mann A, Robertson FP, Cui N, Middleton B, Ackermann K, Kayser M, Thumser AE, Raynaud FI and Skene DJ. Effect of sleep deprivation on the human metabolome. Proc Natl Acad Sci U S A. 2014;111:10761–6.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTExLzI5LzEwNzYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

44. 44.Emerging Risk Factors Collaboration, Kaptoge S, Di Angelantonio E, Lowe G, Pepys MB, Thompson SG, Collins R and Danesh J. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis. Lancet. 2010;375:132–40.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(09)61717-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20031199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273664500030&link_type=ISI) 

45. 45.C. Reactive Protein Coronary Heart Disease Genetics Collaboration, Wensley F, Gao P, Burgess S, Kaptoge S, Di Angelantonio E, Shah T, Engert JC, Clarke R, Davey-Smith G, Nordestgaard BG, Saleheen D, Samani NJ, Sandhu M, Anand S, Pepys MB, Smeeth L, Whittaker J, Casas JP, Thompson SG, Hingorani AD and Danesh J. Association between C reactive protein and coronary heart disease: mendelian randomisation analysis based on individual participant data. BMJ. 2011;342:d548.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE2OiIzNDIvZmViMTVfMi9kNTQ4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDEvMjAyMC4wOC4yNy4yMDE3MzUxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

46. 46.Fuertes-Martin R, Correig X, Vallve JC and Amigo N. Title: Human Serum/Plasma Glycoprotein Analysis by (1)H-NMR, an Emerging Method of Inflammatory Assessment. J Clin Med. 2020;9.
    
    

47. 47.Connelly MA, Otvos JD, Shalaurova I, Playford MP and Mehta NN. GlycA, a novel biomarker of systemic inflammation and cardiovascular disease risk. J Transl Med. 2017;15:219.
    
    

48. 48.Costello LC and Franklin RB. Plasma Citrate Homeostasis: How It Is Regulated; And Its Physiological and Clinical Implications. An Important, But Neglected, Relationship in Medicine. HSOA J Hum Endocrinol. 2016;1.
    
    

49. 49.Kandeger A, Egilmez U, Sayin AA and Selvi Y. The relationship between night eating symptoms and disordered eating attitudes via insomnia and chronotype differences. Psychiatry Res. 2018;268:354–357.
    
    

50. 50.Voight BF, Peloso GM, Orho-Melander M, Frikke-Schmidt R, Barbalic M, Jensen MK, Hindy G, Holm H, Ding EL, Johnson T, Schunkert H, Samani NJ, Clarke R, Hopewell JC, Thompson JF, Li M, Thorleifsson G, Newton-Cheh C, Musunuru K, Pirruccello JP, Saleheen D, Chen L, Stewart A, Schillert A, Thorsteinsdottir U, Thorgeirsson G, Anand S, Engert JC, Morgan T, Spertus J, Stoll M, Berger K, Martinelli N, Girelli D, McKeown PP, Patterson CC, Epstein SE, Devaney J, Burnett MS, Mooser V, Ripatti S, Surakka I, Nieminen MS, Sinisalo J, Lokki ML, Perola M, Havulinna A, de Faire U, Gigante B, Ingelsson E, Zeller T, Wild P, de Bakker PI, Klungel OH, Maitland-van der Zee AH, Peters BJ, de Boer A, Grobbee DE, Kamphuisen PW, Deneer VH, Elbers CC, Onland-Moret NC, Hofker MH, Wijmenga C, Verschuren WM, Boer JM, van der Schouw YT, Rasheed A, Frossard P, Demissie S, Willer C, Do R, Ordovas JM, Abecasis GR, Boehnke M, Mohlke KL, Daly MJ, Guiducci C, Burtt NP, Surti A, Gonzalez E, Purcell S, Gabriel S, Marrugat J, Peden J, Erdmann J, Diemert P, Willenborg C, Konig IR, Fischer M, Hengstenberg C, Ziegler A, Buysschaert I, Lambrechts D, Van de Werf F, Fox KA, El Mokhtari NE, Rubin D, Schrezenmeir J, Schreiber S, Schafer A, Danesh J, Blankenberg S, Roberts R, McPherson R, Watkins H, Hall AS, Overvad K, Rimm E, Boerwinkle E, Tybjaerg-Hansen A, Cupples LA, Reilly MP, Melander O, Mannucci PM, Ardissino D, Siscovick D, Elosua R, Stefansson K, O'Donnell CJ, Salomaa V, Rader DJ, Peltonen L, Schwartz SM, Altshuler D and Kathiresan S. Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study. Lancet. 2012;380:572–80.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(12)60312-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22607825&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307511400030&link_type=ISI) 

51. 51.Schwartz GG, Olsson AG, Abt M, Ballantyne CM, Barter PJ, Brumm J, Chaitman BR, Holme IM, Kallend D, Leiter LA, Leitersdorf E, McMurray JJ, Mundl H, Nicholls SJ, Shah PK, Tardif JC, Wright RS and dal OI. Effects of dalcetrapib in patients with a recent acute coronary syndrome. N Engl J Med. 2012;367:2089–99.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1206797&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23126252&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000311624400006&link_type=ISI) 

52. 52.Barter PJ, Caulfield M, Eriksson M, Grundy SM, Kastelein JJ, Komajda M, Lopez-Sendon J, Mosca L, Tardif JC, Waters DD, Shear CL, Revkin JH, Buhr KA, Fisher MR, Tall AR, Brewer B and Investigators I. Effects of torcetrapib in patients at high risk for coronary events. N Engl J Med. 2007;357:2109–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0706628&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17984165&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251076500003&link_type=ISI) 

53. 53.Ding M and Rexrode KM. A Review of Lipidomics of Cardiovascular Disease Highlights the Importance of Isolating Lipoproteins. Metabolites. 2020;10.
    
    

54. 54.Liu HM, Hu Q, Zhang Q, Su GY, Xiao HM, Li BY, Shen WD, Qiu X, Lv WQ and Deng HW. Causal Effects of Genetically Predicted Cardiovascular Risk Factors on Chronic Kidney Disease: A Two-Sample Mendelian Randomization Study. Front Genet. 2019;10:415.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31130989&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

55. 55.Di Angelantonio E, Danesh J, Eiriksdottir G and Gudnason V. Renal function and risk of coronary heart disease in general populations: new prospective study and systematic review. PLoS Med. 2007;4:e270.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.0040270&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17803353&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 

56. 56.Lotta LA, Scott RA, Sharp SJ, Burgess S, Luan J, Tillin T, Schmidt AF, Imamura F, Stewart ID, Perry JR, Marney L, Koulman A, Karoly ED, Forouhi NG, Sjogren RJ, Naslund E, Zierath JR, Krook A, Savage DB, Griffin JL, Chaturvedi N, Hingorani AD, Khaw KT, Barroso I, McCarthy MI, O'Rahilly S, Wareham NJ and Langenberg C. Genetic Predisposition to an Impaired Metabolism of the Branched-Chain Amino Acids and Risk of Type 2 Diabetes: A Mendelian Randomisation Analysis. PLoS Med. 2016;13:el002179.
    
    

57. 57.VanderWeele TJ, Tchetgen Tchetgen EJ, Cornelis M and Kraft P. Methodological challenges in mendelian randomization. Epidemiology. 2014;25:427–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000081&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24681576&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F01%2F2020.08.27.20173518.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000337313300013&link_type=ISI)