Search

Abstract Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

OBJECTIVES:. We sought to assess whether genetic associations with metabolite concentrations in septic shock patients could be used to identify pathways of potential importance for understanding sepsis pathophysiology. DESIGN:. Retrospective multicenter cohort studies of septic shock patients. SETTING:. All participants who were admitted to 27 participating hospital sites in three countries (Australia, New Zealand, and the United Kingdom) were eligible for inclusion. PATIENTS:. Adult, critically ill, mechanically ventilated patients with septic shock (n = 230) who were a subset of the Adjunctive Corticosteroid Treatment in Critically Ill Patients with Septic Shock trial (ClinicalTrials.gov number: NCT01448109). INTERVENTIONS:. None. MEASUREMENTS AND MAIN RESULTS:. A genome-wide association study was conducted for a range of serum metabolite levels for participants. Genome-wide significant associations (p ≤ 5 × 10–8) were found for the two major ketone bodies (3-hydroxybutyrate [rs2456680] and acetoacetate [rs2213037] and creatinine (rs6851961). One of these single-nucleotide polymorphisms (SNPs) (rs2213037) was located in the alcohol dehydrogenase cluster of genes, which code for enzymes related to the metabolism of acetoacetate and, therefore, presents a plausible association for this metabolite. None of the three SNPs showed strong associations with risk of sepsis, 28- or 90-day mortality, or Acute Physiology and Chronic Health Evaluation score (a measure of sepsis severity). CONCLUSIONS:. We suggest that the genetic associations with metabolites may reflect a starvation response rather than processes involved in sepsis pathophysiology. However, our results require further investigation and replication in both healthy and diseased cohorts including those of different ancestry.

Abstract Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Abstract Parental genes may indirectly influence offspring psychiatric outcomes through the environment that parents create for their children. These indirect genetic effects, also known as genetic nurture, could explain individual differences in common internalising and externalising psychiatric symptoms during childhood. Advanced statistical genetic methods leverage data from families to estimate the overall contribution of parental genetic nurture effects. This study included up to 10,499 children, 5990 mother–child pairs, and 6,222 father–child pairs from the Norwegian Mother Father and Child Study. Genome-based restricted maximum likelihood (GREML) models were applied using software packages GCTA and M-GCTA to estimate variance in maternally reported depressive, disruptive, and attention-deficit hyperactivity disorder (ADHD) symptoms in 8-year-olds that was explained by direct offspring genetic effects and maternal or paternal genetic nurture. There was no strong evidence of genetic nurture in this sample, although a suggestive paternal genetic nurture effect on offspring depressive symptoms (variance explained (V) = 0.098, standard error (SE) = 0.057) and a suggestive maternal genetic nurture effect on ADHD symptoms (V = 0.084, SE = 0.058) was observed. The results indicate that parental genetic nurture effects could be of some relevance in explaining individual differences in childhood psychiatric symptoms. However, robustly estimating their contribution is a challenge for researchers given the current paucity of large-scale samples of genotyped families with information on childhood psychiatric outcomes.

Summary: Background: Hip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest. Methods: GWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA. Findings: 50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82–1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life. Interpretations: One group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism. Funding: Primarily funded by the Medical Research Council and Wellcome Trust.

Abstract Background Observational epidemiological studies suggest a link between several factors related to ovulation and reproductive function and endometrial cancer (EC) risk; however, it is not clear whether these relationships are causal, and whether the risk factors act independently of each other. The aim of this study was to investigate putative causal relationships between the number of live births, age at last live birth, and years ovulating and EC risk. Methods We conducted a series of observational analyses to investigate various risk factors and EC risk in the UK Biobank (UKBB). Additionally, multivariate analysis was performed to elucidate the relationship between the number of live births, age at last live birth, and years ovulating and other related factors such as age at natural menopause, age at menarche, and body mass index (BMI). Secondly, we used Mendelian randomization (MR) to assess if these observed relationships were causal. Genome-wide significant single nucleotide polymorphisms (SNPs) were extracted from previous studies of woman’s number of live births, age at menopause and menarche, and BMI. We conducted a genome-wide association analysis using the UKBB to identify SNPs associated with years ovulating, years using the contraceptive pill, and age at last live birth. Results We found evidence for a causal effect of the number of live births (inverse variance weighted (IVW) odds ratio (OR): 0.537, p = 0.006), the number of years ovulating (IVW OR: 1.051, p = 0.014), in addition to the known risk factors BMI, age at menarche, and age at menopause on EC risk in the univariate MR analyses. Due to the close relationships between these factors, we followed up with multivariable MR (MVMR) analysis. Results from the MVMR analysis showed that number of live births had a causal effect on EC risk (OR: 0.783, p = 0.036) independent of BMI, age at menarche and age at menopause. Conclusions MVMR analysis showed that the number of live births causally reduced the risk of EC.

Abstract Handedness has low heritability and epigenetic mechanisms have been proposed as an etiological mechanism. To examine this hypothesis, we performed an epigenome-wide association study of left-handedness. In a meta-analysis of 3914 adults of whole-blood DNA methylation, we observed that CpG sites located in proximity of handedness-associated genetic variants were more strongly associated with left-handedness than other CpG sites (P = 0.04), but did not identify any differentially methylated positions. In longitudinal analyses of DNA methylation in peripheral blood and buccal cells from children (N = 1737), we observed moderately stable associations across age (correlation range [0.355–0.578]), but inconsistent across tissues (correlation range [− 0.384 to 0.318]). We conclude that DNA methylation in peripheral tissues captures little of the variance in handedness. Future investigations should consider other more targeted sources of tissue, such as the brain.

Summary: Background: Amino acids are key to protein synthesis, energy metabolism, cell signaling and gene expression; however, the contribution of specific maternal amino acids to fetal growth is unclear. Methods: We explored the effect of maternal circulating amino acids on fetal growth, proxied by birthweight, using two-sample Mendelian randomisation (MR) and summary data from a genome-wide association study (GWAS) of serum amino acids levels (sample 1, n = 86,507) and a maternal GWAS of offspring birthweight in UK Biobank and Early Growth Genetics Consortium, adjusting for fetal genotype effects (sample 2, n = 406,063 with maternal and/or fetal genotype effect estimates). A total of 106 independent single nucleotide polymorphisms robustly associated with 19 amino acids (p

Abstract Background Greater maternal adiposity before or during pregnancy is associated with greater offspring adiposity throughout childhood, but the extent to which this is due to causal intrauterine or periconceptional mechanisms remains unclear. Here, we use Mendelian randomisation (MR) with polygenic risk scores (PRS) to investigate whether associations between maternal pre-/early pregnancy body mass index (BMI) and offspring adiposity from birth to adolescence are causal. Methods We undertook confounder adjusted multivariable (MV) regression and MR using mother-offspring pairs from two UK cohorts: Avon Longitudinal Study of Parents and Children (ALSPAC) and Born in Bradford (BiB). In ALSPAC and BiB, the outcomes were birthweight (BW; N = 9339) and BMI at age 1 and 4 years (N = 8659 to 7575). In ALSPAC only we investigated BMI at 10 and 15 years (N = 4476 to 4112) and dual-energy X-ray absorptiometry (DXA) determined fat mass index (FMI) from age 10–18 years (N = 2659 to 3855). We compared MR results from several PRS, calculated from maternal non-transmitted alleles at between 29 and 80,939 single nucleotide polymorphisms (SNPs). Results MV and MR consistently showed a positive association between maternal BMI and BW, supporting a moderate causal effect. For adiposity at most older ages, although MV estimates indicated a strong positive association, MR estimates did not support a causal effect. For the PRS with few SNPs, MR estimates were statistically consistent with the null, but had wide confidence intervals so were often also statistically consistent with the MV estimates. In contrast, the largest PRS yielded MR estimates with narrower confidence intervals, providing strong evidence that the true causal effect on adolescent adiposity is smaller than the MV estimates (P difference = 0.001 for 15-year BMI). This suggests that the MV estimates are affected by residual confounding, therefore do not provide an accurate indication of the causal effect size. Conclusions Our results suggest that higher maternal pre-/early-pregnancy BMI is not a key driver of higher adiposity in the next generation. Thus, they support interventions that target the whole population for reducing overweight and obesity, rather than a specific focus on women of reproductive age.

Estimating direct and indirect effects of genotypes on phenotypes is important for genetic analyses such as Mendelian randomization. Here the authors compare five different methods to estimate direct and indirect genetic effects using summary results statistics that account for sample overlap.

Osteocytes are the master regulatory cells within the skeleton. Here, the authors map the transcriptome of osteocytes from diverse skeletal sites, ages and between sexes and identify an osteocyte transcriptome signature associated with rare skeletal disorders and common complex skeletal diseases.

Observationally, lower birthweight is a risk factor for cardiometabolic disease. Using Mendelian Randomization, the authors investigate whether maternal genetic factors that lower offspring birthweight also increase offspring cardiometabolic risk and show that the observational correlation is unlikely to be due to the intrauterine environment.

Family-based study designs have been applied to resolve confounding by population stratification, dynastic effects and assortative mating in genetic association analyses. Here, Brumpton et al. describe theory and simulations for overcoming such biases in Mendelian randomization through within-family studies.

Genetic variants at the RSPO3 locus are associated with waist-to-hip ratio (adjusted for BMI). Here, Loh et al. describe two independent RSPO3 signals that associate with body fat distribution, perform fine-mapping and explore the function of RSPO3 in human adipocyte biology and body fat distribution in zebrafish

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

Melanocytic nevus count is associated with melanoma risk. In this study, a meta-analysis of 11 nevus GWAS studies identifies novel SNPs in KITLG and 9q32, and bivariate analysis with melanoma GWAS meta-analysis reveals that most nevus genes affect melanoma risk, while melanoma risk loci do not alter the nevus count.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

Abstract Introduction Hearing loss has been identified as the potentially largest modifiable risk factor for Alzheimer's disease (AD), estimated to account for a similar increase in AD risk as the apolipoprotein E (APOE) gene. Methods We investigated the genetic relationship between hearing loss and AD, and sought evidence for a causal relationship. Results We found a significant genetic overlap between hearing impairment and AD and a polygenic risk score for AD was able to significantly predict hearing loss in an independent cohort. Additionally, regions of the genome involved in inflammation were identified to be shared between hearing difficulty and AD. However, causality tests found no significant evidence of a causal relationship between these traits in either direction. Discussion Overall, these results show that the relationship between hearing difficulty and AD may, in part, be due to shared genes and immune response pathways between the traits. However, currently available data do not support a causal relationship.

Epidemiological studies have identified many risk factors for osteoporosis, however it is unclear whether these observational associations reflect true causal effects, or the effects of latent confounding or reverse causality. Mendelian randomization (MR) enables causal relationships to be evaluated, by examining the relationship between genetic susceptibility to the risk factor in question, and the disease outcome of interest. This has been facilitated by the development of two-sample MR analysis, where the exposure and outcome are measured in different studies, and by exploiting summary result statistics from large well-powered genome-wide association studies that are available for thousands of traits. Though MR has several inherent limitations, the field is rapidly evolving and at least 14 methodological extensions have been developed to overcome these. The present paper aims to discuss some of the limitations in the MR analytical framework, and how this method has been applied to the osteoporosis field, helping to reinforce conclusions about causality, and discovering potential new regulatory pathways, exemplified by our recent MR study of sclerostin.

Recall-by-Genotype (RbG) is an approach to recall participants from genetic studies based on their specific genotype for further, more extensive phenotyping. Here, the authors discuss examples of RbG as well as practical and ethical considerations and provide an online tool to aid in designing RbG studies.

The original version of this Article contained errors in the spelling of the authors Fan Liu and M. Arfan Ikram, which were incorrectly given as Fan Lui and Arfan M. Ikram. In addition, the original version of this Article also contained errors in the author affiliations which are detailed in the associated Publisher Correction.

Author summary Bone mineral density and lean skeletal mass are heritable traits. Here, Medina-Gomez and colleagues perform bivariate GWAS analyses of total body lean mass and bone mass density in children, and show genetic loci with pleiotropic effects on both traits.

This report by the Consortium for Refractive Error and Myopia uses gene-environment-wide interaction study (GEWIS) to identify genetic loci that affect environmental influence in myopia development, and identifies ethnic specific genetic loci that attribute to eye refractive errors.

A genome-wide association meta-analysis study here shows novel genetic loci to be associated to body fat percentage, and describes cross-phenotype association that further demonstrate a close relationship between adiposity and cardiovascular disease risk.

Rheumatoid arthritis (RA) is a disease of chronic systemic inflammation (SI). In the present study, we used four datasets to explore whether methylation-derived neutrophil-to-lymphocyte ratio (mdNLR) might be a marker of SI in new onset, untreated, and treated prevalent RA cases and/or a marker of treatment response to the tumour necrosis factor inhibitor (TNFi) etanercept. mdNLR was associated with increased odds of being a new onset RA case (OR = 2.32, 95% CI = 1.95–2.80, P

Alcohol problems represent a classic example of a complex behavioral outcome that is likely influenced by many genes of small effect. A polygenic approach, which examines aggregate measured genetic effects, can have predictive power in cases where individual genes or genetic variants do not. In the current study, we first tested whether polygenic risk for alcohol problems—derived from genome-wide association estimates of an alcohol problems factor score from the age 18 assessment of the Avon Longitudinal Study of Parents and Children (ALSPAC; n = 4304 individuals of European descent; 57% female)—predicted alcohol problems earlier in development (age 14) in an independent sample (FinnTwin12; n = 1162; 53% female). We then tested whether environmental factors (parental knowledge and peer deviance) moderated polygenic risk to predict alcohol problems in the FinnTwin12 sample. We found evidence for both polygenic association and for additive polygene-environment interaction. Higher polygenic scores predicted a greater number of alcohol problems (range of Pearson partial correlations 0.07–0.08, all p-values ≤ 0.01). Moreover, genetic influences were significantly more pronounced under conditions of low parental knowledge or high peer deviance (unstandardized regression coefficients (b), p-values (p), and percent of variance (R2) accounted for by interaction terms: b = 1.54, p = 0.02, R2 = 0.33%; b = 0.94, p = 0.04, R2 = 0.30%, respectively). Supplementary set-based analyses indicated that the individual top single nucleotide polymorphisms (SNPs) contributing to the polygenic scores were not individually enriched for gene-environment interaction. Although the magnitude of the observed effects are small, this study illustrates the usefulness of polygenic approaches for understanding the pathways by which measured genetic predispositions come together with environmental factors to predict complex behavioral outcomes.

Background: Physical activity has been associated with preventing the development of type 2 diabetes and atherosclerotic cardiovascular disease. However, our understanding of the precise molecular mechanisms underlying these effects remains incomplete and good biomarkers to objectively assess physical activity are lacking. Methods: We analyzed 3072 serum proteins in 26 men, normal weight or overweight, undergoing 12 weeks of a combined strength and endurance exercise intervention. We estimated insulin sensitivity with hyperinsulinemic euglycemic clamp, maximum oxygen uptake, muscle strength, and used MRI/MRS to evaluate body composition and organ fat depots. Muscle and subcutaneous adipose tissue biopsies were used for mRNA sequencing. Additional association analyses were performed in samples from up to 47,747 individuals in the UK Biobank, as well as using two-sample Mendelian randomization and mice models. Results: Following 12 weeks of exercise intervention, we observed significant changes in 283 serum proteins. Notably, 66 of these proteins were elevated in overweight men and positively associated with liver fat before the exercise regimen, but were normalized after exercise. Furthermore, for 19.7 and 12.1% of the exercise-responsive proteins, corresponding changes in mRNA expression levels in muscle and fat, respectively, were shown. The protein CD300LG displayed consistent alterations in blood, muscle, and fat. Serum CD300LG exhibited positive associations with insulin sensitivity, and to angiogenesis-related gene expression in both muscle and fat. Furthermore, serum CD300LG was positively associated with physical activity and negatively associated with glucose levels in the UK Biobank. In this sample, the association between serum CD300LG and physical activity was significantly stronger in men than in women. Mendelian randomization analysis suggested potential causal relationships between levels of serum CD300LG and fasting glucose, 2 hr glucose after an oral glucose tolerance test, and HbA1c. Additionally, Cd300lg responded to exercise in a mouse model, and we observed signs of impaired glucose tolerance in male, but not female, Cd300lg knockout mice. Conclusions: Our study identified several novel proteins in serum whose levels change in response to prolonged exercise and were significantly associated with body composition, liver fat, and glucose homeostasis. Serum CD300LG increased with physical activity and is a potential causal link to improved glucose levels. CD300LG may be a promising exercise biomarker and a therapeutic target in type 2 diabetes. Funding: South-Eastern Norway Regional Health Authority, Simon Fougners Fund, Diabetesforbundet, Johan Selmer Kvanes’ legat til forskning og bekjempelse av sukkersyke. The UK Biobank resource reference 53641. Australian National Health and Medical Research Council Investigator Grant (APP2017942). Australian Research Council Discovery Early Career Award (DE220101226). Research Council of Norway (Project grant: 325640 and Mobility grant: 287198). The Medical Student Research Program at the University of Oslo. Novo Nordisk Fonden Excellence Emerging Grant in Endocrinology and Metabolism 2023 (NNF23OC0082123). Clinical trial number: clinicaltrials.gov: NCT01803568.

Disassortative mating refers to the phenomenon in which individuals with dissimilar genotypes and/or phenotypes mate with one another more frequently than would be expected by chance. Although the existence of disassortative mating is well established in plant and animal species, the only documented example of negative assortment in humans involves dissimilarity at the major histocompatibility complex (MHC) locus. Previous studies investigating mating patterns at the MHC have been hampered by limited sample size and contradictory findings. Inspired by the sparse and conflicting evidence, we investigated the role that the MHC region played in human mate selection using genome-wide association data from 872 European American spouses from the Health and Retirement Study (HRS). First, we treated the MHC region as a whole, and investigated genomic similarity between spouses using three levels of genomic variation: single-nucleotide polymorphisms (SNPs), classical human leukocyte antigen (HLA) alleles (both four-digit and two-digit classifications), and amino acid polymorphisms. The extent of MHC dissimilarity between spouses was assessed using a permutation approach. Second, we investigated fine scale mating patterns by testing for deviations from random mating at individual SNPs, HLA genes, and amino acids in HLA molecules. Third, we assessed how extreme the spousal relatedness at the MHC region was compared to the rest of the genome, to distinguish the MHC-specific effects from genome-wide effects. We show that neither the MHC region, nor any single SNPs, classic HLA alleles, or amino acid polymorphisms within the MHC region, were significantly dissimilar between spouses relative to non-spouse pairs. However, dissimilarity in the MHC region was extreme relative to the rest of genome for both spousal and non-spouse pairs. Despite the long-standing controversy, our analyses did not support a significant role of MHC dissimilarity in human mate choice.

Estimating effects of parental and sibling genotypes (indirect genetic effects) can provide insight into how the family environment influences phenotypic variation. There is growing molecular genetic evidence for effects of parental phenotypes on their offspring (e.g. parental educational attainment), but the extent to which siblings affect each other is currently unclear. Here we used data from samples of unrelated individuals, without (singletons) and with biological full-siblings (non-singletons), to investigate and estimate sibling effects. Indirect genetic effects of siblings increase (or decrease) the covariance between genetic variation and a phenotype. It follows that differences in genetic association estimates between singletons and non-singletons could indicate indirect genetic effects of siblings if there is no heterogeneity in other sources of genetic association between singletons and non-singletons. We used UK Biobank data to estimate polygenic score (PGS) associations for height, BMI and educational attainment in self-reported singletons (N = 50,143) and non-singletons (N = 328,549). The educational attainment PGS association estimate was 12% larger (95% C.I. 3%, 21%) in the non-singleton sample than in the singleton sample, but the height and BMI PGS associations were consistent. Birth order data suggested that the difference in educational attainment PGS associations was driven by individuals with older siblings rather than firstborns. The relationship between number of siblings and educational attainment PGS associations was non-linear; PGS associations were 24% smaller in individuals with 6 or more siblings compared to the rest of the sample (95% C.I. 11%, 38%). We estimate that a 1 SD increase in sibling educational attainment PGS corresponds to a 0.025 year increase in the index individual's years in schooling (95% C.I. 0.013, 0.036). Our results suggest that older siblings may influence the educational attainment of younger siblings, adding to the growing evidence that effects of the environment on phenotypic variation partially reflect social effects of germline genetic variation in relatives.