Genetic Polymorphism of miR-196a-2 is Associated with Bone Mineral Density (BMD)

Specific Blood Plasma Circulating miRs Are Associated with the Physiological Impact of Total Fish Meal Replacement with Soybean Meal in Diets for Rainbow Trout (Oncorhynchus mykiss)

High dietary SBM content is known to induce important physiological alterations, hampering its use as a major FM alternative. Rainbow trout (Oncorhynchus mykiss) juveniles were fed two experimental diets during 9 weeks: (i) a FM diet containing 12% FM; and (ii) a vegetable meal (VM) diet totally devoid of FM and based on SBM (26%). Fish fed the VM diet did not show reduced growth performance when compared with fish fed the FM diet. Nevertheless, fish fed the VM diet had an increased viscerosomatic index, lower apparent fat digestibility, higher aminopeptidase enzyme activity and number of villi fusions, and lower α-amylase enzyme activity and brush border integrity. Small RNA-Seq analysis identified six miRs (omy-miR-730a-5p, omy-miR-135c-5p, omy-miR-93a-3p, omy-miR-152-5p, omy-miR-133a-5p, and omy-miR-196a-3p) with higher expression in blood plasma from fish fed the VM diet. Bioinformatic prediction of target mRNAs identified several overrepresented biological processes known to be associated with high dietary SBM content (e.g., lipid metabolism, epithelial integrity disruption, and bile acid status). The present research work increases our understanding of how SBM dietary content has a physiological impact in farmed fish and suggests circulating miRs might be suitable, integrative, and less invasive biomarkers in fish.

Roles of miR-196a and miR-196b in Zebrafish Motor Function

Background: The exertion of motor function depends on various tissues, such as bones and muscles. miR-196 has been widely studied in cancer and other fields, but its effect on bone and skeletal muscle is rarely reported. In order to explore the role of miR-196 family in bone and skeletal muscle, we used the previously successfully constructed miR-196a-1 and miR-196b gene knockout zebrafish animal models for research. Methods: The behavioral trajectories of zebrafish from 4 days post-fertilization (dpf) to 7 dpf were detected to analyze the effect of miR-196a-1 and miR-196b on motor ability. Hematoxylin-eosin (HE) staining and transmission electron microscopy (TEM) were used to detect the dorsal muscle tissue of zebrafish. The bone tissue of zebrafish was detected by microcomputed tomography (micro-CT). Real-time PCR was used to detect the expression levels of related genes, including vcp, dpm1, acta1b, mylpfb, col1a1a, bmp8a, gdf6a, and fgfr3. Results: The behavioral test showed that the total behavioral trajectory, movement time, and movement speed of zebrafish larvae were decreased in the miR-196a-1 and miR-196b gene knockout lines. Muscle tissue analysis showed that the structure of muscle fibers in the zebrafish lacking miR-196a-1 and miR-196b was abnormal and was characterized by vacuolar degeneration of muscle fibers, intranuclear migration, melanin deposition, and inflammatory cell infiltration. Bone CT examination revealed decreased bone mineral density and trabecular bone number. The real-time PCR results showed that the expression levels of vcp, dpm1, gdf6a, fgfr3, and col1a1a were decreased in the miR-196b gene knockout group. The expression levels of dpm1, acta1b, mylpfb, gdf6a, and col1a1a were decreased, and the expression level of fgfr3 was increased in the miR-196b gene knockout group compared with the wild-type group. Conclusions: miR-196a-1 and miR-196b play an important role in muscle fiber structure, bone mineral density, and bone trabecular quantity by affecting the expression of vcp, dpm1, acta1b, mylpfb, gdf6a, fgfr3, and col1a1a and then affect the function of the motor system

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

Critical roles of microRNA-196 in normal physiology and non-malignant diseases: Diagnostic and therapeutic implications

MicroRNAs (miRNAs) have emerged as a critical component of regulatory networks that modulate and fine-tune gene expression in a post-transcriptional manner. The microRNA-196 family is encoded by three loci in the human genome, namely hsa-mir-196a-1, hsa-mir-196a-2, and hsa-mir-196b. Increasing evidence supports the roles of different components of this miRNA family in regulating key cellular processes during differentiation and development, ranging from inflammation and differentiation of stem cells to limb development and remodeling and structure of adipose tissue. This review first discusses about the genomic context and regulation of this miRNA family and then take a bird's eye view on the updated list of its target genes and their biological processes to obtain insights about various functions played by members of the microRNA-196 family. We then describe evidence supporting the involvement of the human microRNA-196 family in regulating critical cellular processes both in physiological and non-malignant inflammatory conditions, highlighting recent seminal findings that carry implications for developing novel therapeutic or diagnostic strategies.

Cancer-associated fibroblasts activated by miR-196a promote the migration and invasion of lung cancer cells

Fibroblasts in the tumor microenvironment, known as cancer-associated fibroblasts (CAFs), promote the migration, invasion, and metastasis of cancer cells when they are activated through diverse processes, including post-transcriptional regulation by microRNAs (miRNAs). To identify the miRNAs that regulate CAF activation, we used NanoString to profile miRNA expression within normal mouse lung fibroblasts (LFs) and CAFs. Based on NanoString profiling, miR-196a was selected as a candidate that was up-regulated in CAFs. miR-196a-overexpressed LFs (LF-196a) promoted the migration and invasion of lung cancer cells in co-culture systems (Transwell migration and spheroid invasion assays). ANXA1 was confirmed as a direct target of miR-196a, and adding back ANXA1 to LF-196a restored the cancer cell invasion promoted by miR-196a. miR-196a increased CCL2 secretion in fibroblasts, and that was suppressed by ANXA1. Furthermore, blocking CCL2 impeded cancer spheroid invasion. In lung adenocarcinoma patients, high miR-196a expression was associated with poor prognosis. Collectively, our results suggest that CAF-specific miR-196a promotes lung cancer progression in the tumor microenvironment via ANXA1 and CCL2 and that miR-196a will be a good therapeutic target or biomarker in lung adenocarcinoma.

Genetic variability in noncoding RNAs: involvement of miRNAs and long noncoding RNAs in osteoporosis pathogenesis

The pathogenesis of osteoporosis is multifactorial and is the consequence of genetic, hormonal and lifestyle factors. Epigenetics, including noncoding RNA (ncRNA) deregulation, represents a link between susceptibility to develop the disease and environmental influences. The majority of studies investigated the expression of ncRNAs in osteoporosis patients; however, very little information is available on their genetic variability. In this review, we focus on two classes of ncRNAs: miRNAs and long noncoding RNAs (lncRNAs). We summarize recent findings on how polymorphisms in miRNAs and lncRNAs can perturb the lncRNA/miRNA/mRNA axis and may be involved in osteoporosis clinical outcome. We also provide a general overview on databases and bioinformatic tools useful for associating miRNAs and lncRNAs variability with complex genetic diseases.

Association between miR-196a2 polymorphism and the development of hepatocellular carcinoma in the Egyptian population

Background

Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide. Circulating microRNAs (miRNAs) are endogenous, small (17–25 nucleotides) non-coding RNAs that are overexpressed in many human cancers including HCC. Single-nucleotide polymorphisms (SNPs) of miRNAs play an important role in the pathogenesis of HCC. In our study, we aimed to evaluate the role of miR-196a2 rs11614913 polymorphism in the development of HCC. A total of 200 subjects, including 80 HCC patients, 60 patients with liver cirrhosis, and 60 healthy controls were selected. The polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was taken to determine miR-196a2 rs11614913 polymorphism.

Results

The genotype distribution of the TC and CC, TC + CC genotypes, and the C allele were significantly higher in HCC patients than control and cirrhotic groups (P = 0.02, P = 0.005, and P = 0.003, respectively). Compared with the wild-type TT genotype, both the variant TC, CC, TC + CC genotypes were associated with an elevated risk of HCC (OR = 2.77, 95% CI = 1.27–6.04), (OR = 4.94, 95% CI = 1.74–14.07), (OR = 3.24, 95% CI = 1.55–6.78) respectively. Moreover, the C allele was correlated with an increased risk of HCC (OR = 2.30, 95% CI = 1.40–3.76) compared to the wide-type T allele. Also, there is no significant correlation between the different miR-196a2 genotypes and either the clinico-pathologic features of HCC or its aggressiveness.

Conclusion

Our results suggest that the miR-196a2 rs11614913 polymorphism is associated with an increased risk of HCC in the Egyptian population.

Evidences from a Systematic Review and Meta-Analysis Unveil the Role of MiRNA Polymorphisms in the Predisposition to Female Neoplasms

Breast (BCa) and gynecological (GCa) cancers constitute a group of female neoplasms that has a worldwide significant contribution to cancer morbidity and mortality. Evidence suggests that polymorphisms influencing miRNA function can provide useful information towards predicting the risk of female neoplasms. Inconsistent findings in the literature should be detected and resolved to facilitate the genetic screening of miRNA polymorphisms, even during childhood or adolescence, and their use as predictors of future malignancies. This study represents a comprehensive systematic review and meta-analysis of the association between miRNA polymorphisms and the risk of female neoplasms. Meta-analysis was performed by pooling odds-ratios (ORs) and generalized ORs while using a random-effects model for 15 miRNA polymorphisms. The results suggest that miR-146a rs2910164 is implicated in the susceptibility to GCa. Moreover, miR-196a2 rs11614913-T had a moderate protective effect against female neoplasms, especially GCa, in Asians but not in Caucasians. MiR-27a rs895819-G might pose a protective effect against BCa among Caucasians. MiR-499 rs3746444-C may slightly increase the risk of female neoplasms, especially BCa. MiR-124 rs531564-G may be associated with a lower risk of female neoplasms. The current evidences do not support the association of the remaining polymorphisms and the risk of female neoplasms.

Overexpressed miR-196a accelerates osteogenic differentiation in osteoporotic mice via GNAS-dependent Hedgehog signaling pathway

Osteoporosis (OP), a common metabolic bone disease, is accompanied by reduced bone mass, bone mineral density (BMD), as well as microstructure destruction of bone. Previously, microRNA-196a-2 (miR-196a-2) and miR-196a-3p were reported for its involvement in BMD. Herein, this study set out to identify the functional relevance of miR-196a in osteogenic differentiation in osteoporotic mice and explore the associated mechanism by establishing an OP mouse model. Guanine nucleotide binding protein, alpha stimulating (GNAS) was verified as a target gene of miR-196a, which was decreased in OP mice. Furthermore, the bone marrow stromal cells (BMSCs) were then extracted from OP mice and treated with miR-196 mimic/inhibitor or small interfering RNA against GNAS to investigate miR-196a interaction with GNAS and the Hedgehog signaling pathway. BMSCs in OP mice transfected with miR-196a mimic or si-GNAS displayed the elevated expression of Smo, ALP, Runx2, and OPN, as well as bone gla protein and tartrate-resistant acid phosphatase, elevated ALP vitality and bone formation ability as well as reduced expression of GNAS and PTCH. Taken conjointly, overexpression of miR-196a repressed GNAS expression by activating the Hedgehog signaling pathway, thus promoting osteogenic differentiation in mice with OP.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Recent Advances in the Genetics of Fractures in Osteoporosis

Genetic susceptibility, together with old age, female sex, and low bone mineral density (BMD) are amongst the strongest determinants of fracture risk. Tmost recent large-scale genome-wide association study (GWAS) meta-analysis has yielded fifteen loci. This review focuses on the advances in the research of genetic determinants of fracture risk. We first discuss the genetic architecture of fracture risk, touching upon different methods and overall findings. We then discuss in a second paragraph the most recent advances in the field and focus on the genetics of fracture risk and also of other endophenotypes closely related to fracture risk such as bone mineral density (BMD). Application of state-of-the-art methodology such as Mendelian randzation in fracture GWAS are reviewed. The final part of this review touches upon potential future directions in genetic research of osteoporotic fractures.

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

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.