A landscape of gene expression regulation for synovium in arthritis.

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

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.

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.

Application of Hi-C technology in three-dimensional genomics research and disease pathogenesis analysis.

3D genomics aims to investigate the spatial structure of chromatin in the nucleus on the basis of genomic sequences, gene structures and relevant regulatory elements. The spatial organization of chromosomes is fundamental for gene expression regulation. Recent advances of high-throughput chromosome conformation capture (Hi-C) technology and its derivatives, has enabled capture of chromatin architecture with high resolution. In this review, we summarize the development and applications of various technologies of 3D genomes in disease research, particularly in the elucidation of pathogenic mechanisms in cancers and other systemic disorders.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis.

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

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.

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).

DDRS: Detection of drug response SNPs specifically in patients receiving drug treatment.

Detecting SNPs associated with drug efficacy or toxicity is helpful to facilitate personalized medicine. Previous studies usually find SNPs associated with clinical outcome only in patients received a specific treatment. However, without information from patients without drug treatment, it is possible that the detected SNPs are associated with patients' clinical outcome even without drug treatment. Here we aimed to detect drug response SNPs based on data from patients with and without drug treatment through combing the cox proportional-hazards model and pairwise Kaplan-Meier survival analysis. A pipeline named Detection of Drug Response SNPs (DDRS) was built and applied to TCGA breast cancer data including 363 patients with doxorubicin treatment and 321 patients without any drug treatment. We identified 548 doxorubicin associated SNPs. Drug response score derived from these SNPs were associated with drug-resistant level (indicated by IC50) of breast cancer cell lines. Enrichment analyses showed that these SNPs were enriched in active epigenetic regulation markers (e.g., H3K27ac). Compared with random genes, the cis-eQTL genes of these SNPs had a shorter protein-protein interaction distance to doxorubicin associated genes. In addition, linear discriminant analysis showed that the eQTL gene expression levels could be used to predict clinical outcome for patients with doxorubicin treatment (AUC = 0.738). Specifically, we identified rs2817101 as a drug response SNP for doxorubicin treatment. Higher expression level of its cis-eQTL gene GSTA1 is associated with poorer survival. This approach can also be applied to identify new drug associated SNPs in other cancers.

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.

An integrative multi-omics network-based approach identifies key regulators for breast cancer.

Although genome-wide association studies (GWASs) have successfully identified thousands of risk variants for human complex diseases, understanding the biological function and molecular mechanisms of the associated SNPs involved in complex diseases is challenging. Here we developed a framework named integrative multi-omics network-based approach (IMNA), aiming to identify potential key genes in regulatory networks by integrating molecular interactions across multiple biological scales, including GWAS signals, gene expression-based signatures, chromatin interactions and protein interactions from the network topology. We applied this approach to breast cancer, and prioritized key genes involved in regulatory networks. We also developed an abnormal gene expression score (AGES) signature based on the gene expression deviation of the top 20 rank-ordered genes in breast cancer. The AGES values are associated with genetic variants, tumor properties and patient survival outcomes. Among the top 20 genes, RNASEH2A was identified as a new candidate gene for breast cancer. Thus, our integrative network-based approach provides a genetic-driven framework to unveil tissue-specific interactions from multiple biological scales and reveal potential key regulatory genes for breast cancer. This approach can also be applied in other complex diseases such as ovarian cancer to unravel underlying mechanisms and help for developing therapeutic targets.

Fibroblast growth factor receptor signaling as therapeutic targets in female reproductive system cancers.

Ovarian cancer, cervical cancer and endometrial cancer are three relatively common malignant cancers of the female reproductive system. Despite improvements in female genital tract cancer detection and development of new therapeutic approaches, there are still poor prognoses and some do not respond to therapeutic patterns, displaying low survival and high frequency of recurrence. In an era of personalized medicine, novel therapeutic approaches with greater efficacy for these cancers represent an unmet need. One of the actionable signaling pathways is the fibroblast growth factor receptor (FGFR) signaling pathway. Several mutations and alterations in FGF/FGFR family members have been reported in human cancers. FGF/FGFR signaling pathway has become a new target for cancer therapy. This review will summarize the role of FGFR pathway and the genetic alterations of the FGF/FGFR related to female reproductive system cancer. We will describe the available inhibitors of FGFR pathway for potential treatment of female reproductive system cancer. Furthermore, we will discuss FGFR-targeted therapies under clinical development for treatment of female reproductive system cancer.

The osteoporosis susceptible SNP rs4325274 remotely regulates the SOX6 gene through enhancers.

Osteoporosis is a typical polygenic disease, and its heritability is as high as 85%. The incidence of osteoporosis has jumped to the fifth among the common diseases. Although a large number of osteoporosis-susceptible SNPs have been identified, most of them are in the non-coding regions of the genome and the functional mechanisms are unknown. The purpose of this study was to explore the function of non-coding osteoporosis-susceptible SNP rs4325274 and dissect the molecular regulatory mechanisms through integrating bioinformatics analysis and functional experiments. Firstly, we found the SNP rs4325274 resided in a putative enhancer element through functional annotation. eQTL and Hi-C analysis found that the SOX6 gene might be a potential distal target of rs4325274. We conducted the motif prediction using multiple databases and verified the result using ChIP-seq data from GEO database. The result showed that the transcription factor HNF1A could preferentially bind to SNP rs4325274-G allele. We further demonstrated that SNP rs4325274 acted as an enhancer regulating SOX6 gene expression by using dual-luciferase reporter assays. Knockdown of HNF1A decreased the SOX6 gene expression. Taken together, our results uncovered a new mechanism of a non-coding functional SNP rs4325274 as a distal enhancer to modulate SOX6 expression, which provides new insights into deciphering molecular regulatory mechanisms underlying non-coding susceptibility SNPs on complex diseases.

An Allele-Specific Functional SNP Associated with Two Systemic Autoimmune Diseases Modulates IRF5 Expression by Long-Range Chromatin Loop Formation.

Both systemic lupus erythematosus (SLE) and systemic sclerosis (SSc) are autoimmune diseases sharing similar genetic backgrounds. Genome-wide association studies have constantly disclosed numerous genetic variants conferring to both disease risks at 7q32.1, but the functional mechanisms underlying them are still largely unknown. Through a series of bioinformatics and functional analyses, we prioritized a potential independent functional single-nucleotide polymorphism (rs13239597) within TNPO3 promoter region, residing in a putative enhancer element and validated that IRF5 is the distal target gene (∼118 kb) of rs13239597, which is a key regulator involved in pathogenic autoantibody dysregulation, increasing risk of both SLE and SSc. We experimentally validated the long-range chromatin interactions between rs13239597 and IRF5 using chromosome conformation capture assay. We further demonstrated that rs13239597-A acted as an allele-specific enhancer regulating IRF5 expression, independently of TNPO3 by using dual-luciferase reporter assays and CRISPR-Cas9. Particularly, the transcription factor EVI1 could preferentially bind to rs13239597-A allele and increase the enhancer activity to regulate IRF5 expression. Taken together, our results uncovered a mechanistic insight of a noncoding functional variant acting as an allele-specific distal enhancer to directly modulate IRF5 expression, which might obligate in understanding of complex genetic architectures of SLE and SSc pathogenesis.

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

Genome-wide association studies (GWASs) have reproducibly associated variants within intergenic regions of 1p36.12 locus with osteoporosis, but the functional roles underlying these noncoding variants are unknown. Through an integrative functional genomic and epigenomic analyses, we prioritized rs6426749 as a potential causal SNP for osteoporosis at 1p36.12. Dual-luciferase assay and CRISPR/Cas9 experiments demonstrate that rs6426749 acts as a distal allele-specific enhancer regulating expression of a lncRNA (LINC00339) (∼360 kb) via long-range chromatin loop formation and that this loop is mediated by CTCF occupied near rs6426749 and LINC00339 promoter region. Specifically, rs6426749-G allele can bind transcription factor TFAP2A, which efficiently elevates the enhancer activity and increases LINC00339 expression. Downregulation of LINC00339 significantly increases the expression of CDC42 in osteoblast cells, which is a pivotal regulator involved in bone metabolism. Our study provides mechanistic insight into how a noncoding SNP affects osteoporosis by long-range interaction, a finding that could indicate promising therapeutic targets for osteoporosis.

Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer.

RANKL is a key regulator involved in bone metabolism, and a drug target for osteoporosis. The clinical diagnosis and assessment of osteoporosis are mainly based on bone mineral density (BMD). Previous powerful genomewide association studies (GWASs) have identified multiple intergenic single-nucleotide polymorphisms (SNPs) located over 100 kb upstream of RANKL and 65 kb downstream of AKAP11 at 13q14.11 for osteoporosis. Whether these SNPs exert their roles on osteoporosis through RANKL is unknown. In this study, we conducted integrative analyses combining expression quantitative trait locus (eQTL), genomic chromatin interaction (high-throughput chromosome conformation capture [Hi-C]), epigenetic annotation, and a series of functional assays. The eQTL analysis identified six potential functional SNPs (rs9533090, rs9594738, r8001611, rs9533094, rs9533095, and rs9594759) exclusively correlated with RANKL gene expression (p < 0.001) at 13q14.11. Co-localization analyses suggested that eQTL signal for RANKL and BMD-GWAS signal shared the same causal variants. Hi-C analysis and functional annotation further validated that the first five osteoporosis SNPs are located in a super-enhancer region to regulate the expression of RANKL via long-range chromosomal interaction. Particularly, dual-luciferase assay showed that the region harboring rs9533090 in the super-enhancer has the strongest enhancer activity, and rs9533090 is an allele-specific regulatory SNP. Furthermore, deletion of the region harboring rs9533090 using CRISPR/Cas9 genome editing significantly reduced RANKL expression in both mRNA level and protein level. Finally, we found that the rs9533090-C robustly recruits transcription factor NFIC, which efficiently elevates the enhancer activity and increases the RANKL expression. In summary, we provided a feasible method to identify regulatory noncoding SNPs to distally regulate their target gene underlying the pathogenesis of osteoporosis by using bioinformatics data analyses and experimental validation. Our findings would be a potential and promising therapeutic target for precision medicine in osteoporosis. © 2018 American Society for Bone and Mineral Research.

Association between fibroblast growth factor 21 and bone mineral density in adults.

PURPOSE: Animal-based studies have reported a decrease in bone mass resulting from high level of fibroblast growth factor 21 (FGF21). However, the correlation between plasma FGF21 levels and bone mineral density (BMD) is paradoxical in previous human-based studies, and the associations between FGF21 gene polymorphisms and BMD haven't been reported yet. Therefore, here, we evaluated plasma FGF21 levels with sufficient study samples, and performed genetic association test to reveal the physiological and genetic role of FGF21 on BMD in adults. METHODS: Plasma and genetic samples containing 168 and 569 Han Chinese subjects, respectively, were employed in this study. Fasting plasma FGF21 levels were determined using enzyme-linked immunosorbent assay (ELISA). Regional BMD values were measured by dual energy X-ray absorptiometry (DXA). Five variants of FGF21 gene were successfully genotyped. RESULTS: Physiological association suggested that plasma FGF21 levels were inversely correlated with BMD in femoral neck (Neck-BMD: P = 0.039) and Ward's triangle (Ward's-BMD: P = 0.002) of hip region. A FGF21 gene variant, rs490942, was significantly associated with the increase of Ward's-BMD in total (P = 0.027) and female (P = 0.016) cohorts, as well as Neck-BMD in female cohort (P = 7.45 × 10-3). Meanwhile, eQTL results indicated that this SNP was related to the decreased level of FGF21 gene expression. CONCLUSIONS: Taking together from both physiological and genetic levels, we suggest that FGF21 is inversely associated with regional BMD. And we haven't observed sex-specific effect in this study.

A functional SNP regulated by miR-196a-3p in the 3'UTR of FGF2 is associated with bone mineral density in the Chinese population.

Previous studies have identified FGF2 as a susceptibility gene for osteoporosis in Caucasians. Evaluating the genetic associations in different ethnicities is necessary. Moreover, elucidating the functional mechanism for the susceptibility loci is important to offer new targets for therapeutic studies. Here, we genotyped 10 SNPs in FGF2 and tested for associations with bone mineral density (BMD) in a discovery sample of 1,300 Chinese subjects. Nominally significant results were subjected to replication in another sample of 1,039 Chinese subjects. We identified one SNP rs1048201:C>T in FGF2 3'untranslated region significantly associated with spine BMD (combined cohorts, P = 1.53×10-3 ). Expression quantitative trait locus analyses revealed that rs1048201 also affected FGF2 gene expression (P = 7.03×10-4 ). Bioinformatics prediction suggested that rs1048201 T allele could disrupt miRNA binding. Luciferase assay validated that the C allele had a repressive effect on FGF2 gene expression. We found that hsa-miR-196a-3p affected expression on both mRNA and protein levels of FGF2. In conclusion, our study provided evidence that a functional SNP rs1048201 was associated with BMD, and SNP rs1048201:C>T variant may act by affecting binding of hsa-miR-196a-3p. The SNP-modified posttranscriptional gene regulation by miRNA could be a potentially pathogenetic mechanism of osteoporosis.

Associations of Plasma FGF2 Levels and Polymorphisms in the FGF2 Gene with Obesity Phenotypes in Han Chinese Population.

Obesity is highly heritable, but the specific genes influencing obesity related traits are largely unknown. Fibroblast growth factor 2 (FGF2) could influence adipocyte differentiation. However, the association of FGF2 polymorphisms and obesity remains unclear. This study aimed to investigate the associations of both the plasma FGF2 levels and SNPs in FGF2 gene with obesity phenotypes in Han Chinese populations. Plasma FGF2 levels were measured and subjected to association analyses in 62 subjects. Eleven SNPs in FGF2 were genotyped and tested for associations in a discovery sample of 1,300 subjects. SNPs significantly associated with obesity were subjected to replication in another independent sample of 1,035 subjects. We found that plasma FGF2 levels were positively correlated with fat mass (P = 0.010). Association analyses in the discovery sample identified three SNPs (rs1449683, rs167428, rs308442) significantly associated with fat mass after multiple testing adjustments (P < 0.0045). Subsequent replication study successfully validated one SNP (rs167428) associated with fat mass (P(combine) = 3.46 × 10(-5)). eQTL analyses revealed that SNPs associated with obesity also affected FGF2 expression. Our findings suggested that high plasma FGF2 level correlated with increased risk of obesity, and FGF2 gene polymorphisms could affect individual variances of obesity in Han Chinese population.