Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

A potential function for MicroRNA-124 in normal and pathological bone conditions

Cells produce short single-stranded non-coding RNAs (ncRNAs) called microRNAs (miRNAs), which actively regulate gene expression at the posttranscriptional level. Several miRNAs have been observed to exert significant impacts on bone health and bone-related disorders. One of these, miR-124, is observed in bone microenvironments and is conserved across species. It affects bone cell growth and differentiation by activating different transcription factors and signaling pathways. In-depth functional analyses of miR-124 have revealed several physiological and pathological roles exerted through interactions with other ncRNAs. Deciphering these RNA-mediated signaling networks and pathways is essential for understanding the potential impacts of dysregulated miRNA functions on bone biology. In this review, we aim to provide a comprehensive analysis of miR-124's involvement in bone physiology and pathology. We highlight the importance of miR-124 in controlling transcription factors and signaling pathways that promote bone growth. This review reveals therapeutic implications for the treatment of bone-related diseases.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

Did giraffe cardiovascular evolution solve the problem of heart failure with preserved ejection fraction?

The evolved adaptations of other species can be a source of insight for novel biomedical innovation. Limitations of traditional animal models for the study of some pathologies are fueling efforts to find new approaches to biomedical investigation. One emerging approach recognizes the evolved adaptations in other species as possible solutions to human pathology. The giraffe heart, for example, appears resistant to pathology related to heart failure with preserved ejection fraction (HFpEF)—a leading form of hypertension-associated cardiovascular disease in humans. Here, we postulate that the physiological pressure-induced left ventricular thickening in giraffes does not result in the pathological cardiovascular changes observed in humans with hypertension. The mechanisms underlying this cardiovascular adaptation to high blood pressure in the giraffe may be a bioinspired roadmap for preventive and therapeutic strategies for human HFpEF.

A towering genome: Experimentally validated adaptations to high blood pressure and extreme stature in the giraffe

The suite of adaptations associated with the extreme stature of the giraffe has long interested biologists and physiologists. By generating a high-quality chromosome-level giraffe genome and a comprehensive comparison with other ruminant genomes, we identified a robust catalog of giraffe-specific mutations. These are primarily related to cardiovascular, bone growth, vision, hearing, and circadian functions. Among them, the giraffe FGFRL1 gene is an outlier with seven unique amino acid substitutions not found in any other ruminant. Gene-edited mice with the giraffe-type FGFRL1 show exceptional hypertension resistance and higher bone mineral density, both of which are tightly connected with giraffe adaptations to high stature. Our results facilitate a deeper understanding of the molecular mechanism underpinning distinct giraffe traits, and may provide insights into the study of hypertension in humans.

Osteoporosis genome-wide association study variant c.3781 C>A is regulated by a novel anti-osteogenic factor miR-345-5p

Upstream transcription factor family member 3 (USF3) c.3781C>A (rs1026364) in the 3'-untranslated region (3'-UTR) has been firmly associated with bone mineral density (BMD) in genome-wide association study (GWAS). However, the molecular mechanism by which it influences BMD and osteoporosis is unknown. Bioinformatics analyses suggested that the risk c.3781A allele creates a target site for hsa-miR-345-5p binding. Luciferase assay validated that the c.3781A allele displayed significantly lower luciferase activities than the c.3781C allele in the human osteoblast cell line hFOB1.19, osteosarcoma cell lines U-2OS and Saos-2, and embryonic kidney cell line 293T. Furthermore, hsa-miR-345-5p regulated USF3 expression on both messenger RNA and protein levels in hFOB1.19 and U937 cells with heterozygous A/C genotype. Transfection of hsa-miR-345-5p antagomiR in heterozygous hFOB1.19 cells significantly increased the expression of osteogenic marker genes RUNX2, OSTERIX, COL1A1, ALP, OPN, OCN, and alkaline phosphatase activity and matrix mineralization level. Importantly, we found that hsa-miR-345-5p also inhibits osteoblast maturation in cell lines U-2OS with hsa-miR-345-5p nonbinding C/C genotype by targeting RUNX3 and SMAD1. Our findings uncovered a novel pathogenetic mechanism of osteoporosis by GWAS variant c.3781C>A-mediated disruption of hsa-miR-345-5p binding at the 3'-UTR of USF3 and the functional role of hsa-miR-345-5p in osteogenic differentiation.

Evaluation of common variants in the CNR2 gene and its interaction with abdominal obesity for osteoporosis susceptibility in Chinese post-menopausal females

Objectives

The objective of this study was to investigate the association of four single-nucleotide polymorphisms (SNPs) of the cannabinoid receptor 2 (CNR2) gene, gene-obesity interaction, and haplotype combination with osteoporosis (OP) susceptibility.

Methods

Chinese patients with OP were recruited between March 2011 and December 2015 from our hospital. In this study, a total of 1267 post-menopausal female patients (631 OP patients and 636 control patients) were selected. The mean age of all subjects was 69.2 years (sd 15.8). A generalized multifactor dimensionality reduction (GMDR) model and logistic regression model were used to examine the interaction between SNP and obesity on OP. For OP patient-control haplotype analyses, the SHEsis online haplotype analysis software (http://analysis.bio-x.cn/) was employed.

Results

The logistic regression model revealed that the C allele of rs2501431 and the G allele of rs3003336 were associated with increased OP risk, compared with those with wild genotype. However, no significant correlations were found when analyzing the association of rs4237 and rs2229579 with OP risk. The GMDR analysis suggested that the interaction model composed of two factors, rs3003336 and abdominal obesity (AO), was the best model with statistical significance (p-value from sign test (P_sign) = 0.012), indicating a potential gene-environment interaction between rs3003336 and AO. Overall, the two-locus models had a cross-validation consistency of 10/10 and had a testing accuracy of 0.641. Abdominally obese subjects with the AG or GG genotype have the highest OP risk, compared with subjects with the AA genotype and normal waist circumference (WC) (odds ratio (OR) 2.23, 95% confidence interval (CI) 1.54 to 3.51). Haplotype analysis also indicated that the haplotype containing the rs3003336-G and rs2501431-C alleles was associated with a statistically increased OP risk.

Conclusion

Our results suggested that the C allele of rs2501431 and the G allele of rs3003336 of the CNR2 gene, interaction between rs3003336 and AO, and the haplotype containing the rs3003336-G and rs2501431-C alleles were all associated with increased OP risk.

Cite this article: Bone Joint Res 2019;8:544–549.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

MicroRNA and Human Bone Health

The small non‐coding microRNAs (miRNAs) are post‐transcription regulators that modulate diverse cellular process in bone cells. Because optimal miRNA targeting is essential for their function, single‐nucleotide polymorphisms (SNPs) within or proximal to the loci of miRNA (miR‐SNPs) or mRNA (PolymiRTS) could potentially disrupt the miRNA‐mRNA interaction, leading to changes in bone metabolism and osteoporosis. Recent human studies of skeletal traits using miRNA profiling, genomewide association studies, and functional studies started to decipher the complex miRNA regulatory network. These studies have indicated that miRNAs may be a promising bone marker. This review focuses on human miRNA studies on bone traits and discusses how genetic variants affect bone metabolic pathways. Major ex vivo investigations using human samples supported with animal and in vitro models have shed light on the mechanistic role of miRNAs. Furthermore, studying the miRNAs’ signatures in secondary osteoporosis and osteoporotic medications such as teriparatide (TPTD) and denosumab (DMab) have provided valuable insight into clinical management of the disease. © 2018 The Authors. JBMR Plus Published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research

Joint study of two genome-wide association meta-analyses identified 20p12.1 and 20q13.33 for bone mineral density

In the present study, aiming to identify loci associated with osteoporosis, we conducted a joint association study of 2 independent genome-wide association meta-analyses of femoral neck and lumbar spine bone mineral densities (BMDs): 1) an in-house study of 6 samples involving 7484 subjects, and 2) the GEFOS-seq study of 7 samples involving 32,965 subjects. The in-house samples were imputed by the 1000 genomes project phase 3 reference panel. SNP-based association test was applied to 7,998,108 autosomal SNPs in each meta-analysis, and for each SNP the 2 association signals were then combined for joint analysis and for mutual replication. Combining the evidence from both studies, we identified 2 novel loci associated with BMDs at the genome-wide significance level (α=5.0×10-8): 20p12.1 (rs73100693 p=2.65×10-8, closest gene MACROD2) and 20q13.33 (rs2380128 p=3.44×10-8, OSBPL2). We also replicated 7 loci that were reported by two recent studies on heel and total body BMD. Our findings provide useful insights that enhance our understanding of bone development, osteoporosis and fracture pathogenesis.

Role of fibroblast growth factor receptors (FGFR) and FGFR like-1 (FGFRL1) in mesenchymal stromal cell differentiation to osteoblasts and adipocytes

Fibroblast growth factors (FGF) and their receptors (FGFRs) regulate many developmental processes including differentiation of mesenchymal stromal cells (MSC). We developed two MSC lines capable of differentiating to osteoblasts and adipocytes and studied the role of FGFRs in this process. We identified FGFR2 and fibroblast growth factor receptor like-1 (FGFRL1) as possible actors in MSC differentiation with gene microarray and qRT-PCR. FGFR2 and FGFRL1 mRNA expression strongly increased during MSC differentiation to osteoblasts. FGF2 treatment, resulting in downregulation of FGFR2, or silencing FGFR2 expression with siRNAs inhibited osteoblast differentiation. During adipocyte differentiation expression of FGFR1 and FGFRL1 increased and was down-regulated by FGF2. FGFR1 knockdown inhibited adipocyte differentiation. Silencing FGFR2 and FGFR1 in MSCs was associated with decreased FGFRL1 expression in osteoblasts and adipocytes, respectively. Our results suggest that FGFR1 and FGFR2 regulate FGFRL1 expression. FGFRL1 may mediate or modulate FGFR regulation of MSC differentiation together with FGFR2 in osteoblastic and FGFR1 in adipocytic lineage.

microRNA Regulation of Skeletal Development

PURPOSE OF REVIEW

Osteogenesis is a complex process involving the specification of multiple progenitor cells and their maturation and differentiation into matrix-secreting osteoblasts. Osteogenesis occurs not only during embryogenesis but also during growth, after an injury, and in normal homeostatic maintenance. While much is known about osteogenesis-associated regulatory genes, the role of microRNAs (miRNAs), which are epigenetic regulators of protein expression, is just beginning to be explored. While miRNAs do not abrogate all protein expression, their purpose is to finely tune it, allowing for a timely and temporary protein down-regulation.

RECENT FINDINGS

The last decade has unveiled a multitude of miRNAs that regulate key proteins within the osteogenic lineage, thus qualifying them as "ostemiRs." These miRNAs may endogenously target an activator or inhibitor of differentiation, and depending on the target, may either lead to the prolongation of a progenitor maintenance state or to early differentiation. Interestingly, cellular identity seems intimately coupled to the expression of miRNAs, which participate in the suppression of previous and subsequent differentiation steps. In such cases where key osteogenic proteins were identified as direct targets of miRNAs in non-bone cell types, or through bioinformatic prediction, future research illuminating the activity of these miRNAs during osteogenesis will be extremely valuable. Many bone-related diseases involve the dysregulation of transcription factors or other proteins found within osteoblasts and their progenitors, and the dysregulation of miRNAs, which target such factors, may play a pivotal role in disease etiology, or even as a possible therapy.

Noncoding RNAs Regulating NF-κB Signaling

As transcription factors that regulate expression of a variety of genes essential for diverse physiological and pathological processes, nuclear factor kappa B (NF-κB) family molecules play important roles in the development and progression of malignant tumor, and constitutive activation of NF-κB has been evidenced in various types of tumor tissues. Underlying its pathologic role, deregulated expression and/or transactivating activity of NF-κB usually involves multiple layers of molecular mechanisms. Noncoding RNAs, including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are known to modulate expression and biological functions of regulatory proteins in a variety of cancer contexts. In this chapter, the regulatory role of miRNAs and lncRNAs in NF-κB signaling in malignant diseases will be discussed.

MicroRNA-140-5p impairs zebrafish embryonic bone development via targeting BMP-2

MicroRNA-140-5p (miRNA-140-5p) is important for embryonic bone development. In this study, we found that miRNA-140-5p and its binding site in the 3'UTR of bone morphogenetic protein 2 (BMP-2) are highly conserved among vertebrates, and miRNA-140-5p negatively regulates both zebrafish and human BMP-2 genes. Microinjection of miRNA-140-5p or BMP-2b morpholino into zebrafish embryos led to a similar phenotype, including shortened tails, curved trunks, and defects in cranial cartilage. Moreover, miRNA-140-5p injection induced zebrafish embryo malformation that could be significantly rescued by microinjection of BMP-2 mRNA. In conclusion, our results indicated that miRNA-140-5p regulates zebrafish embryonic bone development via targeting BMP-2.

Genome-wide association studies identified multiple genetic loci for body size at four growth stages in Chinese Holstein cattle

The growth and maturity of cattle body size affect not only feed efficiency, but also productivity and longevity. Dissecting the genetic architecture of body size is critical for cattle breeding to improve both efficiency and productivity. The volume and weight of body size are indicated by several measurements. Among them, Heart Girth (HG) and Hip Height (HH) are the most important traits. They are widely used as predictors of body weight (BW). Few association studies have been conducted for HG and HH in cattle focusing on single growth stage. In this study, we extended the Genome-wide association studies to a full spectrum of four growth stages (6-, 12-, 18-, and 24-months after birth) in Chinese Holstein heifers. The whole genomic single nucleotide polymorphisms (SNPs) were obtained from the Illumina BovineSNP50 v2 BeadChip genotyped on 3,325 individuals. Estimated breeding values (EBVs) were derived for both HG and HH at the four different ages and analyzed separately for GWAS by using the Fixed and random model Circuitous Probability Unification (FarmCPU) method. In total, 27 SNPs were identified to be significantly associated with HG and HH at different growth stages. We found 66 candidate genes located nearby the associated SNPs, including nine genes that were known as highly related to development and skeletal and muscular growth. In addition, biological function analysis was performed by Ingenuity Pathway Analysis and an interaction network related to development was obtained, which contained 16 genes out of the 66 candidates. The set of putative genes provided valuable resources and can help elucidate the genomic architecture and mechanisms underlying growth traits in dairy cattle.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.

Low density lipoprotein receptor-related protein 5 gene polymorphisms and osteoporosis in Thai menopausal women

Background

Osteoporosis, characterized by low bone mineral density (BMD) and high bone fracture risk, is prevalent in Thai menopausal women. Genetic factors are known to play a key role in BMD. Low density lipoprotein receptor-related protein 5 (LRP5), a co-receptor in the Wnt/beta-catenin pathway, is involved in many aspects of bone biology. As coding single nucleotide polymorphisms (cSNPs) of LRP5, including A1330V (rs3736228), and Asian-related Q89R (rs41494349) and N740N (rs2306862), are associated with lowered BMD, this study aimed to determine the relationship between these LRP5 polymorphisms and BMD in 277 Thai menopausal women.

Results

Only rs3736228 deviated from the Hardy–Weinberg equilibrium of allele frequency (p = 0.022). The median, range and p value for the BMD related to each SNP parameter were compared (Mann–Whitney U test). Significant differences were observed between wild-type and risk alleles for both rs3736228 (total radial, p = 0.011; and radial 33, p = 0.001) and rs2306862 (radial 33: p = 0.015) SNPs, with no significant difference for rs41494349 SNP. Linkage disequilibrium was strong for both rs3736228 and rs2306862 SNPs. Haplotype analysis identified high CC frequency in both normal and osteopenia/osteoporosis groups, with a significant odds ratio for carrying the TT haplotype; however, this was non-significant after adjusting for age. Multivariate binary logistic regression analysis performed for rs3736228 showed that individuals with a body mass index <25 kg/m² had an increased risk of osteoporosis for each decade, but the polymorphism had no effect.

Conclusions

This study did not identify LRP5 polymorphisms as a risk factor for osteoporosis in Thai menopausal women. Further studies with larger sample sizes are needed to further clarify the role of LRP5 as a genetic determinant of osteoporosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12952-016-0059-7) contains supplementary material, which is available to authorized users.

Identification of IDUA and WNT16 Phosphorylation-Related Non-Synonymous Polymorphisms for Bone Mineral Density in Meta-Analyses of Genome-Wide Association Studies

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10(-6) (0.05/9593) and 1.00 × 10(-4), respectively. In stage 2, nine stage 1-discovered phosSNPs (based on α = 1.00 × 10(-4)) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10(-3), 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10(-10), p = 5.26 × 10(-10), and p = 3.01 × 10(-10), respectively) and HIP-BMD (p = 3.26 × 10(-6), p = 1.97 × 10(-6), and p = 1.63 × 10(-12), respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants.