Integrating transcriptome-wide association study and mRNA expression profiling identifies novel genes associated with bone mineral density

Uncovering the candidate genes related to sheep body weight using multi-trait genome-wide association analysis

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74-7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10-17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04-77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Transcriptome-wide association study reveals candidate causal genes for lumbar spinal stenosis

Aims

Lumbar spinal stenosis (LSS) is a common skeletal system disease that has been partly attributed to genetic variation. However, the correlation between genetic variation and pathological changes in LSS is insufficient, and it is difficult to provide a reference for the early diagnosis and treatment of the disease.

Methods

We conducted a transcriptome-wide association study (TWAS) of spinal canal stenosis by integrating genome-wide association study summary statistics (including 661 cases and 178,065 controls) derived from Biobank Japan, and pre-computed gene expression weights of skeletal muscle and whole blood implemented in FUSION software. To verify the TWAS results, the candidate genes were furthered compared with messenger RNA (mRNA) expression profiles of LSS to screen for common genes. Finally, Metascape software was used to perform enrichment analysis of the candidate genes and common genes.

Results

TWAS identified 295 genes with permutation p-values < 0.05 for skeletal muscle and 79 genes associated for the whole blood, such as RCHY1 (P_TWAS = 0.001). Those genes were enriched in 112 gene ontology (GO) terms and five Kyoto Encyclopedia of Genes and Genomes pathways, such as 'chemical carcinogenesis - reactive oxygen species' (LogP value = -2.139). Further comparing the TWAS significant genes with the differentially expressed genes identified by mRNA expression profiles of LSS found 18 overlapped genes, such as interleukin 15 receptor subunit alpha (IL15RA) (P_TWAS = 0.040, P_mRNA = 0.010). Moreover, 71 common GO terms were detected for the enrichment results of TWAS and mRNA expression profiles, such as negative regulation of cell differentiation (LogP value = -2.811).

Conclusion

This study revealed the genetic mechanism behind the pathological changes in LSS, and may provide novel insights for the early diagnosis and intervention of LSS.

Upregulated miR-9-5p inhibits osteogenic differentiation of bone marrow mesenchymal stem cells under high glucose treatment

INTRODUCTION

Diabetic osteoporosis (DOP) is a chronic diabetic complication, which is attributed to high glucose (HG)-induced dysfunction of bone marrow mesenchymal stem cells (BMSCs). Studies have revealed that microRNAs (miRNAs) play critical roles in osteogenic differentiation of BMSCs in DOP. Here, the role of miR-9-5p in DOP progression was explored.

MATERIALS AND METHODS

The rat model of DOP was established by intraperitoneal injection of streptozotocin (STZ). BMSCs were treated with high glucose (HG) to establish in vitro models. Gene expression in BMSCs and bone tissues of rats was tested by RT-qPCR. The degree of osteogenic differentiation of BMSCs was examined by Alizarin Red staining and ALP activity analysis. The protein levels of collagen-I (COL1), osteocalcin (OCN), osteopontin (OPN), runt-related transcription factor-2 (RUNX2), and DEAD-Box Helicase 17 (DDX17) in BMSCs were evaluated by western blotting. The interaction between miR-9-5p and DDX17 was identified by luciferase reporter assay. H&E staining was used to test morphological structure of femurs of rats with STZ treatment.

RESULTS

MiR-9-5p was overexpressed in HG-treated BMSCs, while DDX17 was downregulated. Functionally, miR-9-5p knockdown promoted BMSCs osteogenic differentiation under HG condition. Mechanically, miR-9-5p targeted DDX17. DDX17 knockdown reversed the effect of miR-9-5p silencing on osteogenic differentiation of HG-treated BMSCs. In in vivo studies, miR-9-5p downregulation ameliorated the DOP condition of rats and miR-9-5p expression was negatively correlated with DDX17 expression in bone tissues of rats with STZ treatment.

CONCLUSION

MiR-9-5p knockdown promotes HG-induced osteogenic differentiation BMSCs in vitro and mitigates the DOP condition of rats in vivo by targeting DDX17.

A transcriptome-wide association study to detect novel genes for volumetric bone mineral density

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm2) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm3) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation.

Novel partners of USP6 gene in a spectrum of bone and soft tissue lesions

Nodular fasciitis, primary aneurysmal bone cyst, myositis ossificans, and their related lesions are benign tumors that share common histological features and a chromosomal rearrangement involving the ubiquitin-specific peptidase 6 (USP6) gene. The identification of an increasing number of new partners implicated in USP6 rearrangements demonstrates a complex tumorogenesis of this tumor spectrum. In this study on a series of 77 tumors (28 nodular fasciitis, 42 aneurysmal bone cysts, and 7 myositis ossificans) from the database of the French Sarcoma Group, we describe 7 new partners of the USP6 gene. For this purpose, rearrangements were first researched by multiplexed RT-qPCRs in the entire population. A targeted RNA sequencing was then used on samples selected according to a high USP6-transcription level expression estimated by RT-qPCR. Thanks to this multistep approach, besides the common USP6 fusions observed, we detected novel USP6 partners: PDLIM7 and MYL12A in nodular fasciitis and TPM4, DDX17, GTF2I, KLF3, and MEF2A in aneurysmal bone cysts. In order to try to bring to light the role played by the recently identified USP6 partners in this lesional spectrum, their functions are discussed. Taking into account that a traumatic participation has long been mentioned in the histogenesis of most of these lesions and because of their morphological resemblance to organizing granulation reparative tissue or callus, a focus is placed on their relationship with tissue remodeling and, to a lesser extent, with bone metabolism.

Conserved sequence motifs in human TMTC1, TMTC2, TMTC3, and TMTC4, new O-mannosyltransferases from the GT-C/PMT clan, are rationalized as ligand binding sites

BACKGROUND

The human proteins TMTC1, TMTC2, TMTC3 and TMTC4 have been experimentally shown to be components of a new O-mannosylation pathway. Their own mannosyl-transferase activity has been suspected but their actual enzymatic potential has not been demonstrated yet. So far, sequence analysis of TMTCs has been compromised by evolutionary sequence divergence within their membrane-embedded N-terminal region, sequence inaccuracies in the protein databases and the difficulty to interpret the large functional variety of known homologous proteins (mostly sugar transferases and some with known 3D structure).

RESULTS

Evolutionary conserved molecular function among TMTCs is only possible with conserved membrane topology within their membrane-embedded N-terminal regions leading to the placement of homologous long intermittent loops at the same membrane side. Using this criterion, we demonstrate that all TMTCs have 11 transmembrane regions. The sequence segment homologous to Pfam model DUF1736 is actually just a loop between TM7 and TM8 that is located in the ER lumen and that contains a small hydrophobic, but not membrane-embedded helix. Not only do the membrane-embedded N-terminal regions of TMTCs share a common fold and 3D structural similarity with subgroups of GT-C sugar transferases. The conservation of residues critical for catalysis, for binding of a divalent metal ion and of the phosphate group of a lipid-linked sugar moiety throughout enzymatically and structurally well-studied GT-Cs and sequences of TMTCs indicates that TMTCs are actually sugar-transferring enzymes. We present credible 3D structural models of all four TMTCs (derived from their closest known homologues 5ezm/5f15) and find observed conserved sequence motifs rationalized as binding sites for a metal ion and for a dolichyl-phosphate-mannose moiety.

CONCLUSIONS

With the results from both careful sequence analysis and structural modelling, we can conclusively say that the TMTCs are enzymatically active sugar transferases belonging to the GT-C/PMT superfamily. The DUF1736 segment, the loop between TM7 and TM8, is critical for catalysis and lipid-linked sugar moiety binding. Together with the available indirect experimental data, we conclude that the TMTCs are not only part of an O-mannosylation pathway in the endoplasmic reticulum of upper eukaryotes but, actually, they are the sought mannosyl-transferases.

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

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

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis

CONTEXT

Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms.

OBJECTIVE

Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis.

DESIGN, SETTING, AND PARTICIPANTS

Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings.

RESULTS

In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1.

CONCLUSIONS

Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.