Gene-based association analysis identified novel genes associated with bone mineral density

Variance-components tests for genetic association with multiple interval-censored outcomes

Massive genetic compendiums such as the UK Biobank have become an invaluable resource for identifying genetic variants that are associated with complex diseases. Due to the difficulties of massive data collection, a common practice of these compendiums is to collect interval-censored data. One challenge in analyzing such data is the lack of methodology available for genetic association studies with interval-censored data. Genetic effects are difficult to detect because of their rare and weak nature, and often the time-to-event outcomes are transformed to binary phenotypes for access to more powerful signal detection approaches. However transforming the data to binary outcomes can result in loss of valuable information. To alleviate such challenges, this work develops methodology to associate genetic variant sets with multiple interval-censored outcomes. Testing sets of variants such as genes or pathways is a common approach in genetic association settings to lower the multiple testing burden, aggregate small effects, and improve interpretations of results. Instead of performing inference with only a single outcome, utilizing multiple outcomes can increase statistical power by aggregating information across multiple correlated phenotypes. Simulations show that the proposed strategy can offer significant power gains over a single outcome approach. We apply the proposed test to the investigation that motivated this study, a search for the genes that perturb risks of bone fractures and falls in the UK Biobank.

Genome-wide association study for bone quality of ducks during the laying period

The cage-rearing model of the modern poultry industry makes the bones of birds, especially egg-laying birds, more vulnerable to fracture, which poses serious damage to the health of birds. Research confirms that genetic material plays an important role in regulating bone growth, development, and remodeling. However, the genetic architecture underlying bone traits is not well understood. The objectives of this study are to identify valuable genes and genetic markers through a genome-wide association study (GWAS) for breeding to improve the duck bone quality. First, we quantified the tibia and femur quality traits of 260 laying ducks. Based on GWAS, a total of 75 SNP loci significantly associated with bone quality traits were identified, and 67 potential candidate genes were annotated. According to gene function analysis, genes P4HA2, WNT3A, and BST1 et al may influence bone quality by regulating bone cell activity, calcium and phosphate metabolism, or bone collagen maturation and cross-linking. Meanwhile, combined with the transcriptome results, we found that HOXB cluster genes are also important in bone growth and development. Therefore, our findings were helpful in further understanding the genetic architecture of the duck bone quality and provided a worthy theoretical basis and technological support to improve duck bone quality by breeding.

Integrating genome-wide association study with regulatory SNP annotations identified novel candidate genes for osteoporosis

Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p_{discovery GWAS} = 1.21 × 10^-25, p_{replication GWAS} = 1.80 × 10^-12), CCDC170 (p_{discovery GWAS} = 1.23 × 10^-11, p_{replication GWAS} = 3.22 × 10^-9), and SOX6 (p_{discovery GWAS} = 4.41 × 10^-15, p_{replication GWAS} = 6.57 × 10^-14). Then, for the 36 common target genes, multiple gene ontology (GO) terms were detected for BMD such as positive regulation of cartilage development (p = 9.27 × 10^-3) and positive regulation of chondrocyte differentiation (p = 9.27 × 10^-3). We explored the potential roles of rSNP in the genetic mechanisms of BMD and identified multiple candidate genes. Our study results support the implication of regulatory genetic variants in the development of OP.

microRNAs in newborns with low birth weight: relation to birth size and body composition

BACKGROUND

Children with low birth weight (LBW) have a higher risk of developing endocrine-metabolic disorders later in life. Deregulation of specific microRNAs (miRNAs) could underscore the programming of adult pathologies. We analyzed the miRNA expression pattern in both umbilical cord serum samples from LBW and appropriate-for-gestational-age (AGA) newborns and maternal serum samples in the 3rd trimester of gestation, and delineated the relationships with fetal growth, body composition, and markers of metabolic risk.

METHODS

Serum samples of 12 selected mother-newborn pairs, including 6 LBW and 6 AGA newborns, were used for assessing miRNA profile by RNA-sequencing. The miRNAs with differential expression were validated in a larger cohort [49 maternal samples and 49 umbilical cord samples (24 LBW, 25 AGA)] by RT-qPCR. Anthropometric, endocrine-metabolic markers and body composition (by DXA) in infants were determined longitudinally over 12 months.

RESULTS

LBW newborns presented reduced circulating concentrations of miR-191-3p (P = 0.015). miR-191-3p levels reliably differentiated LBW from AGA individuals (ROC AUC = 0.76) and were positively associated with anthropometric and body composition measures at birth and weight Z-score at 12 months (P < 0.05).

CONCLUSIONS

miR-191-3p was reliably different in LBW individuals, and could be a new player in the epigenetic mechanisms linking LBW and future endocrine-metabolic adverse outcomes.

IMPACT

Children with low birth weight (LBW) have a higher risk of developing endocrine-metabolic disorders. Deregulation of specific microRNAs (miRNAs) could underscore the programming of those pathologies. miR-191-3p is downregulated in serum of LBW newborns, and its concentrations associate positively with neonatal anthropometric measures, with lean mass and bone accretion at age 15 days and with weight Z-score at age 12 months. miR-191-3p was reliably different in individuals with LBW, and could be a new player in the epigenetic mechanisms connecting LBW and future endocrine-metabolic adverse outcomes.

Zuogui Wan slowed senescence of bone marrow mesenchymal stem cells by suppressing Wnt/β-catenin signaling

ETHNOPHARMACOLOGICAL RELEVANCE

In traditional Chinese medicine (TCM), Zuogui Wan (ZGW) is a classical prescription for senile disorders and delay aging. Modern studies show that ZGW promotes central nerve cell regeneration, prevents and cures osteoporosis, enhances the body's antioxidant capacity, regulates the body's immune function, and promotes mesenchymal stem cells (MSCs) proliferation.

AIM OF THE STUDY

It has been shown that MSCs aging is closely associated with organism's aging and age-related disorders. The study aimed to define the effects of ZGW on the aging bone marrow mesenchymal stem cells (BMSCs) and to identify the mechanisms of ZGW delaying BMSCs senescence.

MATERIALS AND METHODS

Network pharmacology analysis combined with GEO data mining, molecular docking and experimental validation were used to evaluate the mechanisms by which ZGW delays MSCs senescence (MSCS). LC-MS was used for quality control analysis of ZGW.

RESULTS

PPI network analysis revealed that EGF, TNF, JUN, MMPs, IL-6, MAPK8, and MYC are components of the core PPI network. GO and KEGG analyses revealed that oxidative stress, regulation of response to DNA damage stimuli, and Wnt signaling were significantly enriched. GEO database validation also indicated that Wnt signaling closely correlated with MSCs aging. Molecular docking analysis of the top-13 active components in the "ZGW-Targets-MSCS" network indicated that most components have strong affinity for key proteins in Wnt signaling, suggesting that modulation of Wnt signaling is an important mechanism of ZGW activity against MSCS. Further experimental validation found that ZGW indeed regulates Wnt signaling and suppresses the expression of age-related factors to enhance cell proliferation, ameliorate DNA damage, and reduce senescence-related secretory phenotype (SASP) secretion, thereby maintaining multidirectional differentiation of rat BMSCs. Similar results were obtained using the Wnt inhibitor, XAV-939.

CONCLUSIONS

Together, our data show that ZGW slows BMSCs aging by suppressing Wnt signaling.

Expression of Wnt signaling agonists and antagonists in periodontitis and healthy subjects, before and after non-surgical periodontal treatment: A systematic review

Periodontitis is a preventable and treatable multifactorial chronic inflammatory disease that can lead to irreversible periodontal destruction and tooth loss. Wnt signaling and its regulators play an important role in periodontal inflammation, destruction, regeneration, and reconstruction. This systematic review aimed at investigating the involvement of Wnt signaling agonists and antagonists in periodontitis and healthy subjects, before and after periodontal treatment. Electronic searches were carried out using MEDLINE/PubMed, EMBASE, and Cochrane Library databases in addition to hand searches. Studies having different designs assessing the levels of Wnt signaling antagonist and agonist levels in gingival crevicular fluid, serum, and tissue in patients diagnosed with periodontitis or gingivitis, compared with healthy individuals were included. In addition, studies compared these levels in periodontitis patients before and after non-surgical periodontal therapy were also eligible. Sixteen studies met the eligibility criteria. Sclerostin (SOST) has been mainly investigated in the literature (8 publications). Sclerostin (5 studies), Wnt-5a (2 studies), secreted frizzled-related protein 1 (SFRP1) (3 studies), and β-catenin (3 studies) show increased levels in periodontitis compared with periodontal health. Strong correlations between marker levels and periodontal clinical parameters were identified for SOST (5 studies), SFRP1 (2 studies), and β-catenin (2 studies). SOST (3 studies) and SFRP1 (1 study) levels significantly decrease following non-surgical periodontal treatment. The present systematic review demonstrated an association between Wnt signaling agonist and antagonist levels and periodontitis. Wnt agonists and antagonists may serve as valuable diagnostic and prognostic markers for periodontitis onset and progression. Further case-control and longitudinal studies should be conducted for different Wnt signaling agonists and antagonists.

Arterial calcification, atherosclerosis and osteoporosis: only clinical associations or a genetic platform?

The review is devoted to the comorbidity of two multifactorial diseases  — atherosclerosis and osteoporosis. Numerous epidemiological, experimental and clinical studies have confirmed the relationship between these diseases based on common risk factors and pathogenetic mechanisms. At the same time, to assess the associations between osteoporosis and atherosclerosis-related cardiovascular diseases, the following surrogate markers are used: vascular calcification, vascular stiffness, bone mineral density. It is known that atherosclerosis and osteoporosis depend on the human genotype, and they are caused by the interaction between the environment and genes. The modifiable risk factors for these diseases are largely similar, and the common features of atherosclerosis and osteoporosis pathogenesis make it possible to formulate the concept of a unified genetic basis of their development. Advances in molecular technology have made it possible to conduct a genome-wide association study (GWAS) and successfully identify genetic markers associated with both atherosclerosis and osteoporosis. The review aim was to describe the genes associated with developing atherosclerosis, arterial calcification and osteoporosis, as well as to provide information on the current understanding of the general genetic basis for plaque formation, vascular calcium deposition, and a decrease in bone mass. The analysis of publications from the PubMed, Medline, Web of Science and Cochrane Library databases since 2000 have been carried out. The article describes the genetic markers associated with atherosclerosis and osteoporosis, as well as considers the achievements in studying genetics of osteoporosis and atherosclerosis-related cardiovascular diseases. In addition, modern approaches and directions for further research of these diseases was established. The review can be useful for medical practitioners to clarify various genetic associations and mechanisms that lead to this comorbidity.

Exploring the Epigenetic Regulatory Role of m6A-Associated SNPs in Type 2 Diabetes Pathogenesis

Purpose

Genetic factors in type 2 diabetes (T2D) pathogenesis have been widely explored by the genome-wide association studies (GWAS), identifying a great amount of susceptibility loci. With the development of high-resolution sequencing, the N(6)-methyladenosine (m6A) RNA modification has been proved to be affected by genetic variation. In this study, we identified the T2D-associated m6A-SNPs from T2D GWAS data and explored the underlying mechanism of the pathogenesis of T2D.

Methods

We examined the association of m6A-SNPs with T2D among large-scale T2D GWAS summary statistics and further performed multi-omics integrated analysis to explore the potential role of the identified m6A-SNPs in T2D pathogenesis.

Results

Among the 15,124 T2D-associated m6A-SNPs, 71 of them reach the genome-wide significant threshold (5.0e-05). The leading SNP rs4993986 (C>G), which is located near the m6A modification site at the 3' end of the HLA-DQB1 transcript, is expected to participate in the pathogenesis of T2D by influencing m6A modification to regulate the HLA-DQB1 expression.

Conclusion

The current study has suggested a potential correlation between m6A-SNPs and T2D pathogenesis and also provided new insights into the pathogenic mechanism of the T2D susceptibility loci identified by GWAS.

Gene-based association analysis identified a novel gene associated with systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with strong genetic predisposition. Genome-wide association studies (GWAS) of SLE have identified more than 50 robust susceptibility loci. However, traditional individual SNP-based GWAS have made it difficult to identify variants with small effects. Moreover, variants revealed by GWAS only explain a limited disease heritability, suggesting that many susceptibility genes remain uncovered.

METHODS

We first curated the published SLE GWAS data from 1047 SLE patients and 1205 healthy controls of Chinese ancestry and performed a gene-based association study. Then quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was conducted to verify novel genes identified above.

RESULTS

Gene-based association study identified 10 SLE-associated genes, nine of which were reported by previous GWAS, the other one, ILRUN, is a newly identified gene and was further validated by qRT-PCR. Gene expression analysis of Gene Expression Omnibus (GEO) datasets also showed that the expression of ILRUN in patients with SLE was lower than that in normal subjects.

CONCLUSION

In this study, gene-based association study and qRT-PCR identified that ILRUN is a novel susceptibility gene of SLE. ILRUN may regulate inflammation and antiviral response through its effect on the transcription of type I interferons )I-IFN, and participate in the pathogenesis of SLE.

Identifying Susceptibility Loci for Cutaneous Squamous Cell Carcinoma Using a Fast Sequence Kernel Association Test

Cutaneous squamous cell carcinoma (cSCC) accounts for about 20% of all skin cancers, the most common type of malignancy in the United States. Genome-wide association studies (GWAS) have successfully identified multiple genetic variants associated with the risk of cSCC. Most of these studies were single-locus-based, testing genetic variants one-at-a-time. In this article, we performed gene-based association tests to evaluate the joint effect of multiple variants, especially rare variants, on the risk of cSCC by using a fast sequence kernel association test (fastSKAT). The study included 1,710 cSCC cases and 24,304 cancer-free controls from the Nurses' Health Study, the Nurses' Health Study II and the Health Professionals Follow-up Study. We used UCSC Genome Browser to define gene units as candidate loci, and further evaluated the association between all variants within each gene unit and disease outcome. Four genes HP1BP3, DAG1, SEPT7P2, and SLFN12 were identified using Bonferroni adjusted significance level. Our study is complementary to the existing GWASs, and our findings may provide additional insights into the etiology of cSCC. Further studies are needed to validate these findings.

Three functional polymorphisms in CCDC170 were associated with osteoporosis phenotype

MicroRNAs (miRNAs) play essential roles in regulating bone formation and homeostasis. Genomic variations within miRNA target sites may therefore be important sources of genetic differences in osteoporosis risk. The function of CCDC170 in bone biology is still unclear. To verify the function of CCDC170, we knocked down CCDC170 in cells and mice and searched for miRNA recognition sites within CCDC170 using the TargetScan, miRNASNP, and miRBase databases. In this study, our results demonstrated that CCDC170 plays an important role in the positive regulation of bone formation. MiR-153-3p, miR-374b-3p, miR-4274, miR-572 and miR-2964a-5p inhibited CCDC170 expression in an allele-specific manner by binding GWAS lead SNPs rs6932603, rs3757322 and rs3734806. These findings may improve our understanding of the association between CCDC170, miRNAs, GWAS lead SNPs, and osteoporosis pathogenesis and may provide a potential therapeutic target for osteoporosis therapy.

Summary: CCDC170 may provide a potential therapeutic target for osteoporosis therapy.

A novel computational strategy for DNA methylation imputation using mixture regression model (MRM)

BACKGROUND

DNA methylation is an important heritable epigenetic mark that plays a crucial role in transcriptional regulation and the pathogenesis of various human disorders. The commonly used DNA methylation measurement approaches, e.g., Illumina Infinium HumanMethylation-27 and -450 BeadChip arrays (27 K and 450 K arrays) and reduced representation bisulfite sequencing (RRBS), only cover a small proportion of the total CpG sites in the human genome, which considerably limited the scope of the DNA methylation analysis in those studies.

RESULTS

We proposed a new computational strategy to impute the methylation value at the unmeasured CpG sites using the mixture of regression model (MRM) of radial basis functions, integrating information of neighboring CpGs and the similarities in local methylation patterns across subjects and across multiple genomic regions. Our method achieved a better imputation accuracy over a set of competing methods on both simulated and empirical data, particularly when the missing rate is high. By applying MRM to an RRBS dataset from subjects with low versus high bone mineral density (BMD), we recovered methylation values of ~ 300 K CpGs in the promoter regions of chromosome 17 and identified some novel differentially methylated CpGs that are significantly associated with BMD.

CONCLUSIONS

Our method is well applicable to the numerous methylation studies. By expanding the coverage of the methylation dataset to unmeasured sites, it can significantly enhance the discovery of novel differential methylation signals and thus reveal the mechanisms underlying various human disorders/traits.

Association of HIVEP3 Gene and Lnc RNA with Femoral Neck Bone Mineral Content and Hip Geometry by Genome-Wide Association Analysis in Chinese People

PURPOSE

GWAS has successfully located and analyzed the pathogenic genes of osteoporosis. Genetic studies have found that heritability of BMD is 50%-85% while the other half is caused by hip geometric parameters and tissue horizontal characteristics. This study was designed to study the GWAS of osteoporosis in Shanghai Han population.

METHODS

We collected 1224 unrelated healthy young men (20-40 years old), young women (20-40 years old), and postmenopausal women (over 50 years old) who lived in Shanghai. BMD and hip geometric parameters were measured by dual-energy X-ray absorptiometry. The genomic DNA of peripheral blood was extracted and analyzed by using Illumina Human Asian Screening Array-24 + v1.0 (ASA) gene chip. Statistical analysis was carried out to evaluate the relationship between these SNPs and BMD and hip geometric parameters.

RESULTS

A total of 1155 subjects were included. We found that one SNP rs35282355 located in the human immunodeficiency virus type 1 enhancer-binding protein 3 gene (HIVEP3) and another 25 SNPs located in LINC RNA were significantly correlated with bone mineral content (BMC) in the femoral neck (P= 2.30 × 10^-9, P < 5 × 10^-8). We also found that the correlation between SNP rs35282355 and cross-sectional area (CSA) of hip geometry was a significant marginal statistical difference (P = 5.95 × 10^-8).

CONCLUSIONS

Through this study, we found that HIVEP3 gene and LINC RNA are potentially correlated with femoral neck BMC. These results provide important information for us to further understand the etiology and genetic pathogenesis of osteoporosis. In the future, we will expand the sample size to verify these loci and carry out molecular research.

WNT gene polymorphisms and predisposition to apical periodontitis

Single nucleotide polymorphisms (SNPs) in WNT genes may impact gene/protein function and contribute to individual predisposition to apical periodontitis (AP). Here, we investigated the association of SNPs in/nearby WNT3, WNT3A, WNT5A, WNT8A, WNT9B and WNT11 genes with AP using a case-control dataset. Cases were defined as individuals with deep caries and AP (n = 188); controls had deep caries and no AP (n = 230). Genotyping was performed using Taqman chemistry in real time PCR. Data analyses was performed using Fisher Exact tests assuming a Bonferroni correction threshold value of 0.005. Single-SNP association analysis revealed a trend for association with WNT3 rs9890413 genotypes (P = 0.009) under a dominant model and allelic association for WNT3A rs1745420 (P = 0.009). Haplotypes involving WNT3-WNT9B-WNT3A alleles were also significantly associated with AP (P ≤ 0.003). Luciferase reporter assays showed higher transcriptional activity (1.4-fold) with the alternate G allele in rs1745420. Expression of WNT3, WNT3A and WNT5A in AP tissues was significantly higher than in control tissues, and inversely correlated with the expression of SERPINB1, COL1A1 and TIMP1 (P < 0.05). Our results suggest that WNT genes have a role in modulating AP and polymorphisms in these genes may increase susceptibility to AP.

Genome-wide association study of lncRNA polymorphisms with bone mineral density

Recent studies suggested that long noncoding RNAs (lncRNAs) were widely transcribed in the genome, but their potential roles in the genetic complexity of human disorders required further exploration. The purpose of the present study was to explore genetic polymorphisms of lncRNAs associated with bone mineral density (BMD) and its potential value. Based on the lncRNASNP database, 55,906 lncSNPs were selected to conduct a genome-wide association study meta-analysis among 11,140 individuals of seven independent studies for BMDs at femoral neck (FN), lumbar spine, and total hip (HIP). Promising results were replicated in Genetic Factors for Osteoporosis Consortium (GEFOS Sequencing, n = 32,965). We found two lncRNA loci that were significantly associated with BMD. MEF2C antisense RNA 1 (MEF2C-AS1) located at 5q14.3 was significantly associated with FN-BMD after Bonferroni correction, and the strongest association signal was detected at rs6894139 (P = 3.03 × 10-9 ). LOC100506136 rs6465531 located at 7q21.3 showed significant association with HIP-BMD (P = 7.43 × 10-7 ). MEF2C-AS1 rs6894139 was replicated in GEFOS Sequencing with P-value of 1.43 × 10-23 . Our results illustrated the important role of polymorphisms in lncRNAs in determining variations of BMD and provided justification and evidence for subsequent functional studies.

Genome-wide identification of m6A-associated SNPs as potential functional variants for bone mineral density

This study investigated the effect of the N⁶-methyladenosine (m⁶A)-associated SNPs on bone mineral density (BMD) and found plenty of m⁶A-SNPs that were associated with BMD. This study increases our understanding on the regulation patterns of SNP and may provide new clues for further detection of functional mechanism underlying the associations between SNPs and osteoporosis.

INTRODUCTION

m⁶A plays critical roles in many fundamental biological processes and a variety of diseases. The m⁶A-associated SNPs may be potential functional variants for BMD. The aim of this study was to investigate the effect of the genome-wide m⁶A-SNPs on BMD.

METHODS

We examined the association of m⁶A-SNPs with femoral neck (FN) and lumbar spine (LS) BMD in 32,961 individuals and quantitative heel ultrasounds (eBMD) in 142,487 individuals. Furthermore, we performed expression quantitative trait locus (eQTL) analyses for the m⁶A-SNPs using whole genome data of about 10.5 million SNPs and 21,323 mRNAs from 43 Chinese individuals, as well as public available data. Differential expression analyses were also performed to support the identified genes.

RESULTS

We found 138, 125, and 993 m⁶A-SNPs which were associated with FN-BMD, LS-BMD, and eBMD (P < 0.05), respectively. The associations of rs11614913 (P = 8.92 × 10^-10) in MIR196A2 and rs1110720 (P = 2.05 × 10^-10) in ESPL1 with LS-BMD reached the genome-wide significance level. In addition, a total of 24 m⁶A-SNPs were significantly associated with eBMD (P < 5.0 × 10^-8). Further eQTL analyses showed that 47 of these BMD-associated m⁶A-SNPs were associated with expressions of the 46 corresponding local genes. Moreover, the expressions of 26 of these genes were associated with BMD.

CONCLUSION

The present study represents the first effort of investigating the associations and the mechanisms underlying the link between m⁶A-SNPs and BMD. The results suggested that m⁶A-SNP may play important roles in the pathology of osteoporosis.

Homozygous Dkk1 Knockout Mice Exhibit High Bone Mass Phenotype Due to Increased Bone Formation

Wnt antagonist Dkk1 is a negative regulator of bone formation and Dkk1 ^+/- heterozygous mice display a high bone mass phenotype. Complete loss of Dkk1 function disrupts embryonic head development. Homozygous Dkk1 ^-/- mice that were heterozygous for Wnt3 loss of function mutation (termed Dkk1 KO) are viable and allowed studying the effects of homozygous inactivation of Dkk1 on bone formation. Dkk1 KO mice showed a high bone mass phenotype exceeding that of heterozygous mice as well as a high incidence of polydactyly and kinky tails. Whole body bone density was increased in the Dkk1 KO mice as shown by longitudinal dual-energy X-ray absorptiometry. MicroCT analysis of the distal femur revealed up to 3-fold increases in trabecular bone volume and up to 2-fold increases in the vertebrae, compared to wild type controls. Cortical bone was increased in both the tibiae and vertebrae, which correlated with increased strength in tibial 4-point bending and vertebral compression tests. Dynamic histomorphometry identified increased bone formation as the mechanism underlying the high bone mass phenotype in Dkk1 KO mice, with no changes in bone resorption. Mice featuring only Wnt3 heterozygosity showed no evident bone phenotype. Our findings highlight a critical role for Dkk1 in the regulation of bone formation and a gene dose-dependent response to loss of DKK1 function. Targeting Dkk1 to enhance bone formation offers therapeutic potential for osteoporosis.

Functional relevance for associations between osteoporosis and genetic variants

Osteoporosis is characterized by increased bone loss and deterioration of bone microarchitecture, which will lead to reduced bone strength and increased risk of fragility fractures. Previous studies have identified many genetic loci associated with osteoporosis, but functional mechanisms underlying the associations have rarely been explored. In order to explore the potential molecular functional mechanisms underlying the associations for osteoporosis, we performed integrative analyses by using the publically available datasets and resources. We searched 128 identified osteoporosis associated SNPs (P<10-6), and 8 SNPs exert cis-regulation effects on 11 eQTL target genes. Among the 8 SNPs, 2 SNPs (RPL31 rs2278729 and LRP5 rs3736228) were confirmed to impact the expression of 3 genes (RPL31, CPT1A and MTL5) that were differentially expressed between human subjects of high BMD group and low BMD group. All of the functional evidence suggested the important functional mechanisms underlying the associations of the 2 SNPs (rs2278729 and rs3736228) and 3 genes (RPL31, CPT1A and MTL5) with osteoporosis. This study may provide novel insights into the functional mechanisms underlying the osteoporosis associated genetic variants, which will help us to comprehend the potential mechanisms underlying the genetic association for osteoporosis.

Gene-Based Genome-Wide Association Analysis in European and Asian Populations Identified Novel Genes for Rheumatoid Arthritis

Objective

Rheumatoid arthritis (RA) is a complex autoimmune disease. Using a gene-based association research strategy, the present study aims to detect unknown susceptibility to RA and to address the ethnic differences in genetic susceptibility to RA between European and Asian populations.

Methods

Gene-based association analyses were performed with KGG 2.5 by using publicly available large RA datasets (14,361 RA cases and 43,923 controls of European subjects, 4,873 RA cases and 17,642 controls of Asian Subjects). For the newly identified RA-associated genes, gene set enrichment analyses and protein-protein interactions analyses were carried out with DAVID and STRING version 10.0, respectively. Differential expression verification was conducted using 4 GEO datasets. The expression levels of three selected ‘highly verified’ genes were measured by ELISA among our in-house RA cases and controls.

Results

A total of 221 RA-associated genes were newly identified by gene-based association study, including 71‘overlapped’, 76 ‘European-specific’ and 74 ‘Asian-specific’ genes. Among them, 105 genes had significant differential expressions between RA patients and health controls at least in one dataset, especially for 20 genes including 11 ‘overlapped’ (ABCF1, FLOT1, HLA-F, IER3, TUBB, ZKSCAN4, BTN3A3, HSP90AB1, CUTA, BRD2, HLA-DMA), 5 ‘European-specific’ (PHTF1, RPS18, BAK1, TNFRSF14, SUOX) and 4 ‘Asian-specific’ (RNASET2, HFE, BTN2A2, MAPK13) genes whose differential expressions were significant at least in three datasets. The protein expressions of two selected genes FLOT1 (P value = 1.70E-02) and HLA-DMA (P value = 4.70E-02) in plasma were significantly different in our in-house samples.

Conclusion

Our study identified 221 novel RA-associated genes and especially highlighted the importance of 20 candidate genes on RA. The results addressed ethnic genetic background differences for RA susceptibility between European and Asian populations and detected a long list of overlapped or ethnic specific RA genes. The study not only greatly increases our understanding of genetic susceptibility to RA, but also provides important insights into the ethno-genetic homogeneity and heterogeneity of RA in both ethnicities.

Genetic control of bone mass

Bone mineral density (BMD) is a quantitative traits used as a surrogate phenotype for the diagnosis of osteoporosis, a common metabolic disorder characterized by increased fracture risk as a result of a decreased bone mass and deterioration of the microarchitecture of the bone. Normal variation in BMD is determined by both environmental and genetic factors. According to heritability studies, 50-85% of the variance in BMD is controlled by genetic factors which are mostly polygenic. In contrast to the complex etiology of osteoporosis, there are disorders with deviating BMD values caused by one mutation with a large impact. These mutations can result in monogenic bone disorders with either an extreme high (sclerosteosis, Van Buchem disease, osteopetrosis, high bone mass phenotype) or low BMD (osteogenesis imperfecta, juvenile osteoporosis, primary osteoporosis). Identification of the disease causing genes, increased the knowledge on the regulation of BMD and highlighted important signaling pathways and novel therapeutic targets such as sclerostin, RANKL and cathepsin K. Genetic variation in genes involved in these pathways are often also involved in the regulation of normal variation in BMD and osteoporosis susceptibility. In the last decades, identification of genetic factors regulating BMD has proven to be a challenge. Several approaches have been tested such as linkage studies and candidate and genome wide association studies. Although, throughout the years, technological developments made it possible to study increasing numbers of genetic variants in populations with increasing sample sizes at the same time, only a small fraction of the genetic impact can yet be explained. In order to elucidate the missing heritability, the focus shifted to studying the role of rare variants, copy number variations and epigenetic influences. This review summarizes the genetic cause of different monogenic bone disorders with deviating BMD and the knowledge on genetic factors explaining normal variation in BMD and osteoporosis risk.

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

Osteoporosis is a medical condition of global concern, with increasing incidence in both sexes. Bone mineral density (BMD), a highly heritable trait, has been proven a useful diagnostic factor in predicting fracture. Because medical information is lacking about male osteoporotic genetics, we conducted a genome-wide association study of BMD in Korean men. With 1,176 participants, we analyzed 4,414,664 single nucleotide polymorphisms (SNPs) after genomic imputation, and identified five SNPs and three loci correlated with bone density and strength. Multivariate linear regression models were applied to adjust for age and body mass index interference. Rs17124500 (p = 6.42 × 10^-7), rs34594869 (p = 6.53 × 10^-7) and rs17124504 (p = 6.53 × 10^-7) in 14q31.3 and rs140155614 (p = 8.64 × 10^-7) in 15q25.1 were significantly associated with lumbar spine BMD (LS-BMD), while rs111822233 (p = 6.35 × 10^-7) was linked with the femur total BMD (FT-BMD). Additionally, we analyzed the relationship between BMD and five genes previously identified in Korean men. Rs61382873 (p = 0.0009) in LRP5, rs9567003 (p = 0.0033) in TNFSF11 and rs9935828 (p = 0.0248) in FOXL1 were observed for LS-BMD. Furthermore, rs33997547 (p = 0.0057) in ZBTB and rs1664496 (p = 0.0012) in MEF2C were found to influence FT-BMD and rs61769193 (p = 0.0114) in ZBTB to influence femur neck BMD. We identified five SNPs and three genomic regions, associated with BMD. The significance of our results lies in the discovery of new loci, while also affirming a previously significant locus, as potential osteoporotic factors in the Korean male population.

Diametrical diseases reflect evolutionary-genetic tradeoffs: Evidence from psychiatry, neurology, rheumatology, oncology and immunology

Tradeoffs centrally mediate the expression of human adaptations. We propose that tradeoffs also influence the prevalence and forms of human maladaptation manifest in disease. By this logic, increased risk for one set of diseases commonly engenders decreased risk for another, diametric, set of diseases. We describe evidence for such diametric sets of diseases from epidemiological, genetic and molecular studies in four clinical domains: (i) psychiatry (autism vs psychotic-affective conditions), (ii) rheumatology (osteoarthritis vs osteoporosis), (iii) oncology and neurology (cancer vs neurodegenerative disorders) and (iv) immunology (autoimmunity vs infectious disease). Diametric disorders are important to recognize because genotypes or environmental factors that increase risk for one set of disorders protect from opposite disorders, thereby providing novel and direct insights into disease causes, prevention and therapy. Ascertaining the mechanisms that underlie disease-related tradeoffs should also indicate means of circumventing or alleviating them, and thus reducing the incidence and impacts of human disease in a more general way.