Genome wide analysis for mouth ulcers identifies associations at immune regulatory loci

Mouth ulcers are the most common ulcerative condition and encompass several clinical diagnoses, including recurrent aphthous stomatitis (RAS). Despite previous evidence for heritability, it is not clear which specific genetic loci are implicated in RAS. In this genome-wide association study (n = 461,106) heritability is estimated at 8.2% (95% CI: 6.4%, 9.9%). This study finds 97 variants which alter the odds of developing non-specific mouth ulcers and replicate these in an independent cohort (n = 355,744) (lead variant after meta-analysis: rs76830965, near IL12A, OR 0.72 (95% CI: 0.71, 0.73); P = 4.4e−483). Additional effect estimates from three independent cohorts with more specific phenotyping and specific study characteristics support many of these findings. In silico functional analyses provide evidence for a role of T cell regulation in the aetiology of mouth ulcers. These results provide novel insight into the pathogenesis of a common, important condition.

Oral ulcerations are sores of the mucous membrane of the mouth and highly prevalent in the population. Here, in a genome-wide association study, the authors identify 97 loci associated with mouth ulcers highlighting genes involved in T cell-mediated immunity and T_H1 responses.

Coronary artery disease, genetic risk and the metabolome in young individuals

Background: Genome-wide association studies have identified genetic variants associated with coronary artery disease (CAD) in adults - the leading cause of death worldwide. It often occurs later in life, but variants may impact CAD-relevant phenotypes early and throughout the life-course. Cohorts with longitudinal and genetic data on thousands of individuals are letting us explore the antecedents of this adult disease. Methods: 148 metabolites, with a focus on the lipidome, measured using nuclear magnetic resonance ( 1H-NMR) spectroscopy, and genotype data were available from 5,907 individuals at ages 7, 15, and 17 years from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Linear regression was used to assess the association between the metabolites and an adult-derived genetic risk score (GRS) of CAD comprising 146 variants. Individual variant-metabolite associations were also examined. Results: The CAD-GRS associated with 118 of 148 metabolites (false discovery rate [FDR] < 0.05), the strongest associations being with low-density lipoprotein (LDL) and atherogenic non-LDL subgroups. Nine of 146 variants in the GRS associated with one or more metabolites (FDR < 0.05). Seven of these are within lipid loci: rs11591147 PCSK9, rs12149545 HERPUD1-CETP, rs17091891 LPL, rs515135 APOB, rs602633 CELSR2-PSRC1, rs651821 APOA5, rs7412 APOE-APOC1. All associated with metabolites in the LDL or atherogenic non-LDL subgroups or both including aggregate cholesterol measures. The other two variants identified were rs112635299 SERPINA1 and rs2519093 ABO. Conclusions: Genetic variants that influence CAD risk in adults are associated with large perturbations in metabolite levels in individuals as young as seven. The variants identified are mostly within lipid-related loci and the metabolites they associated with are primarily linked to lipoproteins. Along with further research, this knowledge could allow for preventative measures, such as increased monitoring of at-risk individuals and perhaps treatment earlier in life, to be taken years before any symptoms of the disease arise.

Apparent latent structure within the UK Biobank sample has implications for epidemiological analysis

Large studies use genotype data to discover genetic contributions to complex traits and infer relationships between those traits. Co-incident geographical variation in genotypes and health traits can bias these analyses. Here we show that single genetic variants and genetic scores composed of multiple variants are associated with birth location within UK Biobank and that geographic structure in genotype data cannot be accounted for using routine adjustment for study centre and principal components derived from genotype data. We find that major health outcomes appear geographically structured and that coincident structure in health outcomes and genotype data can yield biased associations. Understanding and accounting for this phenomenon will be important when making inference from genotype data in large studies.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Investigating the impact of cigarette smoking behaviours on DNA methylation patterns in adolescence

Smoking usually begins in adolescence, and early onset of smoking has been linked to increased risk of later life disease. There is a need to better understand the biological impact of cigarette smoking behaviours in adolescence. DNA methylation profiles related to smoking behaviours and cessation in adulthood have been previously identified, but alterations arising from smoking initiation have not been thoroughly investigated. We aimed to investigate DNA methylation in the Avon Longitudinal Study of Parents and Children in relation to (1) different smoking measures, (2) time since smoking initiation and frequency of smoke exposure and (3) latent classes of smoking behaviour. Using 2620 CpG sites previously associated with cigarette smoking, we investigated DNA methylation change in relation to own smoking measures, smoke exposure duration and frequency, and using longitudinal latent class analysis of different smoking behaviour patterns in 968 adolescents. Eleven CpG sites located in seven gene regions were differentially methylated in relation to smoking in adolescence. While only AHRR (cg05575921) showed a robust pattern of methylation in relation to weekly smoking, several CpGs showed differences in methylation among individuals who had tried smoking compared with non-smokers. In relation to smoke exposure duration and frequency, cg05575921 showed a strong dose-response relationship, while there was evidence for more immediate methylation change at other sites. Our findings illustrate the impact of cigarette smoking behaviours on DNA methylation at some smoking-responsive CpG sites, even among individuals with a short smoking history.

The ‘ALSPAC in London’ dataset: adiposity, cardiometabolic risk profiles, and the emerging arterial phenotype in young adulthood

Rising rates of adiposity in the young pose one of the greatest threats to future population burden of cardiovascular disease. Understanding the contribution of genetic and early-life influences to adiposity profiles in young adulthood – when the first signs of subclinical cardiovascular disease commonly appear – are vital if effective lifetime prevention strategies are to be developed. This data note documents the extensive range of genotypic and phenotypic data available from a London-based sub-study of the long-running Avon Longitudinal Study of Parents and Children (ALSPAC)—the ‘ALSPAC in London’ Study—in which extensive adipose and cardiovascular phenotyping was carried out in participants recruited based on a genetic predisposition to obesity.

Associations of Body Mass and Fat Indexes With Cardiometabolic Traits

BACKGROUND

Body mass index (BMI) is criticized for not distinguishing fat from lean mass and ignoring fat distribution, leaving its ability to detect health effects unclear.

OBJECTIVES

The aim of this study was to compare BMI with total and regional fat indexes from dual-energy x-ray absorptiometry in their associations with cardiometabolic traits. Duration of exposure to and change in each index across adolescence were examined in relation to detailed traits in young adulthood.

METHODS

BMI was examined alongside total, trunk, arm, and leg fat indexes (each in kilograms per square meter) from dual-energy x-ray absorptiometry at ages 10 and 18 years in relation to 230 traits from targeted metabolomics at age 18 years in 2,840 offspring from the Avon Longitudinal Study of Parents and Children.

RESULTS

Higher total fat mass index and BMI at age 10 years were similarly associated with cardiometabolic traits at age 18 years, including higher systolic and diastolic blood pressure, higher very low-density lipoprotein and low-density lipoprotein cholesterol, lower high-density lipoprotein cholesterol, higher triglycerides, and higher insulin and glycoprotein acetyls. Associations were stronger for both indexes measured at age 18 years and for gains in each index from age 10 to 18 years (e.g., 0.45 SDs [95% confidence interval: 0.38 to 0.53] in glycoprotein acetyls per SD unit gain in fat mass index vs. 0.38 SDs [95% confidence interval: 0.27 to 0.48] per SD unit gain in BMI). Associations resembled those for trunk fat index. Higher lean mass index was weakly associated with traits and was not protective against higher fat mass index.

CONCLUSIONS

The results of this study support abdominal fatness as a primary driver of cardiometabolic dysfunction and BMI as a useful tool for detecting its effects.

Reassessing the Association between Circulating Vitamin D and IGFBP-3: Observational and Mendelian Randomization Estimates from Independent Sources

BACKGROUND

Circulating insulin-like growth factor binding protein 3 (IGFBP-3) has been associated with prostate cancer. Preclinical studies found that vitamin D regulates IGFBP-3 expression, although evidence from epidemiologic studies is conflicting.

METHODS

Mendelian randomization analyses (MR) were conducted to reassess associations between IGFBP-3 and prostate cancer risk and advanced prostate cancer using summary statistics from the PRACTICAL consortium (44,825 cases; 27,904 controls). Observational and MR analyses were conducted to assess the relationship between inactive vitamin D [25(OH)D] and IGFBP-3 using data from the ProtecT study (1,366 cases;1,071 controls) and summary statistics from the CHARGE consortium (n = 18,995).

RESULTS

The OR for prostate cancer per SD unit increase in circulating IGFBP-3 was 1.14 [95% confidence interval (CI), 1.02-1.28]. The OR for advanced prostate cancer per SD unit increase in IGFBP-3 was 1.22 (95% CI, 1.07-1.40). Observationally, a SD increase in 25(OH)D was associated with a 0.1SD (95% CI, 0.05-0.14) increase in IGFBP-3. MR analyses found little evidence for a causal relationship between circulating 25(OH)D and IGFBP-3 in the circulation.

CONCLUSIONS

This study provided confirmatory evidence that IGFBP-3 is a risk factor for prostate cancer risk and progression. Observationally, there was evidence that 25(OH)D is associated with IGFBP-3, but MR analyses suggested that these findings were unlikely to be causal. Findings may be limited by the nature of instrumentation of 25(OH)D and IGFBP-3 and the utility of circulating measures. 25(OH)D appears unlikely to be causally related to IGFBP-3 in the circulation, however, our findings do not preclude causal associations at the tissue level.

IMPACT

IGFBP-3 is a prostate cancer risk factor but 25(OH)D are unlikely to be causally related to IGFBP-3 in the circulation.

Tooth loss is a complex measure of oral disease: Determinants and methodological considerations

OBJECTIVES

Counts of missing teeth or measures of incident tooth loss are gaining attention as a simple way to measure dental status in large population studies. We explore the meaning of these metrics and how missing teeth might influence other measures of dental status.

METHODS

An observational study was performed in 2 contrasting adult populations. In total, 62 522 adult participants were available with clinically assessed caries and periodontal indices from the Swedish arm of the Gene-Lifestyle Interactions and Dental Endpoints Study (GLIDE) and the Korea National Health and Nutrition Examination Survey (KNHANES) in the Republic of Korea. Longitudinal measures of tooth loss were available for 28 244 participants in GLIDE with median follow-up of 10.6 years.

RESULTS

In longitudinal analysis, hazard for tooth loss was associated with baseline dental status (previous tooth loss, periodontal status and caries status) and socio-demographic variables (age, smoking status and highest educational level). Analysis of cross-sectional data suggested that indices of caries exposure were not independent of periodontal status. The strength and direction of association varied between groups, even for measures specifically intended to avoid measuring tooth loss. Individuals with impaired periodontal health (community periodontal index [CPI] 3 or higher in any sextant) had higher standardized decayed and filled surfaces (DFS; number of DFS divided by total number of tooth surfaces) in GLIDE (incidence risk ratio [IRR] 1.05 [95% CI: 1.04, 1.07], but lower standardized DFS in KNHANES (IRR: 0.95 [0.92, 0.98]) than individuals with better periodontal health (CPI <3 in all sextants).

CONCLUSIONS

Incident tooth loss is a complex measure of dental disease, with multiple determinants. The relative importance of dental caries and periodontal disease as drivers of tooth loss differs between age groups. Measures of dental caries exposure are associated with periodontal status in the studied populations, and these associations can be population-specific. Consideration of the study-specific properties of these metrics may be required for valid inference in large population studies.

Assessing the causal role of body mass index on cardiovascular health in young adults: Mendelian randomization and recall-by-genotype analyses

Background

Body mass index (BMI) has been suggested to be causally related to cardiovascular health in mid-to-late life, but this has not been explored systematically at younger ages - nor with detailed cardiovascular phenotyping. Recall-by-Genotype (RbG) is an approach that enables the collection of precise phenotypic measures in smaller studies, whilst maintaining statistical power and ability for causal inference.

Methods

In this study, we used a combination of conventional multivariable regression analysis, Mendelian randomization (MR) and sub-sample RbG methodologies to estimate the causal effect of BMI on gross-level and detailed cardiovascular health in healthy participants from the Avon Longitudinal Study of Parents and Children at age 17 (N=1420-3108 for different outcomes) and an independent sample from the same cohort (for RbG) study at age 21 (N=386-418).

Results

In both MR and RbG analyses, results suggested that higher BMI causes higher blood pressure (BP) and left ventricular mass index (LVMI) in young adults (e.g., difference in LVMI per kg/m2 using MR: 1.07g/m2.7; 95% CI: 0.62, 1.52; P=3.87x10-06 and per 3.58kg/m2 using RbG: 1.65g/m2.7 95% CI: 0.83, 2.47; P=0.0001). Additionally, RbG results suggested a causal role of higher BMI on higher stroke volume (SV: difference per 3.58kg/m2: 1.49ml/m2.04; 95% CI: 0.62, 2.35; P=0.001) and cardiac output (CO: difference per 3.58kg/m2: 0.11l/min/m1.83; 95% CI: 0.03, 0.19; P=0.01) but no strong evidence for a causal role on systemic vascular resistance or total arterial compliance. Neither analysis supported a causal role of higher BMI on heart rate.

Conclusions

Complementary MR and RbG causal methodologies, together with a range of sensitivity analyses, suggest that higher BMI is likely to cause worse cardiovascular health, specifically higher BP and LVMI, even in youth. Higher BMI also resulted in increased CO in the RbG study, which appeared to be solely driven by SV, as neither MR nor RbG analyses suggested a causal effect of BMI on heart rate. These consistent results support efforts to reduce BMI from a young age to prevent later adverse cardiovascular health and illustrate the potential for phenotypic resolution with maintained analytical power using RbG.

Author Correction: Discovery and refinement of genetic loci associated with cardiometabolic risk using dense imputation maps

In the version of the article published, the surname of author Aaron Isaacs is misspelled as Issacs.

BMI and Mortality in UK Biobank: Revised Estimates Using Mendelian Randomization

OBJECTIVE

The aim of this study was to obtain estimates of the causal relationship between BMI and mortality.

METHODS

Mendelian randomization (MR) with BMI-associated genotypic variation was used to test the causal effect of BMI on all-cause and cause-specific mortality in UK Biobank participants of White British ancestry.

RESULTS

MR analyses supported a causal association between higher BMI and greater risk of all-cause mortality (hazard ratio [HR] per 1 kg/m2 : 1.03; 95% CI: 0.99-1.07) and mortality from cardiovascular diseases (HR: 1.10; 95% CI: 1.01-1.19), specifically coronary heart disease (HR: 1.12; 95% CI: 1.00-1.25) and those excluding coronary heart disease/stroke/aortic aneurysm (HR: 1.24; 95% CI: 1.03-1.48), stomach cancer (HR: 1.18; 95% CI: 0.87-1.62), and esophageal cancer (HR: 1.22; 95% CI: 0.98-1.53), and a decreased risk of lung cancer mortality (HR: 0.96; 95% CI: 0.85-1.08). Sex stratification supported the causal role of higher BMI increasing bladder cancer mortality risk (males) but decreasing respiratory disease mortality risk (males). The J-shaped observational association between BMI and mortality was visible with MR analyses, but the BMI at which mortality was minimized was lower and the association was flatter over a larger BMI range.

CONCLUSIONS

Results support a causal role of higher BMI in increasing the risk of all-cause mortality and mortality from several specific causes.

Evaluation of the causal effects between subjective wellbeing and cardiometabolic health: mendelian randomisation study

OBJECTIVES

To investigate whether the association between subjective wellbeing (subjective happiness and life satisfaction) and cardiometabolic health is causal.

DESIGN

Two sample, bidirectional mendelian randomisation study.

SETTING

Genetic data taken from various cohorts comprised of the general population (mostly individuals of European ancestry, plus a small proportion of other ancestries); follow-up analysis included individuals from the United Kingdom.

PARTICIPANTS

Summary data were used from previous genome wide association studies (number of participants ranging from 83 198 to 339 224), which investigated traits related to cardiovascular or metabolic health, had the largest sample sizes, and consisted of the most similar populations while minimising sample overlap. A follow-up analysis included 337 112 individuals from the UK Biobank (54% female (n=181 363), mean age 56.87 years (standard deviation 8.00) at recruitment).

MAIN OUTCOME MEASURES

Subjective wellbeing and 11 measures of cardiometabolic health (coronary artery disease; myocardial infarction; total, high density lipoprotein, and low density lipoprotein cholesterol; diastolic and systolic blood pressure; body fat; waist to hip ratio; waist circumference; and body mass index).

RESULTS

Evidence of a causal effect of body mass index on subjective wellbeing was seen; each 1 kg/m2 increase in body mass index caused a -0.045 (95% confidence interval -0.084 to -0.006, P=0.02) standard deviation reduction in subjective wellbeing. Follow-up analysis of this association in an independent sample from the UK Biobank provided strong evidence of an effect of body mass index on satisfaction with health (β=-0.035 unit decrease in health satisfaction (95% confidence interval -0.043 to -0.027) per standard deviation increase in body mass index, P<0.001). No clear evidence of a causal effect was seen between subjective wellbeing and the other cardiometabolic health measures, in either direction.

CONCLUSIONS

These results suggest that a higher body mass index is associated with a lower subjective wellbeing. A follow-up analysis confirmed this finding, suggesting that the effect in middle aged people could be driven by satisfaction with health. Body mass index is a modifiable determinant, and therefore, this study provides further motivation to tackle the obesity epidemic because of the knock-on effects of higher body mass index on subjective wellbeing.

FUT2 secretor genotype and susceptibility to infections and chronic conditions in the ALSPAC cohort

Background: The FUT2 (fucosyltransferase-2) gene determines blood group secretor status. Being homozygous for the inactive "non-secretor" rs601338(A) allele confers resistance to certain infections (e.g. Norovirus, Rotavirus) and susceptibility to others (e.g. Haemophilus influenza, Streptococcus pneumonia). Non-secretors also have an increased risk of type 1 diabetes and inflammatory bowel disease. We examined FUT2 genotype, infections and chronic conditions in a population-based cohort. Methods: We studied 7,582 pregnant women from the ALSPAC pregnancy cohort. Infections (measles, mumps, chicken pox, whooping cough, meningitis, herpes, gonorrhea and urinary infections) and chronic conditions (kidney disease, hypertension, diabetes, rheumatism, arthritis, psoriasis, hay fever, asthma, eczema and allergies) were self-reported. FUT2 secretor status was determined from the rs601338 genotype. ABO blood type was obtained from clinical records. Results: Overall, 1920 women (25.3%) were homozygous for the non-secretor allele (AA). Secretor status was associated with mumps, with 68% of non-secretors experiencing this infection, compared to 48% of secretors (RR, 1.40; 95% CI, 1.34-1.46). A weaker association was observed for measles infection (76% vs. 72%; RR, 1.05; 95% CI, 1.02-1.09). Non-secretors also experienced an increased risk of kidney disease (5.4% vs. 3.9%; RR, 1.39; 95% CI, 1.11-1.75). Independent of secretor status, AB blood type was a risk factor for mumps (RR 1.15; 95%CI, 1.03, 1.28 compared to type O). We found no evidence of interaction between secretor status and blood type.  For some conditions, including asthma and arthritis, FUT2 heterozygosity (GA) appeared to confer an intermediate phenotype. There was no strong evidence of association between secretor status and other infections or chronic conditions, although statistical power was limited for rare outcomes. Conclusion: Our results identify an association between FUT2 secretor status and self-reported kidney disease, and confirm a recently reported association with susceptibility to mumps infection. The clinical implications of these associations warrant further investigation.

Systematic Mendelian randomization framework elucidates hundreds of CpG sites which may mediate the influence of genetic variants on disease

We have undertaken a systematic Mendelian randomization (MR) study using methylation quantitative trait loci (meQTL) as genetic instruments to assess the relationship between genetic variation, DNA methylation and 139 complex traits. Using two-sample MR, we identified 1148 associations across 61 traits where genetic variants were associated with both proximal DNA methylation (i.e. cis-meQTL) and complex trait variation (P < 1.39 × 10-08). Joint likelihood mapping provided evidence that the genetic variant which influenced DNA methylation levels for 348 of these associations across 47 traits was also responsible for variation in complex traits. These associations showed a high rate of replication in the BIOS QTL and UK Biobank datasets for 14 selected traits, as 101 of the attempted 128 associations survived multiple testing corrections (P < 3.91 × 10-04). Integrating expression quantitative trait loci (eQTL) data suggested that genetic variants responsible for 306 of the 348 refined meQTL associations also influence gene expression, which indicates a coordinated system of effects that are consistent with causality. CpG sites were enriched for histone mark peaks in tissue types relevant to their associated trait and implicated genes were enriched across relevant biological pathways. Though we are unable to distinguish mediation from horizontal pleiotropy in these analyses, our findings should prove valuable in prioritizing candidate loci where DNA methylation may influence traits and help develop mechanistic insight into the aetiology of complex disease.

Coronary artery disease, genetic risk and the metabolome in young individuals

Background: Genome-wide association studies have identified genetic variants associated with coronary artery disease (CAD) in adults - the leading cause of death worldwide. It often occurs later in life, but variants may impact CAD-relevant phenotypes early and throughout the life-course. Cohorts with longitudinal and genetic data on thousands of individuals are letting us explore the antecedents of this adult disease. Methods: 149 metabolites, with a focus on the lipidome, measured using nuclear magnetic resonance ( 1H-NMR) spectroscopy, and genotype data were available from 5,905 individuals at ages 7, 15, and 17 years from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Linear regression was used to assess the association between the metabolites and an adult-derived genetic risk score (GRS) of CAD comprising 146 variants. Individual variant-metabolite associations were also examined. Results: The CAD-GRS associated with 118 of 149 metabolites (false discovery rate [FDR] < 0.05), the strongest associations being with low-density lipoprotein (LDL) and atherogenic non-LDL subgroups. Nine of 146 variants in the GRS associated with one or more metabolites (FDR < 0.05). Seven of these are within lipid loci: rs11591147 PCSK9, rs12149545 HERPUD1-CETP, rs17091891 LPL, rs515135 APOB, rs602633 CELSR2-PSRC1, rs651821 APOA5, rs7412 APOE-APOC1. All associated with metabolites in the LDL or atherogenic non-LDL subgroups or both including aggregate cholesterol measures. The other two variants identified were rs112635299 SERPINA1 and rs2519093 ABO. Conclusions: Genetic variants that influence CAD risk in adults are associated with large perturbations in metabolite levels in individuals as young as seven. The variants identified are mostly within lipid-related loci and the metabolites they associated with are primarily linked to lipoproteins. This knowledge could allow for preventative measures, such as increased monitoring of at-risk individuals and perhaps treatment earlier in life, to be taken years before any symptoms of the disease arise.

Associations of device-measured physical activity across adolescence with metabolic traits: Prospective cohort study

BACKGROUND

Multiple occasions of device-measured physical activity have not been previously examined in relation to metabolic traits. We described associations of total activity, moderate-to-vigorous physical activity (MVPA), and sedentary time from three accelerometry measures taken across adolescence with detailed traits related to systemic metabolism.

METHODS AND FINDINGS

There were 1,826 male and female participants recruited at birth in 1991-1992 via mothers into the Avon Longitudinal Study of Parents and Children offspring cohort who attended clinics in 2003-2005, 2005-2006, and 2006-2008 who were included in ≥1 analysis. Waist-worn uniaxial accelerometers measured total activity (counts/min), MVPA (min/d), and sedentary time (min/d) over ≥3 d at mean age 12y, 14y, and 15y. Current activity (at age 15y), mean activity across occasions, interaction by previous activity, and change in activity were examined in relation to systolic and diastolic blood pressure, insulin, C-reactive protein, and 230 traits from targeted metabolomics (nuclear magnetic resonance spectroscopy), including lipoprotein cholesterol and triglycerides, amino and fatty acids, glycoprotein acetyls, and others, at age 15y. Mean current total activity was 477.5 counts/min (SD = 164.0) while mean MVPA and sedentary time durations were 23.6 min/d (SD = 17.9) and 522.1 min/d (SD = 66.0), respectively. Mean body mass index at age 15y was 21.4 kg/m2 (SD = 3.5). Correlations between first and last activity measurement occasions were low (e.g., r = 0.40 for counts/min). Current activity was most strongly associated with cholesterol and triglycerides in high-density lipoprotein (HDL) and very low-density lipoprotein (VLDL) particles (e.g., -0.002 mmol/l or -0.18 SD units; 95% CI -0.24--0.11 for triglycerides in chylomicrons and extremely large very low-density lipoprotein [XL VLDL]) and with glycoprotein acetyls (-0.02 mmol/l or -0.16 SD units; 95% CI -0.22--0.10), among others. Associations were similar for mean activity across 3 occasions. Attenuations were modest with adjustment for fat mass index based on dual-energy X-ray absorptiometry (DXA). In mutually adjusted models, higher MVPA and sedentary time were oppositely associated with cholesterol and triglycerides in VLDL and HDL particles (MVPA more strongly with glycoprotein acetyls and sedentary time more strongly with amino acids). Associations appeared less consistent for sedentary time than for MVPA based on longer-term measures and were weak for change in all activity types from age 12y-15y. Evidence was also weak for interaction between activity types at age 15y and previous activity measures in relation to most traits (minimum P = 0.003; median P = 0.26 for counts/min) with interaction coefficients mostly positive. Study limitations include modest sample sizes and relatively short durations of accelerometry measurement on each occasion (3-7 d) and of time lengths between first and last accelerometry occasions (<4 years), which can obscure patterns from chance variation and limit description of activity trajectories. Activity was also recorded using uniaxial accelerometers which predated more sensitive triaxial devices.

CONCLUSIONS

Our results support associations of physical activity with metabolic traits that are small in magnitude and more robust for higher MVPA than lower sedentary time. Activity fluctuates over time, but associations of current activity with most metabolic traits do not differ by previous activity. This suggests that the metabolic effects of physical activity, if causal, depend on most recent engagement.

Assessing the causal association between 25-hydroxyvitamin D and the risk of oral and oropharyngeal cancer using Mendelian randomization

Circulating 25-hydroxyvitamin D (25OHD) is an appealing potential intervention for cancer risk and has been associated with oral and oropharyngeal cancer risk but evidence is inconsistent. The availability of genetic variants, uncorrelated with known confounders, but predictive of 25OHD and genetic data in a large oral and oropharyngeal cancer collaboration aids causal inference when assessing this association. A total of 5,133 oral and oropharyngeal cancer cases and 5,984 controls with genetic data were included in the study. Participants were based in Europe, North America and South America and were part of the Genetic Associations and Mechanisms in Oncology (GAME-ON) Network. Five genetic variants reliably associated with circulating 25OHD were used to create a relative genetic measure of 25OHD. In the absence of measured 25OHD, two-sample Mendelian randomization using individual level outcome data were used to estimate causal odds ratios (OR) for cancer case status per standard deviation increase in log25OHD. Analyses were replicated in an independent population-based cohort (UK Biobank). In the GAME-ON study, there was little evidence of a causal association between circulating 25OHD and oral cancer (OR = 0.86 [0.68;1.09], p = 0.22), oropharyngeal cancer (OR = 1.28 [0.72;2.26], p = 0.40) or when sites were combined (OR = 1.01 [0.74;1.40], p = 0.93). Replication in UK Biobank and pooled estimates produced similar results. Our study suggests that a clinically relevant protective effect of 25OHD on oral and oropharyngeal cancer risk is unlikely and supplementation of the general population with 25OHD is unlikely to be beneficial in preventing these cancers.

Consortium-based genome-wide meta-analysis for childhood dental caries traits

Prior studies suggest dental caries traits in children and adolescents are partially heritable, but there has been no large-scale consortium genome-wide association study (GWAS) to date. We therefore performed GWAS for caries in participants aged 2.5-18.0 years from nine contributing centres. Phenotype definitions were created for the presence or absence of treated or untreated caries, stratified by primary and permanent dentition. All studies tested for association between caries and genotype dosage and the results were combined using fixed-effects meta-analysis. Analysis included up to 19 003 individuals (7530 affected) for primary teeth and 13 353 individuals (5875 affected) for permanent teeth. Evidence for association with caries status was observed at rs1594318-C for primary teeth [intronic within ALLC, odds ratio (OR) 0.85, effect allele frequency (EAF) 0.60, P 4.13e-8] and rs7738851-A (intronic within NEDD9, OR 1.28, EAF 0.85, P 1.63e-8) for permanent teeth. Consortium-wide estimated heritability of caries was low [h2 of 1% (95% CI: 0%: 7%) and 6% (95% CI 0%: 13%) for primary and permanent dentitions, respectively] compared with corresponding within-study estimates [h2 of 28% (95% CI: 9%: 48%) and 17% (95% CI: 2%: 31%)] or previously published estimates. This study was designed to identify common genetic variants with modest effects which are consistent across different populations. We found few single variants associated with caries status under these assumptions. Phenotypic heterogeneity between cohorts and limited statistical power will have contributed; these findings could also reflect complexity not captured by our study design, such as genetic effects which are conditional on environmental exposure.

Genotype-Based Recall Studies in Complex Cardiometabolic Traits

In genotype-based recall (GBR) studies, people (or their biological samples) who carry genotypes of special interest for a given hypothesis test are recalled from a larger cohort (or biobank) for more detailed investigations. There are several GBR study designs that offer a range of powerful options to elucidate (1) genotype-phenotype associations (by increasing the efficiency of genetic association studies, thereby allowing bespoke phenotyping in relatively small cohorts), (2) the effects of environmental exposures (within the Mendelian randomization framework), and (3) gene-treatment interactions (within the setting of GBR interventional trials). In this review, we overview the literature on GBR studies as applied to cardiometabolic health outcomes. We also review the GBR approaches used to date and outline new methods and study designs that might enhance the utility of GBR-focused studies. Specifically, we highlight how GBR methods have the potential to augment randomized controlled trials, providing an alternative application for the now increasingly accepted Mendelian randomization methods usually applied to large-scale population-based data sets. Further to this, we consider how functional and basic science approaches alongside GBR designs offer intellectually intriguing and potentially powerful ways to explore the implications of alterations to specific (and potentially druggable) biological pathways.

Influence of puberty timing on adiposity and cardiometabolic traits: A Mendelian randomisation study

BACKGROUND

Earlier puberty is widely linked with future obesity and cardiometabolic disease. We examined whether age at puberty onset likely influences adiposity and cardiometabolic traits independent of childhood adiposity.

METHODS AND FINDINGS

One-sample Mendelian randomisation (MR) analyses were conducted on up to 3,611 white-European female and male offspring from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort recruited at birth via mothers between 1 April 1991 and 31 December 1992. Time-sensitive exposures were age at menarche and age at voice breaking. Outcomes measured at age 18 y were body mass index (BMI), dual-energy X-ray absorptiometry-based fat and lean mass indices, blood pressure, and 230 cardiometabolic traits derived from targeted metabolomics (150 concentrations plus 80 ratios from nuclear magnetic resonance [NMR] spectroscopy covering lipoprotein subclasses of cholesterol and triglycerides, amino acids, inflammatory glycoproteins, and others). Adjustment was made for pre-pubertal BMI measured at age 8 y. For negative control MR analyses, BMI and cardiometabolic trait measures taken at age 8 y (before puberty, and which therefore cannot be an outcome of puberty itself) were used. For replication analyses, 2-sample MR was conducted using summary genome-wide association study data on up to 322,154 adults for post-pubertal BMI, 24,925 adults for post-pubertal NMR cardiometabolic traits, and 13,848 children for pre-pubertal obesity (negative control). Like observational estimates, 1-sample MR estimates in ALSPAC using 351 polymorphisms for age at menarche (explaining 10.6% of variance) among 2,053 females suggested that later age at menarche (per year) was associated with -1.38 kg/m2 of BMI at age 18 y (or -0.34 SD units, 95% CI -0.46, -0.23; P = 9.77 × 10-09). This coefficient attenuated 10-fold upon adjustment for BMI at age 8 y, to -0.12 kg/m2 (or -0.03 SDs, 95% CI -0.13, 0.07; P = 0.55). Associations with blood pressure were similar, but associations across other traits were small and inconsistent. In negative control MR analyses, later age at menarche was associated with -0.77 kg/m2 of pre-pubertal BMI measured at age 8 y (or -0.39 SDs, 95% CI -0.50, -0.29; P = 6.28 × 10-13), indicating that variants influencing menarche also influence BMI before menarche. Cardiometabolic trait associations were weaker and less consistent among males and both sexes combined. Higher BMI at age 8 y (per 1 kg/m2 using 95 polymorphisms for BMI explaining 3.4% of variance) was associated with earlier menarche among 2,648 females (by -0.26 y, 95% CI -0.37, -0.16; P = 1.16 × 10-06), likewise among males and both sexes combined. In 2-sample MR analyses using 234 polymorphisms and inverse variance weighted (IVW) regression, each year later age at menarche was associated with -0.81 kg/m2 of adult BMI (or -0.17 SD units, 95% CI -0.21, -0.12; P = 4.00 × 10-15). Associations were weaker with cardiometabolic traits. Using 202 polymorphisms, later menarche was associated with lower odds of childhood obesity (IVW-based odds ratio = 0.52 per year later, 95% CI 0.48, 0.57; P = 6.64 × 10-15). Study limitations include modest sample sizes for 1-sample MR, lack of inference to non-white-European populations, potential selection bias through modest completion rates of puberty questionnaires, and likely disproportionate measurement error of exposures by sex. The cardiometabolic traits examined were heavily lipid-focused and did not include hormone-related traits such as insulin and insulin-like growth factors.

CONCLUSIONS

Our results suggest that puberty timing has a small influence on adiposity and cardiometabolic traits and that preventive interventions should instead focus on reducing childhood adiposity.

Longitudinal serological measures of common infection in the Avon Longitudinal Study of Parents and Children cohort

Antibodies against pathogens provide information on exposure to infectious agents and are meaningful measures of past and present infection. Antibodies were measured in the plasma of children that are the offspring in a population-based birth cohort, the Avon Longitudinal Study of Parents and Children (ALSPAC). Plasma was collected during clinics at age 5, 7, 11 and 15 years. The antigens examined include: fungal ( Saccharomyces cerevisiae); protozoan ( Toxoplasma gondii and surface antigen 1 of  T. gondii); herpes viruses (cytomegalovirus, Epstein-Barr virus, herpes simplex virus type 1); common colds (influenza virus subtypes H1N1 and H3N2); other antigens (measles); animal (feline herpes virus, Theiler's virus); bacteria ( Helicobacter pylori); dietary antigens (bovine casein alpha protein, bovine casein beta protein). Alongside the depth of data available within the ALSPAC cohort, this longitudinal resource will enable the investigation of the association between infections and a wide variety of outcomes.

Longitudinal serological measures of common infection in the Avon Longitudinal Study of Parents and Children cohort

Antibodies against pathogens provide information on exposure to infectious agents and are meaningful measures of past and present infection. Antibodies were measured in the plasma of children that are the offspring in a population-based birth cohort, the Avon Longitudinal Study of Parents and Children (ALSPAC). Plasma was collected during clinics at age 5, 7, 11 and 15 years. The antigens examined include: fungal ( Saccharomyces cerevisiae); protozoan ( Toxoplasma gondii and surface antigen 1 of  T. gondii); herpes viruses (cytomegalovirus, Epstein-Barr virus, herpes simplex virus type 1); common colds (influenza virus subtypes H1N1 and H3N2); other antigens (measles); animal (feline herpes virus, Theiler’s virus); bacteria ( Helicobacter pylori); dietary antigens (bovine casein alpha protein, bovine casein beta protein). Alongside the depth of data available within the ALSPAC cohort, this longitudinal resource will enable the investigation of the association between infections and a wide variety of outcomes.

The MR-Base platform supports systematic causal inference across the human phenome

Results from genome-wide association studies (GWAS) can be used to infer causal relationships between phenotypes, using a strategy known as 2-sample Mendelian randomization (2SMR) and bypassing the need for individual-level data. However, 2SMR methods are evolving rapidly and GWAS results are often insufficiently curated, undermining efficient implementation of the approach. We therefore developed MR-Base (<ext-link ext-link-type="uri" xlink:href="http://www.mrbase.org">http://www.mrbase.org</ext-link>): a platform that integrates a curated database of complete GWAS results (no restrictions according to statistical significance) with an application programming interface, web app and R packages that automate 2SMR. The software includes several sensitivity analyses for assessing the impact of horizontal pleiotropy and other violations of assumptions. The database currently comprises 11 billion single nucleotide polymorphism-trait associations from 1673 GWAS and is updated on a regular basis. Integrating data with software ensures more rigorous application of hypothesis-driven analyses and allows millions of potential causal relationships to be efficiently evaluated in phenome-wide association studies.