Cross-trait multivariate GWAS confirms health implications of pubertal timing

Pubertal timing is highly variable and is associated with long-term health outcomes. Phenotypes associated with pubertal timing include age at menarche, age at voice break, age at first facial hair and growth spurt, and pubertal timing seems to have a shared genetic architecture between the sexes. However, puberty phenotypes have primarily been assessed separately, failing to account for shared genetics, which limits the reliability of the purported health implications. Here, we model the common genetic architecture for puberty timing using a multivariate GWAS, with an effective population of 514,750 European participants. We find 266 independent variants in 197 loci, including 18 novel variants. Transcriptomic, proteome imputation and fine-mapping analyses reveal genes causal for pubertal timing, including KDM4C, LEPR, CCNC, ACP1, and PCSK1. Linkage disequilibrium score regression and Mendelian randomisation analysis establish causal associations between earlier puberty and both accelerated ageing and the risk of developing cardiovascular disease and osteoporosis. We find that alanine aminotransferase, glycated haemoglobin, high-density lipoprotein cholesterol and Parabacteroides levels are mediators of these relationships, and establish that controlling oily fish and retinol intake may be beneficial for promoting healthy pubertal development.