Unraveling the relationship between high-sensitivity C-reactive protein and frailty: evidence from longitudinal cohort study and genetic analysis.

BACKGROUND: This study aimed to investigate the association of high-sensitivity C-reactive protein (hs-CRP) with incident frailty as well as its effects on pre-frailty progression and regression among middle-aged and older adults. METHODS: Based on the frailty index (FI) calculated with 41 items, 6890 eligible participants without frailty at baseline from China Health and Retirement Longitudinal Study (CHARLS) were categorized into health, pre-frailty, and frailty groups. Logistic regression models were used to estimate the longitudinal association between baseline hs-CRP and incident frailty. Furthermore, a series of genetic approaches were conducted to confirm the causal relationship between CRP and frailty, including Linkage disequilibrium score regression (LDSC), pleiotropic analysis, and Mendelian randomization (MR). Finally, we evaluated the association of hs-CRP with pre-frailty progression and regression. RESULTS: The risk of developing frailty was 1.18 times (95% CI: 1.03-1.34) higher in participants with high levels of hs-CRP at baseline than low levels of hs-CRP participants during the 3-year follow-up. MR analysis suggested that genetically determined hs-CRP was potentially positively associated with the risk of frailty (OR: 1.06, 95% CI: 1.03-1.08). Among 5241 participants with pre-frailty at baseline, we found pre-frailty participants with high levels of hs-CRP exhibit increased odds of progression to frailty (OR: 1.39, 95% CI: 1.09-1.79) and decreased odds of regression to health (OR: 0.84, 95% CI: 0.72-0.98) when compared with participants with low levels of hs-CRP. CONCLUSIONS: Our results suggest that reducing systemic inflammation is significant for developing strategies for frailty prevention and pre-frailty reversion in the middle-aged and elderly population.

PKM2 is a Novel Osteoporosis-Associated Protein in Chinese.

Purpose:Osteoporosis is characterized by low bone mineral density (BMD) and high risk of osteoporotic fracture (OF). Peripheral blood monocytes (PBM) can differentiate into osteoclasts to resorb bone. This study was to identify PBM-expressed proteins significant for osteoporosis in Chinese Han elderly population (>65 years), and focused on two phenotypes of osteoporosis: low BMD and OF. METHODS: Label-free quantitative proteomics was employed to profile PBM proteome and to identify differentially expressed proteins (DEPs) between OF (N=27) vs. non-fractured (NF, N=24) subjects and between low BMD (N=12) vs. high BMD (N=12) subjects in women. Western blotting (WB) was conducted to validate differential expression, and ELISA to evaluate translational value for secretory protein of interest. RESULTS: We discovered 59 DEPs with fold change (FC)>1.3 (P<1×10-5), and validated the significant up-regulation of pyruvate kinase isozyme 2 (PKM2) with osteoporosis (P<0.001). PKM2 protein upregulation with OF was replicated with PBM in men (P=0.04). Plasma PKM2 protein level was significantly elevated with OF in an independent sample (N=100, FC=1.68, P=0.01). Pursuant functional assays showed that extracellular PKM2 protein supplement not only promoted monocyte trans-endothelial migration, growth, and osteoclast differentiation (marker gene expression), but also inhibited osteoblast growth, differentiation (ALP gene expression), and activity. CONCLUSION: The above findings suggest that PKM2 protein is a novel osteoporosis-associated functional protein in Chinese Han elderly population. It may serve as a risk biomarker and drug target for osteoporosis.

Rheumatoid arthritis and osteoporosis: shared genetic effect, pleiotropy and causality.

Rheumatoid arthritis (RA) is associated with increased localized and generalized bone loss, but the complex genetic mechanism between them is still unknown. By leveraging large-scale genome-wide association studies summary statistics and individual-level datasets (i.e. UK Biobank), a series of genetic approaches were conducted. Linkage disequilibrium score regression reveals a shared genetic correlation between RA and estimated bone mineral density (eBMD) (rg = -0.059, P = 0.005). The PLACO analysis has identified 74 lead (8 novel) pleiotropic loci that could be mapped to 99 genes, the genetic functions of which reveal the possible mechanism underlying RA and osteoporosis. In European, genetic risk score (GRS) and comprehensive Mendelian randomization (MR) were utilized to evaluate the causal association between RA and osteoporosis in European and Asian. The increase in GRS of RA could lead to a decrease of eBMD (beta = -0.008, P = 3.77E-6) and a higher risk of facture [odds ratio (OR) = 1.012, P = 0.044]. MR analysis identified that genetically determined RA was causally associated with eBMD (beta = -0.021, P = 4.14E-05) and fracture risk (OR = 1.036, P = 0.004). Similar results were also observed in Asian that osteoporosis risk could be causally increased by RA (OR = 1.130, P = 1.04E-03) as well as antibodies against citrullinated proteins-positive RA (OR = 1.083, P = 0.015). Overall, our study reveals complex genetic mechanism between RA and osteoporosis and provides strong evidence for crucial role of RA in pathogenesis of osteoporosis.

A novel long non-coding RNA, lnc-RNU12, influences the T-cell cycle via c-JUN and CCNL2 in rheumatoid arthritis.

OBJECTIVES: Long non-coding RNAs (lncRNAs) play important roles in RA pathogenesis. However, specific lncRNAs that regulate gene expression in RA pathogenesis are poorly known. This study was undertaken to characterize a novel lncRNA (lnc-RNU12) that has a lower-than-normal expression level in RA patients. METHODS: We performed initial genome-wide lncRNA microarray screening in peripheral blood mononuclear cells from 28 RA cases and 18 controls. Multiple methods were used to validate the detected associations between lncRNAs and RA. Furthermore, we identified the source and characteristics of the highlighted lncRNAs, detected the target genes, and determined the functional effect on immune cells through lncRNA knock-down in Jurkat T cell lines. RESULTS: lnc-RNU12 was downregulated in peripheral blood mononuclear cells and T cell subtypes of RA patients and was genetically associated with RA risk. lnc-RNU12 mediates the effect of microbiome alterations on RA risk. Activation of T cells caused low expression of lnc-RNU12. Knock-down of lnc-RNU12 in Jurkat T cells caused cell cycle S-phase arrest and altered the expression of protein-coding genes related to the cell cycle and apoptosis (e.g. c-JUN, CCNL2, CDK6, MYC, RNF40, PKM, VPS35, DNAJB6 and FLCN). Finally, c-JUN and CCNL2 were identified as target genes of lnc-RNU12 at the mRNA and protein expression levels. RNA-binding protein immunoprecipitation assays verified the interaction between lnc-RNU12 and the two proteins (c-Jun and cyclin L2) in Jurkat cells. CONCLUSIONS: Our study suggested that lnc-RNU12 was involved in the pathogenesis of RA by influencing the T cell cycle by targeting c-JUN and CCNL2.

Establishment of reference intervals for bone turnover markers in healthy Chinese older adults.

BACKGROUND: Reference ranges for bone turnover markers (BTMs) are still lacking in the healthy Chinese population. AIM: To establish reference intervals for BTMs and to investigate the correlations between BTMs and bone mineral density (BMD) in Chinese older adults. SUBJECTS AND METHODS: A community-based cross-sectional study was conducted among 2511 Chinese subjects aged over 50 yrs residing in Zhenjiang, Southeast China. Reference intervals for BTMs (i.e. procollagen type I N-terminal propeptide, P1NP; ß cross-linked C-terminal telopeptide of type I collagen, ß-CTX) were calculated as the central 95% range of all measurements in Chinese older adults. RESULTS: The reference intervals of P1NP, ß-CTX and P1NP/ß-CTX were 15.8-119.9 ng/mL, 0.041-0.675 ng/mL and 49.9-1261.5 for females and 13.6-111.4 ng/mL, 0.038-0.627 ng/mL and 41.0-1269.1 for males, respectively. In the multiple linear regression analysis, only ß-CTX was negatively associated with BMD after adjusting for age and body mass index (BMI) in both sex-stratified groups (all p < .05). CONCLUSION: This study established age- and sex-specific reference intervals for BTMs in a large sample of healthy Chinese participants ≥ 50 and < 80 years of age and explored the correlations between BTMs and BMD, which provides an effective reference for the assessment of bone turnover in the clinical practice of osteoporosis.

Why SNP rs3755955 is associated with human bone mineral density? A molecular and cellular study in bone cells.

SNP rs3755955 (major/minor allele: G/A) located in Iduronidase-Alpha-L- (IDUA) gene was reported to be significant for human bone mineral density (BMD). This follow-up study was to uncover the underlying association mechanism through molecular and cellular functional assays relevant to bone. We tested the effects of single nucleotide polymorphisms (SNP) rs3755955 (defined allele G as wild-type and allele A as variant-type) on osteoblastic and osteoclastic functions, as well as protein phosphorylation in stably transfected human fetal osteoblast (hFOB) cell and mononuclear-macrophage (RAW264.7) cell. In hFOB cells, transfection with variant-type IDUA significantly decreased osteoblastic gene expression (OPN, COL1A1 and RANKL) (p < 0.01), impeded cell proliferation (p < 0.05), stimulated cell apoptosis (p < 0.001) and decreased ALP enzyme activity, as compared with that of wild-type IDUA transfection. In RAW264.7 cells, transfection with variant-type IDUA significantly inhibited cell apoptosis (p < 0.01), promoted osteoclastic precursor cell migration (p < 0.0001), growth (p < 0.01), osteoclastic gene expression (TRAP, RANK, Inte-αv and Cath-K) (p < 0.05) and TRAP enzyme activity (p < 0.001), as compared with that of wild-type IDUA transfection. In both hFOB and RAW264.7 cells, the total protein and IDUA protein-specific phosphorylation levels were significantly reduced by variant IDUA transfection, as compared with that of wild-type IDUA transfection (p < 0.05). Variant allele A of phosSNP rs3755955 in IDUA gene regulates protein phosphorylation, inhibits osteoblast function and promotes osteoclastic activity. The SNP rs3755955 could alter IDUA protein phosphorylation, significantly regulates human osteoblastic and osteoclastic gene expression, and influences the growth, differentiation and activity of osteoblast and osteoclast, hence to affect BMD.

Body Surface Area (BSA) is a Better Osteoporosis Associated Anthropometric Parameter Than Other Anthropometric Parameters in Elderly Population.

Body surface area (BSA) is widely used for adjusting drug dose, while few studies have yet systematically evaluated its association with osteoporosis and compared its advantage with other anthropometric parameters in osteoporotic risk prediction. A total of 10,021 Chinese individuals aged over 65 years were enrolled in our study. Bone mineral density (BMD) was measured, and demographic information was also collected. Pearson correlation analysis, receiver operating characteristic (ROC) curves and predictive analysis were performed to assess the clinical practice of BSA for osteoporosis. BSA had the strongest correlation with BMD (0.544, p < 0.001) compared with conventional anthropometric indices. Besides, BSA had the highest power in osteoporosis prediction, with an area under the curve (AUC) reaching 0.81. After incorporating BSA into the osteoporosis risk prediction model, the AUC improved from 0.82 to 0.83 (p < 0.01). We found BSA provided additional diagnostic value beyond conventional anthropometric information with continuous and category NRIs were 30.40% (p < 0.01) and 3.29% (p < 0.01), respectively, and the IDI was 1.85% (p < 0.01). BSA was positively associated with osteoporosis and showed superior discriminative ability for osteoporosis risk prediction compared with other anthropometric parameters in the Chinese elderly population.

Does obesity mediate the relationship between diabetes and osteoporosis in Chinese elderly population?

OBJECTIVES: Diabetes mellitus (DM), osteoporosis (OP), and obesity (OB) are three complex diseases. OB is associated with both DM and OP, but it is unclear whether OB mediates association between DM and OP. The study aimed to investigate the potential mediation effects of OB on association between DM and bone mineral density (BMD) by the causal inference tests (CIT). METHODS: A total of 5682 Chinese aged over 65 years were enrolled in an ongoing cohort: Osteoporosis Preventive Project (OPP). Obesity-related indexes, including body mass index (BMI), waist circumference, and waist circumference-hip circumference-ratio (WHR), and BMD at total hip (TH) and femur neck (FN) were measured. RESULTS: Subjects with DM had significant greater values of age, weight, BMI, waist circumference, WHR, and BMD than non-DM subjects. BMD at TH and FN was significantly associated with DM (p < 0.05) with adjustment of age both in males and females. Further CIT showed that OB-related indexes (BMI, waist circumference, and WHR) are significantly mediators in the associations between DM and BMD in females, but not in males. Furthermore, the mediation effects of waist circumference were detected on DM and TH BMD in the females of normal-weight group. CONCLUSIONS: Obesity-related indexes, especially waist circumference, serve as significant mediator(s) between DM and OP in Chinese female elderly. Diabetes increases BMD by increasing obesity-related indexes. The findings established the intermediate role of OB underlying the association between DM and OP in human population.

PhosSNPs-Regulated Gene Network and Pathway Significant for Rheumatoid Arthritis.

OBJECTIVES: Peripheral blood mononuclear cells (PBMCs) are critical for immunity and participate in multiple human diseases, including rheumatoid arthritis (RA). PhosSNPs are nonsynonymous SNPs influencing protein phosphorylation, thus probably modulate cell signaling and gene expression. We aimed to identify phosSNPs-regulated gene network/pathway potentially significant for RA. METHODS: We collected genome-wide phosSNP genotyping data and transcriptome-wide mRNA expression data from PBMCs of a Chinese sample. We discovered and verified with public datasets differentially expressed genes (DEGs) associated with RA, and replicated RA-associated SNPs in our study sample. We performed a targeted expression quantitative trait locus (eQTL) study on significant phosSNPs and DEGs. RESULTS: We identified 29 nominally significant eQTL phosSNPs and 83 target genes, and constructed comprehensive regulatory/interaction networks, highlighting the vital effects of two eQTL phosSNPs (rs371513 and rs4824675, FDR <0.05) and four critical node genes (HSPA4, NDUFA2, MRPL15, and ATP5O). Besides, two node/key genes NDUFA2 and ATP5O, regulated by rs371513, were significantly enriched in mitochondrial oxidative phosphorylation pathway. Besides, four pairs of eQTL effects were replicated independently in whole blood and/or transformed fibroblasts. CONCLUSIONS: The findings delineated a potential role of protein phosphorylation and genetic variations in RA and warranted the significant roles of phosSNPs in regulating RA-associated genes expression in PBMCs. The results pointed out the relevance and significance of oxidative phosphorylation pathway to RA.

Identification of Potential Pleiotropic Genes for Immune and Skeletal Diseases Using Multivariate MetaCCA Analysis.

Background: Immune and skeletal systems physiologically and pathologically interact with each other. Immune and skeletal diseases may share potential pleiotropic genetics factors, but the shared specific genes are largely unknown. Objective: This study aimed to investigate the overlapping genetic factors between multiple diseases (including rheumatoid arthritis (RA), psoriasis, osteoporosis, osteoarthritis, sarcopenia, and fracture). Methods: The canonical correlation analysis (metaCCA) approach was used to identify the shared genes for six diseases by integrating genome-wide association study (GWAS)-derived summary statistics. The versatile Gene-based Association Study (VEGAS2) method was further applied to refine and validate the putative pleiotropic genes identified by metaCCA. Results: About 157 (p<8.19E-6), 319 (p<3.90E-6), and 77 (p<9.72E-6) potential pleiotropic genes were identified shared by two immune diseases, four skeletal diseases, and all of the six diseases, respectively. The top three significant putative pleiotropic genes shared by both immune and skeletal diseases, including HLA-B, TSBP1, and TSBP1-AS1 (p

The different correlations between obesity and osteoporosis after adjustment of static mechanical loading from weight and fat free mass.

OBJECTIVES: To explore complex correlations between obesity (OB) and osteoporosis (OP) after adjustment of static mechanical loading from weight and fat free mass (FFM). METHODS: A total of 3749 Chinese aged ≥65 years were selected from our ongoing cohort study. OB indices and bone mineral density (BMD) were measured for each subject. Linear regression analyses were performed to explore the correlations between OB indices and OP under three adjustment models (unadjusted, adjusted with weight and adjusted with FFM). RESULTS: Under no adjustment, three general obesity indices (body mass index: BMI, fat mass: FM, and percentage FM: PFM) were positively associated with BMD at three skeletal sites (P<0.001) in the regression analyses. However, after the adjustment with weight, these associations were mostly significant but reverse i.e., negatively in direction. After adjustment with FFM, the three indices were still positively and significantly (P<0.001) associated with BMD but regression coefficients were smaller compared to the unadjusted associations. Similar associations were observed for central adiposity and lower limb adiposity indices. CONCLUSIONS: The combined relation of OB to OP due to the physiological factors secreted from adipose tissues and the static mechanical loading from FM is positive in direction.

Global Public Interests and Dynamic Trends in Osteoporosis From 2004 to 2019: Infodemiology Study.

BACKGROUND: With the prolonging of human life expectancy and subsequent population aging, osteoporosis (OP) has become an important public health issue. OBJECTIVE: This study aimed to understand the global public search interests and dynamic trends in "osteoporosis" using the data derived from Google Trends. METHODS: An online search was performed using the term "osteoporosis" in Google Trends from January 1, 2004, to December 31, 2019, under the category "Health." Cosinor analysis was used to test the seasonality of relative search volume (RSV) for "osteoporosis." An analysis was conducted to investigate the public search topic rising in RSV for "osteoporosis." RESULTS: There was a descending trend of global RSV for "osteoporosis" from January 2004 to December 2014, and a slowly increasing trend from January 2015 to December 2019. Cosinor analysis showed significant seasonal variations in global RSV for "osteoporosis" (P=.01), with a peak in March and a trough in September. In addition, similar decreasing trends of RSV for "osteoporosis" were found in Australia, New Zealand, Ireland, and Canada from January 2004 to December 2019. Cosinor test revealed significant seasonal variations in RSV for "osteoporosis" in Australia, New Zealand, Canada, Ireland, UK, and USA (all P<.001). Furthermore, public search rising topics related to "osteoporosis" included denosumab, fracture risk assessment tool, bone density, osteopenia, osteoarthritis, and risk factor. CONCLUSIONS: Our study provided evidence about the public search interest and dynamic trends in OP using web-based data, which would be helpful for public health and policy making.

Potential Biomolecules Capable of Assessing Risk of Osteoporotic Fracture: A Scoping Review.

Osteoporosis as a systemic chronic skeletal disease is characterized by low bone mineral density and increased risk to osteoporotic fractures. Osteoporosis is prevalent in the middle-aged and elderly population, especially in the postmenopausal women. With population aging, osteoporosis has become a world-wide serious public health problem. Early recognition of the high-risk population followed by timely and efficient intervention and/or treatment is important for preventing osteoporotic fractures. In light of the high heritability and complex pathogenesis of osteoporosis, comprehensive consideration of vital biological/biochemical factors is necessary for accurate risk evaluation of fractures. For this purpose, we review recent research progress on molecules which can be applied to assess risk for osteoporotic fractures. Future integrative analyses and systematic evaluation of these molecules may facilitate developing novel methodologies and/or test strategies, i.e., biochips, for early recognition of osteoporosis, hence contributing to preventing osteoporotic fractures.

Protein array test detected three osteoporosis related plasma inflammatory cytokines in Chinese postmenopausal women.

Inflammatory cytokines were involved in pathological conditions of osteoporosis (OP). However, the specific OP-associated inflammatory cytokines are still awaiting to be detected by using a systemic method. Herein, we adopted an extreme sampling scheme and examined inflammatory cytokines between subjects with low and high bone mineral density (BMD) through protein microarray. First, 8 candidate cytokines including B lymphocyte chemoattractant (BLC), osteopontin (OPN) and insulin-like growth factor-binding protein 4 (IGFBP4) were identified in the discovery extreme sampling subgroup. Then, the different expressions for BLC, OPN and IGFBP4 were validated and replicated in two independent extreme sampling subgroups. Further functional experiments showed that the cytokine BLC was involved in bone metabolism by inhibiting bone formation and promoting bone resorption. Together, this study further revealed that inflammatory cytokines were closely related with OP, and that they highlighted critical roles of BLC in the pathogenesis of OP.

Global correlation analysis for miRNA and protein expression profiles in human peripheral blood mononuclear cells.

Micro-RNAs (miRNAs) are small noncoding RNAs that negatively regulate gene expression at protein level by protein translation inhibition or mRNA degradation. However, the global correlation patterns between miRNA and protein have not been studied yet. To establish the global correlation patterns in human peripheral blood mononuclear cells (PBMCs), this study conducted multiple types of miRNA-protein correlation analyses in 28 Chinese subjects. Pearson correlation analysis showed a negative but relatively small global correlation in each subject. Among the 371 constructed miRNA-protein pairs (60 unique miRNAs, and 150 unique proteins), 10.5% of pairs have significant correlations (P < 0.05). Some highlighted miRNAs (e.g., hsa-miR-590-3p, hsa-miR-520d-3p) exerted significant regulation on multiple genes. Simultaneously, some genes (e.g., HSP90B1) were targeted by multiple miRNAs. The target genes associated with miRNAs tend to enrich in some important GO terms: biological processes (e.g., gene expression, protein binding and RNA binding), and molecular functions (protein binding: GO:0005515; RNA binding: GO:0003723). The results provided a global view of the miRNA-protein expression correlation profile in human PBMCs, which would facilitate in-depth investigation of biological functions of key miRNAs/proteins and better understanding of the pathogenesis underlying PBMC related diseases.

Why SNP rs227584 is associated with human BMD and fracture risk? A molecular and cellular study in bone cells.

A large number of SNPs significant for osteoporosis (OP) had been identified by genome-wide association studies. However, the underlying association mechanisms were largely unknown. From the perspective of protein phosphorylation, gene expression regulation, and bone cell activity, this study aims to illustrate association mechanisms for representative SNPs of interest. We utilized public databases and bioinformatics tool to identify OP-associated SNPs which potentially influence protein phosphorylation (phosSNPs). Associations with hip/spine BMD, as well as fracture risk, in human populations for one significant phosSNP, that is, rs227584 (major/minor allele: C/A, EAS population) located in C17orf53 gene, were suggested in prior meta-analyses. Specifically, carriers of allele C had significant higher BMD and lower risk of low-trauma fractures than carriers of A. We pursued to test the molecular and cellular functions of rs227584 in bone through osteoblastic cell culture and multiple assays. We identified five phosSNPs significant for OP (P < 0.01). The osteoblastic cells, which was transfected with wild-type C17orf53 (allele C at rs227584, P126), demonstrated specific interaction with NEK2 kinase, increased expression levels of osteoblastic genes significantly (OPN, OCN, COL1A1, P < 0.05), and promoted osteoblast growth and ALP activity, in contrast to those transfected with mutant C17orf53 (allele A at rs227584, T126). In the light of the consistent evidences between the present functional study in human bone cells and the prior association studies in human populations, we conclude that the SNP rs227584, via altering protein-kinase interaction, regulates osteoblastic gene expression, influences osteoblast growth and activity, hence to affect BMD and fracture risk in humans.

Multi-omics integrative analysis identified SNP-methylation-mRNA: Interaction in peripheral blood mononuclear cells.

Genetic variants have potential influence on DNA methylation and thereby regulate mRNA expression. This study aimed to comprehensively reveal the relationships among SNP, methylation and mRNA, and identify methylation-mediated regulation patterns in human peripheral blood mononuclear cells (PBMCs). Based on in-house multi-omics datasets from 43 Chinese Han female subjects, genome-wide association trios were constructed by simultaneously testing the following three association pairs: SNP-methylation, methylation-mRNA and SNP-mRNA. Causal inference test (CIT) was used to identify methylation-mediated genetic effects on mRNA. A total of 64,184 significant cis-methylation quantitative trait loci (meQTLs) were identified (FDR < 0.05). Among the 745 constructed trios, 464 trios formed SNP-methylation-mRNA regulation chains (CIT). Network analysis (Cytoscape 3.3.0) constructed multiple complex regulation networks among SNP, methylation and mRNA (eg a total of 43 SNPs simultaneously connected to cg22517527 and further to PRMT2, DIP2A and YBEY). The regulation chains were supported by the evidence from 4DGenome database, relevant to immune or inflammatory related diseases/traits, and overlapped with previous eQTLs from dbGaP and GTEx. The results provide new insights into the regulation patterns among SNP, DNA methylation and mRNA expression, especially for the methylation-mediated effects, and also increase our understanding of functional mechanisms underlying the established associations.

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.

Rheumatoid arthritis-associated DNA methylation sites in peripheral blood mononuclear cells.

OBJECTIVES: To identify novel DNA methylation sites significant for rheumatoid arthritis (RA) and comprehensively understand their underlying pathological mechanism. METHODS: We performed (1) genome-wide DNA methylation and mRNA expression profiling in peripheral blood mononuclear cells from RA patients and health controls; (2) correlation analysis and causal inference tests for DNA methylation and mRNA expression data; (3) differential methylation genes regulatory network construction; (4) validation tests of 10 differential methylation positions (DMPs) of interest and corresponding gene expressions; (5) correlation between PARP9 methylation and its mRNA expression level in Jurkat cells and T cells from patients with RA; (6) testing the pathological functions of PARP9 in Jurkat cells. RESULTS: A total of 1046 DNA methylation positions were associated with RA. The identified DMPs have regulatory effects on mRNA expressions. Causal inference tests identified six DNA methylation-mRNA-RA regulatory chains (eg, cg00959259-PARP9-RA). The identified DMPs and genes formed an interferon-inducible gene interaction network (eg, MX1, IFI44L, DTX3L and PARP9). Key DMPs and corresponding genes were validated their differences in additional samples. Methylation of PARP9 was correlated with mRNA level in Jurkat cells and T lymphocytes isolated from patients with RA. The PARP9 gene exerted significant effects on Jurkat cells (eg, cell cycle, cell proliferation, cell activation and expression of inflammatory factor IL-2). CONCLUSIONS: This multistage study identified an interferon-inducible gene interaction network associated with RA and highlighted the importance of PARP9 gene in RA pathogenesis. The results enhanced our understanding of the important role of DNA methylation in pathology of RA.

SAMD9 is a (epi-) genetically regulated anti-inflammatory factor activated in RA patients.

To identify PBMC-expressed genes significant for RA, and to ascertain their upstream regulatory factors, as well as downstream functional effects relevant to RA pathogenesis. We performed peripheral blood mononuclear cells (PBMCs) transcriptome-wide mRNA expression profiling in a case-control discovery sample. Differentially expressed genes (DEGs) were identified and validated in PBMCs in independent samples. We also generated genome-wide SNP genotyping data, and collected miRNA expression data and DNA methylation data from PBMCs of the discovery sample. Pearson correlation analyses were conducted to identify miRNAs/DNA methylations influencing DEG expression. Association analyses were conducted to identify expression-regulating SNPs. The key DEG, SAMD9, which was reported to function as a tumor suppressor gene, was assessed for its effects on T cell proliferation, apoptosis, and inflammatory cytokine expression. A total of 181 DEGs (Fold Change ≥ 2.0, Bonferroni adjusted p ≤ 0.05) were discovered in PBMCs. Four DEGs (SAMD9, CKLF, PARP9, and GUSB), upregulated with RA, were validated independently in PBMCs. Specifically, SAMD9 mRNA expression level was significantly upregulated in PHA-activated Jurkat T cells in vitro, and correlated with 8 miRNAs and associated with 22 SNPs in PBMCs in vivo. Knockdown of SAMD9 could transiently promote Jurkat T cell proliferation within 48 h and significantly induce TNF-α and IL-8 expression in T cells. SAMD9 expression is (epi-) genetically regulated, and significantly upregulated in PBMCs in RA patients and in activated T cells in vitro. SAMD9 might serve as a T cell activation marker but act as an anti-inflammatory factor.