Abstract The aim of this international meta-analysis was to quantify the predictive value of BMI for incident fracture and relationship of this risk with age, sex, follow-up time, and BMD. A total of 1 667 922 men and women from 32 countries (63 cohorts), followed for a total of 16.0 million person-years were studied. 293 325 had FN BMD measured (2.2 million person-years follow-up). An extended Poisson model in each cohort was used to investigate relationships between WHO-defined BMI categories (Underweight: <18.5 kg/m2; Normal: 18.5-24.9 kg/m2; Overweight: 25.0-29.9 kg/m2; Obese I: 30.0-34.9 kg/m2; Obese II: ≥35.0 kg/m2) and risk of incident osteoporotic, major osteoporotic and hip fracture (HF). Inverse-variance weighted β-coefficients were used to merge the cohort-specific results. For the subset with BMD available, in models adjusted for age and follow-up time, the hazard ratio (95% CI) for HF comparing underweight with normal weight was 2.35 (2.10-2.60) in women and for men was 2.45 (1.90-3.17). Hip fracture risk was lower in overweight and obese categories compared to normal weight [obese II vs normal: women 0.66 (0.55-0.80); men 0.91 (0.66-1.26)]. Further adjustment for FN BMD T-score attenuated the increased risk associated with underweight [underweight vs normal: women 1.69 (1.47-1.96); men 1.46 (1.00-2.13)]. In these models, the protective effects of overweight and obesity were attenuated, and in both sexes, the direction of association reversed to higher fracture risk in Obese II category [Obese II vs Normal: women 1.24 (0.97-1.58); men 1.70 (1.06-2.75)]. Results were similar for other fracture outcomes. Underweight is a risk factor for fracture in both men and women regardless of adjustment for BMD. However, while overweight/obesity appeared protective in base models, they became risk factors after additional adjustment for FN BMD, particularly in the Obese II category. This effect in the highest BMI categories was of greater magnitude in men than women. These results will inform the second iteration of FRAX®.
Novel Genetic Determinants of Dental Maturation in Children
Dental occlusion requires harmonious development of teeth, jaws, and other elements of the craniofacial complex, which are regulated by environmental and genetic factors. We performed the first genome-wide association study (GWAS) on dental development (DD) using the Demirjian radiographic method. Radiographic assessments from participants of the Generation R Study (primary study population, N 1 = 2,793; mean age of 9.8 y) were correlated with ~30 million genetic variants while adjusting for age, sex, and genomic principal components (proxy for population stratification). Variants associated with DD at genome-wide significant level ( P < 5 × 10 −8 ) mapped to 16q12.2 ( IRX5) (lead variant rs3922616, B = 0.16; P = 2.2 × 10 −8 ). We used Fisher’s combined probability tests weighted by sample size to perform a meta-analysis ( N = 14,805) combining radiographic DD at a mean age of 9.8 y from Generation R with data from a previous GWAS ( N 2 = 12,012) on number of teeth (NT) in infants used as proxy of DD at a mean age of 9.8 y (including the ALSPAC and NFBC1966). This GWAS meta-analysis revealed 3 novel loci mapping to 7p15.3 ( IGF2BP3: P = 3.2 × 10 −8 ), 14q13.3 ( PAX9: P = 1.9 × 10 −8 ), and 16q12.2 ( IRX5: P = 1.2 × 10 −9 ) and validated 8 previously reported NT loci. A polygenic allele score constructed from these 11 loci was associated with radiographic DD in an independent Generation R set of children ( N = 703; B = 0.05, P = 0.004). Furthermore, profiling of the identified genes across an atlas of murine and human stem cells observed expression in the cells involved in the formation of bone and/or dental tissues (>0.3 frequency per kilobase of transcript per million mapped reads), likely reflecting functional specialization. Our findings provide biological insight into the polygenic architecture of the pediatric dental maturation process.
Gender and Geographic Origin as Determinants of Manuscript Publication Outcomes: JBMR® Bibliometric Analysis from 2017 to 2019
Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium ( n Stage1 : 111,666; n Stage2 : 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea ( PPM1J , EDEM3, ACP1, SPEG, EYA4, CYP1A1 , and ATXN2L ; P Stage1 <3.7×10 −7 ), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 ( P =5.4×10 −8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2 -knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.
Shared Genetic Risk Factors of Intracranial, Abdominal, and Thoracic Aneurysms
Background Intracranial aneurysms ( IAs ), abdominal aortic aneurysms ( AAAs ), and thoracic aortic aneurysms ( TAAs ) all have a familial predisposition. Given that aneurysm types are known to co‐occur, we hypothesized that there may be shared genetic risk factors for IAs , AAAs, and TAAs . Methods and Results We performed a mega‐analysis of 1000 Genomes Project‐imputed genome‐wide association study ( GWAS ) data of 4 previously published aneurysm cohorts: 2 IA cohorts (in total 1516 cases, 4305 controls), 1 AAA cohort (818 cases, 3004 controls), and 1 TAA cohort (760 cases, 2212 controls), and observed associations of 4 known IA , AAA , and/or TAA risk loci (9p21, 18q11, 15q21, and 2q33) with consistent effect directions in all 4 cohorts. We calculated polygenic scores based on IA ‐, AAA ‐, and TAA ‐associated SNP s and tested these scores for association to case‐control status in the other aneurysm cohorts; this revealed no shared polygenic effects. Similarly, linkage disequilibrium–score regression analyses did not show significant correlations between any pair of aneurysm subtypes. Last, we evaluated the evidence for 14 previously published aneurysm risk single‐nucleotide polymorphisms through collaboration in extended aneurysm cohorts, with a total of 6548 cases and 16 843 controls ( IA ) and 4391 cases and 37 904 controls ( AAA ), and found nominally significant associations for IA risk locus 18q11 near RBBP 8 to AAA (odds ratio [ OR ]=1.11; P =4.1×10 −5 ) and for TAA risk locus 15q21 near FBN 1 to AAA ( OR =1.07; P =1.1×10 −3 ). Conclusions Although there was no evidence for polygenic overlap between IAs , AAAs , and TAAs , we found nominally significant effects of two established risk loci for IAs and TAAs in AAAs . These two loci will require further replication.
A Meta-Analysis of the Association of Fracture Risk and Body Mass Index in Women
Genetic College Many genomic elements in humans are associated with behavior, including educational attainment. In a genome-wide association study including more than 100,000 samples, Rietveld et al. (p. 1467 , published online 30 May; see the Perspective by Flint and Munafò ) looked for genes related to educational attainment in Caucasians. Small genetic effects at three loci appeared to impact educational attainment.
A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11