Studies populations
Lipid GWAS bottom line analytics was in fact obtained from the fresh new Mil Experienced System (MVP) (doing 215,551 Western european some one) and Globally Lipids Family genes Consortium (GLGC) (to 188,577 genotyped some body) . Since more exposures inside multivariable MR analyses, we put Bmi conclusion analytics out of a good meta-investigation out-of GWASs inside the doing 795,640 anybody and you can many years at menarche realization statistics from a great meta-research out-of GWASs inside the around 329,345 females regarding Eu origins [17,23]. This new MVP gotten ethical and read process acceptance in the Veteran Affair Central Institutional Review Panel in accordance with the principles detail by detail on the Statement off Helsinki, and you may written concur was taken from all professionals. Towards the Willer and you will associates and you will BCAC data establishes, i send your reader towards the number one GWAS manuscripts as well as their supplementary issue getting all about agree standards for each and every of their particular cohorts. Addiitional information in these cohorts have the fresh new S1 Text message.
Lipid meta-research
I did a fixed-consequences meta-data between for every lipid attribute (Complete cholesterol levels [TC], LDL, HDL, and you can triglycerides [TGs]) inside the GLGC together with related lipid characteristic on the MVP cohort [twelve,22] making use of the default configurations inside the PLINK . There can be certain genomic rising cost of living within these meta-analysis association statistics, however, linkage disequilibrium (LD)-score regression intercepts demonstrate that that it rising prices is in higher part due to polygenicity rather than population stratification (S1 Fig).
MR analyses
MR analyses were performed using the TwoSampleMR R package version 0.4.13 ( . For all analyses, we used a 2-sample MR framework, with exposure(s) (lipids, BMI, age at menarche) and outcome (BC) genetic associations from separate cohorts. Unless otherwise noted, MR results reported in this manuscript used inverse-variance weighting assuming a multiplicative random effects model. For single-trait MR analyses, we additionally employed Egger regression , weighted median , and mode-based estimates. SNPs associated with each lipid trait were filtered for genome-wide significance (P < 5 ? 10 ?8 ) from the MVP lipid study , and then we removed SNPs in LD (r 2 < 0.001 in UK10K consortium) in order to obtain independent variants. All genetic variants were harmonized using the TwoSampleMR harmonization function with default parameters. Each of these independent, genome-wide significant SNPs was termed a genetic instrument. We estimated that these single-trait MR genetic instruments had 80% power to reject the null hypothesis, with a 1% error rate, for the following odds ratio (OR) increases in BC risk due to a standard deviation increase in lipid levels: HDL, 1.057; LDL, 1.058; TGs, 1.055; TC, 1.060 [30,31]. We tested for directional pleiotropy using the MR-Egger regression test . To reduce heterogeneity in our genetic instruments for single-trait MR, we employed a pruning procedure (S1 Text). Genetic instruments used in single-trait MR are listed in S1 Table. For multivariable MR experiments [32,33], we generated genetic instruments by first filtering the genotyped variants for those present across all data sets. For each trait and data set combination (Yengo and colleagues for BMI; Day and colleagues for age at menarche ; MVP and GLGC for HDL, LDL, and TGs), we then filtered for genome-wide significance (P < 5 ? 10 ?8 ) and for linkage disequilibrium (r 2 < 0.001 in UK10K consortium) . We performed tests for instrument strength and validity , and each multivariable MR experiment had sufficient instrument strength. We removed variants driving heterogeneity in the ratio of outcome/exposure effects causing instrument invalidity (S1 Text). Genetic instruments used in multivariable MR are listed in S2 Table. Because the MR methods and tests we employed are highly correlated, we did not apply a multiple testing correction to the reported P-values.
