Osteoporosis- and obesity-risk interrelationships: an epigenetic analysis of GWAS-derived SNPs at the developmental gene TBX15

A genome-wide association study of morphometric traits in dromedaries

BACKGROUND

Investigating genomic regions associated with morphometric traits in camels is valuable, because it allows a better understanding of adaptive and productive features to implement a sustainable management and a customised breeding program for dromedaries.

OBJECTIVES

With a genome-wide association study (GWAS) including 96 Iranian dromedaries phenotyped for 12 morphometric traits and genotyped-by-sequencing (GBS) with 14,522 SNPs, we aimed at identifying associated candidate genes.

METHODS

The association between SNPs and morphometric traits was investigated using a linear mixed model with principal component analysis (PCA) and kinship matrix.

RESULTS

With this approach, we detected 59 SNPs located in 37 candidate genes potentially associated to morphometric traits in dromedaries. The top associated SNPs were related to pin width, whither to pin length, height at whither, muzzle girth, and tail length. Interestingly, the results highlight the association between whither height, muzzle circumference, tail length, whither to pin length. The identified candidate genes were associated with growth, body size, and immune system in other species.

CONCLUSIONS

We identified three key hub genes in the gene network analysis including ACTB, SOCS1 and ARFGEF1. In the central position of gene network, ACTB was detected as the most important gene related to muscle function. With this initial GWAS using GBS on dromedary camels for morphometric traits, we show that this SNP panel can be effective for genetic evaluation of growth in dromedaries. However, we suggest a higher-density SNP array may greatly improve the reliability of the results.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Association of air pollutants and osteoporosis risk: The modifying effect of genetic predisposition

BACKGROUND

Limited studies have examined the association between air pollutants and osteoporosis incidence; however, the results are conflicting. We aimed to quantify the effects of selected air pollutants on osteoporosis risk and explore the modifying effect of genetic predisposition.

METHODS

A total of 422,955 subjects who did not have osteoporosis at baseline in the UK Biobank were included from 2006 to 2010. We conducted a Cox proportional hazards model with adjustment for covariates to examine the association between air pollutant scores and individual air pollutants and incident osteoporosis. Furthermore, a polygenic risk score (PRS) of osteoporosis was built and examined to determine whether genetic susceptibility modified the effect of air pollutants on osteoporosis. The relationship between air pollutants and osteoporosis was examined by using a restricted cubic spline (RCS) method.

RESULTS

After confounder adjustment, the results showed a remarkable increase in the risk of osteoporosis with each 10 unit increase in exposure to air pollution (hazard ratio: 1.06, 95 % confidence interval: 1.03-1.08), PM_2.5 (1.94, 1.52-2.48), NO₂ (1.06, 1.02-1.10), and NO_X (1.03, 1.01-1.04). However, no significant association was observed between PM₁₀ or PM_2.5-10 exposure and osteoporosis. Subjects with high air pollutant exposure levels and a high PRS had a noteworthy increase in osteoporosis risk compared to those with low air pollutant exposure levels and a low PRS. Air pollutants and genetic variants exerted additive effects on the risk of osteoporosis. Positive correlations were observed between osteoporosis and PM_2.5 (P < 0.001), NO₂ (P = 0.001), and NOx (P = 0.002) exposure.

CONCLUSIONS

Exposure to PM_2.5, NO₂ and NOx was associated with an increase in osteoporosis risk, and this effect was more pronounced in populations with high genetic risk. The association between PM_2.5, NO₂ and NOx exposure and osteoporosis is modified by genetic variations.

Epigenomic and Transcriptomic Prioritization of Candidate Obesity-Risk Regulatory GWAS SNPs

Concern about rising rates of obesity has prompted searches for obesity-related single nucleotide polymorphisms (SNPs) in genome-wide association studies (GWAS). Identifying plausible regulatory SNPs is very difficult partially because of linkage disequilibrium. We used an unusual epigenomic and transcriptomic analysis of obesity GWAS-derived SNPs in adipose versus heterologous tissues. From 50 GWAS and 121,064 expanded SNPs, we prioritized 47 potential causal regulatory SNPs (Tier-1 SNPs) for 14 gene loci. A detailed examination of seven loci revealed that four (CABLES1, PC, PEMT, and FAM13A) had Tier-1 SNPs positioned so that they could regulate use of alternative transcription start sites, resulting in different polypeptides being generated or different amounts of an intronic microRNA gene being expressed. HOXA11 and long noncoding RNA gene RP11-392O17.1 had Tier-1 SNPs in their 3' or promoter region, respectively, and strong preferences for expression in subcutaneous versus visceral adipose tissue. ZBED3-AS1 had two intragenic Tier-1 SNPs, each of which could contribute to mediating obesity risk through modulating long-distance chromatin interactions. Our approach not only revealed especially credible novel regulatory SNPs, but also helped evaluate previously highlighted obesity GWAS SNPs that were candidates for transcription regulation.

BTOB: Extending the Biased GWAS to Bivariate GWAS

In recent years, a number of literatures published large-scale genome-wide association studies (GWASs) for human diseases or traits while adjusting for other heritable covariate. However, it is known that these GWASs are biased, which may lead to biased genetic estimates or even false positives. In this study, we provide a method called "BTOB" which extends the biased GWAS to bivariate GWAS by integrating the summary association statistics from the biased GWAS and the GWAS for the adjusted heritable covariate. We employ the proposed BTOB method to analyze the summary association statistics from the large scale meta-GWASs for waist-to-hip ratio (WHR) and body mass index (BMI), and show that the proposed approach can help identify more susceptible genes compared with the corresponding univariate GWASs. Theoretical results and simulations also confirm the validity and efficiency of the proposed BTOB method.

Prioritization of Osteoporosis-Associated Genome-wide Association Study (GWAS) Single-Nucleotide Polymorphisms (SNPs) Using Epigenomics and Transcriptomics

Genetic risk factors for osteoporosis, a prevalent disease associated with aging, have been examined in many genome-wide association studies (GWASs). A major challenge is to prioritize transcription-regulatory GWAS-derived variants that are likely to be functional. Given the critical role of epigenetics in gene regulation, we have used an unusual epigenetics-based and transcription-based approach to identify some of the credible regulatory single-nucleotide polymorphisms (SNPs) relevant to osteoporosis from 38 reported bone mineral density (BMD) GWASs. Using Roadmap databases, we prioritized SNPs based upon their overlap with strong enhancer or promoter chromatin preferentially in osteoblasts relative to 12 heterologous cell culture types. We also required that these SNPs overlap open chromatin (Deoxyribonuclease I [DNaseI]-hypersensitive sites) and DNA sequences predicted to bind to osteoblast-relevant transcription factors in an allele-specific manner. From >50,000 GWAS-derived SNPs, we identified 14 novel and credible regulatory SNPs (Tier-1 SNPs) for osteoporosis risk. Their associated genes, BICC1, LGR4, DAAM2, NPR3, or HMGA2, are involved in osteoblastogenesis or bone homeostasis and regulate cell signaling or enhancer function. Four of these genes are preferentially expressed in osteoblasts. BICC1, LGR4, and DAAM2 play important roles in canonical Wnt signaling, a pathway critical for bone formation and repair. The transcription factors predicted to bind to the Tier-1 SNP-containing DNA sequences also have bone-related functions. We present evidence that some of the Tier-1 SNPs exert their effects on BMD risk indirectly through little-studied long noncoding RNA (lncRNA) genes, which may, in turn, control the nearby bone-related protein-encoding gene. Our study illustrates a method to identify novel BMD-related causal regulatory SNPs for future study and to prioritize candidate regulatory GWAS-derived SNPs, in general. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Epigenetic therapies of osteoporosis

The study of epigenetics reaches its 50th anniversary, however, its clinical application is gradually coming into the clinical setting. Osteoporosis is one of the major and widely diffused bone diseases. Pathogenic mechanisms at the epigenetic level may interfere with bone remodeling occurring during osteoporosis. Preclinical models were used to understand whether such events may interfere with the disease. Besides, observational clinical trials investigated epigenetic-related biomarkers. This effort leads to some epigenetic-related therapies in clinical trials for the treatment of osteoporosis. Bisphosphonates (BPs), target therapy blocking RANK/RANKL pathway, and anti-sclerostin antibody (SOST) are the main therapeutic approaches. However, future large trials will reveal whether epigenetic therapies of osteoporosis will remain a work in progress or data will become more robust in the real-world management of these frailty patients.