Increased serum neopterin had been described in older age two decades ago. Neopterin is a biomarker of systemic adaptive immune activation that could be potentially implicated in metabolic syndrome (MetS). Measurements of waist circumference, triglycerides, high-density lipoprotein cholesterol (HDLC), systolic and diastolic blood pressure, glycated hemoglobin as components of MetS definition, and plasma total neopterin concentrations were performed in 594 participants recruited in the European Prospective Investigation into Cancer and Nutrition (EPIC).

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

For complex traits, most associated single nucleotide variants (SNV) discovered to date have a small effect, and detection of association is only possible with large sample sizes. Because of patient confidentiality concerns, it is often not possible to pool genetic data from multiple cohorts, and meta-analysis has emerged as the method of choice to combine results from multiple studies. Many meta-analysis methods are available for single SNV analyses. As new approaches allow the capture of low frequency and rare genetic variation, it is of interest to jointly consider multiple variants to improve power. However, for the analysis of haplotypes formed by multiple SNVs, meta-analysis remains a challenge, because different haplotypes may be observed across studies. We propose a two-stage meta-analysis approach to combine haplotype analysis results. In the first stage, each cohort estimate haplotype effect sizes in a regression framework, accounting for relatedness among observations if appropriate. For the second stage, we use a multivariate generalized least square meta-analysis approach to combine haplotype effect estimates from multiple cohorts. Haplotype-specific association tests and a global test of independence between haplotypes and traits are obtained within our framework. We demonstrate through simulation studies that we control the type-I error rate, and our approach is more powerful than inverse variance weighted meta-analysis of single SNV analysis when haplotype effects are present. We replicate a published haplotype association between fasting glucose-associated locus (G6PC2) and fasting glucose in seven studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and we provide more precise haplotype effect estimates.

The bone-sparing effect of estrogens is mediated primarily via estrogen receptor (ER)α, which stimulates gene transcription through activation function (AF)-1 and AF-2. The role of ERαAF-1 for the estradiol (E2) effects is tissue specific. The selective ER modulators (SERMs) raloxifene (Ral), lasofoxifene (Las), and bazedoxifene (Bza) can be used to treat postmenopausal osteoporosis. They all reduce the risk for vertebral fractures, whereas Las and partly Bza, but not Ral, reduce the risk for nonvertebral fractures. Here, we have compared the tissue specificity of Ral, Las, and Bza and evaluated the role of ERαAF-1 for the effects of these SERMs, with an emphasis on bone parameters. We treated ovariectomized (OVX) wild-type (WT) mice and OVX mice lacking ERαAF-1 (ERαAF-1(0)) with E2, Ral, Las, or Bza. All three SERMs increased trabecular bone mass in the axial skeleton. In the appendicular skeleton, only Las increased the trabecular bone volume/tissue volume and trabecular number, whereas both Ral and Las increased the cortical bone thickness and strength. However, Ral also increased cortical porosity. The three SERMs had only a minor effect on uterine weight. Notably, all evaluated effects of these SERMs were absent in ovx ERαAF-1(0) mice. In conclusion, all SERMs had similar effects on axial bone mass. However, the SERMs had slightly different effects on the appendicular skeleton since only Las increased the trabecular bone mass and only Ral increased the cortical porosity. Importantly, all SERM effects require a functional ERαAF-1 in female mice. These results could lead to development of more specific treatments for osteoporosis.

Screening gene function in vivo is a powerful approach to discover novel drug targets. We present high-throughput screening (HTS) data for 3 762 distinct global gene knockout (KO) mouse lines with viable adult homozygous mice generated using either gene-trap or homologous recombination technologies. Bone mass was determined from DEXA scans of male and female mice at 14 weeks of age and by microCT analyses of bones from male mice at 16 weeks of age. Wild-type (WT) cagemates/littermates were examined for each gene KO. Lethality was observed in an additional 850 KO lines. Since primary HTS are susceptible to false positive findings, additional cohorts of mice from KO lines with intriguing HTS bone data were examined. Aging, ovariectomy, histomorphometry and bone strength studies were performed and possible non-skeletal phenotypes were explored. Together, these screens identified multiple genes affecting bone mass: 23 previously reported genes (Calcr, Cebpb, Crtap, Dcstamp, Dkk1, Duoxa2, Enpp1, Fgf23, Kiss1/Kiss1r, Kl (Klotho), Lrp5, Mstn, Neo1, Npr2, Ostm1, Postn, Sfrp4, Slc30a5, Slc39a13, Sost, Sumf1, Src, Wnt10b), five novel genes extensively characterized (Cldn18, Fam20c, Lrrk1, Sgpl1, Wnt16), five novel genes with preliminary characterization (Agpat2, Rassf5, Slc10a7, Slc26a7, Slc30a10) and three novel undisclosed genes coding for potential osteoporosis drug targets.

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

The third generation selective estrogen receptor modulators lasofoxifene (las) and bazedoxifene (bza) are indicated for treatment of postmenopausal osteoporosis. 17β-Estradiol (E2) and the second generation SERM raloxifene (ral) have major effects on the immune system, particularly on B cells. Treatment with E2 or ral inhibits B lymphopoiesis and treatment with E2, but not ral, stimulates antibody production. The effects of las and bza on the immune system have not been studied. Therefore, the aim of this study was to investigate their role in B cell development, maturation, and function. C57BL/6 mice were sham-operated or ovariectomized (ovx) and treated with vehicle, E2, ral, las, or bza. All substances increased total bone mineral density in ovx mice, as measured by peripheral quantitative computed tomography. In uterus, bza alone lacked agonistic effect in ovx mice and even acted as an antagonist in sham mice. As expected, E2 decreased B cell numbers at all developmental stages from pre-BI cells (in bone marrow) to transitional 1 (T1) B cells (in spleen) and increased marginal zone (MZ) B cells as determined by flow cytometry. However, treatment with las or bza only decreased the last stages of bone marrow B cell development and splenic T1 B cells, but had no effect MZ B cells. E2 increased antibody-producing cells quantified by ELISPOT, but las or bza did not. In conclusion, las and bza differ from E2 by retaining normal number of cells at most B cell stages during B lymphopoiesis and maturation and by not increasing antibody-producing cells.

Treatment with anti-inflammatory glucocorticoids is associated with osteoporosis. Many of the treated patients are postmenopausal women, who even without treatment have an increased risk of osteoporosis. Lymphocytes have been shown to play a role in postmenopausal and arthritis-induced osteoporosis, and they are targeted by glucocorticoids. The aim of this study was to investigate the mechanisms behind effects of glucocorticoids on bone during health and menopause, focusing on lymphocytes. Female C57BL/6 or SCID mice were therefore sham-operated or ovariectomized and 2 weeks later treatment with dexamethasone (dex), the nonsteroidal anti-inflammatory drug carprofen, or vehicle was started and continued for 2.5 weeks. At the termination of experiments, femurs were phenotyped using peripheral quantitative computed tomography and high-resolution micro-computed tomography, and markers of bone turnover were analyzed in serum. T and B lymphocyte populations in bone marrow and spleen were analyzed by flow cytometry. Dex-treated C57BL/6 mice had increased trabecular bone mineral density, but lower cortical content and thickness compared with vehicle-treated mice. The dex-treated mice also had lower levels of bone turnover markers and markedly decreased numbers of spleen T and B lymphocytes. In contrast, these effects could not be repeated when mice were treated with the nonsteroidal anti-inflammatory drug carprofen. In addition, dex did not increase trabecular bone in ovariectomized SCID mice lacking functional T and B lymphocytes. In contrast to most literature, the results from this study indicate that treatment with dex increased trabecular bone density, which may indicate that this effect is associated with corticosteroid-induced alterations of the lymphocyte populations.

Many genetic loci affect circulating lipid levels, but it remains unknown whether lifestyle factors, such as physical activity, modify these genetic effects. To identify lipid loci interacting with physical activity, we performed genome-wide analyses of circulating HDL cholesterol, LDL cholesterol, and triglyceride levels in up to 120,979 individuals of European, African, Asian, Hispanic, and Brazilian ancestry, with follow-up of suggestive associations in an additional 131,012 individuals. We find four loci, in/near CLASP1, LHX1, SNTA1, and CNTNAP2, that are associated with circulating lipid levels through interaction with physical activity; higher levels of physical activity enhance the HDL cholesterol-increasing effects of the CLASP1, LHX1, and SNTA1 loci and attenuate the LDL cholesterol-increasing effect of the CNTNAP2 locus. The CLASP1, LHX1, and SNTA1 regions harbor genes linked to muscle function and lipid metabolism. Our results elucidate the role of physical activity interactions in the genetic contribution to blood lipid levels.

Both fibroblast growth factors (FGFs), by binding to FGF receptors (FGFRs), and activation of the gravitostat, by artificial loading, decrease the body weight (BW). Previous studies demonstrate that both the FGF system and loading have the capacity to regulate BW independently of leptin. The aim of the current study was to determine the possible interactions between the effect of increased loading and the FGF system for the regulation of BW. We observed that the BW-reducing effect of increased loading was abolished in mice treated with a monoclonal antibody directed against FGFR1c, suggesting interactions between the two systems. As serum levels of endocrine FGF21 and hepatic FGF21 mRNA were increased in the loaded mice compared with the control mice, we first evaluated the loading response in FGF21 over expressing mice with constant high FGF21 levels. Leptin treatment, but not increased loading, decreased the BW in the FGF21-overexpressing mice, demonstrating that specifically the loading effect is attenuated in the presence of high activity in the FGF system. However, as FGF21 knockout mice displayed a normal loading response on BW, FGF21 is neither mediating nor essential for the loading response. In conclusion, the BW-reducing effect of increased loading but not of leptin treatment is blocked by high activity in the FGF system. We propose that both the gravitostat and the FGF system regulate BW independently of leptin and that pharmacologically enhanced activity in the FGF system reduces the sensitivity of the gravitostat.

Glaucoma is the leading cause of irreversible blindness globally 1 . Despite its gravity, the disease is frequently undiagnosed in the community 2 . Raised intraocular pressure (IOP) is the most important risk factor for primary open-angle glaucoma (POAG)3,4. Here we present a meta-analysis of 139,555 European participants, which identified 112 genomic loci associated with IOP, 68 of which are novel. These loci suggest a strong role for angiopoietin-receptor tyrosine kinase signaling, lipid metabolism, mitochondrial function and developmental processes underlying risk for elevated IOP. In addition, 48 of these loci were nominally associated with glaucoma in an independent cohort, 14 of which were significant at a Bonferroni-corrected threshold. Regression-based glaucoma-prediction models had an area under the receiver operating characteristic curve (AUROC) of 0.76 in US NEIGHBORHOOD study participants and 0.74 in independent glaucoma cases from the UK Biobank. Genetic-prediction models for POAG offer an opportunity to target screening and timely therapy to individuals most at risk.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Age at adiposity rebound (AR) is associated with obesity and Type 2 Diabetes in adults. The aim of the present study was to investigate the role of age at AR in adult fat mass, fat distribution and pubertal timing for a Swedish cohort.

Excess body weight, physical activity, smoking, alcohol consumption and certain dietary factors are individually related to colorectal cancer (CRC) risk; however, little is known about their joint effects. The aim of this study was to develop a healthy lifestyle index (HLI) composed of five potentially modifiable lifestyle factors--healthy weight, physical activity, non-smoking, limited alcohol consumption and a healthy diet, and to explore the association of this index with CRC incidence using data collected within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged > 55 years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey–Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey–Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han–Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p < 5 × 10− 8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p = 4.6 × 10− 8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98–1.14; p = 0.17), displaying high degree of heterogeneity (I2 = 57%; Qhet p = 0.0006). Under Han–Eskin alternative random effects model the summary effect was significant (p = 0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size > 1.25) may still be consistent with an effect size < 1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Statin treatment and variants in the gene encoding HMG-CoA reductase are associated with reductions in both the concentration of LDL cholesterol and the risk of coronary heart disease, but also with modest hyperglycaemia, increased bodyweight, and modestly increased risk of type 2 diabetes, which in no way offsets their substantial benefits. We sought to investigate the associations of LDL cholesterol-lowering PCSK9 variants with type 2 diabetes and related biomarkers to gauge the likely effects of PCSK9 inhibitors on diabetes risk.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Vitamin D insufficiency is common, correctable, and influenced by genetic factors, and it has been associated with risk of several diseases. We sought to identify low-frequency genetic variants that strongly increase the risk of vitamin D insufficiency and tested their effect on risk of multiple sclerosis, a disease influenced by low vitamin D concentrations. We used whole-genome sequencing data from 2,619 individuals through the UK10K program and deep-imputation data from 39,655 individuals genotyped genome-wide. Meta-analysis of the summary statistics from 19 cohorts identified in CYP2R1 the low-frequency (minor allele frequency = 2.5%) synonymous coding variant g.14900931G>A (p.Asp120Asp) (rs117913124[A]), which conferred a large effect on 25-hydroxyvitamin D (25OHD) levels (-0.43 SD of standardized natural log-transformed 25OHD per A allele; p value = 1.5 × 10-88). The effect on 25OHD was four times larger and independent of the effect of a previously described common variant near CYP2R1. By analyzing 8,711 individuals, we showed that heterozygote carriers of this low-frequency variant have an increased risk of vitamin D insufficiency (odds ratio [OR] = 2.2, 95% confidence interval [CI] = 1.78-2.78, p = 1.26 × 10-12). Individuals carrying one copy of this variant also had increased odds of multiple sclerosis (OR = 1.4, 95% CI = 1.19-1.64, p = 2.63 × 10-5) in a sample of 5,927 case and 5,599 control subjects. In conclusion, we describe a low-frequency CYP2R1 coding variant that exerts the largest effect upon 25OHD levels identified to date in the general European population and implicates vitamin D in the etiology of multiple sclerosis.

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.