Gene-gene and gene-lifestyle interactions of AKAP11, KCNMA1, PUM1, SPTBN1, and EPDR1 on osteoporosis risk in middle-aged adults

Utilize polygenic risk score to enhance fracture risk estimation and improve the performance of FRAX in patients with osteoporosis

This study examined the use of polygenic risk scores (PGS) in combination with the Fracture Risk Assessment Tool (FRAX) to enhance fragility fractures risk estimation in osteoporosis patients. Analyzing data from over 57,000 participants, PGS improved fracture risk estimation, especially for individuals with intermediate to low risks, allowing personalized preventive strategies.

INTRODUCTION

Osteoporosis and fragility fractures are multifactorial, with contributions from both clinical and genetic determinants. However, whether using polygenic risk scores (PGS) may enhance the risk estimation of osteoporotic fracture in addition to Fracture Risk Assessment Tool (FRAX) remains unknown. This study investigated the collective association of PGS and FRAX with fragility fracture.

METHODS

We conducted a cohort study from the Taiwan Precision Medicine Initiative (TPMI) at Taichung Veterans General Hospital, Taiwan. Genotyping was performed to compute PGS associated with bone mineral density (BMD). Phenome-wide association studies were executed to pinpoint phenotypes correlated with the PGS. Logistic regression analysis was conducted to ascertain factors associated with osteoporotic fractures.

RESULTS

Among all 57,257 TPMI participants, 3744 (904 men and 2840 women, with a mean age of 66.7) individuals had BMD testing, with 540 (14.42%) presenting with fractures. The 3744 individuals who underwent BMD testing were categorized into four quartiles (Q1-Q4) based on PGS; 540 (14.42%) presented with fractures. Individuals with PGS-Q1 exhibited lower BMD, a higher prevalence of major fractures, and elevated FRAX-major and FRAX-hip than those with PGS-Q4. PGS was associated with major fractures after adjusting age, sex, and FRAX scores. Notably, the risk of major fractures (PGS-Q1 vs. Q4) was significantly higher in the subgroups of FRAX-major scores < 10% and 10-20%, but not in participants with a FRAX-major score ≧ 20%.

CONCLUSIONS

Our study highlights the potential of PGS to augment fracture risk estimation in conjunction with FRAX, particularly in individuals with middle to low risks. Incorporating genetic testing could empower physicians to tailor personalized preventive strategies for osteoporosis.

Impact of Lung-Related Polygenic Risk Scores on Chronic Obstructive Pulmonary Disease Risk and Their Interaction with w-3 Fatty Acid Intake in Middle-Aged and Elderly Individuals

Chronic obstructive pulmonary disease (COPD) is a complex, progressive respiratory disorder with persistent airflow limitation and tissue destruction. We aimed to explore the genetic impact of COPD and its interaction with nutrient intake in 8840 middle-aged and elderly individuals from the Ansan/Ansung cohorts. Participants were diagnosed with COPD if the ratio of forced expiratory volume in 1 s (FEV1) to forced vital capacity (FVC) was less than 0.7 using spirometry, and if they were previously diagnosed with COPD by a physician. Genome-wide association studies (GWAS) were performed to screen for genetic variants associated with COPD risk. Among them, we selected the genetic variants that exhibited interactions using the generalized multifactor dimensionality reduction (GMDR) method. The polygenic risk score (PRS) was computed by summing the number of risk alleles in the SNP-SNP interaction models that adhered to specific rules. Subsequently, participants were categorized into low-PRS, medium-PRS, and high-PRS groups. The participants with COPD exhibited significantly lower FEV1/FVC ratios (0.64) than those without COPD (0.82). It was positively associated with inflammation markers (serum C-reactive protein and white blood cell levels). A higher proportion of COPD participants were smokers and engaged in regular exercise. The 5-SNP model consisted of FAM13A_rs1585258, CAV1_rs1997571, CPD_rs719601, PEPD_rs10405598, and ITGA1_rs889294, and showed a significant association with COPD risk (p < 0.001). Participants in the high-PRS group of this model had a 2.2-fold higher risk of COPD than those in the low-PRS group after adjusting for covariates. The PRS interacted with w-3 fatty acid intake and exercise, thus influencing the risk of COPD. There was an increase in COPD incidence among individuals with a higher PRS, particularly those with low consumption of w-3 fatty acid and engaged in high levels of exercise. In conclusion, adults with a high-PRS are susceptible to COPD risk, and w-3 fatty acid intake and exercise may impact the risk of developing COPD, potentially applying to formulate precision medicines to prevent COPD.

A Prediction Model for Osteoporosis Risk Using a Machine-Learning Approach and Its Validation in a Large Cohort

BACKGROUND

Osteoporosis develops in the elderly due to decreased bone mineral density (BMD), potentially increasing bone fracture risk. However, the BMD is not regularly measured in a clinical setting. This study aimed to develop a good prediction model for the osteoporosis risk using a machine learning (ML) approach in adults over 40 years in the Ansan/Anseong cohort and the association of predicted osteoporosis risk with a fracture in the Health Examinees (HEXA) cohort.

METHODS

The 109 demographic, anthropometric, biochemical, genetic, nutrient, and lifestyle variables of 8,842 participants were manually selected in an Ansan/Anseong cohort and included in the ML algorithm. The polygenic risk score (PRS) of osteoporosis was generated with a genome-wide association study and added for the genetic impact of osteoporosis. Osteoporosis was defined with < -2.5 T scores of the tibia or radius compared to people in their 20s-30s. They were divided randomly into the training (n = 7,074) and test (n = 1,768) sets-Pearson's correlation between the predicted osteoporosis risk and fracture in the HEXA cohort.

RESULTS

XGBoost, deep neural network, and random forest generated the prediction model with a high area under the curve (AUC, 0.86) of the receiver operating characteristic (ROC) with 10, 15, and 20 features; the prediction model by XGBoost had the highest AUC of ROC, high accuracy and k-fold values (> 0.85) in 15 features among seven ML approaches. The model included the genetic factor, genders, number of children and breastfed children, age, residence area, education, seasons to measure, height, smoking status, hormone replacement therapy, serum albumin, hip circumferences, vitamin B6 intake, and body weight. The prediction models for women alone were similar to those for both genders, with lower accuracy. When the prediction model was applied to the HEXA study, the correlation between the fracture incidence and predicted osteoporosis risk was significant but weak (r = 0.173, P < 0.001).

CONCLUSION

The prediction model for osteoporosis risk generated by XGBoost can be applied to estimate osteoporosis risk. The biomarkers can be considered for enhancing the prevention, detection, and early therapy of osteoporosis risk in Asians.

Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity

Obesity has become a severe public health challenge globally. The present study aimed to identify separate and interactive dietary, genetic, and other factors that increase the risk of obesity as measured by body fat (BF) mass. We utilized a genome-wide association study to identify genetic variants associated with high fat mass (obesity; n = 10,502) and combined them to generate polygenic risk scores (PRS) of genetic variants interacting with each other in adults aged over 40 while excluding body-fat-related diseases in a city-hospital-based cohort (n = 53,828). It was validated in Ansan/Ansung plus rural cohorts (n = 13,007). We then evaluated dietary and lifestyle factors in subjects to assess what factors might help overcome a genetic propensity for higher BF. The three-SNP model included brain-derived neurotrophic factor (BDNF)_rs6265, fat-mass- and obesity-associated protein (FTO)_rs1421085, and SEC16B_rs509325. The genes with the minor alleles of ADCY3_rs6545790 and BAIAP2_rs35867081 increased their gene expression in the visceral and subcutaneous adipocytes, but their gene expression decreased in the hypothalamus in eQTL analysis. In the three-SNP model, the PRS was associated with BF mass by 1.408 and 1.396 times after adjusting covariates 1 (age, gender, survey year, residence area, education, and income) and 2 (covariates in model 1 plus energy intake, alcohol intake, regular exercise, and smoking status), respectively. However, when separating subjects by PRS of the three-SNP model, a plant-based diet was the most significant factor associated with low BF, followed by high-protein diets and lower energy intakes. They could offset the effects of high genetic risk for high BF. In conclusion, modulating nutrient intakes might overcome a high genetic risk for obesity. Dietary choices favoring more plant-based and higher-protein foods might help prevent increased BF in Asians and potentially people of other ethnicities with high polygenetic risk scores.

Identification of recurrent variants implicated in disease in bicuspid aortic valve patients through whole-exome sequencing

Bicuspid aortic valve (BAV) is the most common congenital heart defect in human beings, with an estimated prevalence in the general population of between 0.5 and 2%. Moreover, BAV is the most common cause of aortic stenosis in the pediatric population. Patients with BAV may have no symptoms for life, and some of them may progress to aortic stenosis. Genetic factors increase the susceptibility and development of BAV. However, the pathogenesis and BAV are still unclear, and more genetic variants are still needed for elucidating the molecular mechanism and stratification of patients. The present study carried out screening of variants implicated in disease in BAV patients. The whole-exome sequencing (WES) was performed in 20 BAV patients and identified 40 different heterozygous missense mutations in 36 genes (MIB2, FAAH, S100A1, RGS16, MAP3K19, NEB, TTN, TNS1, CAND2, CCK, KALRN, ATP10D, SLIT3, ROS1, FABP7, NUP205, IL11RA, NPR2, COL5A1, CUBN, JMJD1C, ANXA7, TRIM8, LGR4, TPCN2, APOA5, GPR84, LRP1, NCOR2, AKAP11, ESRRB, NGB, AKAP13, WWOX, KCNJ12, ARHGEF1). The mutations in these genes were identified as recurrent variants implicated in disease by in silico prediction tool analysis. Nine genes (MIB2, S100A1, TTN, CCK, NUP205, LGR4, NCOR2, ESRRB, and WWOX) among the 36 genes were identified as variants implicated in disease via unanimous agreement of in silico prediction tool analysis and sequenced in an independent cohort of 137 BAV patients to validate the results of WES. BAV patients carrying these variants demonstrated reduced left ventricular ejection fractions (LVEF) (63.8 ± 7.5% vs. 58.4 ± 5.2%, P < 0.001) and larger calcification volume [(1129.3 ± 154) mm³ vs. (1261.8 ± 123) mm³, P < 0.001]. The variants in TTN, NUP205 and NCOR2 genes are significantly associated with reduced LVEF, and the variants in S100A1, LGR4, ESRRB, and WWOX genes are significantly associated with larger calcification volume. We identified a panel of recurrent variants implicated in disease in genes related to the pathogenesis of BAV. Our data speculate that these variants are promising markers for risk stratification of BAV patients with increased susceptibility to aortic stenosis.

Protection against Osteoporosis by Fermented Mulberry Vinegar Supplementation via Inhibiting Osteoclastic Activity in Ovariectomized Rats and Osteoclastic Cells

Menopause increases the osteoporosis risk, to which phytoestrogen intake can be beneficial. This study hypothesized that mulberry vinegar had a preventive effect on osteoporosis by decreasing osteoclastic activity. The hypothesis was tested in ovariectomized (OVX) rats and RANKL-differentiated osteoclast cells. OVX rats were given 0(OVX-CON), 0.5(OVX-MVL), 1(OVX-MVM), and 2(OVX-MVH) fermented mulberry vinegar (MV) mL/kg body weight (BW) daily for 12 weeks. Sham-operated rats had no MV supplementation (Normal-CON). The osteoporosis-related biomarkers were measured, and Micro-CT determined the bone mass of the femur. RANKL-differentiated Raw 264.7 cells were treated with MV (0–100 μg/mL). The cell viability, osteoporosis-related mRNA expression, and protein contents were measured. MV contained Acetobacter pasteurianus (7.31 log CFU/mL), citric acid (106 mg/mL), lactic acid (19.2 mg/mL), acetic acid (15.0 mg/mL), and rutin (0.36 mg/mL). OVX-MVM elevated the serum 17β-estradiol concentration similar to the Normal-CON group, but it did not prevent the decrease in uterine weight. OVX-MVM prevented the increase in osteoclastic-related parameters, including cathepsin K(CtsK), receptor activator of NF-κB ligand (RANKL), and tartrate-resistant acid phosphatase (TRAP) in the circulation. OVX-MVH also lowered C-telopeptide of type Ⅰ collagen as much as the Normal-CON group (p < 0.05). By contrast, OVX-MVH increased the serum osteoprotegerin concentration, an inhibitor of osteoclasts, better than the Normal-CON group (p < 0.05). These changes were integrated to alter the bone mineral density (BMD) in Micro-CT analysis: OVX-MVM and OVX-MVH prevented BMD decrease after OVX as much as the Normal-CON. In RANKL-differentiated osteoclast cells, the MV treatment for 24 and 48 h decreased RANKL-induced differentiation in osteoclast cells dose-dependently up to 100 µg/mL. Its decrease was related to inhibiting the TRAP activity and reducing TRAP-positive multinucleated cells during the five-day administration of RANKL. MV treatments also decreased mRNA expression of osteoclast-related genes (TRAP, Ctsk, OSCAR, and NFATc1). MV suppressed the protein contents of NFATc1 and c-FOS-related osteoclast. In conclusion, MV intake (1 mg/kg bw) protected against BMD loss mainly by inhibiting the osteoclastic activity (RANKL/RANK/TRAP) in OVX rats. MV may develop as a functional food for anti-osteoporosis in menopausal women.

Polymorphisms in Genes Involved in Osteoblast Differentiation and Function Are Associated with Anthropometric Phenotypes in Spanish Women

Much of the genetic variance associated with osteoporosis is still unknown. Bone mineral density (BMD) is the main predictor of osteoporosis risk, although other anthropometric phenotypes have recently gained importance. The aim of this study was to analyze the association of SNPs in genes involved in osteoblast differentiation and function with BMD, body mass index (BMI), and waist (WC) and hip (HC) circumferences. Four genes that affect osteoblast differentiation and/or function were selected from among the differentially expressed genes in fragility hip fracture (FOXC1, CTNNB1, MEF2C, and EBF2), and an association study of four single-nucleotide polymorphisms (SNPs) was conducted in a cohort of 1001 women. Possible allelic imbalance was also studied for SNP rs87939 of the CTNNB1 gene. We found significant associations of SNP rs87939 of the CTNNB1 gene with LS-sBMD, and of SNP rs1366594 of the MEF2C gene with BMI, after adjustment for confounding variables. The SNP of the MEF2C gene also showed a significant trend to association with FN-sBMD (p = 0.009). A possible allelic imbalance was ruled out as no differences for each allele were detected in CTNNB1 expression in primary osteoblasts obtained from homozygous women. In conclusion, we demonstrated that two SNPs in the MEF2C and CTNNB1 genes, both implicated in osteoblast differentiation and/or function, are associated with BMI and LS-sBMD, respectively.

Interactions between Polygenic Risk Scores, Dietary Pattern, and Menarche Age with the Obesity Risk in a Large Hospital-Based Cohort

Obese Asians are more susceptible to metabolic diseases than obese Caucasians of the same body mass index (BMI). We hypothesized that the genetic variants associated with obesity risk interact with the lifestyles of middle-aged and elderly adults, possibly allowing the development of personalized interventions based on genotype. We aimed to examine this hypothesis in a large city hospital-based cohort in Korea. The participants with cancers, thyroid diseases, chronic kidney disease, or brain-related diseases were excluded. The participants were divided into case and control according to their BMI: ≥25 kg/m² (case; n = 17,545) and <25 kg/m² (control; n = 36,283). The genetic variants that affected obesity risk were selected using a genome-wide association study, and the genetic variants that interacted with each other were identified by generalized multifactor dimensionality reduction analysis. The selected genetic variants were confirmed in the Ansan/Ansung cohort, and polygenetic risk scores (PRS)−nutrient interactions for obesity risk were determined. A high BMI was associated with a high-fat mass (odds ratio (OR) = 20.71) and a high skeletal muscle-mass index (OR = 3.38). A high BMI was positively related to metabolic syndrome and its components, including lipid profiles, whereas the initial menstruation age was inversely associated with a high BMI (OR = 0.78). The best model with 5-SNPs included SEC16B_rs543874, DNAJC27_rs713586, BDNF_rs6265, MC4R_rs6567160, and GIPR_rs1444988703. The high PRS with the 5-SNP model was positively associated with an obesity risk of 1.629 (1.475–1.798) after adjusting for the covariates. The 5-SNP model interacted with the initial menstruation age, fried foods, and plant-based diet for BMI risk. The participants with a high PRS also had a higher obesity risk when combined with early menarche, low plant-based diet, and a high fried-food intake than in participants with late menarche, high plant-based diet, and low fried-food intake. In conclusion, people with a high PRS and earlier menarche age are recommended to consume fewer fried foods and a more plant-based diet to decrease obesity risk. This result can be applied to personalized nutrition for preventing obesity.

An emerging spectrum of variants and clinical features in KCNMA1-linked channelopathy

KCNMA1-linked channelopathy is an emerging neurological disorder characterized by heterogeneous and overlapping combinations of movement disorder, seizure, developmental delay, and intellectual disability. KCNMA1 encodes the BK K⁺ channel, which contributes to both excitatory and inhibitory neuronal and muscle activity. Understanding the basis of the disorder is an important area of active investigation; however, the rare prevalence has hampered the development of large patient cohorts necessary to establish genotype-phenotype correlations. In this review, we summarize 37 KCNMA1 alleles from 69 patients currently defining the channelopathy and assess key diagnostic and clinical hallmarks. At present, 3 variants are classified as gain-of-function with respect to BK channel activity, 14 loss-of-function, 15 variants of uncertain significance, and putative benign/VUS. Symptoms associated with these variants were curated from patient-provided information and prior publications to define the spectrum of clinical phenotypes. In this newly expanded cohort, seizures showed no differential distribution between patients harboring GOF and LOF variants, while movement disorders segregated by mutation type. Paroxysmal non-kinesigenic dyskinesia was predominantly observed among patients with GOF alleles of the BK channel, although not exclusively so, while additional movement disorders were observed in patients with LOF variants. Neurodevelopmental and structural brain abnormalities were prevalent in patients with LOF mutations. In contrast to mutations, disease-associated KCNMA1 single nucleotide polymorphisms were not predominantly related to neurological phenotypes but covered a wider set of peripheral physiological functions. Together, this review provides additional evidence exploring the genetic and biochemical basis for KCNMA1-linked channelopathy and summarizes the clinical repository of patient symptoms across multiple types of KCNMA1 gene variants.

βII spectrin (SPTBN1): biological function and clinical potential in cancer and other diseases

βII spectrin, the most common isoform of non-erythrocyte spectrin, is a cytoskeleton protein present in all nucleated cells. Interestingly, βII spectrin is essential for the development of various organs such as nerve, epithelium, inner ear, liver and heart. The functions of βII spectrin include not only establishing and maintaining the cell structure but also regulating a variety of cellular functions, such as cell apoptosis, cell adhesion, cell spreading and cell cycle regulation. Notably, βII spectrin dysfunction is associated with embryonic lethality and the DNA damage response. More recently, the detection of altered βII spectrin expression in tumors indicated that βII spectrin might be involved in the development and progression of cancer. Its mutations and disorders could result in developmental disabilities and various diseases. The versatile roles of βII spectrin in disease have been examined in an increasing number of studies; nonetheless, the exact mechanisms of βII spectrin are still poorly understood. Thus, we summarize the structural features and biological roles of βII spectrin and discuss its molecular mechanisms and functions in development, homeostasis, regeneration and differentiation. This review highlight the potential effects of βII spectrin dysfunction in cancer and other diseases, outstanding questions for the future investigation of therapeutic targets. The investigation of the regulatory mechanism of βII spectrin signal inactivation and recovery may bring hope for future therapy of related diseases.

Polygenetic-Risk Scores Related to Crystallin Metabolism Are Associated with Age-Related Cataract Formation and Interact with Hyperglycemia, Hypertension, Western-Style Diet, and Na Intake

Age-related cataract (ARC) development is associated with loss of crystalline lens transparency related to interactions between genetic and environmental factors. We hypothesized that polygenetic risk scores (PRS) of the selected genetic variants among the ARC-related genes might reveal significant genetic impacts on ARC risk, and the PRS might have gene–gene and gene–lifestyle interactions. We examined the hypothesis in 1972 and 39,095 subjects aged ≥50 years with and without ARC, respectively, in a large-scale hospital-based cohort study conducted from 2004 to 2013. Single nucleotide polymorphisms (SNPs) of the genes related to ARC risk were identified, and polygenetic risk scores (PRS) were generated based on the results of a generalized multifactor dimensionality reduction analysis. Lifestyle interactions with PRS were evaluated. The PRS derived from the best model included the following six SNPs related to crystallin metabolism: ULK4_rs1417380362, CRYAB_rs2070894, ACCN1_rs55785344, SSTR2_rs879419608, PTN_rs322348, and ICA1_rs200053781. The risk of ARC in the high-PRS group was 2.47-fold higher than in the low-PRS group after adjusting for confounders. Age, blood pressure, and glycemia interacted with PRS to influence the risk of ARC: the incidence of ARC was much higher in the elderly (≥65 years) and individuals with hypertension or hyperglycemia. The impact of PRS on ARC risk was greatest in middle-aged individuals with hypertension or hyperglycemia. Na, coffee, and a Western-style diet intake also interacted with PRS to influence ARC risk. ARC risk was higher in the high-PRS group than in the low-PRS group, and high Na intake, Western-style diet, and low coffee intake elevated its risk. In conclusion, ARC risk had a positive association with PRS related to crystallin metabolism. The genetic impact was greatest among those with high Na intake or hypertension. These results can be applied to precision nutrition interventions to prevent ARC.