Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk.

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.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Dnmt3a-mutated clonal hematopoiesis promotes osteoporosis.

Osteoporosis is caused by an imbalance of osteoclasts and osteoblasts, occurring in close proximity to hematopoietic cells in the bone marrow. Recurrent somatic mutations that lead to an expanded population of mutant blood cells is termed clonal hematopoiesis of indeterminate potential (CHIP). Analyzing exome sequencing data from the UK Biobank, we found CHIP to be associated with increased incident osteoporosis diagnoses and decreased bone mineral density. In murine models, hematopoietic-specific mutations in Dnmt3a, the most commonly mutated gene in CHIP, decreased bone mass via increased osteoclastogenesis. Dnmt3a-/- demethylation opened chromatin and altered activity of inflammatory transcription factors. Bone loss was driven by proinflammatory cytokines, including Irf3-NF-κB-mediated IL-20 expression from Dnmt3a mutant macrophages. Increased osteoclastogenesis due to the Dnmt3a mutations was ameliorated by alendronate or IL-20 neutralization. These results demonstrate a novel source of osteoporosis-inducing inflammation.

Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

Development of a polygenic risk score to improve screening for fracture risk: A genetic risk prediction study.

BACKGROUND: Since screening programs identify only a small proportion of the population as eligible for an intervention, genomic prediction of heritable risk factors could decrease the number needing to be screened by removing individuals at low genetic risk. We therefore tested whether a polygenic risk score for heel quantitative ultrasound speed of sound (SOS)-a heritable risk factor for osteoporotic fracture-can identify low-risk individuals who can safely be excluded from a fracture risk screening program. METHODS AND FINDINGS: A polygenic risk score for SOS was trained and selected in 2 separate subsets of UK Biobank (comprising 341,449 and 5,335 individuals). The top-performing prediction model was termed "gSOS", and its utility in fracture risk screening was tested in 5 validation cohorts using the National Osteoporosis Guideline Group clinical guidelines (N = 10,522 eligible participants). All individuals were genome-wide genotyped and had measured fracture risk factors. Across the 5 cohorts, the average age ranged from 57 to 75 years, and 54% of studied individuals were women. The main outcomes were the sensitivity and specificity to correctly identify individuals requiring treatment with and without genetic prescreening. The reference standard was a bone mineral density (BMD)-based Fracture Risk Assessment Tool (FRAX) score. The secondary outcomes were the proportions of the screened population requiring clinical-risk-factor-based FRAX (CRF-FRAX) screening and BMD-based FRAX (BMD-FRAX) screening. gSOS was strongly correlated with measured SOS (r2 = 23.2%, 95% CI 22.7% to 23.7%). Without genetic prescreening, guideline recommendations achieved a sensitivity and specificity for correct treatment assignment of 99.6% and 97.1%, respectively, in the validation cohorts. However, 81% of the population required CRF-FRAX tests, and 37% required BMD-FRAX tests to achieve this accuracy. Using gSOS in prescreening and limiting further assessment to those with a low gSOS resulted in small changes to the sensitivity and specificity (93.4% and 98.5%, respectively), but the proportions of individuals requiring CRF-FRAX tests and BMD-FRAX tests were reduced by 37% and 41%, respectively. Study limitations include a reliance on cohorts of predominantly European ethnicity and use of a proxy of fracture risk. CONCLUSIONS: Our results suggest that the use of a polygenic risk score in fracture risk screening could decrease the number of individuals requiring screening tests, including BMD measurement, while maintaining a high sensitivity and specificity to identify individuals who should be recommended an intervention.

A Metabolic Screen in Adolescents Reveals an Association Between Circulating Citrate and Cortical Bone Mineral Density.

Observations that insulin and adiponectin levels are related to cortical bone size in adolescents, independently of body composition, suggest factors related to fat metabolism directly influence skeletal development. To explore this question, we examined associations between a metabolic screen focusing on fat metabolism, and peripheral quantitative computed tomography (pQCT) measures of the mid-tibia, in 15-year-olds from the Avon Longitudinal Study of Parents and Children. Metabolic profiles were generated by proton nuclear magnetic resonance spectroscopy, from blood samples obtained at the same time as pQCT scans. Ordinary least squares linear regression was used to investigate relationships between metabolic measures and periosteal circumference (PC), cortical thickness (CT), and cortical bone mineral density (BMDC ). Metabolic profiles yielded 22 independent components following principal component analysis (PCA), giving a Bonferroni-adjusted threshold for statistical significance of p = 0.002. Data were available in 1121 subjects (487 males), mean age 15 years. Several metabolites related to lipid and cholesterol metabolism were associated with PC, CT, and BMDC after adjustment for age, sex, and Tanner stage. After additional adjustment for height, fat, and lean mass, only the association between citrate and BMDC remained below the Bonferroni-significant threshold (ß = -0.14 [-0.18, -0.09]) (ß represents a standardized coefficient). Citrate also showed evidence of association with PC (ß = 0.06 [0.03, 0.10]) and strength strain index (SSI; ß = 0.04 [0.01, 0.08]). Subsequently, we investigated whether these relationships were explained by increased bone resorption. Citrate was strongly related to serum ß-C-telopeptides of type I collagen (ß-CTX) (ß = 0.20 [0.16, 0.23]). After additional adjustment for ß-CTX the above associations between citrate and BMDC (ß = -0.04 [-0.08, 0.01]), PC (ß = 0.03 [-0.01, 0.07]) and SSI (ß = 0.03 [-0.01, 0.07]) were no longer observed. We conclude that in adolescents, circulating levels of citrate are inversely related to BMDC and positively related to PC, reflecting associations with higher bone turnover. Further studies are justified to elucidate possible contributions of citrate, a constituent of bone matrix, to bone resorption and cortical density. © 2019 American Society for Bone and Mineral Research.

Genome-wide association study of extreme high bone mass: Contribution of common genetic variation to extreme BMD phenotypes and potential novel BMD-associated genes.

BACKGROUND: Generalised high bone mass (HBM), associated with features of a mild skeletal dysplasia, has a prevalence of 0.18% in a UK DXA-scanned adult population. We hypothesized that the genetic component of extreme HBM includes contributions from common variants of small effect and rarer variants of large effect, both enriched in an extreme phenotype cohort. METHODS: We performed a genome-wide association study (GWAS) of adults with either extreme high or low BMD. Adults included individuals with unexplained extreme HBM (n = 240) from the UK with BMD Z-scores ≥+3.2, high BMD females from the Anglo-Australasian Osteoporosis Genetics Consortium (AOGC) (n = 1055) with Z-scores +1.5 to +4.0 and low BMD females also part of AOGC (n = 900), with Z-scores -1.5 to -4.0. Following imputation, we tested association between 6,379,332 SNPs and total hip and lumbar spine BMD Z-scores. For potential target genes, we assessed expression in human osteoblasts and murine osteocytes. RESULTS: We observed significant enrichment for associations with established BMD-associated loci, particularly those known to regulate endochondral ossification and Wnt signalling, suggesting that part of the genetic contribution to unexplained HBM is polygenic. Further, we identified associations exceeding genome-wide significance between BMD and four loci: two established BMD-associated loci (5q14.3 containing MEF2C and 1p36.12 containing WNT4) and two novel loci: 5p13.3 containing NPR3 (rs9292469; minor allele frequency [MAF] = 0.33%) associated with lumbar spine BMD and 11p15.2 containing SPON1 (rs2697825; MAF = 0.17%) associated with total hip BMD. Mouse models with mutations in either Npr3 or Spon1 have been reported, both have altered skeletal phenotypes, providing in vivo validation that these genes are physiologically important in bone. NRP3 regulates endochondral ossification and skeletal growth, whilst SPON1 modulates TGF-ß regulated BMP-driven osteoblast differentiation. Rs9292469 (downstream of NPR3) also showed some evidence for association with forearm BMD in the independent GEFOS sample (n = 32,965). We found Spon1 was highly expressed in murine osteocytes from the tibiae, femora, humeri and calvaria, whereas Npr3 expression was more variable. CONCLUSION: We report the most extreme-truncate GWAS of BMD performed to date. Our findings, suggest potentially new anabolic bone regulatory pathways that warrant further study.

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis.

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

Variants in the fetal genome near FLT1 are associated with risk of preeclampsia.

Preeclampsia, which affects approximately 5% of pregnancies, is a leading cause of maternal and perinatal death. The causes of preeclampsia remain unclear, but there is evidence for inherited susceptibility. Genome-wide association studies (GWAS) have not identified maternal sequence variants of genome-wide significance that replicate in independent data sets. We report the first GWAS of offspring from preeclamptic pregnancies and discovery of the first genome-wide significant susceptibility locus (rs4769613; P = 5.4 × 10-11) in 4,380 cases and 310,238 controls. This locus is near the FLT1 gene encoding Fms-like tyrosine kinase 1, providing biological support, as a placental isoform of this protein (sFlt-1) is implicated in the pathology of preeclampsia. The association was strongest in offspring from pregnancies in which preeclampsia developed during late gestation and offspring birth weights exceeded the tenth centile. An additional nearby variant, rs12050029, associated with preeclampsia independently of rs4769613. The newly discovered locus may enhance understanding of the pathophysiology of preeclampsia and its subtypes.

Incorporating Known Genetic Variants Does Not Improve the Accuracy of PSA Testing to Identify High Risk Prostate Cancer on Biopsy.

INTRODUCTION: Prostate-specific antigen (PSA) testing is a widely accepted screening method for prostate cancer, but with low specificity at thresholds giving good sensitivity. Previous research identified four single nucleotide polymorphisms (SNPs) principally associated with circulating PSA levels rather than with prostate cancer risk (TERT rs2736098, FGFR2 rs10788160, TBX3 rs11067228, KLK3 rs17632542). Removing the genetic contribution to PSA levels may improve the ability of the remaining biologically-determined variation in PSA to discriminate between high and low risk of progression within men with identified prostate cancer. We investigate whether incorporating information on the PSA-SNPs improves the discrimination achieved by a single PSA threshold in men with raised PSA levels. MATERIALS AND METHODS: Men with PSA between 3-10 ng/mL and histologically-confirmed prostate cancer were categorised as high or low risk of progression (Low risk: Gleason score≤6 and stage T1-T2a; High risk: Gleason score 7-10 or stage T2C). We used the combined genetic effect of the four PSA-SNPs to calculate a genetically corrected PSA risk score. We calculated the Area under the Curve (AUC) to determine how well genetically corrected PSA risk scores distinguished men at high risk of progression from low risk men. RESULTS: The analysis includes 868 men with prostate cancer (Low risk: 684 (78.8%); High risk: 184 (21.2%)). Receiver operating characteristic (ROC) curves indicate that including the 4 PSA-SNPs does not improve the performance of measured PSA as a screening tool for high/low risk prostate cancer (measured PSA level AUC = 59.5% (95% CI: 54.7,64.2) vs additionally including information from the 4 PSA-SNPs AUC = 59.8% (95% CI: 55.2,64.5) (p-value = 0.40)). CONCLUSION: We demonstrate that genetically correcting PSA for the combined genetic effect of four PSA-SNPs, did not improve discrimination between high and low risk prostate cancer in men with raised PSA levels (3-10 ng/mL). Replication and gaining more accurate estimates of the effects of the 4 PSA-SNPs and additional variants associated with PSA levels and not prostate cancer could be obtained from subsequent GWAS from larger prospective studies.

The effects of height and BMI on prostate cancer incidence and mortality: a Mendelian randomization study in 20,848 cases and 20,214 controls from the PRACTICAL consortium.

BACKGROUND: Epidemiological studies suggest a potential role for obesity and determinants of adult stature in prostate cancer risk and mortality, but the relationships described in the literature are complex. To address uncertainty over the causal nature of previous observational findings, we investigated associations of height- and adiposity-related genetic variants with prostate cancer risk and mortality. METHODS: We conducted a case-control study based on 20,848 prostate cancers and 20,214 controls of European ancestry from 22 studies in the PRACTICAL consortium. We constructed genetic risk scores that summed each man's number of height and BMI increasing alleles across multiple single nucleotide polymorphisms robustly associated with each phenotype from published genome-wide association studies. RESULTS: The genetic risk scores explained 6.31 and 1.46% of the variability in height and BMI, respectively. There was only weak evidence that genetic variants previously associated with increased BMI were associated with a lower prostate cancer risk (odds ratio per standard deviation increase in BMI genetic score 0.98; 95% CI 0.96, 1.00; p = 0.07). Genetic variants associated with increased height were not associated with prostate cancer incidence (OR 0.99; 95% CI 0.97, 1.01; p = 0.23), but were associated with an increase (OR 1.13; 95 % CI 1.08, 1.20) in prostate cancer mortality among low-grade disease (p heterogeneity, low vs. high grade <0.001). Genetic variants associated with increased BMI were associated with an increase (OR 1.08; 95 % CI 1.03, 1.14) in all-cause mortality among men with low-grade disease (p heterogeneity = 0.03). CONCLUSIONS: We found little evidence of a substantial effect of genetically elevated height or BMI on prostate cancer risk, suggesting that previously reported observational associations may reflect common environmental determinants of height or BMI and prostate cancer risk. Genetically elevated height and BMI were associated with increased mortality (prostate cancer-specific and all-cause, respectively) in men with low-grade disease, a potentially informative but novel finding that requires replication.

Associations of vitamin D pathway genes with circulating 25-hydroxyvitamin-D, 1,25-dihydroxyvitamin-D, and prostate cancer: a nested case-control study.

PURPOSE: Vitamin D pathway single nucleotide polymorphisms (SNPs) are potentially useful proxies for investigating whether circulating vitamin D metabolites [total 25-hydroxyvitamin-D, 25(OH)D; 1,25-dihydroxyvitamin, 1,25(OH)2D] are causally related to prostate cancer. We investigated associations of sixteen SNPs across seven genes with prostate-specific antigen-detected prostate cancer. METHODS: In a nested case-control study (within the ProtecT trial), we estimated odds ratios and 95 % confidence intervals (CIs) quantifying associations between SNPs and prostate cancer. Subgroup analyses investigated whether associations were stronger in men who had high/low sun exposure [a proxy for 25(OH)D]. We quantified associations of SNPs with stage (T1-T2/T3-T4) and grade (<7/≥7). Multiple variant scores included SNPs encoding proteins involved in 25(OH)D synthesis and metabolism. RESULTS: We included 1,275 prostate cancer cases (141 locally advanced, 385 high grades) and 2,062 healthy controls. Vitamin D-binding protein SNPs were associated with prostate cancer (rs4588-A: OR 1.20, CI 1.01, 1.41, p = 0.04; rs7041-T: OR 1.19, CI 1.02, 1.38, p = 0.03). Low 25(OH)D metabolism score was associated with high (vs low) grade (OR 0.76, CI 0.63, 0.93, p = 0.01); there was a similar association of its component variants: rs6013897-A in CYP24A1 (OR 0.78, CI 0.60, 1.01, p = 0.06) and rs10877012-T in CYP27B1 (OR 0.80, CI 0.63, 1.02, p = 0.07). There was no evidence that associations differed by level of sun exposure. CONCLUSION: We found some evidence that vitamin D pathway SNPs were associated with prostate cancer risk and grade, but not stage. There was no evidence of an association in men with deficient vitamin D (measured by having low sun exposure).

Cis and trans effects of human genomic variants on gene expression.

Gene expression is a heritable cellular phenotype that defines the function of a cell and can lead to diseases in case of misregulation. In order to detect genetic variations affecting gene expression, we performed association analysis of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with gene expression measured in 869 lymphoblastoid cell lines of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort in cis and in trans. We discovered that 3,534 genes (false discovery rate (FDR) = 5%) are affected by an expression quantitative trait locus (eQTL) in cis and 48 genes are affected in trans. We observed that CNVs are more likely to be eQTLs than SNPs. In addition, we found that variants associated to complex traits and diseases are enriched for trans-eQTLs and that trans-eQTLs are enriched for cis-eQTLs. As a variant affecting both a gene in cis and in trans suggests that the cis gene is functionally linked to the trans gene expression, we looked specifically for trans effects of cis-eQTLs. We discovered that 26 cis-eQTLs are associated to 92 genes in trans with the cis-eQTLs of the transcriptions factors BATF3 and HMX2 affecting the most genes. We then explored if the variation of the level of expression of the cis genes were causally affecting the level of expression of the trans genes and discovered several causal relationships between variation in the level of expression of the cis gene and variation of the level of expression of the trans gene. This analysis shows that a large sample size allows the discovery of secondary effects of human variations on gene expression that can be used to construct short directed gene regulatory networks.

Genetic variation in prostate-specific antigen-detected prostate cancer and the effect of control selection on genetic association studies.

BACKGROUND: Only a minority of the genetic components of prostate cancer risk have been explained. Some observed associations of SNPs with prostate cancer might arise from associations of these SNPs with circulating prostate-specific antigen (PSA) because PSA values are used to select controls. METHODS: We undertook a genome-wide association study (GWAS) of screen-detected prostate cancer (ProtecT: 1,146 cases and 1,804 controls); meta-analyzed the results with those from the previously published UK Genetic Prostate Cancer Study (1,854 cases and 1,437 controls); investigated associations of SNPs with prostate cancer using either "low" (PSA < 0.5 ng/mL) or "high" (PSA ≥ 3 ng/mL, biopsy negative) PSA controls; and investigated associations of SNPs with PSA. RESULTS: The ProtecT GWAS confirmed previously reported associations of prostate cancer at three loci: 10q11.23, 17q24.3, and 19q13.33. The meta-analysis confirmed associations of prostate cancer with SNPs near four previously identified loci (8q24.21,10q11.23, 17q24.3, and 19q13.33). When comparing prostate cancer cases with low PSA controls, alleles at genetic markers rs1512268, rs445114, rs10788160, rs11199874, rs17632542, rs266849, and rs2735839 were associated with an increased risk of prostate cancer, but the effect-estimates were attenuated to the null when using high PSA controls (Pheterogeneity in effect-estimates < 0.04). We found a novel inverse association of rs9311171-T with circulating PSA. CONCLUSIONS: Differences in effect-estimates for prostate cancer observed when comparing low versus high PSA controls may be explained by associations of these SNPs with PSA. IMPACT: These findings highlight the need for inferences from genetic studies of prostate cancer risk to carefully consider the influence of control selection criteria.

Does bone resorption stimulate periosteal expansion? A cross-sectional analysis of ß-C-telopeptides of type I collagen (CTX), genetic markers of the RANKL pathway, and periosteal circumference as measured by pQCT.

We hypothesized that bone resorption acts to increase bone strength through stimulation of periosteal expansion. Hence, we examined whether bone resorption, as reflected by serum ß-C-telopeptides of type I collagen (CTX), is positively associated with periosteal circumference (PC), in contrast to inverse associations with parameters related to bone remodeling such as cortical bone mineral density (BMDC ). CTX and mid-tibial peripheral quantitative computed tomography (pQCT) scans were available in 1130 adolescents (mean age 15.5 years) from the Avon Longitudinal Study of Parents and Children (ALSPAC). Analyses were adjusted for age, gender, time of sampling, tanner stage, lean mass, fat mass, and height. CTX was positively related to PC (ß=0.19 [0.13, 0.24]) (coefficient=SD change per SD increase in CTX, 95% confidence interval)] but inversely associated with BMDC (ß=-0.46 [-0.52,-0.40]) and cortical thickness [ß=-0.11 (-0.18, -0.03)]. CTX was positively related to bone strength as reflected by the strength-strain index (SSI) (ß=0.09 [0.03, 0.14]). To examine the causal nature of this relationship, we then analyzed whether single-nucleotide polymorphisms (SNPs) within key osteoclast regulatory genes, known to reduce areal/cortical BMD, conversely increase PC. Fifteen such genetic variants within or proximal to genes encoding receptor activator of NF-κB (RANK), RANK ligand (RANKL), and osteoprotegerin (OPG) were identified by literature search. Six of the 15 alleles that were inversely related to BMD were positively related to CTX (p<0.05 cut-off) (n=2379). Subsequently, we performed a meta-analysis of associations between these SNPs and PC in ALSPAC (n=3382), Gothenburg Osteoporosis and Obesity Determinants (GOOD) (n=938), and the Young Finns Study (YFS) (n=1558). Five of the 15 alleles that were inversely related to BMD were positively related to PC (p<0.05 cut-off). We conclude that despite having lower BMD, individuals with a genetic predisposition to higher bone resorption have greater bone size, suggesting that higher bone resorption is permissive for greater periosteal expansion.

Fraction of exhaled nitric oxide values in childhood are associated with 17q11.2-q12 and 17q12-q21 variants.

BACKGROUND: The fraction of exhaled nitric oxide (Feno) value is a biomarker of eosinophilic airway inflammation and is associated with childhood asthma. Identification of common genetic variants associated with childhood Feno values might help to define biological mechanisms related to specific asthma phenotypes. OBJECTIVE: We sought to identify the genetic variants associated with childhood Feno values and their relation with asthma. METHODS: Feno values were measured in children age 5 to 15 years. In 14 genome-wide association studies (N = 8,858), we examined the associations of approximately 2.5 million single nucleotide polymorphisms (SNPs) with Feno values. Subsequently, we assessed whether significant SNPs were expression quantitative trait loci in genome-wide expression data sets of lymphoblastoid cell lines (n = 1,830) and were related to asthma in a previously published genome-wide association data set (cases, n = 10,365; control subjects: n = 16,110). RESULTS: We identified 3 SNPs associated with Feno values: rs3751972 in LYR motif containing 9 (LYRM9; P = 1.97 × 10(-10)) and rs944722 in inducible nitric oxide synthase 2 (NOS2; P = 1.28 × 10(-9)), both of which are located at 17q11.2-q12, and rs8069176 near gasdermin B (GSDMB; P = 1.88 × 10(-8)) at 17q12-q21. We found a cis expression quantitative trait locus for the transcript soluble galactoside-binding lectin 9 (LGALS9) that is in linkage disequilibrium with rs944722. rs8069176 was associated with GSDMB and ORM1-like 3 (ORMDL3) expression. rs8069176 at 17q12-q21, but not rs3751972 and rs944722 at 17q11.2-q12, were associated with physician-diagnosed asthma. CONCLUSION: This study identified 3 variants associated with Feno values, explaining 0.95% of the variance. Identification of functional SNPs and haplotypes in these regions might provide novel insight into the regulation of Feno values. This study highlights that both shared and distinct genetic factors affect Feno values and childhood asthma.

Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity.

The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10(-8)) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure.

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10⁻¹4; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹°; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻9). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10⁻9). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60-0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.