Long-Fei Wu

Functional relevance for multiple sclerosis-associated genetic variants

Multiple sclerosis (MS) is an inflammatory and demyelinating disease of central nervous system. Many genetic variants associated with MS have been identified by genome-wide association studies, but functional mechanism underlying the associations is largely unclear. Utilizing the publically available datasets, we carried out gene relationships among implicated loci (GRAIL) analyses to search for MS-associated SNPs/genes. Expression quantitative trait loci (eQTLs) analyses were conducted to identify eQTL SNPs/target genes. Further, functional prediction for SNP, differential gene expression, and functional annotation clustering analyses for gene were conducted to explore their functional relevance to MS. Among the 284 identified MS-associated SNPs (P < 10−4), eQTL analysis showed that 45 SNPs act as cis-effect regulators on 19 MS-associated genes. Among the 19 eQTL target genes, 14 showed significantly differential expressions in MS-related cells. Among the 45 SNPs, 15 were predicted most likely located in transcription factor (TF) binding sites, and five predicted SNPs (rs3095329 of TUBB, rs9469220/rs2647046 of HLA-DQB1, rs11154801 of AHI1, and rs1062158 of NDFIP1) have corresponding target genes with significantly differential expressions in multiple cell groups, while rs7194 of HLA-DRA was predicted in the has-miR-6507-3p binding site. The functional evidence, taken together, highlighted the functional relevance of the six SNPs to MS. The present findings provide novel insights into the functional mechanisms underlying the MS-associated genetic variants, which improve our understanding of the genetic association for MS.

Integrated Analyses of Gene Expression Profiles Digs out Common Markers for Rheumatic Diseases

Objective Rheumatic diseases have some common symptoms. Extensive gene expression studies, accumulated thus far, have successfully identified signature molecules for each rheumatic disease, individually. However, whether there exist shared factors across rheumatic diseases has yet to be tested. Methods We collected and utilized 6 public microarray datasets covering 4 types of representative rheumatic diseases including rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis, and osteoarthritis. Then we detected overlaps of differentially expressed genes across datasets and performed a meta-analysis aiming at identifying common differentially expressed genes that discriminate between pathological cases and normal controls. To further gain insights into the functions of the identified common differentially expressed genes, we conducted gene ontology enrichment analysis and protein-protein interaction analysis. Results We identified a total of eight differentially expressed genes (TNFSF10, CX3CR1, LY96, TLR5, TXN, TIA1, PRKCH, PRF1), each associated with at least 3 of the 4 studied rheumatic diseases. Meta-analysis warranted the significance of the eight genes and highlighted the general significance of four genes (CX3CR1, LY96, TLR5, and PRF1). Protein-protein interaction and gene ontology enrichment analyses indicated that the eight genes interact with each other to exert functions related to immune response and immune regulation. Conclusion The findings support that there exist common factors underlying rheumatic diseases. For rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and osteoarthritis diseases, those common factors include TNFSF10, CX3CR1, LY96, TLR5, TXN, TIA1, PRKCH, and PRF1. In-depth studies on these common factors may provide keys to understanding the pathogenesis and developing intervention strategies for rheumatic diseases.

Identification and evaluation of lncRNA and mRNA integrative modules in human peripheral blood mononuclear cells.

AIM To identify sets of functionally related long noncoding RNAs (lncRNAs) and mRNAs and to evaluate the importance of lncRNAs in an lncRNA-mRNA network. METHODS We carried out weighted gene co-expression network analysis and enrichment analyses to identify functional modules of co-expressed lncRNAs and mRNAs in human peripheral blood mononuclear cells of 43 females. RESULTS We identified seven modules and found hub lncRNAs in each module. Four of the seven modules had significant gene ontology enrichments. Some of the hub lncRNAs (e.g., SSX8, UCA1, HOXA-AS2, STARD4-AS1 and PCBP1-AS1) have known functions related with diseases such as cancers. CONCLUSION We identified seven biologically important lncRNA and mRNA integrative modules in females and showed that lncRNAs might play important roles in lncRNA-mRNA co-expression modules.

Relative abundance of mature myostatin rather than total myostatin is negatively associated with bone mineral density in Chinese

Myostatin is mainly secreted by skeletal muscle and negatively regulates skeletal muscle growth. However, the roles of myostatin on bone metabolism are still largely unknown. Here, we recruited two large populations containing 6308 elderly Chinese and conducted comprehensive statistical analyses to evaluate the associations among lean body mass (LBM), plasma myostatin, and bone mineral density (BMD). Our data revealed that total myostatin in plasma was mainly determined by LBM. The relative abundance of mature myostatin (mature/total) was significantly lower in high versus low BMD subjects. Moreover, the relative abundance of mature myostatin was positively correlated with bone resorption marker. Finally, we carried out in vitro experiments and found that myostatin has inhibitory effects on the proliferation and differentiation of human osteoprogenitor cells. Taken together, our results have demonstrated that the relative abundance of mature myostatin in plasma is negatively associated with BMD, and the underlying functional mechanism for the association is most likely through inhibiting osteoblastogenesis and promoting osteoclastogenesis.

Functional mechanisms for type 2 diabetes-associated genetic variants

AIMS Type 2 diabetes (T2D) is a complex endocrine and metabolic disorder, characterized by hyperglycemia due to insulin resistance and relative lack of insulin. Several recent studies have identified a large number of genetic loci associated with T2D without exploring functional mechanisms underlying the associations. This study established integrative analyses to detect the functional mechanisms for T2D-related associations. METHODS Based on the public available datasets and resources, this study performed integrative analyses (gene relationships among implicated loci (GRAIL), expression quantitative trait loci (eQTL) analysis, differential gene expression analysis and functional prediction analysis) to detect the molecular functional mechanisms underlying the associations. RESULTS Two single nucleotide polymorphisms (SNPs) (rs7593730, rs2439312) have been found to act as cis-effect regulators of two corresponding eQTL genes (RBMS1, NRG1) among 252 selected (P<E-4) genetic associations that were archived in the public databases. These two non-HLA genes were also differentially expressed in T2D-related cell groups. The two SNPs were predicted as regulatory sites by utilizing online prediction tools. CONCLUSIONS This study detected potential regulatory mechanisms underlying the associations between T2D and two identified SNPs. Integrative analysis can be used to provide suggestive clues for the molecular functional mechanisms in T2D.

Integrative analysis identified mediation effects of lncRNAs on the correlations between methylation and mRNA

The aim of this study was to construct DNA methylation-lncRNA-mRNA interaction trios in peripheral blood mononuclear cells. We first conducted eQTL analyses using genome-wide methylation, lncRNA and mRNA expression data from 43 Chinese females. Next, causal inference test (CIT) was used to detect the lncRNA mediation effects on methylation and mRNA. Methylation-lncRNA cis-eQTL analysis identified 11 significant cis-methylation-lncRNA pairs. Combined with the results from the next lncRNA-mRNA eQTL and methylation-mRNA eQTL analyses, the 11 significant pairs and their corresponding 11,204 target e-mRNAs formed 12,245 trios. Further CIT identified six lncRNAs as mediators in regulating the corresponding pairs between methylation and mRNA. This study detected lncRNAs with mediation effects on the correlations between DNA methylations and a large number of mRNAs.

Anxa2 attenuates osteoblast growth and is associated with hip BMD and osteoporotic fracture in Chinese elderly

Low bone mineral density (BMD) is a risk factor of osteoporotic fracture (OF). Peripheral blood monocytes (PBM) can differentiate into osteoclasts to resorb bone. It was known that PBM-expressed Anxa2 protein is associated with BMD, and extracellular Anxa2 protein promotes osteoclastogenesis. This study aimed to test 1) whether Anxa2 protein level in PBM differs significantly between subjects with OF and without fracture history (NF); 2) whether Anxa2 level in plasma is associated with BMD; 3) how Anxa2 protein at various concentrations would affect osteoblastic activity in vitro. All the study subjects were Chinese Han elderly. Firstly, Anxa2 protein in PBM was identified and quantitated by LC-MS/MS and compared between 45 OF cases and 42 healthy controls. Secondly, plasma Anxa2 protein level was quantitated by ELISA and compared between unrelated subjects with extremely low vs. high hip BMD (0.63±0.10 vs. 1.05±0.10 g/cm2, n = 75). Furthermore, in vitro functional assay was utilized to test the effects of extracellular Anxa2 protein on osteoblastic growth. We found that Anxa2 protein expression in PBM was significantly up-regulated in OF vs. NF subjects (fold change [FC)] = 1.16, P<0.05). Plasma Anxa2 protein concentration (range: 31.69–227.35ng/ml) was significantly elevated in low vs. high BMD subjects (84.85 vs. 66.15ng/ml, FC = 1.28, P<0.05). Cellular dynamical monitoring demonstrated that the general shape of dose-response relationship is the inverse U-shaped curve. Specifically, lower dose of Anxa2 protein may promote osteoblast growth and the optimal concentration for osteoblastic growth was around 50ng/ml, but even higher concentration could attenuate hFOB1.19 osteoprogenitor cell growth. We concluded that Anxa2 protein could attenuate osteoblast growth and be associated with hip BMD and OF in Chinese elderly.

Plasma gelsolin is associated with hip BMD in Chinese postmenopausal women

Gelsolin (GSN) protein, expressed in circulating monocytes, was previously reported to be associated with osteoporosis in both Chinese and Caucasian women. This study aims to test if plasma GSN protein level is associated with hip bone mineral density (BMD) in Chinese population. Based on two study Groups containing 6,308 old Chinese, we adopted extreme sampling scheme and selected 3 independent samples (Subgroups 1–3) for discovery, replication, and validation purposes. We tested plasma GSN concentration, and analyzed whether plasma GSN level differs between subjects with extremely low vs. high hip BMD. In Group 1 (N = 1,860), the plasma GSN level increased in the female with low BMD, which was discovered in the Subgroup 1 (N = 42, p = 0.093) and replicated in the Subgroup 2 (N = 39, p = 0.095). With more extreme sampling for the Subgroup 3 from the Group 2 (N = 4,448), the difference of plasma GSN level in the female with low BMD vs. high BMD is more significant (N = 45, p = 0.037). After the subjects were pooled from Subgroups 2 and 3, the difference in plasma GSN between low and high BMD subjects became even more significant (p = 0.016). The plasma GSN level was negatively correlated with total hip BMD (r = -0.26, p = 0.033). We concluded that plasma GSN was associated with hip BMD in Chinese postmenopausal women and plasma GSN might be a potential risk biomarker for osteoporosis.

Detection of lncRNA-mRNA interaction modules by integrating eQTL with weighted gene co-expression network analysis

One major function of lncRNA is to regulate the expression of mRNA, but the patterns of their interactions were largely unknown. We attempted to construct lncRNA-mRNA interaction modules at a genome-wide scale. We performed a genome-wide lncRNA-mRNA eQTL analysis in peripheral blood mononuclear cells of 43 individuals, followed by weighted gene co-expression network analysis and functional enrichment analysis which sought to detect functional modules. There were 4627 significant cis lnc-eQTL pairs (P < 1.4 × 10−6) and 1,587,128 significant trans lnc-eQTL pairs (P < 3.46 × 10−9). We detected 11 eQTL modules for the lnc-eQTL networks. Among them, five modules showed significant enrichments in GO terms, and three modules showed significant enrichments in specific KEGG pathways (e.g., Toll-like receptor, PI3K-Akt, NF-kappa B, and TNF signaling pathways). lncRNA-protein interaction analysis showed that some well-known functional lncRNAs (HOTAIR, CCDC26, RHPN1-AS1, WT1-AS, and TCL6) in the eQTL module interacted with genes in focal adhesion and PI3K-Akt signaling pathway. We identified biologically functional lncRNA-mRNA interaction modules by integrating eQTL and weighted gene co-expression network analysis. Integrative analysis of lncRNA and mRNA data by applying eQTL analysis and weighted gene co-expression network analysis methods could be helpful for functional annotation of lncRNAs.

Association of Plasma Irisin with Bone Mineral Density in a Large Chinese Population Using an Extreme Sampling Design

Irisin, a myokine produced by skeletal muscle in response to physical exercise, promotes trans-differentiation of white adipose tissue into brown adipose tissue. Recent evidences suggested that irisin also plays an important role in the control of bone metabolism. This study aimed to ascertain the relationship between plasma irisin and bone mineral density (BMD) in Chinese population by adoption of an extreme sampling method. Based on a large and screened Chinese elderly population (N = 6308), two subgroups with extremely high and low hip BMD were selected for discovery (N = 80, high vs. low BMD = 44:36) and validation (N = 60, high vs. low BMD = 30:30), respectively. Plasma irisin, P1NP, and β-CTx were measured using commercially available ELISA kits. Other metabolic parameters (e.g., blood glucose, total cholesterol and triglycerides) were collected. Student’s t test and Spearman correlation analyses were conducted in SPSS. Significant difference was discovered for plasma irisin between females and age-matched males (N = 80, male vs. female = 42:38, P = 0.002). The plasma irisin levels were significantly higher in high BMD subjects than in low BMD subjects, which was observed in both discovery (P = 0.012) and validation samples (P = 0.022). However, such observation was limited to males only. Further correlation analyses in males showed that plasma irisin was correlated with BMD (r = 0.362, P = 0.025) and triglyceride (r = − 0.354, P = 0.032). Plasma irisin levels were associated with hip BMD in Chinese elderly men. This study represented the first effort of investigating the relationship of plasma irisin and BMD in elderly population. The positive correlation between plasma irisin and BMD hints intrinsic communication between muscle and bone.