Genome-Wide Association Study of Bone Mineral Density in Korean Men

Multiple polygenic scores improve bone mineral density prediction in an independent sample of Caucasian women

PURPOSE OF THE STUDY

To determine if multiple Genetic Risk Scores (GRSs) improve bone mineral density (BMD) prediction over single GRS in an independent sample of Caucasian women.

STUDY DESIGN

Based on summary statistics of four genome-wide association studies related to two osteoporosis-associated traits, namely BMD and heel quantitative ultrasound derived estimated BMD (eBMD), four GRSs were derived for 1205 individuals in the Genome-Wide Scan for Female Osteoporosis Gene Study. The effect of each GRS on BMD variation was assessed using multivariable linear regression, with conventional risk factors adjusted for. Next, the eBMD-related GRS that explained the most variance in BMD was selected to be entered into a multi-score model, along with the BMD-related GRS. Elastic net regularised regression was used to develop the multiscore model, which estimated the joint effect of two GRSs (GRS_BMD and GRS_eBMD) on BMD variation, after being adjusted for conventional risk factors.

RESULTS

With the same clinical risk factors having been adjusted for, the model that included GRS_BMD performed best by explaining 32.53% of the variance in BMD; the single-score model that included GRS_eBMD explained 34.03% of BMD variance. The model that includes both GRS_BMD and GRS_ eBMD, as well as the clinical risk factors, aggregately explained 35.05% in BMD variation. Compared with the single GRS models, the multiscore model explained significantly more variance in BMD.

CONCLUSIONS

The multipolygenic score model explained a considerable amount of BMD variation. Compared with single score models, multipolygenic score model provided significant improvement in explaining BMD variation.

A population-specific low-frequency variant of SLC22A12 (p.W258*) explains nearby genome-wide association signals for serum uric acid concentrations among Koreans

Prolonged hyperuricemia is a cause of gout and an independent risk factor for chronic health conditions including diabetes and chronic kidney diseases. Genome-wide association studies (GWASs) for serum uric acid (SUA) concentrations have repeatedly confirmed genetic loci including those encoding uric acid transporters such as ABCG2 and SLC9A2. However, many single nucleotide polymorphisms (SNPs) found in GWASs have been common variants with small effects and unknown functions. In addition, there is still much heritability to be explained. To identify the causative genetic variants for SUA concentrations in Korean subjects, we conducted a GWAS (1902 males) and validation study (2912 males and females) and found four genetic loci reaching genome-wide significance on chromosomes 4 (ABCG2) and 11 (FRMD8, EIF1AD and SLC22A12-NRXN2). Three loci on chromosome 11 were distributed within a distance of 1.3 megabases and they were in weak or moderate linkage disequilibrium (LD) states (r² = 0.02–0.68). In a subsequent association analysis on the GWAS loci of chromosome 11 using closely positioned markers derived from whole genome sequencing data, the most significant variant to be linked with the nearby GWAS signal was rs121907892 (c.774G>A, p.W258*) of the SLC22A12 gene. This variant, and each of the three GWAS SNPs, were in LD (r² = 0.06–0.32). The strength of association of SNPs with SUA concentration (negative logarithm of P-values) and the genetic distance (r² of LD) between rs121907892 and the surrounding SNPs showed a quantitative correlation. This variant has been found only in Korean and Japanese subjects and is known to lower the SUA concentration in the general population. Thus, this low-frequency variant, rs121907892, known to regulate SUA concentrations in previous studies, is responsible for the nearby GWAS signals.

Obesity and Bone Health Revisited: A Mendelian Randomization Study for Koreans

Epidemiologic studies have replicated positive associations between obesity and bone health, but their mechanisms are still debatable. We aimed to scrutinize an association between bone health and obesity using genetic instrumental variables (IVs) with the distinction of general versus abdominal obesity. We selected independent IVs of body mass index (BMI) and BMI-adjusted waist circumference (aWC, a proxy of a central fat distribution) by combining novel genomewide searches from the Korean Genome Epidemiology Study (KoGES) consortium and existing reports. We evaluated the associations of obesity indices with bone health measures for weight-bearing and non-weight-bearing bones, applying standard Mendelian randomization analyses. The IVs for BMI and aWC selected from KoGES cohort studies (n = 14,389) explained its own trait only, negating the mutual correlation at the phenotypic level. Two-stage least squares analyses using an independent cohort study (n = 2507, mean age = 44.4 years, men = 44.3%) showed that BMI but not aWC was positively associated with bone mineral density (BMD for weight-bearing bones: 0.063 ± 0.016 g/cm² per one standard deviation increase in BMI), implying the fat distribution might be neutral. The association was weaker for non-weight-bearing bones (BMI on BMD: 0.034 ± 0.011 g/cm² ), and for postmenopausal women the association was absent. Obesity increased both bone area and bone mineral content (BMC) to a lesser degree, but the increase in BMC was not evident for menopausal women. When we stratified the weight into lean body mass and fat mass, the increase in BMD was more evident for lean body mass, and fat mass showed a beneficial role only for men and premenopausal women. Our findings suggest that bone health might gain little from obesity, if any, through its added weight, and other means to prevent bone loss would be essential for postmenopausal women. © 2019 American Society for Bone and Mineral Research.

CYP1A1 gene polymorphisms modify the association between PM10 exposure and lung function

Genetic epidemiological studies have provided evidence that several genes modify the link between air pollution and lung function. We assessed whether the adverse impacts of particulate matter with an aerodynamic diameter ≤10 μm (PM₁₀) on lung function are modified by CYP1A1 gene polymorphisms in Korean adults. We used health check-up data from 1817 men, and the annual mean concentrations of ambient PM₁₀ estimated from the ambient data. Three single nucleotide polymorphisms (SNPs) of CYP1A1 were selected for our study. We identified significant CYP1A1 SNPs-by-PM₁₀ interactions for forced expiratory volume 1 s (FEV₁) and forced vital capacity (FVC) (all p_int < 0.05). Minor allele carriers of the SNPs were more susceptible to PM₁₀-induced FEV₁ and FVC reduction. The subgroup analysis of SNP genotypes showed that no significant association between PM₁₀ and FEV₁ or FVC was observed in homozygous reference genotype groups of all SNPs (all p_assoc > 0.05), whereas in heterozygous or homozygous alternate genotype groups, PM₁₀ was significantly associated with decreased FEV₁ (all p_assoc for FEV₁ < 0.05). The association between persistent exposure to PM₁₀ and lung function decline in Korean men may be determined in part by several functional variants of the CYP1A1 gene.

RANKL/RANK/OPG Polymorphisms and Heel Quantitative Ultrasound in Young Adults

BACKGROUND

The receptor activator of the nuclear factor-kappa B ligand (RANKL), the receptor activator of nuclear factor-kappa B (RANK), and the osteoprotegerin (OPG) signaling pathway play an important role in the regulation of bone remodeling and osteoclast differentiation. Quantitative ultrasound (QUS) is a relatively recent and noninvasive method providing structural information on microstructure, bone elasticity, and connectivity. However, in contrast to bone mineral density measurements, the possible association of the RANKL/RANK/OPG pathway with heel QUS has not been analyzed.

OBJECTIVES

The aim of this study was to assess, for the first time, the contribution of the RANKL/RANK/OPG pathway genes in the genetic background of heel QUS parameters.

METHODS

Ten single-nucleotide polymorphisms (SNPs) of RANKL (rs9594759, rs12585014, rs7988338, rs2148073), RANK (rs1805034, rs12458117, rs3018362), and OPG (rs4355801, rs3102735, rs2073618) were selected as genetic markers and genotyped using Open Array technology in 575 self-reported Caucasian individuals aged 18-25. Bone mass in the right calcaneus was estimated with QUS to obtain the broadband ultrasound attenuation (BUA) measurement (dB/MHz). Linear regression analyses were performed to test the possible association between the SNPs and BUA.

RESULTS

Linear regression analysis of all the tested SNPs revealed no significant association with the BUA parameter after adjusting for age, gender, weight, height, physical activity, and calcium intake. The lowest p-value was observed for the rs9594759 RANKL polymorphism and heel QUS (p = .06; b* = -.075, 95% CI [-0.960, 0.028]).

CONCLUSION

Our results suggest that the polymorphism of the RANKL, RANK, and OPG genes does not make a significant genetic contribution to heel ultrasound measurements in a population of young Caucasian adults. Further studies replicating the results in independent populations are needed to support these initial findings.