Genome-wide identification of RNA modification-related single nucleotide polymorphisms associated with rheumatoid arthritis.

BACKGROUND: RNA modification plays important roles in many biological processes, such as gene expression control. The aim of this study was to identify single nucleotide polymorphisms related to RNA modification (RNAm-SNPs) for rheumatoid arthritis (RA) as putative functional variants. METHODS: We examined the association of RNAm-SNPs with RA in summary data from a genome-wide association study of 19,234 RA cases and 61,565 controls. We performed eQTL and pQTL analyses for the RNAm-SNPs to find associated gene expression and protein levels. Furthermore, we examined the associations of gene expression and circulating protein levels with RA using two-sample Mendelian randomization analysis methods. RESULTS: A total of 160 RNAm-SNPs related to m6A, m1A, A-to-I, m7G, m5C, m5U and m6Am modifications were identified to be significantly associated with RA. These RNAm-SNPs were located in 62 protein-coding genes, which were significantly enriched in immune-related pathways. RNAm-SNPs in important RA susceptibility genes, such as PADI2, SPRED2, PLCL2, HLA-A, HLA-B, HLA-DRB1, HLA-DPB1, TRAF1 and TXNDC11, were identified. Most of these RNAm-SNPs showed eQTL effects, and the expression levels of 26 of the modifiable genes (e.g., PADI2, TRAF1, HLA-A, HLA-DRB1, HLA-DPB1 and HLA-B) in blood cells were associated with RA. Circulating protein levels, such as CFB, GZMA, HLA-DQA2, IL21, LRPAP1 and TFF3, were affected by RNAm-SNPs and were associated with RA. CONCLUSION: The present study identified RNAm-SNPs in the reported RA susceptibility genes and suggested that RNAm-SNPs may affect RA risk by affecting the expression levels of corresponding genes and proteins.

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.

RNA modification-related variants in genomic loci associated with body mass index.

BACKGROUND: Genome-wide association studies (GWASs) have identified hundreds of loci for body mass index (BMI), but functional variants in these loci are less known. The purpose of this study was to identify RNA modification-related SNPs (RNAm-SNPs) for BMI in GWAS loci. BMI-associated RNAm-SNPs were identified in a GWAS of approximately 700,000 individuals. Gene expression and circulating protein levels affected by the RNAm-SNPs were identified by QTL analyses. Mendelian randomization (MR) methods were applied to test whether the gene expression and protein levels were associated with BMI. RESULTS: A total of 78 RNAm-SNPs associated with BMI (P < 5.0 × 10-8) were identified, including 65 m6A-, 10 m1A-, 3 m7G- and 1 A-to-I-related SNPs. Two functional loss, high confidence level m6A-SNPs, rs6713978 (P = 6.4 × 10-60) and rs13410999 (P = 8.2 × 10-59), in the intron of ADCY3 were the top significant SNPs. These two RNAm-SNPs were associated with ADCY3 gene expression in adipose tissues, whole blood cells, the tibial nerve, the tibial artery and lymphocytes, and the expression levels in these tissues were associated with BMI. Proteins enriched in specific KEGG pathways, such as natural killer cell-mediated cytotoxicity, the Rap1 signaling pathway and the Ras signaling pathway, were affected by the RNAm-SNPs, and circulating levels of some of these proteins (ADH1B, DOCK9, MICB, PRDM1, STOM, TMPRSS11D and TXNDC12) were associated with BMI in MR analyses. CONCLUSIONS: Our study identified RNAm-SNPs in BMI-related genomic loci and suggested that RNA modification may affect BMI by affecting the expression levels of corresponding genes and proteins.

A polygenic risk score improves risk stratification of coronary artery disease: a large-scale prospective Chinese cohort study.

AIMS: To construct a polygenic risk score (PRS) for coronary artery disease (CAD) and comprehensively evaluate its potential in clinical utility for primary prevention in Chinese populations. METHODS AND RESULTS: Using meta-analytic approach and large genome-wide association results for CAD and CAD-related traits in East Asians, a PRS comprising 540 genetic variants was developed in a training set of 2800 patients with CAD and 2055 controls, and was further assessed for risk stratification for CAD integrating with the guideline-recommended clinical risk score in large prospective cohorts comprising 41 271 individuals. During a mean follow-up of 13.0 years, 1303 incident CAD cases were identified. Individuals with high PRS (the highest 20%) had about three-fold higher risk of CAD than the lowest 20% (hazard ratio 2.91, 95% confidence interval 2.43-3.49), with the lifetime risk of 15.9 and 5.8%, respectively. The addition of PRS to the clinical risk score yielded a modest yet significant improvement in C-statistic (1%) and net reclassification improvement (3.5%). We observed significant gradients in both 10-year and lifetime risk of CAD according to the PRS within each clinical risk strata. Particularly, when integrating high PRS, intermediate clinical risk individuals with uncertain clinical decision for intervention would reach the risk levels (10-year of 4.6 vs. 4.8%, lifetime of 17.9 vs. 16.6%) of high clinical risk individuals with intermediate (20-80%) PRS. CONCLUSION: The PRS could stratify individuals into different trajectories of CAD risk, and further refine risk stratification for CAD within each clinical risk strata, demonstrating a great potential to identify high-risk individuals for targeted intervention in clinical utility.

Mendelian randomization analysis revealed potential causal factors for systemic lupus erythematosus.

Genome-wide association studies (GWAS) have identified many loci for systemic lupus erythematosus (SLE). However, identification of functionally relevant genes remains a challenge. The aim of this study was to highlight potential causal genes for SLE in the GWAS loci. By applying Mendelian randomization (MR) methods, such as summary data-based MR (SMR), generalized SMR and MR pleiotropy residual sum and outlier, we identified DNA methylations in 15 loci and mRNA expression of 21 genes that were causally associated with SLE. The identified genes enriched in 14 specific KEGG pathways (e.g. SLE, viral carcinogenesis) and two GO terms (interferon-γ-mediated signaling pathway and innate immune response). Among the identified genes, UBE2L3 and BLK variants were significantly associated with UBE2L3 and BLK methylations and gene expressions, respectively. UBE2L3 was up-regulated in SLE patients in several types of immune cells. Methylations (e.g. cg06850285) and mRNA expression of UBE2L3 were causally associated with SLE. Methylation site cg09528494 and mRNA expression of BLK were causally associated with SLE. BLK single nucleotide polymorphisms that were significantly associated with SLE were strongly associated with plasma cathepsin B level. Deep analysis identified that plasma cathepsin B level was causally associated with SLE. In summary, this study identified hundreds of DNA methylations and genes as potential risk factors for SLE. Genetic variants in UBE2L3 gene might affect SLE by influencing gene expression. Genetic variants in BLK gene might affect SLE by influencing BLK gene expression and plasma cathepsin B protein level.

Phosphorylation-related SNPs influence lipid levels and rheumatoid arthritis risk by altering gene expression and plasma protein levels.

OBJECTIVES: Phosphorylation-related single-nucleotide polymorphisms (phosSNPs) are missense SNPs that may influence protein phosphorylation. The aim of this study was to evaluate the effect of phosSNPs on lipid levels and RA. METHODS: We examined the association of phosSNPs with lipid levels and RA in large-scale genome-wide association studies (GWAS) and performed random sampling and fgwas analyses to determine whether the phosSNPs associated with lipid levels and RA were significantly enriched. Furthermore, we performed QTL analysis and Mendelian randomization analysis to obtain additional evidence to be associated with the identified phosSNPs and genes. RESULTS: We found 483 phosSNPs for lipid levels and 243 phosSNPs for RA in the GWAS loci (P < 1.0 × 10-5). SNPs associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, Total cholesterol (TC) and RA were significantly enriched with phosSNPs. Almost all of the identified phosSNPs showed expression quantitative trait loci (eQTL) effects. A total of 48 protein QTLs and 9 metabolite QTLs were found. The phosSNP rs3184504 (p.Trp262Arg) at SH2B3 was significantly associated with RA, SH2B3 expression level, and plasma levels of high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, TC, hypoxanthine and 80 proteins, including beta-2-microglobulin. SH2B3 was differentially expressed between RA cases and controls in peripheral blood mononuclear cells and synovial tissues. Mendelian randomization analysis showed that SH2B3 expression level was significantly associated with TC level and RA. Plasma beta-2-microglobulin level was causally associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, TC levels and RA. CONCLUSION: The findings suggested that phosSNPs may play important roles in lipid metabolism and the pathological mechanisms of RA. PhosSNPs may influence lipid levels and RA risk by altering gene expression and plasma protein levels.

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.

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.

Genome-wide enrichment of m6A-associated single-nucleotide polymorphisms in the lipid loci.

N6-methyladenosine (m6A) plays critical roles in many fundamental biological processes and a variety of diseases. The aim of this study was to investigate the effect of m6A-SNPs on lipid levels. We examined the association of m6A-SNPs with lipid levels in a genome-wide association studies (GWAS) of 188,578 individuals. Furthermore, we performed expression quantitative trait loci and differential expression analyses to add additional information for the identified m6A-SNPs. We found 1,655 m6A-SNPs in the GWAS dataset. Among them, 395 (23.9%) were nominally (P < 0.05) associated with lipid levels, and 22 reached the genome-wide significance level (P < 5.0 × 10-8). Using the fgwas method we found that SNPs, which influence high-density lipoprotein cholesterol (log2 enrichment of 3.35, 95% CI: (0.92, 4.48)) and TG (log enrichment of 3.22, 95% CI: (1.18, 4.44)), were enriched in m6A methylation. The high confidence (determined by miCLIP experiment) m6A-SNP rs6859 at the 3'-untranslated region of PVRL2 was associated with high-density lipoprotein cholesterol (P = 1.21 × 10-15), low-density lipoprotein cholesterol (P = 1.77 × 10-106), total cholesterol (P = 4.82 × 10-82), and triglycerides (P = 8.10 × 10-5) levels, coronary artery disease (P = 0.01), as well as PVRL2 mRNA expression in artery tibial (P = 2.38 × 10-6) and whole blood (P = 5.59 × 10-19). Moreover, PVRL2 was differentially expressed in adipose tissue of familial combined hyperlipidemia (P = 9.27 × 10-4). The present study found plenty of lipid-associated m6A-SNPs and demonstrated that m6A-SNPs may play important roles in lipid metabolisms. Further studies were needed to elucidate the mechanisms.

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.

Associations between potentially functional CORIN SNPs and serum corin levels in the Chinese Han population.

BACKGROUND: Corin is an important convertase involved in the natriuretic peptide system and may indirectly regulate blood pressure. Genetic factors relate to corin remain unclear. The purpose of the current study was to comprehensively examine the associations among CORIN SNPs, methylations, serum corin levels and hypertension. RESULTS: We genotyped 9 tag-SNPs in the CORIN gene and measured serum corin levels in 731 new-onset hypertensive cases and 731 age- and sex-matched controls. DNA methylations were tested in 43 individuals. Mendelian randomization was used to investigate the causal associations. Under additive models, we observed associations of rs2289433 (p.Cys13Tyr), rs6823184, rs10517195, rs2271037 and rs12509275 with serum corin levels after adjustment for covariates (P = 0.0399, 0.0238, 0.0016, 0.0148 and 0.0038, respectively). The tag-SNP rs6823184 and SNPs that are in strong linkage disequilibrium with it, i.e., rs10049713, rs6823698 and rs1866689, were associated with CORIN gene expression (P = 2.38 × 10- 24, 5.94 × 10- 27, 6.31 × 10- 27 and 6.30 × 10- 27, respectively). Neither SNPs nor corin levels was found to be associated with hypertension. SNP rs6823184, which is located in a DNase hypersensitivity cluster, a CpG island and transcription factor binding sites, was significantly associated with cg02955940 methylation levels (P = 1.54 × 10- 7). A putative causal association between cg02955940 methylation and corin levels was detected (P = 0.0011). CONCLUSION: This study identified potentially functional CORIN SNPs that were associated with serum corin level in the Chinese Han population. The effect of CORIN SNPs on corin level may be mediated by DNA methylation.

Detection of Putative Functional Single Nucleotide Polymorphisms in Blood Pressure Loci and Validation of Association Between Single Nucleotide Polymorphism in WBP1L and Hypertension in the Chinese Han Population.

We have performed a gene-based association study and detected several important blood pressure (BP)-associated genes. In this study, we explored functional variants in these genes by bioinformatics analysis and validated the associations between the functional single nucleotide polymorphisms (SNPs) and hypertension with public data and our in-house data of 857 cases and 927 controls. We found various functional variants in the BP-associated genes, including missense mutations and phosphorylation-related SNPs. Most of these SNPs were associated with expressions of the local genes. Some of these SNPs were associated with coronary artery disease or ischemic stroke. The associations between 12 functional SNPs in 7 genes and BP were validated (P < 5 × 10). The intronic SNP rs176185, which may influence promoter histone, enhancer histone, DNase and regulatory motifs and showed cis-eQTL effect on WBP1L, was associated with hypertension in the Chinese Han population (P = 0.0119). Our study detected plenty of potential functional SNPs in the BP-associated genes and demonstrated that rs176185 may be associated with hypertension in the Chinese Han population.

SNPs rs10224002 in PRKAG2 may disturb gene expression and consequently affect hypertension.

Genome-wide association studies have identified a large number of genetic loci for blood pressure in European populations. The associations in other populations are needed to determine. The purpose of this study was to examine the associations between the single nucleotide polymorphisms (SNPs) identified in European populations and hypertension in the Chinese Han population, and highlight the potential roles. Seven tag-SNPs were genotyped in 857 hypertension cases and 927 controls to test the associations. The intronic SNP rs10224002 (PRKAG2) which could affect DNase and regulatory motif was associated with hypertension (P = 0.024). This SNP was also found to be associated with coronary artery disease and stroke. We searched for potential functional variants by bioinformatics analysis and found various kinds of variants such as missense mutations, phosphorylation-related SNPs and SNPs that have regulatory potentials in the blood pressure loci. We performed expression quantitative trait locus (eQTL) and differential expression analyses for the identified variants and genes. eQTL analysis found that rs10224002 was associated with PRKAG2 gene expression in peripheral blood (P = 0.0016). PRKAG2 was differentially expressed between hypertension cases and controls (P = 0.0133), coronary artery disease cases and controls (P = 0.02112) and stroke cases and controls (P = 0.0059). Our study demonstrated that SNPs rs10224002 may be associated with hypertension in the Chinese Han population and PRKAG2 may play a role in the etiology of hypertension and cardiovascular diseases.

Integrative Analysis Confirmed the Association between Osteoprotegerin and Osteoporosis.

Objective This study aimed to verify the association between osteoprotegerin gene (OPG) and its variants with osteoporosis (OP) by performing integrative analysis.Methods We used the KGG software to perform gene-based association analysis, which integrated all publicly available single-nucleotide polymorphism (SNP)-based P values and obtained an overall P value for the OPG. The significant SNPs were screened for expression quantitative trait loci (eQTLs). Meta-analysis was used to combine the associations between the variants of OPG and bone mineral density (BMD) reported in the literatures. Then we performed dual-luciferase reporter gene systems for the functional verification of the variants of OPG in vitro.Results In the gene-based association analysis, the over all P value of OPG was 6.24×10 -13for BMD at femoral neck (FN) and 7.37×10 -17 for BMD at lumbar spine (LS), indicating the importance of OPG for OP. The publicly available eQTL database identified 5 eQTLs which exert cis-regulation effects on OPG at FN and LS. Literature searching found that rs2073617 (known as T950C) was the hot spot SNP. There were 13 relevant studies on rs2073617 besides the GEFOS-2 study identified from the PubMed. Significant differences among TT, TC and CC genotypes at FN (P= 0.047) and LS (P= 0.025) were shown by meta-analysis, demonstrating the associations between T950C polymorphism and BMD. Luciferase gene expression was significantly higher at the presence of allele C than allele T in the 293T cells (t=-9.47, P<0.01). Conclusion The integrative analysis further confirmed the importance of OPG in OP and the correlation of T950C polymorphism with BMD of OP. The strategy can be used as a reference for functional interpretation of other disease-related genes.

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.

Multiple correlation analyses revealed complex relationship between DNA methylation and mRNA expression in human peripheral blood mononuclear cells.

DNA methylation is an important regulator on the mRNA expression. However, a genome-wide correlation pattern between DNA methylation and mRNA expression in human peripheral blood mononuclear cells (PBMCs) is largely unknown. The comprehensive relationship between mRNA and DNA methylation was explored by using four types of correlation analyses and a genome-wide methylation-mRNA expression quantitative trait locus (eQTL) analysis in PBMCs in 46 unrelated female subjects. An enrichment analysis was performed to detect biological function for the detected genes. Single pair correlation coefficient (r T1) between methylation level and mRNA is moderate (-0.63-0.62) in intensity, and the negative and positive correlations are nearly equal in quantity. Correlation analysis on each gene (T4) found 60.1% genes showed correlations between mRNA and gene-based methylation at P < 0.05 and more than 5.96% genes presented very strong correlation (R T4 > 0.8). Methylation sites have regulation effects on mRNA expression in eQTL analysis, with more often observations in region of transcription start site (TSS). The genes under significant methylation regulation both in correlation analysis and eQTL analysis tend to cluster to the categories (e.g., transcription, translation, regulation of transcription) that are essential for maintaining the basic life activities of cells. Our findings indicated that DNA methylation has predictive regulation effect on mRNA with a very complex pattern in PBMCs. The results increased our understanding on correlation of methylation and mRNA and also provided useful clues for future epigenetic studies in exploring biological and disease-related regulatory mechanisms in PBMC.

Correlation analyses revealed global microRNA-mRNA expression associations in human peripheral blood mononuclear cells.

MicroRNAs (miRNAs) can regulate gene expression through binding to complementary sites in the 3'-untranslated regions of target mRNAs, which will lead to existence of correlation in expression between miRNA and mRNA. However, the miRNA-mRNA correlation patterns are complex and remain largely unclear yet. To establish the global correlation patterns in human peripheral blood mononuclear cells (PBMCs), multiple miRNA-mRNA correlation analyses and expression quantitative trait locus (eQTL) analysis were conducted in this study. We predicted and achieved 861 miRNA-mRNA pairs (65 miRNAs, 412 mRNAs) using multiple bioinformatics programs, and found global negative miRNA-mRNA correlations in PBMC from all 46 study subjects. Among the 861 pairs of correlations, 19.5% were significant (P < 0.05) and ~70% were negative. The correlation network was complex and highlighted key miRNAs/genes in PBMC. Some miRNAs, such as hsa-miR-29a, hsa-miR-148a, regulate a cluster of target genes. Some genes, e.g., TNRC6A, are regulated by multiple miRNAs. The identified genes tend to be enriched in molecular functions of DNA and RNA binding, and biological processes such as protein transport, regulation of translation and chromatin modification. The results provided a global view of the miRNA-mRNA expression correlation profile in human PBMCs, which would facilitate in-depth investigation of biological functions of key miRNAs/mRNAs and better understanding of the pathogenesis underlying PBMC-related diseases.

The distribution and functional relevance analysis of runs of homozygosity (ROHs) in Chinese Han female population.

Extended homozygosity is a genomic region in which the copies inherited from parents are identical, and has obvious inter-individual differences in length and frequency. Runs of homozygosity (ROHs), regarded as a type of structure variations, may have potential capacity in regulating gene transcription. To learn more about the genome-wide distribution of ROH regions in humans and understand the potential roles, this study applied ROH-based approach to quantify and characterize ROHs in 41 Chinese Han female subjects, and test potential associations between ROHs and mRNA expressions by eQTL analysis to ascertain whether ROHs are relevant to gene transcription in peripheral blood mononuclear cells (PBMCs). 10,884 ROH regions were identified in human genome. The average cumulative length of ROH regions was 217,250 ± 20,241 kb. The number of core segments in each chromosome generally matched the total length of corresponding chromosome, i.e., the longer the chromosome, the more the core segments. Genes located in the core regions of ROH were significantly enriched in multiple basic metabolism pathways. A total of 226 cis-eQTLs and 178 trans-eQTLs were identified. The cis-effect size was mainly concentrated at ± 0.5; and the trans-effect size was mainly concentrated at -1.5 and 1.0. Genes with eQTL effects were significantly enriched in functions related to protein binding, cytosol, nucleoplasm, nuclear membrane, protein binding and citrate metabolic process. This study described comprehensive distributions and characteristics of ROH in Han female Chinese, and recognized the significant role of ROH associated with gene transcription in human PBMC.

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.