An Osteoporosis Risk SNP at 1p36.12 Acts as an Allele-Specific Enhancer to Modulate LINC00339 Expression via Long-Range Loop Formation

Oncogenes and tumor suppressor genes: functions and roles in cancers

Cancer, being the most formidable ailment, has had a profound impact on the human health. The disease is primarily associated with genetic mutations that impact oncogenes and tumor suppressor genes (TSGs). Recently, growing evidence have shown that X-linked TSGs have specific role in cancer progression and metastasis as well. Interestingly, our genome harbors around substantial portion of genes that function as tumor suppressors, and the X chromosome alone harbors a considerable number of TSGs. The scenario becomes even more compelling as X-linked TSGs are adaptive to key epigenetic processes such as X chromosome inactivation. Therefore, delineating the new paradigm related to X-linked TSGs, for instance, their crosstalk with autosome and involvement in cancer initiation, progression, and metastasis becomes utmost importance. Considering this, herein, we present a comprehensive discussion of X-linked TSG dysregulation in various cancers as a consequence of genetic variations and epigenetic alterations. In addition, the dynamic role of X-linked TSGs in sex chromosome-autosome crosstalk in cancer genome remodeling is being explored thoroughly. Besides, the functional roles of ncRNAs, role of X-linked TSG in immunomodulation and in gender-based cancer disparities has also been highlighted. Overall, the focal idea of the present article is to recapitulate the findings on X-linked TSG regulation in the cancer landscape and to redefine their role toward improving cancer treatment strategies.

A regulatory variant impacting TBX1 expression contributes to basicranial morphology in Homo sapiens

Changes in gene regulatory elements play critical roles in human phenotypic divergence. However, identifying the base-pair changes responsible for the distinctive morphology of Homo sapiens remains challenging. Here, we report a noncoding single-nucleotide polymorphism (SNP), rs41298798, as a potential causal variant contributing to the morphology of the skull base and vertebral structures found in Homo sapiens. Screening for differentially regulated genes between Homo sapiens and extinct relatives revealed 13 candidate genes associated with basicranial development, with TBX1, implicated in DiGeorge syndrome, playing a pivotal role. Epigenetic markers and in silico analyses prioritized rs41298798 within a TBX1 intron for functional validation. CRISPR editing revealed that the 41-base-pair region surrounding rs41298798 modulates gene expression at 22q11.21. The derived allele of rs41298798 acts as an allele-specific enhancer mediated by E2F1, resulting in increased TBX1 expression levels compared to the ancestral allele. Tbx1-knockout mice exhibited skull base and vertebral abnormalities similar to those seen in DiGeorge syndrome. Phenotypic differences associated with TBX1 deficiency are observed between Homo sapiens and Neanderthals (Homo neanderthalensis). In conclusion, the regulatory divergence of TBX1 contributes to the formation of skull base and vertebral structures found in Homo sapiens.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

METTL3-stabilized super enhancers-lncRNA SUCLG2-AS1 mediates the formation of a long-range chromatin loop between enhancers and promoters of SOX2 in metastasis and radiosensitivity of nasopharyngeal carcinoma

BACKGROUND

Super enhancers (SE) play pivotal roles in cell identity and diseases occur including tumorigenesis. The depletion of SE-associated lncRNA transcripts, also known as super-lncRNA, causes the activity of SE to be dysregulated.

METHODS

We screened and identified an elevated metastasis-associated SE-lncRNA SUCLG2-AS1 in nasopharyngeal carcinoma (NPC) using RNA-sequencing, real-time quantitative polymerase chain reaction (RT-qPCR) and bioinformatics. Western blotting, RT-qPCR, methylated RNA immunoprecipitation (MeRIP), RNA immunoprecipitation, chromatin immunoprecipitation, RNA pull-down and 3C (chromosome conformation capture assays) were used for mechanistic studies.

RESULTS

SUCLG2-AS1 was correlated with a poor prognosis. SUCLG2-AS1 promotes NPC cell invasion and metastasis while repressing apoptosis and radiosensitivity in vitro and in vivo. Mechanistically, high SUCLG2-AS1 expression occurred in an m6A-dependent manner. SUCLG2-AS1 was found to be located in the SE region of SOX2, and it regulated the expression of SOX2 via long-range chromatin loop formation, which via mediating CTCF (transcription factor) occupied the SE and promoter region of SOX2, thus regulating the metastasis and radiosensitivity of NPC.

CONCLUSIONS

Taken together, our data suggest that SUCLG2-AS1 may serve as a novel intervention target for the clinical treatment of NPC.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

lncRNA MIR31HG Regulates Proliferation and Migration by Targeting Matrix Gla Protein in Nonsyndromic Cleft Lip With or Without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common craniofacial birth defect with complex etiologies. Recently, the dysregulation of long noncoding RNAs (lncRNAs) has been implicated in many developmental diseases, including NSCL/P. However, the functions and mechanisms of lncRNAs in NSCL/P have not been fully elucidated. In this study, we found that lncRNA MIR31HG in NSCL/P patients was significantly downregulated than that in healthy individuals (GSE42589, GSE183527). In addition, single nucleotide polymorphism rs58751040 in MIR31HG was nominally associated with NSCL/P susceptibility (odds ratio: 1.29, 95% confidence interval: 1.03-1.54, p = 4.93 × 10^-2) through a case-control study (504 NSCL/P cases and 455 controls). Luciferase activity assay showed that the C allele of rs58751040 revealed a decreased transcription activity of MIR31HG than the G allele. Moreover, knockdown of MIR31HG promoted cell proliferation and migration in human oral keratinocytes and human embryonic palate mesenchyme. Bioinformatic analysis and cellular studies suggested that MIR31HG may confer susceptibility to risk of NSCL/P through matrix Gla protein (MGP) signaling. In summary, we identified a novel lncRNA involved in the development of NSCL/P.

New insights into the role of long non-coding RNAs in osteoporosis

Osteoporosis is a common disease in elderly individuals, and osteoporosis can easily lead to bone and hip fractures that seriously endanger the health of elderly individuals. At present, the treatment of osteoporosis is mainly anti-osteoporosis drugs, but there are side effects associated with anti-osteoporosis drugs. Therefore, it is very important to develop early diagnostic indicators and new therapeutic drugs for the prevention and treatment of osteoporosis. Long noncoding RNAs (lncRNAs), noncoding RNAs longer than 200 nucleotides, can be used as diagnostic markers for osteoporosis, and lncRNAs play an important role in the progression of osteoporosis. Many studies have shown that lncRNAs can be the target of osteoporosis. Therefore, herein, the role of lncRNAs in osteoporosis is summarized, aiming to provide some information for the prevention and treatment of osteoporosis.

Differential Expression of Non-Coding RNAs in Stem Cell Development and Therapeutics of Bone Disorders

Stem cells' self-renewal and multi-lineage differentiation are regulated by a complex network consisting of signaling factors, chromatin regulators, transcription factors, and non-coding RNAs (ncRNAs). Diverse role of ncRNAs in stem cell development and maintenance of bone homeostasis have been discovered recently. The ncRNAs, such as long non-coding RNAs, micro RNAs, circular RNAs, small interfering RNA, Piwi-interacting RNAs, etc., are not translated into proteins but act as essential epigenetic regulators in stem cells' self-renewal and differentiation. Different signaling pathways are monitored efficiently by the differential expression of ncRNAs, which function as regulatory elements in determining the fate of stem cells. In addition, several species of ncRNAs could serve as potential molecular biomarkers in early diagnosis of bone diseases, including osteoporosis, osteoarthritis, and bone cancers, ultimately leading to the development of new therapeutic strategies. This review aims to explore the specific roles of ncRNAs and their effective molecular mechanisms in the growth and development of stem cells, and in the regulation of osteoblast and osteoclast activities. Furthermore, we focus on and explore the association of altered ncRNA expression with stem cells and bone turnover.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

Osteoarthritis year in review: genetics, genomics and epigenetics

This "year in review" provides a summary of the research findings on the topic of genetics, genomics and epigenetics for osteoarthritis (OA) between Mar 2021-Apr 2022. A search routine of the literature in PubMed for the keyword, osteoarthritis, together with topics on genetics, genomics, epigenetics, polymorphism, DNA methylation, noncoding RNA, lncRNA, proteomics, and single cell RNA sequencing, returned key research articles and relevant reviews. Following filtering of duplicates across search routines, 695 unique research articles and 112 reviews were identified. We manually curated these articles and selected 90 as references for this review. However, we were unable to refer to all these articles, and only used selected articles to highlight key outcomes and trends. The trend in genetics is on the meta-analysis of existing cohorts with comparable genetic and phenotype characterisation of OA; in particular, clear definition of endophenotypes to enhance the genetic power. Further, many researchers are realizing the power of big data and multi-omics approaches to gain molecular insights for OA, and this has opened innovative approaches to include transcriptomics and epigenetics data as quantitative trait loci (QTLs). Given that most of the genetic loci for OA are not located within coding regions of genes, implying the impact is likely to be on gene regulation, epigenetics is a hot topic, and there is a surge in studies relating to the role of miRNA and long non-coding RNA on cartilage biology and pathology. The findings are exciting and new insights are provided in this review to summarize a year of research and the road map to capture all new innovations to achieve the desired goal in OA prevention and treatment.

Nucleus pulposus related lncRNA and mRNA expression profiles in intervertebral disc degeneration

In the present study, we aimed to have a comprehensive understanding of nucleus pulposus related long noncoding RNA (lncRNA) and mRNA expression profiles in intervertebral disc degeneration (IDD). In total, 2418 mRNAs and 528 lncRNAs were found to be differentially expressed in the IDD group compared with the Control group. Combining microarray datasets and sequencing data, 5 overlapping DEMs and 7 overlapping DELs were identified. NF-κB signaling pathway, PI3K-Akt signaling pathway and Wnt/β-catenin signaling pathway were strongly linked with enriched GO terms and KEGG pathways. The ceRNA network suggested that lnc-TMEM44-AS1-hsa-miR-206-HDAC4 may be one crucial axis in IDD. PPI network analysis was constructed with 309 nodes and 129 edges. And the highest connectivity degrees were ALB, APOB and CCL2. This study suggested that specific lncRNAs and ceRNA axes may be crucial in the development of IDD. It provides a new perspective for delaying IDD process and enhancing intervertebral disc repair.

Role of a small GTPase Cdc42 in aging and age-related diseases

A small GTPase, Cdc42 is evolutionarily one of the most ancient members of the Rho family, which is ubiquitously expressed and involved in a wide range of fundamental cellular functions. The crucial role of Cdc42 includes regulation of the actin cytoskeleton, cell polarity, morphology and migration, endocytosis and exocytosis, cell cycle, and proliferation in many different cell types. Many studies have provided compelling yet contradicting evidence that Cdc42 dysregulation plays an important role in cellular and tissue aging. Furthermore, Cdc42 is a critical factor in the development and progression of aging-related pathologies, such as neurodegenerative and cardiovascular disorders, diabetes type 2, and aging-related disorders of the joints and bones, and the inhibition of the Cdc42 demonstrates potentially significant therapeutic and anti-aging effects in animal models of aging and disease. However, regulation of Cdc42 expression and activity is very complex and depends on many factors, such as the origin and complexity of the tissues, hormonal status, etc. Therefore, this review is focused on current advances in understanding the underlying cellular and molecular mechanisms associated with Cdc42 activity and regulation of senescence in different cell types since they may provide a foundation for novel therapeutic strategies and targeted drugs to reverse the aging process and treat aging-associated disorders.

Identification of a Functional Susceptibility Variant for Adolescent Idiopathic Scoliosis that Upregulates Early Growth Response 1 (EGR1)-Mediated UNCX Expression

Adolescent idiopathic scoliosis (AIS) is a serious health problem affecting 3% of live births all over the world. Many loci associated with AIS have been identified by previous genome wide association studies, but their biological implication remains mostly unclear. In this study, we evaluated the AIS-associated variants in the 7p22.3 locus by combining in silico, in vitro, and in vivo analyses. rs78148157 was located in an enhancer of UNCX, a homeobox gene and its risk allele upregulated the UNCX expression. A transcription factor, early growth response 1 (EGR1), transactivated the rs78148157-located enhancer and showed a higher binding affinity for the risk allele of rs78148157. Furthermore, zebrafish larvae with UNCX messenger RNA (mRNA) injection developed body curvature and defective neurogenesis in a dose-dependent manner. rs78148157 confers the genetic susceptibility to AIS by enhancing the EGR1-regulated UNCX expression. © 2022 American Society for Bone and Mineral Research (ASBMR).

A novel lnc-LAMC2-1:1 SNP promotes colon adenocarcinoma progression by targeting miR-216a-3p/HMGB3

Single nucleotide polymorphisms (SNPs) was associated with altering the secondary structure of long non-coding RNA (lncRNA). Increasing reports showed that lnc-LAMC2-1:1 SNP played an important role in cancer development and invasion. This study is to elucidate the molecular function of lnc-LAMC2-1:1 SNP rs2147578 promoting tumor progression in colon adenocarcinoma (COAD). In this study, we found that the lnc-LAMC2-1:1 SNP rs2147578 was upregulated in COAD cell lines. Furthermore, lnc-LAMC2-1:1 SNP rs2147578 promoted colon cancer migration, invasion, and proliferation. Interestingly, lnc-LAMC2-1:1 SNP rs2147578 positively regulated HMGB3 expression via miR-216a-3p in colon cancer cells. Functional enrichment analysis showed that targeting genes of miR-216a-3p were enriched in regulating the pluripotency of stem cells, MAPK signaling pathway, TNF signaling pathway, neurotrophin signaling pathway, relaxin signaling pathway, and FoxO signaling pathway. Tumor Immune Estimation Resource (TIMER) database revealed that there was a significantly positive correlation between HMGB3 expression and the infiltration of CD8⁺ T cells, B cells, neutrophils, macrophages, and CD4⁺ T cells. Finally, HMGB3 overexpression was validated in external data. In conclusions, lnc-LAMC2-1:1 SNP rs2147578 was involved in promoting COAD progression by targeting miR-216a-3p/HMGB3, and this study will provide a novel molecular target for COAD.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

Lineage-selective super enhancers mediate core regulatory circuitry during adipogenic and osteogenic differentiation of human mesenchymal stem cells

Human mesenchymal stem cells (hMSCs) can be differentiated into osteoblasts and adipocytes. During these processes, super enhancers (SEs) play important roles. Here, we performed comprehensive characterization of the SEs changes associated with adipogenic and osteogenic differentiation of hMSCs, and revealed that SEs changed more dramatically compared with typical enhancers. We identified a set of lineage-selective SEs, whose target genes were enriched with cell type-specific functions. Functional experiments in lineage-selective SEs demonstrated their specific roles in directed differentiation of hMSCs. We also found that some key transcription factors regulated by lineage-selective SEs could form core regulatory circuitry (CRC) to regulate each other's expression and control the hMSCs fate determination. In addition, we found that GWAS SNPs of osteoporosis and obesity were significantly enriched in osteoblasts-selective SEs or adipocytes-selective SEs, respectively. Taken together, our studies unveiled important roles of lineage-selective SEs in hMSCs differentiation into osteoblasts and adipocytes.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Identification of Known and Novel Long Noncoding RNAs Potentially Responsible for the Effects of Bone Mineral Density (BMD) Genomewide Association Study (GWAS) Loci

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 noncoding 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 single-nucleotide polymorphisms (SNPs) demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (false discovery rate [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. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Risk SNP-mediated LINC01614 upregulation drives head and neck squamous cell carcinoma progression via PI3K/AKT signaling pathway

As potential biomarkers and therapeutic targets, long noncoding RNAs (lncRNAs) are involved in the tumorigenesis of various tumors. Genetic variation in long noncoding regions can lead to lncRNA dysfunction and even cancer. Nevertheless, studies on the association between lncRNA-associated single-nucleotide polymorphisms (SNPs) and the risk of head and neck squamous cell carcinoma (HNSCC) remain inadequate. Here, we aimed to explore the association between SNPs in LINC01614 and HNSCC risk, and the potential role of LINC01614 in tumorigenesis. In this study, we found that rs16854802 A > G (odds ratio [OR] = 1.42, 95% confidence interval [CI]: 1.22-1.77, p < 0.001) and rs3113503 G > C (OR = 1.38, 95% CI: 1.15-1.64, p < 0.001) in LINC01614 increased the risk of HNSCC in the Chinese population. Functional bioinformatic analysis and luciferase reporter assay revealed that rs3113503 G > C variant disrupted the binding of miRNA-616-3p to LINC01614, which resulted in the increased expression of LINC01614. Further analysis of the TCGA database demonstrated that the upregulated LINC01614 in HNSCC cancer tissues was associated with poor prognostic in HNSCC patients. In vitro experiments showed that knockdown of LINC01614 inhibited the proliferation, invasion, and migration ability of HNSCC cells. Mechanistically, allele C of rs3113503 in LINC01614 was more effective than allele G in activating the PI3K/AKT signaling pathway. Moreover, the reduced expression of LINC01614 also inhibited the activation of the PI3K/AKT signaling pathway. In summary, our findings revealed that the risk SNP rs3113503 G > C in LINC01614 altered the binding to miR-616-3p, which led to increased LINC01614 expression and promoted HNSCC progression by activating the PI3K/AKT signaling pathway.

An Osteoporosis Susceptibility Allele at 11p15 Regulates SOX6 Expression by Modulating TCF4 Chromatin Binding

Osteoporosis is an age-related complex disease clinically diagnosed with bone mineral density (BMD). Although several genomewide association studies (GWASs) have discovered multiple noncoding genetic variants at 11p15 influencing osteoporosis risk, the functional mechanisms of these variants remain unknown. Through integrating bioinformatics and functional experiments, a potential functional single-nucleotide polymorphism (SNP; rs1440702) located in an enhancer element was identified and the A allele of rs1440702 acted as an allelic specificities enhancer to increase its distal target gene SOX6 (~600 Kb upstream) expression, which plays a key role in bone formation. We also validated this long-range regulation via conducting chromosome conformation capture (3C) assay. Furthermore, we demonstrated that SNP rs1440702 with a risk allele (rs1440702-A) could increase the activity of the enhancer element by altering the binding affinity of the transcription factor TCF4, resulting in the upregulation expression of SOX6 gene. Collectively, our integrated analyses revealed how the noncoding genetic variants (rs1440702) affect osteoporosis predisposition via long-range gene regulatory mechanisms and identified its target gene SOX6 for downstream biomarker and drug development. © 2022 American Society for Bone and Mineral Research (ASBMR).

LINC00339: An emerging major player in cancer and metabolic diseases

Long noncoding RNAs (lncRNAs) are a series of RNA molecules without ability to code proteins. LncRNAs have emerged as significant players in almost all aspects of gene function and regulation and play crucial roles in many human diseases. Particular lncRNAs are deemed to be promising molecular biomarkers used for diagnosing diseases and determining patient prognoses and treatment efficacies. LINC00339 is a new budding lncRNA and much of evidence shows that it is abnormally expressed in multifarious diseases, including endometriosis, cardiomyocyte apoptosis, osteoporosis, digestive-system tumors, respiratory-system tumors, nervous-system tumors, and diseases involving other systems. Additionally, LINC00339 is remarkably associated with different clinical features, such as tumor size, TNM stage, and pathological grade. LINC00339 expression has been proved to upregulate in the aforementioned diseases and has been identified to promote disease occurrence and development. It is also reported that LINC00339 is associated with various cellular events, such as tumor cell proliferation, motility and invasiveness, the expression has also been proved that it is closely related to clinical symptoms in cancer patients. This review summarizes the relationships among expression levels, biological features, clinical symptoms, and regulatory mechanisms of LINC00339 in several diseases and discusses the clinical applications of LINC00339 as a cancer diagnostic, prognostic and treatment efficacy biomarker.

Relationship between rs7586085, GALNT3 and CCDC170 gene polymorphisms and the risk of osteoporosis among the Chinese Han population

Osteoporosis (OP) has plagued many women for years, and bone density loss is an indicator of OP. The purpose of this study was to evaluate the relationship between the polymorphism of the rs7586085, CCDC170 and GALNT3 gene polymorphisms and the risk of OP in the Chinese Han population. Using the Agena MassArray method, we identified six candidate SNPs on chromosomes 2 and 6 in 515 patients with OP and 511 healthy controls. Genetic model analysis was performed to evaluate the significant association between variation and OP risk, and meanwhile, the multiple tests were corrected by false discovery rate (FDR). Haploview 4.2 was used for haplotype analysis. In stratified analysis of BMI ˃ 24, rs7586085, rs6726821, rs6710518, rs1346004, and rs1038304 were associated with the risk of OP based on the results of genetic models among females even after the correction of FDR (qd < 0.05). In people at age ≤ 60 years, rs1038304 was associated with an increased risk of OP under genetic models after the correction of FDR (qd < 0.05). Our study reported that GALNT3 and CCDC170 gene polymorphisms and rs7586085 are the effective risk factors for OP in the Chinese Han population.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.

In the latest statistics, the incidence and mortality rate of colorectal cancer (CRC) remains high. Genome-wide association studies (GWAS) have become a powerful tool for identifying genetic susceptibility loci that confer significant risk on disease, and have identified more than 160 risk loci associated with CRC. However, it has proven quite difficult to identify the regulatory variants and target genes involved behind these GWAS signals. Here, we take advantage of multi-omics data and multiple biological experiments to reveal new biological pathways affecting susceptibility to CRC. We show that a specific genetic variant, rs7229639, and three other high linked functional variants (rs77544449, rs60385309 and rs72917785) at 18q21.1 might regulate the expression of LIPG, a gene that was shown to exhibit an oncogenic function by our in-vitro experiments and clinical data analysis. The link between genetic variants, gene expression and CRC phenotype established by us could provide references for follow-up basic and clinical studies.

The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation

Previous powerful genome-wide association studies (GWASs) and whole-genome sequencing have identified multiple single-nucleotide polymorphisms (SNPs) located over 69 kb upstream of CTNNB1 at 3p22.1 locus associated with osteoporosis. The CTNNB1 gene encodes β-catenin that is an integral part of adherens junctions and the primary mediator of the canonical Wnt signaling pathway. The causal variants and underlying molecular mechanisms of the osteoporosis susceptibility locus 3p22.1 remains unknown. Through comprehensive computational analyses, including expression quantitative trait locus (eQTL), high-throughput chromatin interaction (Hi-C), epigenomic and functional annotation, four enhancer SNPs (rs9820407, rs9878224, rs454690 and rs9832204) were prioritized as potential causal SNPs at 3p22.1 for osteoporosis. Rs9820407 displayed the strongest enhancer activity in dual-luciferase assays. Specifically, the minor rs9820407-A can preferentially bind transcription factor FOXC1, elevate the enhancer activity and increase CTNNB1 expression. The architectural protein CTCF was presumably involved in long-range chromatin interaction between rs9820407 and CTNNB1. Our study provided a mechanistic insight into how noncoding enhancer SNP rs9820407 distally regulates CTNNB1 expression and modulates osteoporosis risk.

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.

Single Nucleotide Polymorphisms of Porcine lncMGPF Regulate Meat Production Traits by Affecting RNA Stability

lncMGPF is a novel positive regulator of myogenic differentiation, muscle growth and regeneration in mouse, pig, and human. But whether natural mutations within lncMGPF gene regulate animal meat production traits is unclear. In this study, ten single nucleotide polymorphisms (SNPs) of pig lncMGPF (plncMGPF) gene were identified among commercial pig breeds and Chinese local pig breeds. These SNPs are highly linked and constructed into multiple haplotypes, and haplotype ATTCATGTTC (H1) mainly exists in commercial pig breeds while haplotype GCCTGCACCT (H3) is more frequent in Chinese local pig breeds. Association analysis indicated that all SNPs are significantly associated with the backfat thickness and loin muscle area (P < 0.05), respectively, and homologous H1 individuals have higher loin muscle area and lower backfat thickness than H3 pigs. Bioinformatics and functional analysis showed that haplotype H1 has a longer half-life and more stable RNA secondary structure than haplotype H3. plncMGPF haplotype H1 has stronger effects on pig primary myogenic progenitor cells differentiation and muscle growth than haplotype H3. Further experiments showed that two SNPs (rs81403974 and rs325492834) function together to confer plncMGPF stability and function. Our observation suggested that the SNPs in lncMGPF can change the RNA stabilities and lncMGPF function, thereby affecting the porcine meat production traits.

Recent advances in the epigenetics of bone metabolism

Osteoporosis is a common form of metabolic bone disease that is costly to treat and is primarily diagnosed on the basis of bone mineral density. As the influences of genetic lesions and environmental factors are increasingly studied in the pathological development of osteoporosis, regulated epigenetics are emerging as the important pathogenesis mechanisms in osteoporosis. Recently, osteoporosis genome-wide association studies and multi-omics technologies have revealed that susceptibility loci and the misregulation of epigenetic modifiers are key factors in osteoporosis. Over the past decade, extensive studies have demonstrated epigenetic mechanisms, such as DNA methylation, histone/chromatin modifications, and non-coding RNAs, as potential contributing factors in osteoporosis that affect disease initiation and progression. Herein, we review recent advances in epigenetics in osteoporosis, with a focus on exploring the underlying mechanisms and potential diagnostic/prognostic biomarker applications for osteoporosis.

Long non-coding RNA RP11-84C13.1 promotes osteogenic differentiation of bone mesenchymal stem cells and alleviates osteoporosis progression via the miR-23b-3p/RUNX2 axis

The objective of the present study was to determine the role of RP11-84C13.1 in osteoporosis (OP) and its molecular mechanism. First, clinical samples were collected from OP patients and normal control patients. Human bone marrow stromal cells (hBMSCs) were extracted from femoral head tissues. Runt-related transcription factor 2 (RUNX2) and RP11-84C13.1 serum levels were assessed by reverse transcription-quantitative (RT-q)PCR. Following transfection of pcDNA-RP11-84C13.1, si-RP11-84C13.1, microRNA (miRNA)-23b-3p mimic and miRNA-23b-3p inhibitor, the expression levels of RUNX2 and RP11-84C13.1 were determined by RT-qPCR. In addition, the osteogenic ability of hBMSCs was assessed by Alizarin Red staining. The binding of RP11-84C13.1 to miRNA-23b-3p and the binding of miRNA-23b-3p to RUNX2 was confirmed by dual-luciferase reporter gene assay. Long non-coding RNA (lncRNA) RP11-84C13.1 was significantly downregulated in the serum of OP patients. The osteogenic differentiation-related genes RUNX2 and RP11-84C13.1 were markedly upregulated in a time-dependent manner, while the miRNA-23b-3p level gradually decreased in hBMSCs with the prolongation of osteogenesis. RP11-84C13.1 knockdown inhibited the osteogenic differentiation of hBMSCs. Furthermore, RP11-84C13.1 regulated RUNX2 expression by targeting miRNA-23b-3p. Overexpression of miRNA-23b-3p partially reversed the promoting effect of RP11-84C13.1 on the osteogenesis of hBMSCs. In conclusion, lncRNA RP11-84C13.1 upregulated RUNX2 by absorbing miRNA-23b-3p, and thus induced hBMSC osteogenesis to alleviate osteoporosis.

CRISPR-Cas9-Mediated Genome Editing Confirms EPDR1 as an Effector Gene at the BMD GWAS-Implicated 'STARD3NL' Locus

Genome-wide-association studies (GWASs) have discovered genetic signals robustly associated with BMD, but typically not the precise localization of effector genes. By intersecting genome-wide promoter-focused Capture C and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data generated in human mesenchymal progenitor cell (hMSC)-derived osteoblasts, consistent contacts were previously predicted between the EPDR1 promoter and multiple BMD-associated candidate causal variants at the 'STARD3NL' locus. RNAi knockdown of EPDR1 expression in hMSC-derived osteoblasts was shown to lead to inhibition of osteoblastogenesis. To fully characterize the physical connection between these putative noncoding causal variants at this locus and the EPDR1 gene, clustered regularly interspaced short-palindromic repeat Cas9 endonuclease (CRISPR-Cas9) genome editing was conducted in hFOB1.19 cells across the single open-chromatin region harboring candidates for the underlying causal variant, rs1524068, rs6975644, and rs940347, all in close proximity to each other. RT-qPCR and immunoblotting revealed dramatic and consistent downregulation of EPDR1 specifically in the edited differentiated osteoblast cells. Consistent with EPDR1 expression changes, alkaline phosphatase staining was also markedly reduced in the edited differentiated cells. Collectively, CRISPR-Cas9 genome editing in the hFOB1.19 cell model supports previous observations, where this regulatory region harboring GWAS-implicated variation operates through direct long-distance physical contact, further implicating a key role for EPDR1 in osteoblastogenesis and BMD determination. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

STAT3-induced ZBED3-AS1 promotes the malignant phenotypes of melanoma cells by activating PI3K/AKT signaling pathway

Melanoma is considered as the most frequent primary malignancy occurring in skin. Accumulating studies have suggested that long non-coding RNAs (lncRNAs) play critical parts in multiple cancers. In this study, we explored the molecular mechanism of ZBED3 antisense RNA 1 (ZBED3-AS1) in melanoma. We observed that ZBED3-AS1 expression was remarkably up-regulated in melanoma tissues, and high ZBED3-AS1 level was linked to unsatisfactory survival of melanoma patients. Then, we discovered that ZBED3-AS1 was overexpressed in melanoma cells compared with human epidermal melanocytes. In addition, loss-of-function assays verified that ZBED3-AS1 knockdown restrained cell proliferation, migration, epithelial-mesenchymal transition (EMT), and stemness in melanoma. In addition, signal transducer and activator of transcription 3 (STAT3), which also showed tumour-facilitating functions in melanoma, was confirmed as a transcriptional activator of ZBED3-AS1. Moreover, ZBED3-AS1 enhanced the expression of AT-rich interaction domain 4B (ARID4B) through sequestering miR-381-3p. Importantly, we further confirmed that ZBED3-AS1 promoted the malignant progression of melanoma by regulating miR-381-3p/ARID4B axis to activate the phosphatidylinositol 3-kinase/AKT serine/threonine kinase (PI3K/AKT) signalling pathway. In a word, our research might provide a novel therapeutic target for melanoma.

An Intronic Risk SNP rs12454712 for Central Obesity Acts As an Allele-Specific Enhancer To Regulate BCL2 Expression

Genome-wide association studies (GWAS) have reproducibly associated the single nucleotide polymorphism (SNP) rs12454712 with waist-to-hip ratio adjusted for BMI (WHRadjBMI), but the functional role underlying this intronic variant is unknown. Integrative genomic and epigenomic analyses supported rs12454712 as a functional independent variant. We further demonstrated that rs12454712 acted as an allele-specific enhancer regulating expression of its located gene BCL2 by using dual-luciferase reporter assays and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9. Specifically, the rs12454712-C allele can bind transcription factor ZNF329, which efficiently elevates the enhancer activity and increases BCL2 expression. Knocking down Bcl2 in 3T3-L1 cells led to the downregulation of adipogenic differentiation marker genes and increased cell apoptosis. A significant negative correlation between BCL2 expression in subcutaneous adipose tissues and obesity was observed. Our findings illustrate the molecular mechanisms behind the intronic SNP rs12454712 for central obesity, which would be a potential and promising target for developing appropriate therapies.

Polymorphism of rs6426749 at 1p36.12 is associated with the risk of osteoarthritis in Taiwanese female population

BACKGROUND

Osteoarthritis (OA) is a multifactorial disease that is associated with several genetic factors. TFAP2A with a motif of C allele at rs6426749 demonstrates a higher binding ability, thereby increasing CDC42 expression, potentially affecting OA occurrence. In this study, we evaluated the role of rs6426749 polymorphisms on knee OA in a female Taiwanese population.

METHODS

We performed a case-control study of 368 OA cases and 379 controls between March 2017 and October 2018. Knee OA was defined using the Kellgren-Lawrence grading system, and genotypes were determined using the Sequenom MassArray iPLEX Gold assay. Stratified sex and body mass index (BMI) analyses were performed using logistic regression to explore interactions between genes and the environment. We also used expression quantitative trait loci data from the genotype-tissue expression project to conduct functional analyses.

RESULTS

The C allele of rs6426749 was associated with the risk of knee OA (odds ratio [OR] = 1.31, 95% confidence interval [CI], 1.01-1.71; p = 0.042), after adjusting for gender, age, and BMI. In addition, subgroup analyses indicated that females expressing C alleles showed an increased risk for knee OA (OR = 1.56; 95% CI, 1.12-2.18; p = 0.009). Females with a normal BMI and the C allele had the highest OA risk (OR = 1.73; 95% CI, 1.08-2.76; p = 0.022).

CONCLUSION

Our findings indicated that rs6426749 may be related to OA susceptibility in the Taiwanese population. This was particularly true for women with normal BMI.

Editing GWAS: experimental approaches to dissect and exploit disease-associated genetic variation

Genome-wide association studies (GWAS) have uncovered thousands of genetic variants that influence risk for human diseases and traits. Yet understanding the mechanisms by which these genetic variants, mainly noncoding, have an impact on associated diseases and traits remains a significant hurdle. In this review, we discuss emerging experimental approaches that are being applied for functional studies of causal variants and translational advances from GWAS findings to disease prevention and treatment. We highlight the use of genome editing technologies in GWAS functional studies to modify genomic sequences, with proof-of-principle examples. We discuss the challenges in interrogating causal variants, points for consideration in experimental design and interpretation of GWAS locus mechanisms, and the potential for novel therapeutic opportunities. With the accumulation of knowledge of functional genetics, therapeutic genome editing based on GWAS discoveries will become increasingly feasible.

Opportunities and Challenges in Functional Genomics Research in Osteoporosis: Report From a Workshop Held by the Causes Working Group of the Osteoporosis and Bone Research Academy of the Royal Osteoporosis Society on October 5th 2020

The discovery that sclerostin is the defective protein underlying the rare heritable bone mass disorder, sclerosteosis, ultimately led to development of anti-sclerostin antibodies as a new treatment for osteoporosis. In the era of large scale GWAS, many additional genetic signals associated with bone mass and related traits have since been reported. However, how best to interrogate these signals in order to identify the underlying gene responsible for these genetic associations, a prerequisite for identifying drug targets for further treatments, remains a challenge. The resources available for supporting functional genomics research continues to expand, exemplified by "multi-omics" database resources, with improved availability of datasets derived from bone tissues. These databases provide information about potential molecular mediators such as mRNA expression, protein expression, and DNA methylation levels, which can be interrogated to map genetic signals to specific genes based on identification of causal pathways between the genetic signal and the phenotype being studied. Functional evaluation of potential causative genes has been facilitated by characterization of the "osteocyte signature", by broad phenotyping of knockout mice with deletions of over 7,000 genes, in which more detailed skeletal phenotyping is currently being undertaken, and by development of zebrafish as a highly efficient additional in vivo model for functional studies of the skeleton. Looking to the future, this expanding repertoire of tools offers the hope of accurately defining the major genetic signals which contribute to osteoporosis. This may in turn lead to the identification of additional therapeutic targets, and ultimately new treatments for osteoporosis.

Elucidating the role of long intergenic non-coding RNA 339 in human endometrium and endometriosis

Endometriosis is a complex disease, influenced by genetic factors. Genetic markers associated with endometriosis exist at chromosome 1p36.12 and lead to altered expression of the long intergenic non-coding RNA 339 (LINC00339), however, the role of LINC00339 in endometriosis pathophysiology remains unknown. The aim of this work was to characterize the expression patterns of LINC00339 mRNA in endometrium and endometriotic lesions in situ and to determine the functional role of LINC00339 in human endometrium. We employed RNA-sequencing (RNA-seq), quantitative RT-PCR and in situ hybridization to investigate the abundance of LINC00339 transcripts in endometrium and endometrial cell lines and to describe the pattern and localization of LINC00339 expression in endometrium and endometriotic lesions. LINC00339 mRNA expression was manipulated (overexpressed and silenced) in endometrial stromal cell lines and RNA-seq data from overexpression models were analysed using online bioinformatics platforms (STRING and Ingenuity Pathway Analysis) to determine functional processes. We demonstrated the expression of LINC00339 in endometriotic lesions for the first time; we found LINC00339 expression was restricted to the lesion foci and absent in surrounding non-lesion tissue. Furthermore, manipulation of LINC00339 expression in endometrial stromal cell lines significantly impacted the expression of genes involved in immune defence pathways. These studies identify a novel mechanism for LINC00339 activity in endometrium and endometriosis, paving the way for future work, which is essential for understanding the pathogenesis of endometriosis.

Functional genomics of autoimmune diseases

For more than a decade, genome-wide association studies have been applied to autoimmune diseases and have expanded our understanding on the pathogeneses. Genetic risk factors associated with diseases and traits are essentially causative. However, elucidation of the biological mechanism of disease from genetic factors is challenging. In fact, it is difficult to identify the causal variant among multiple variants located on the same haplotype or linkage disequilibrium block and thus the responsible biological genes remain elusive. Recently, multiple studies have revealed that the majority of risk variants locate in the non-coding region of the genome and they are the most likely to regulate gene expression such as quantitative trait loci. Enhancer, promoter and long non-coding RNA appear to be the main target mechanisms of the risk variants. In this review, we discuss functional genetics to challenge these puzzles.

Identification of Critical Functional Modules and Signaling Pathways in Osteoporosis

Background:

Osteoporosis is the most common bone metabolic disease. Abnormal osteoclast formation and resorption play a fundamental role in osteoporosis pathogenesis. Recent researches have greatly broadened our understanding of molecular mechanisms of osteoporosis. However, the molecular mechanisms leading to osteoporosis are still not entirely clear.

Objective:

The purpose of this work is to study the critical regulatory genes, functional modules, and signaling pathways.

Methods:

Differential expression analysis, network topology-based analysis, and overrepresentation enrichment analysis (ORA) were used to identify differentially expressed genes (DEGs), gene subnetworks, and signaling pathways related to osteoporosis, respectively.

Results:

Differential expression analysis identified DEGs, such as POGLUT1, DAPK3 and NFKBIA, associated with osteoclastogenesis, which highlighted Notch, apoptosis and NF-kB signaling pathways. Network topology-based analysis identified the upregulated subnetwork characterized by EXOSC8 and DIS3L from the RNA exosome complex, and the downregulated subnetwork composed of histone deacetylases and the cofactors, MORF4L1 and JDP2. Furthermore, the overrepresentation enrichment analysis highlighted that corticotrophin-releasing hormone signaling pathway might affect osteoclastogenesis through its component NR4A1, and suppressing osteoclast differentiation and osteoclast bone resorption with urocortin (UCN).

Conclusion:

Our systematic analysis not only discovered novel molecular mechanisms but also proposed potential drug targets for osteoporosis.

SNP-adjacent super enhancer network mediates enhanced osteogenic differentiation of MSCs in ankylosing spondylitis

Ankylosing spondylitis (AS) is a rheumatic disease with pathological osteogenesis that causes bony ankylosis and even deformity over time. Mesenchymal stem cells (MSCs) are multipotent stem cells that are the main source of osteoblasts. We previously demonstrated that enhanced osteogenic differentiation of MSCs from AS patients (ASMSCs) is related to pathological osteogenesis in AS. However, the more concrete mechanism needs further exploration. Super enhancers (SEs) are dense clusters of stitched enhancers that control cell identity determination and disease development. Single-nucleotide polymorphisms (SNPs) regulate the formation and interaction of SEs and denote genes accounting for AS susceptibility. Via integrative analysis of multiomic data, including histone 3 lysine 27 acetylation (H3K27ac), chromatin immunoprecipitation sequencing (ChIP-seq), SNPs and RNA sequencing (RNA-seq) data, we discovered a transcription network mediated by AS SNP-adjacent SEs (SASEs) in ASMSCs and identified key genes, such as Toll-like receptor 4 (TLR4), interleukin 18 receptor 1 (IL18R1), insulin-like growth factor binding protein 4 (IGFBP4), transportin 1 (TNPO1) and proprotein convertase subtilisin/kexin type 5 (PCSK5), which are pivotal in osteogenesis and AS pathogenesis. The SASE-regulated network modulates the enhanced osteogenic differentiation of ASMSCs by synergistically activating the PI3K-Akt, NF-kappaB and Hippo signaling pathways. Our results emphasize the crucial role of the SASE-regulated network in pathological osteogenesis in AS, and the preferential inhibition of ASMSC osteogenic differentiation by JQ1 indicates that SEs may be attractive targets in future treatment for new bone formation in AS.

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.

LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60 000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

Wnt-regulated lncRNA discovery enhanced by in vivo identification and CRISPRi functional validation

BACKGROUND

Wnt signaling is an evolutionarily conserved developmental pathway that is frequently hyperactivated in cancer. While multiple protein-coding genes regulated by Wnt signaling are known, the functional lncRNAs regulated by Wnt signaling have not been systematically characterized.

METHODS

We comprehensively mapped Wnt-regulated lncRNAs from an orthotopic Wnt-addicted pancreatic cancer model and examined the response of lncRNAs to Wnt inhibition between in vivo and in vitro cancer models. We further annotated and characterized these Wnt-regulated lncRNAs using existing genomic classifications (using data from FANTOM5) in the context of Wnt signaling and inferred their role in cancer pathogenesis (using GWAS and expression data from the TCGA). To functionally validate Wnt-regulated lncRNAs, we performed CRISPRi screens to assess their role in cancer cell proliferation both in vivo and in vitro.

RESULTS

We identified 3633 lncRNAs, of which 1503 were regulated by Wnt signaling in an orthotopic Wnt-addicted pancreatic cancer model. These lncRNAs were much more sensitive to changes in Wnt signaling in xenografts than in cultured cells. Our analysis suggested that Wnt signaling inhibition could influence the co-expression relationship of Wnt-regulated lncRNAs and their eQTL-linked protein-coding genes. Wnt-regulated lncRNAs were also implicated in specific gene networks involved in distinct biological processes that contribute to the pathogenesis of cancers. Consistent with previous genome-wide lncRNA CRISPRi screens, around 1% (13/1503) of the Wnt-regulated lncRNAs were found to modify cancer cell growth in vitro. This included CCAT1 and LINC00263, previously reported to regulate cancer growth. Using an in vivo CRISPRi screen, we doubled the discovery rate, identifying twice as many Wnt-regulated lncRNAs (25/1503) that had a functional effect on cancer cell growth.

CONCLUSIONS

Our study demonstrates the value of studying lncRNA functions in vivo, provides a valuable resource of lncRNAs regulated by Wnt signaling, and establishes a framework for systematic discovery of functional lncRNAs.

SNPs in Sites for DNA Methylation, Transcription Factor Binding, and miRNA Targets Leading to Allele-Specific Gene Expression and Contributing to Complex Disease Risk: A Systematic Review

INTRODUCTION

The complex genetic diversity among human populations results from an assortment of factors acting at various sequential levels, including mutations, population migrations, genetic drift, and selection. Although there are a plethora of DNA sequence variations identified through genome-wide association studies (GWAS), the challenge remains to explain the mechanisms underlying interindividual phenotypic disparity accounting for disease susceptibility. Single nucleotide polymorphisms (SNPs) present in the sites for DNA methylation, transcription factor (TF) binding, or miRNA targets can alter the gene expression. The systematic review aimed to evaluate the complex crosstalk among SNPs, miRNAs, DNA methylation, and TFs for complex multifactorial disease risk.

METHODS

PubMed and Scopus databases were used from inception until May 15, 2019. Initially, screening of articles involved studies assessing the interaction of SNPs with TFs, DNA methylation, or miRNAs resulting in allele-specific gene expression in complex multifactorial diseases. We also included the studies which provided experimental validation of the interaction of SNPs with each of these factors. The results from various studies on multifactorial diseases were assessed.

RESULTS

A total of 11 articles for SNPs interacting with DNA methylation, 30 articles for SNPs interacting with TFs, and 11 articles for SNPs in miRNA binding sites were selected. The interactions of SNPs with epigenetic factors were found to be implicated in different types of cancers, autoimmune diseases, cardiovascular diseases, diabetes, and asthma.

CONCLUSION

The systematic review provides evidence for the interplay between genetic and epigenetic risk factors through allele-specific gene expression in various complex multifactorial diseases.

Multiomics dissection of molecular regulatory mechanisms underlying autoimmune-associated noncoding SNPs

More than 90% of autoimmune-associated variants are located in noncoding regions, leading to challenges in deciphering the underlying causal roles of functional variants and genes and biological mechanisms. Therefore, to reduce the gap between traditional genetic findings and mechanistic understanding of disease etiologies and clinical drug development, it is important to translate systematically the regulatory mechanisms underlying noncoding variants. Here, we prioritized functional noncoding SNPs with regulatory gene targets associated with 19 autoimmune diseases by incorporating hundreds of immune cell-specific multiomics data. The prioritized SNPs are associated with transcription factor (TF) binding, histone modification, or chromatin accessibility, indicating their allele-specific regulatory roles. Their target genes are significantly enriched in immunologically related pathways and other known immunologically related functions. We found that 90.1% of target genes are regulated by distal SNPs involving several TFs (e.g., the DNA-binding protein CCCTC-binding factor [CTCF]), suggesting the importance of long-range chromatin interaction in autoimmune diseases. Moreover, we predicted potential drug targets for autoimmune diseases, including 2 genes (NFKB1 and SH2B3) with known drug indications on other diseases, highlighting their potential drug repurposing opportunities. Taken together, these findings may provide useful information for future experimental follow-up and drug applications on autoimmune diseases.

lncRNA Neat1 Stimulates Osteoclastogenesis Via Sponging miR-7

Increasing evidence uncover the essential role of long noncoding RNA (lncRNAs) in bone metabolism and the association of lncRNA with genetic risk of osteoporosis. However, whether lncRNA nuclear paraspeckle assembly transcript 1 (Neat1) is involved remains largely unknown. In the present study, we found that Neat1 is induced by osteoclastic differentiation stimuli. Knockdown of Neat1 attenuates osteoclast formation whereas overexpression of Neat1 accelerates osteoclast formation. In vivo evidence showed that enhanced Neat1 expression stimulates osteoclastogenesis and reduces bone mass in mice. Mechanically, Neat1 competitively binds with microRNA 7 (miR-7) and blocks its function for regulating protein tyrosine kinase 2 (PTK2). Intergenic SNP rs12789028 acts as allele-specific long-range enhancer for NEAT1 via chromatin interactions. We establish for the first time that Neat1 plays an essential role in osteoclast differentiation, and provide genetic mechanism underlying the association of NEAT1 locus with osteoporosis risk. These results enrich the current knowledge of NEAT1 function, and uncover the potential of NEAT1 as a therapeutic target for osteoporosis. © 2020 American Society for Bone and Mineral Research.

Significance of Single-Nucleotide Variants in Long Intergenic Non-protein Coding RNAs

Single-nucleotide variants (SNVs) are the most common genetic variants and universally present in the human genome. Genome-wide association studies (GWASs) have identified a great number of disease or trait-associated variants, many of which are located in non-coding regions. Long intergenic non-protein coding RNAs (lincRNAs) are the major subtype of long non-coding RNAs; lincRNAs play crucial roles in various disorders and cellular models via multiple mechanisms. With rapid growth in the number of the identified lincRNAs and genetic variants, there is great demand for an investigation of SNVs in lincRNAs. Hence, in this article, we mainly summarize the significant role of SNVs within human lincRNA regions. Some pivotal variants may serve as risk factors for the development of various disorders, especially cancer. They may also act as important regulatory signatures involved in the modulation of lincRNAs in a tissue- or disorder-specific manner. An increasing number of researches indicate that lincRNA variants would potentially provide additional options for genetic testing and disease risk assessment in the personalized medicine era.

A Non-Coding Disease Modifier of Pancreatic Agenesis Identified by Genetic Correction in a Patient-Derived iPSC Line

GATA6 is a critical regulator of pancreatic development, with heterozygous mutations in this transcription factor being the most common cause of pancreatic agenesis. To study the variability in disease phenotype among individuals harboring these mutations, a patient-induced pluripotent stem cell model was used. Interestingly, GATA6 protein expression remained depressed in pancreatic progenitor cells even after correction of the coding mutation. Screening the regulatory regions of the GATA6 gene in these patient cells and 32 additional agenesis patients revealed a higher minor allele frequency of a SNP 3' of the GATA6 coding sequence. Introduction of this minor allele SNP by genome editing confirmed its functionality in depressing GATA6 expression and the efficiency of pancreas differentiation. This work highlights a possible genetic modifier contributing to pancreatic agenesis and demonstrates the usefulness of using patient-induced pluripotent stem cells for targeted discovery and validation of non-coding gene variants affecting gene expression and disease penetrance.

Long Non-coding RNAs Mechanisms of Action in HIV-1 Modulation and the Identification of Novel Therapeutic Targets

This review aims to highlight the role of long non-coding RNAs in mediating human immunodeficiency virus (HIV-1) viral replication, latency, disease susceptibility and progression. In particular, we focus on identifying possible lncRNA targets and their purported mechanisms of action for future drug design or gene therapeutics.