Dataset Information

Identification of known and novel long non-coding RNAs potentially responsible for the effects of BMD GWAS loci

ABSTRACT: Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long non-coding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor SNPs demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (FDR correction<0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability (RCP)>0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis.

ORGANISM(S): Homo sapiens

PROVIDER: GSE186922 | GEO | 2021/11/04

REPOSITORIES: GEO

ACCESS DATA

Similar Datasets

Project description:Genome-wide association studies (GWASs) have revealed thousands of associations in many complex traits and diseases. Previous studies suggest that a subset of associations are due to alterations in splicing; however, interpreting the effects of splicing on protein isoforms is hindered by limitations in defining full-length transcript isoforms using short-read RNA-seq data. Long-read RNA-seq represents a powerful approach to define and quantify transcript isoforms. In this study, we developed a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease relevant model to infer the effects of sQTL on the ultimate protein isoform products they encode. Such information enables identification of genes potentially responsible for GWAS associations. As a proof-of-concept, we generated deep coverage (N=~22 million full-length reads) PacBio long-read RNAseq data on human fetal osteoblasts (hFOBs), a cell-line of relevance to the regulation of bone mineral density (BMD). We identified 68,326 protein-coding isoforms, including 17,375 (25%) which were novel. Next, we used Bayesian colocalization to identify 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H4PP > 0.75). A total of 836 junctions with colocalizing sQTLs in 459 (of the 732) genes were expressed in hFOB long-read RNA-seq data. With these data, we formulated hypotheses regarding the potential mechanism of action of each sQTL. For example, we identified 7 junctions with colocalizing sQTLs (maximum H4PP = 0.98-0.99) in TPM2 for splice junctions between two nearly mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in hFOBs showed two TPM2 isoforms with opposing effects on mineralization. Our results suggest that splicing is a major mechanism underlying GWAS associations and long-read proteogenomics data is critical to precisely define the protein isoforms that are produced from splicing alterations.

Dataset Information

Identification of known and novel long non-coding RNAs potentially responsible for the effects of BMD GWAS loci

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure