Identification of novel rheumatoid arthritis-associated MiRNA-204-5p from plasma exosomes.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by infiltration of immune cells in the synovium. However, the crosstalk of immune cells and synovial fibroblasts is still largely unknown. Here, global miRNA screening in plasma exosomes was carried out with a custom microarray (RA patients vs. healthy controls = 9:9). A total of 14 exosomal miRNAs were abnormally expressed in the RA patients. Then, downregulated expression of exosomal miR-204-5p was confirmed in both the replication (RA patients vs. healthy controls = 30:30) and validation groups (RA patients vs. healthy controls = 56:60). Similar to the findings obtained in humans, a decreased abundance of exosomal miR-204-5p was observed in mice with collagen-induced arthritis (CIA). Furthermore, Spearman correlation analysis indicated that plasma exosomal miR-204-5p expression was inversely correlated with disease parameters of RA patients, such as rheumatoid factor, erythrocyte sedimentation rate, and C-reactive protein. In vitro, our data showed that human T lymphocytes released exosomes containing large amounts of miR-204-5p, which can be transferred into synovial fibroblasts, inhibiting cell proliferation. Overexpression of miR-204-5p in synovial fibroblasts suppressed synovial fibroblast activation by targeting genes related to cell proliferation and invasion. In vivo assays found that administration of lentiviruses expressing miR-204-5p markedly alleviated the disease progression of the mice with CIA. Collectively, this study identified a novel RA-associated plasma exosomal miRNA-204-5p that mediates the communication between immune cells and synovial fibroblasts and can be used as a potential biomarker for RA diagnosis and treatment.

Phosphorylation of transcription factor bZIP21 by MAP kinase MPK6-3 enhances banana fruit ripening.

Ripening of fleshy fruits involves both diverse post-translational modifications (PTMs) and dynamic transcriptional reprogramming, but the interconnection between PTMs, such as protein phosphorylation and transcriptional regulation, in fruit ripening remains to be deciphered. Here, we conducted a phosphoproteomic analysis during banana (Musa acuminata) ripening and identified 63 unique phosphopeptides corresponding to 49 proteins. Among them, a Musa acuminata basic leucine zipper transcription factor21 (MabZIP21) displayed elevated phosphorylation level in the ripening stage. MabZIP21 transcript and phosphorylation abundance increased during banana ripening. Genome-wide MabZIP21 DNA binding assays revealed MabZIP21-regulated functional genes contributing to banana ripening, and electrophoretic mobility shift assay, chromatin immunoprecipitation coupled with quantitative polymerase chain reaction, and dual-luciferase reporter analyses demonstrated that MabZIP21 stimulates the transcription of a subset of ripening-related genes via directly binding to their promoters. Moreover, MabZIP21 can be phosphorylated by MaMPK6-3, which plays a role in banana ripening, and T318 and S436 are important phosphorylation sites. Protein phosphorylation enhanced MabZIP21-mediated transcriptional activation ability, and transient overexpression of the phosphomimetic form of MabZIP21 accelerated banana fruit ripening. Additionally, MabZIP21 enlarges its role in transcriptional regulation by activating the transcription of both MaMPK6-3 and itself. Taken together, this study reveals an important machinery of protein phosphorylation in banana fruit ripening in which MabZIP21 is a component of the complex phosphorylation pathway linking the upstream signal mediated by MaMPK6-3 with transcriptional controlling of a subset of ripening-associated genes.

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.

Deciphering transcriptional regulators of banana fruit ripening by regulatory network analysis.

Fruit ripening is a critical phase in the production and marketing of fruits. Previous studies have indicated that fruit ripening is a highly coordinated process, mainly regulated at the transcriptional level, in which transcription factors play essential roles. Thus, identifying key transcription factors regulating fruit ripening as well as their associated regulatory networks promises to contribute to a better understanding of fruit ripening. In this study, temporal gene expression analyses were performed to investigate banana fruit ripening with the aim to discern the global architecture of gene regulatory networks underlying fruit ripening. Eight time points were profiled covering dynamic changes of phenotypes, the associated physiology and levels of known ripening marker genes. Combining results from a weighted gene co-expression network analysis (WGCNA) as well as cis-motif analysis and supported by EMSA, Y1H, tobacco-, banana-transactivation experimental results, the regulatory network of banana fruit ripening was constructed, from which 25 transcription factors were identified as prime candidates to regulate the ripening process by modulating different ripening-related pathways. Our study presents the first global view of the gene regulatory network involved in banana fruit ripening, which may provide the basis for a targeted manipulation of fruit ripening to attain higher banana and loss-reduced banana commercialization.

Framework Phylogeny, Evolution and Complex Diversification of Chinese Oaks.

Oaks (Quercus L.) are ideal models to assess patterns of plant diversity. We integrated the sequence data of five chloroplast and two nuclear loci from 50 Chinese oaks to explore the phylogenetic framework, evolution and diversification patterns of the Chinese oak's lineage. The framework phylogeny strongly supports two subgenera Quercus and Cerris comprising four infrageneric sections Quercus, Cerris, Ilex and Cyclobalanopsis for the Chinese oaks. An evolutionary analysis suggests that the two subgenera probably split during the mid-Eocene, followed by intergroup divergence within the subgenus Cerris around the late Eocene. The initial diversification of sections in the subgenus Cerris was dated between the mid-Oligocene and the Oligocene-Miocene boundary, while a rapid species radiation in section Quercus started in the late Miocene. Diversification simulations indicate a potential evolutionary shift on section Quercus, while several phenotypic shifts likely occur among all sections. We found significant negative correlations between rates of the lineage diversification and phenotypic turnover, suggesting a complex interaction between the species evolution and morphological divergence in Chinese oaks. Our infrageneric phylogeny of Chinese oaks accords with the recently proposed classification of the genus Quercus. The results point to tectonic activity and climatic change during the Tertiary as possible drivers of evolution and diversification in the Chinese oak's lineage.

Spatial Genetic Structure and Demographic History of the Dominant Forest Oak Quercus fabri Hance in Subtropical China.

Oak trees (Quercus L.) are important models for estimating abiotic impacts on the population structure and demography of long life span tree species. In this study, we generated genetic data for 17 nuclear microsatellite loci in 29 natural populations of Quercus fabri to estimate the population genetic structure. We also integrated approximate Bayesian computation (ABC) and ecological niche analysis to infer the population differentiation processes and demographic history of this oak species. The genetic analyses indicated two genetic clusters across the 29 populations collected, where most approximately corresponded to the intraspecific differentiation among populations from western and eastern China, whereas admixed populations were mainly found in central mountains of China. The best model obtained from hierarchical ABC simulations suggested that the initial intraspecific divergence of Q. fabri potentially occurred during the late Pliocene (ca. 3.99 Ma) to form the two genetic clusters, and the admixed population group might have been generated by genetic admixture of the two differentiated groups at ca. 53.76 ka. Ecological analyses demonstrated clear differentiation among the Q. fabri population structures, and association estimations also indicated significant correlations between geography and climate with the genetic variation in this oak species. Our results suggest abiotic influences, including past climatic changes and ecological factors, might have affected the genetic differentiation and demographic history of Q. fabri in subtropical China.

Banana MaBZR1/2 associate with MaMPK14 to modulate cell wall modifying genes during fruit ripening.

KEY MESSAGE: Banana MaBZR1/2 interact with MaMPK14 to enhance the transcriptional inhibition of cell wall modifying genes including MaEXP2, MaPL2 and MaXET5. Fruit ripening and softening, the major attributes to perishability in fleshy fruits, are modulated by various plant hormones and gene expression. Banana MaBZR1/2, the central transcription factors of brassinosteroid (BR) signaling, mediate fruit ripening through regulation of ethylene biosynthesis, but their possible roles in fruit softening as well as the underlying mechanisms remain to be determined. In this work, we found that MaBZR1/2 directly bound to and repressed the promoters of several cell wall modifying genes such as MaEXP2, MaPL2 and MaXET5, whose transcripts were elevated concomitant with fruit ripening. Moreover, yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays indicated that MaBZR1/2 physically interacted with a mitogen-activated protein kinase MaMPK14, and this interaction strengthened the MaBZR1/2-mediated transcriptional inhibitory abilities. Collectively, our study provides insight into the mechanism of MaBZR1/2 contributing to fruit ripening and softening, which may have potential for banana molecular improvement.

Integrative analysis revealed potential causal genetic and epigenetic factors for multiple sclerosis.

OBJECTIVE: Many genomic loci have been identified for multiple sclerosis (MS) by genome-wide association studies (GWAS). Discrimination of the most functionally relevant genes in these loci remains challenging. The aim of this study was to highlight potential causal genes for MS. METHODS: We detected potential causal DNA methylations and gene expressions for MS by integrating data from large scale GWAS and quantitative trait locus (QTL) studies using the summary data-based Mendelian randomization method. Potential functional SNPs in the identified genes were searched. RESULTS: We found 178 DNA methylation sites and mRNA expressions of 29 genes that were causally associated with MS. The identified genes enriched in 21 specific KEGG pathways and 80 GO terms (e.g., antigen processing and presentation, interferon gamma mediated signaling pathway). Among the identified non-MHC genes, METTL21B, METTL1 and TSFM were strongly connected. MS-associated SNPs in DDR1 were strongly associated with plasma MHC class I polypeptide-related sequence B (MICB) and Granzyme A levels. And plasma MICB and Granzyme A levels were causally associated with MS. Many SNPs in the causal genes showed QTL effects. The association between m6A-SNPs rs923829 and METTL21B expression level was validated in 40 unrelated Chinese Han individuals. CONCLUSIONS: This study identified many DNA methylations and genes as important risk factors for MS and provided novel evidence on the association between circulating MICB and Granzyme A and MS. We also showed that the interaction among DDR1, MICB and GZMA and interaction among METTL21B, METTL1 and TSFM may participate in the pathogenesis of MS.

MaBZR1/2 act as transcriptional repressors of ethylene biosynthetic genes in banana fruit.

Banana fruit (Musa acuminate L.) ripening is a complex genetical process affected by multiple phytohormones and expression of various genes. However, whether plant hormone brassinosteroid (BR) is involved in this process remains obscure. In this work, three genes that encode BR core signaling components brassinazole resistant (BZR) proteins, namely MaBZR1 to MaBZR3, were characterized from banana fruit. MaBZR1-MaBZR3 exhibited both nuclear and cytoplasmic localization and behaved as transcription inhibitors. Expression analysis showed that MaBZR1/2/3 were continuously decreased as fruit ripening proceeded, indicating their negative roles in banana ripening. Moreover, gel shift and transient expression assays demonstrated that MaBZR1/2 could suppress the transcription of ethylene biosynthetic genes, including MaACS1, MaACO13 and MaACO14, which increased gradually during the banana ripening, via specifically binding to CGTGT/CG sequence in their promoters. Importantly, exogenous application of BRs promotes banana ripening, which is presumably due to the accelerated expression of MaACS1 and MaACO13/14, and consequently the ethylene production. Our study indicates that MaBZR1/2 act as transcriptional repressors of ethylene biosynthetic genes during banana fruit ripening.

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.

Identification of Two Transcriptional Activators MabZIP4/5 in Controlling Aroma Biosynthetic Genes during Banana Ripening.

The transcriptional regulation of aroma formation genes remains poorly understood in the banana. In this work, we found that the expressions of a subset of aroma biosynthetic genes including MaOMT1, MaMT1, MaGT1, MaBCAT1, MaACY1, MaAGT1, and BanAAT, as well as two bZIP genes, MabZIP4 and MabZIP5, were down-regulated when prestored at 7 °C compared to those prestored at 22 °C during the ripening process of banana. Furthermore, MabZIP4 and MabZIP5 were shown to be able to activate the transcription of these aroma biosynthetic genes. Importantly, MabZIP4 directly binds to BanAAT promoter, while MabZIP5 binds to the promoters of MaMT1, MaACY1, MaAGT1, and BanAAT via the G-box motif, implicating the diverse functional significances of MabZIPs in controlling aroma biosynthesis in banana. Overall, this work sheds new insights on the understanding of transcriptional regulatory mechanisms associated with aroma formation during banana ripening.

Identification of expression quantitative trait loci (eQTLs) in human peripheral blood mononuclear cells (PBMCs) and shared with liver and brain.

PBMCs are essential for immunity and involved in various diseases. To identify genetic variations contributing to PBMCs transcriptome-wide gene expression, we performed a genome-wide eQTL analysis by using genome-wide SNPs data and transcriptome-wide mRNA expression data. To assess whether there are common regulation patterns shared among different tissues/organs, public datasets were utilized to identify common eQTLs shared with PBMCs in lymphoblastoid, monocytes, liver, and brain. Allelic expression imbalance (AEI) assay was employed to validate representative eQTLs identified. We identified 443 cis- and 2386 trans-eSNPs (FDR <0.05), which regulated 128 and 635 target genes, respectively. A transcriptome-wide expression regulation network was constructed, highlighting the importance of 28 pleiotropic eSNPs and 18 dually (cis- and trans-) regulated genes. Three genes, that is, TIPRL, HSPB8, and EGLN3, were commonly regulated by hundreds of eSNPs and constituted a very complex interaction network. Strikingly, the missense SNP rs371513 trans- regulated 25 target genes, which were functionally related to poly(A) RNA binding. Among 8904 eQTLs (P < 0.001) identified herein in PBMCs, a minority (163) was overlapped with lymphoblastoid, monocytes, liver, and/or brain. Besides, two cis-eSNPs in PBMC were confirmed by AEI. The present results demonstrated a comprehensive expression regulation network for human PBMCs and may provide novel insights into the pathogenesis of immunological diseases related to PBMCs.

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.

Integrative multi-omics analysis revealed SNP-lncRNA-mRNA (SLM) networks in human peripheral blood mononuclear cells.

Long non-coding RNAs (lncRNAs) serve as important controller of cellular functions via regulating RNA transcription, degradation and translation. However, what are the regulation patterns of lncRNAs on downstream mRNA and how the upstream genetic variants regulate lncRNAs are largely unknown. We first performed a comprehensive expression quantitative trait locus (eQTL) analysis (MatrixeQTL package, R) using genome-wide lncRNA expression and SNP genotype data from human peripheral blood mononuclear cells (PBMCs) of 43 unrelated individuals. Subsequently, multi-omics integrative network analysis was applied to construct SNP-lncRNA-mRNA (SLM) interaction networks. The causal inference test (CIT) was used to identify lncRNA-mediated (epi-) genetic regulation on mRNA expressions. Our eQTL analysis detected 707 pairs of cis-effect associations (p < 5.64E-06) and 6657 trans-effect associations (p < 3.51E-08), respectively. We also found that top significant cis-eSNPs were enriched around the lncRNA transcription start site regions, and that enrichment patterns of cis-eSNPs differs among different lncRNA sizes (small, medium and large).The constructed SLM interaction networks (1 primary networks and four small separate networks) showed various complex interaction patterns. Especially, the in-depth CIT detected 50 significant lncRNA-mediated SLM trios, and some hotspots (e.g., SNPs: rs926370, rs7716167 and rs16880521; lncRNAs: HIT000061975 and ENST00000579057.1). This study represents the first effort of dissecting the SLM interaction patterns in PBMCs by multi-omics integrative network analysis and causal inference test for clearing the regulation chain. The results provide novel insights into the regulation patterns of lncRNA, and may facilitate investigations of PBMC-related immune physiological process and immunological diseases in the future.

Integrative analysis for identification of shared markers from various functional cells/tissues for rheumatoid arthritis.

Rheumatoid arthritis (RA) is a systemic autoimmune disease. So far, it is unclear whether there exist common RA-related genes shared in different tissues/cells. In this study, we conducted an integrative analysis on multiple datasets to identify potential shared genes that are significant in multiple tissues/cells for RA. Seven microarray gene expression datasets representing various RA-related tissues/cells were downloaded from the Gene Expression Omnibus (GEO). Statistical analyses, testing both marginal and joint effects, were conducted to identify significant genes shared in various samples. Followed-up analyses were conducted on functional annotation clustering analysis, protein-protein interaction (PPI) analysis, gene-based association analysis, and ELISA validation analysis in in-house samples. We identified 18 shared significant genes, which were mainly involved in the immune response and chemokine signaling pathway. Among the 18 genes, eight genes (PPBP, PF4, HLA-F, S100A8, RNASEH2A, P2RY6, JAG2, and PCBP1) interact with known RA genes. Two genes (HLA-F and PCBP1) are significant in gene-based association analysis (P = 1.03E-31, P = 1.30E-2, respectively). Additionally, PCBP1 also showed differential protein expression levels in in-house case-control plasma samples (P = 2.60E-2). This study represented the first effort to identify shared RA markers from different functional cells or tissues. The results suggested that one of the shared genes, i.e., PCBP1, is a promising biomarker for RA.

Elevated Plasma P-Selectin Autoantibodies in Primary Sjögren Syndrome Patients with Thrombocytopenia.

BACKGROUND Primary Sjögren's syndrome (pSS) is one of the most common chronic systemic autoimmune diseases, and thrombocytopenia is one of the hematological manifestations of pSS. When platelet and endothelial cells are activated, P-selectin is expressed on the cell surface. This study aimed to investigate the role of P-selectin autoantibodies in the pathogenesis of thrombocytopenia in pSS. MATERIAL AND METHODS P-selectin autoantibodies were measured by enzyme-linked immunosorbent assay (ELISA) in 38 pSS patients without thrombocytopenia and 32 pSS patients with thrombocytopenia, 32 idiopathic thrombocytopenic purpura (ITP) patients, and 35 healthy controls. RESULTS The plasma P-selectin autoantibodies (A490) in ITP patients and pSS patients with/without thrombocytopenia were significantly higher than those in healthy controls, but there were no significant differences between ITP patients and pSS patients with thrombocytopenia. The positive rate of P-selectin autoantibodies in pSS patients with thrombocytopenia was significantly higher than that in ITP patients. The platelet count was lower in P-selectin autoantibodies-positive patients, while among pSS patients with thrombocytopenia, the platelet count was lower in P-selectin autoantibodies-positive patients than in P-selectin autoantibodies-negative patients. In ITP patients and pSS patients with thrombocytopenia, the platelet count was lower in P-selectin autoantibodies-positive patients. CONCLUSIONS Elevated plasma P-selectin autoantibodies may play a role in the pathogenesis of thrombocytopenia in pSS patients.

Functional relevance for associations between genetic variants and systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a serious prototype autoimmune disease characterized by chronic inflammation, auto-antibody production and multi-organ damage. Recent association studies have identified a long list of loci that were associated with SLE with relatively high statistical power. However, most of them only established the statistical associations of genetic markers and SLE at the DNA level without supporting evidence of functional relevance. Here, using publically available datasets, we performed integrative analyses (gene relationship across implicated loci analysis, differential gene expression analysis and functional annotation clustering analysis) and combined with expression quantitative trait loci (eQTLs) results to dissect functional mechanisms underlying the associations for SLE. We found that 14 SNPs, which were significantly associated with SLE in previous studies, have cis-regulation effects on four eQTL genes (HLA-DQA1, HLA-DQB1, HLA-DQB2, and IRF5) that were also differentially expressed in SLE-related cell groups. The functional evidence, taken together, suggested the functional mechanisms underlying the associations of 14 SNPs and SLE. The study may serve as an example of mining publically available datasets and results in validation of significant disease-association results. Utilization of public data resources for integrative analyses may provide novel insights into the molecular genetic mechanisms underlying human diseases.