The Crosstalk of Pathways Involved in Immune Response Maybe the Shared Molecular Basis of Rheumatoid Arthritis and Type 2 Diabetes

Rheumatoid arthritis (RA) and Type 2 diabetes (T2D) are both systemic diseases linked with altered immune response, moderate mortality when present together. The treatment for both RA and T2D are not satisfied, partly because of the linkage between them has not yet been appreciated. A comprehensive study for the potential associations between the two disorders is needed. In this study, we used RNA sequencing to explore the differently expressed genes (DEGs) in peripheral blood mononuclear cells (PBMC) of 10 RA and 10 T2D patients comparing with 10 healthy volunteers (control). We used bioinformatics analysis and the Ingenuity Pathways Analysis (IPA) to predict the commonalities on signaling pathways and molecular networks between those two diseases. 212 DEGs in RA and 114 DEGs in T2D patients were identified compared with healthy controls, respectively. 32 DEGs were shared between the two comparisons. The top 10 shared pathways interacted in cross-talking networks, regulated by 5 shared predicted upstream regulators, leading to the activated immune response were explored, which was considered as partly of the association mechanism of this two disorders. These discoveries would be considered as new understanding on the associations between RA and T2D, and provide novel treatment or prevention strategy.     

Knowledge-based analysis of genetic associations of rheumatoid arthritis to inform studies searching for pleiotropic genes: a literature review and network analysis

IntroductionPleiotropy describes the genetic effect of a single gene on multiple phenotypic traits. Gene variants directly affect the normal processes of a series of physiological and biochemical reactions, and therefore cause a variety of diseases traits to be changed accordingly. Moreover, a shared genetic susceptibility mechanism may exist between different diseases. Therefore, shared genes, with pleiotropic effects, are important to understand the sharing pathogenesis and hence the mechanisms underlying comorbidity.MethodsIn this study, we proposed combining genome-wide association studies (GWAS) and public knowledge databases to search for potential pleiotropic genes associated with rheumatoid arthritis (RA) and eight other related diseases. Here, a GWAS-based network analysis is used to recognize risk genes significantly associated with RA. These RA risk genes are re-extracted as potential pleiotropic genes if they have been proved to be susceptible genes for at least one of eight other diseases in the OMIM or PubMed databases.ResultsIn total, we extracted 116 potential functional pleiotropic genes for RA and eight other diseases, including five hub pleiotropic genes, BTNL2, HLA-DRA, NOTCH4, TNXB, and C6orf10, where BTNL2, NOTCH4, and C6orf10 are novel pleiotropic genes identified by our analysis.ConclusionsThis study demonstrates that pleiotropy is a common property of genes associated with disease traits. Our results ascertained the shared genetic risk profiles that predisposed individuals to RA and other diseases, which could have implications for identification of molecular targets for drug development, and classification of diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0715-1) contains supplementary material, which is available to authorized users.     

The Paradox of High Availability and Low Recognition of Soluble HLA-G by LILRB1 Receptor in Rheumatoid Arthritis Patients

HLA-G is a regulatory molecule involved in immunologic tolerance. Growing evidence indicates that HLA-G plays a role in the regulation of inflammatory processes and autoimmune diseases. This study aimed at a systematic evaluation of soluble HLA-G (sHLA-G) in plasma of rheumatoid arthritis (RA) patients with long-lasting chronic inflammation. RA patients (n=68) and healthy controls (n=26) had their plasmatic sHLA-G measured by ELISA whereas the binding capability of sHLA-G to its cognate LILRB1 receptor was measured by a Luminex-based assay. All subjects were PCR-genotyped for HLA-G 14bp polymorphism (rs66554220). Significantly higher sHLA-G levels were observed in patients (p<0.001), however no significant differences were observed in LILRB1 binding capacity between RA patients and controls. Remarkably, the proportion of patients presenting specific binding of sHLA-G to LILRB1 was significantly decreased as compared to controls (56% vs. 81%, p=0.027). Patients without rheumatoid factor (RF-) were significantly overrepresented in the group of patients positive for LILRB1 binding as compared to patients without LILRB1 binding (31% vs 10%, p=0.033). Furthermore, methotrexate treated patients (n=58) revealed significantly lower LILRB1 binding to sHLA-G molecules than non-treated patients (medians: 12.2 vs. 67.7 units/ml, p=0.031). Unlike in controls, no significant differences in sHLA-G levels were observed among patients grouped by 14pb genotype. Thus, in a substantial number of late RA patients, the circulating sHLA-G molecules are impaired regarding LILRB1 recognition, meaning that although increased levels are observed; these molecules are not qualified to exert their protective functions against inflammation. Our findings offer new insights into the immunopathology of RA patients with long-lasting anti-RA-treatment and highlight the importance to also measure the binding capability of sHLA-G to LILRB1.     

Determination of oxidant stress in plasma of rheumatoid arthritis and primary osteoarthritis patients

Determination of oxidant stress in plasma of rheumatoid arthritis (RA) and primary osteoarthritis (POA) patients is important in understanding the pathogenesis of these diseases. In this study, we examined the relationship between oxidant stress and inflammation by measuring protein carbonyl content, thiol levels and plasma protein fractions as the oxidation markers and erythrocyte sedimentation rate (ESR), C-reactive protein (CRP) tests as inflammation markers. Protein carbonyls content was higher in RA and POA patients, as compared to controls (p<0.0001), while the plasma thiol levels in both groups of patients were significantly lower than controls (p<0.0001). Increased levels of proteins under 40 kDa molecular mass were detected in the RA and POA patients compared to that of controls (p<0.0001) both in HPLC and SDS-PAGE analysis. Total protein concentration in plasma of RA patients was higher than the controls (p<0.001), while in POA patients was lower than that of controls (p<0.001). ESR and CRP levels were higher in both the patient groups than the normal group (p<0.001). These results suggested that alterations in the oxidant stress markers could be the cause of inflammation in these diseases. Thus, while working for RA/POA treatment strategies, consideration of the relationship between oxidant stress and inflammation would be worth evaluating.     

Integrated microRNA–mRNA analysis of coronary artery disease

Although patients with coronary artery disease (CAD) have a high mortality rate, the pathogenesis of CAD is still poorly understood. The purpose of this study was to explore the underlying molecular mechanisms and potential target molecules for CAD. The platelet miRNA (GSE28858) and blood mRNA (GSE42148) expression profiles of patients with CAD and healthy controls were downloaded from Gene Expression Omnibus. Differentially expressed miRNAs and genes (DEGs) were identified by significant analysis of microarray algorithm after data preprocessing. Furthermore, the miRNA-target gene regulatory network was constructed based on miRecords database. The spearman correlation coefficients (ρ) between miRNAs and their target genes were calculated. Six up- (miR-340, miR-545, miR-451, miR454-5p, miR-624 and miR-585) and four down-regulated (miR-199a, miR-17-3p, miR-154 and miR-339) miRNAs were screened. Total 295 target genes of miR-545, miR-451, miR-585 and miR-154 were predicted. Among these 295 target genes, 7 genes were DEGs. Further analysis showed miR-545-TFEC and miR-585-SPOCK1 were highly positively correlated (ρ = 0.808091264; ρ = 0.874680776) in CAD samples. Therefore, differentially expressed miRNAs might participate in the pathogenesis of CAD by regulating their target genes.     

Common Genetic Determinants of Lung Function,Subclinical Atherosclerosis and Risk of Coronary Artery Disease

Chronic obstructive pulmonary disease (COPD) independently associates with an increased risk of coronary artery disease (CAD), but it has not been fully investigated whether this co-morbidity involves shared pathophysiological mechanisms. To identify potential common pathways across the two diseases, we tested all recently published single nucleotide polymorphisms (SNPs) associated with human lung function (spirometry) for association with carotid intima-media thickness (cIMT) in 3,378 subjects with multiple CAD risk factors, and for association with CAD in a case-control study of 5,775 CAD cases and 7,265 controls. SNPs rs2865531, located in the CFDP1 gene, and rs9978142, located in the KCNE2 gene, were significantly associated with CAD. In addition, SNP rs9978142 and SNP rs3995090 located in the HTR4 gene, were associated with average and maximal cIMT measures. Genetic risk scores combining the most robustly spirometry–associated SNPs from the literature were modestly associated with CAD, (odds ratio (OR) (95% confidence interval (CI₉₅) = 1.06 (1.03, 1.09); P-value = 1.5×10⁻⁴, per allele). In conclusion, our study suggests that some genetic loci implicated in determining human lung function also influence cIMT and susceptibility to CAD. The present results should help elucidate the molecular underpinnings of the co-morbidity observed across COPD and CAD.     

Identification of differentially expressed genes and the role of PDK4 in CD14+ monocytes of coronary artery disease

Background. Coronary artery disease (CAD) is a chronic inflammatory disease caused by development of atherosclerosis (AS), which is the leading cause of mortality and disability. Our study aimed to identify the differentially expressed genes (DEGs) in CD14⁺ monocytes from CAD patients compared with those from non-CAD controls, which might pave the way to diagnosis and treatment for CAD. Methods. The RNA-sequencing (RNA-seq) was performed by BGISEQ-500, followed by analyzing with R package to screening DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by R package. In addition, we validated the results of RNA-seq using real-time quantitative polymerase chain reaction (RT-qPCR). Furthermore, we explored the function of selected ten genes in LDL-treated CD14⁺ monocytes by RT-qPCR. Results. a total of 2897 DEGs were identified, including 753 up- and 2144 down-regulated genes in CD14⁺ monocytes from CAD patients. These DEGs were mainly enriched in plasma membrane and cell periphery of cell component, immune system process of biological process, NF-κB signaling pathway, cell adhesion molecules signaling pathway and cytokine–cytokine receptor interaction signaling pathway. In LDL-treated CD14⁺ monocytes, the mRNA expression of pyruvate dehydrogenase kinase 4 (PDK4) was significantly up-regulated. Conclusion. In the present study, we suggested that PDK4 might play a role in progression of CAD. The study will provide some pieces of evidence to investigate the role and mechanism of key genes in the pathogenesis of CAD.     

Identification of novel biomarkers and candidate small molecule drugs in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis of high‐throughput data

Background: Rheumatoid arthritis (RA) and osteoarthritis (OA) are two major types of joint diseases. The present study aimed to identify hub genes involved in the pathogenesis and further explore the potential treatment targets of RA and OA.Methods: The gene expression profile of {"type":"entrez-geo","attrs":{"text":"GSE12021","term_id":"12021"}}GSE12021 was downloaded from Gene Expression Omnibus (GEO). Total 31 samples (12 RA, 10 OA and 9 NC samples) were used. The differentially expressed genes (DEGs) in RA versus NC, OA versus NC and RA versus OA groups were screened using limma package. We also verified the DEGs in {"type":"entrez-geo","attrs":{"text":"GSE55235","term_id":"55235"}}GSE55235 and {"type":"entrez-geo","attrs":{"text":"GSE100786","term_id":"100786"}}GSE100786. Functional annotation and protein–protein interaction (PPI) network construction of OA‐ and RA‐specific DEGs were performed. Finally, the candidate small molecules as potential drugs to treat RA and OA were predicted in CMap database.Results: 165 up-regulated and 163 down-regulated DEGs between RA and NC samples, 73 up-regulated and 293 down-regulated DEGs between OA and NC samples, 92 up-regulated and 98 down-regulated DEGs between RA and OA samples were identified. Immune response and TNF signaling pathway were significantly enriched pathways for RA‐ and OA‐specific DEGs, respectively. The hub genes were mainly associated with ‘Primary immunodeficiency’ (RA vs. NC group), ‘Ribosome’ (OA vs. NC group), and ‘Chemokine signaling pathway’ (RA vs. OA group). Arecoline and Cefamandole were the most promising small molecule to reverse the RA and OA gene expression.Conclusion: Our findings suggest new insights into the underlying pathogenesis of RA and OA, which may improve the diagnosis and treatment of these intractable chronic diseases.     

A meta-analysis of genome-wide association scans identifies IL18RAP, PTPN2, TAGAP, and PUS10 as shared risk loci for Crohn's disease and celiac disease

Crohn''s disease (CD) and celiac disease (CelD) are chronic intestinal inflammatory diseases, involving genetic and environmental factors in their pathogenesis. The two diseases can co-occur within families, and studies suggest that CelD patients have a higher risk to develop CD than the general population. These observations suggest that CD and CelD may share common genetic risk loci. Two such shared loci, IL18RAP and PTPN2, have already been identified independently in these two diseases. The aim of our study was to explicitly identify shared risk loci for these diseases by combining results from genome-wide association study (GWAS) datasets of CD and CelD. Specifically, GWAS results from CelD (768 cases, 1,422 controls) and CD (3,230 cases, 4,829 controls) were combined in a meta-analysis. Nine independent regions had nominal association p-value <1.0×10⁻⁵ in this meta-analysis and showed evidence of association to the individual diseases in the original scans (p-value <1×10⁻² in CelD and <1×10⁻³ in CD). These include the two previously reported shared loci, IL18RAP and PTPN2, with p-values of 3.37×10⁻⁸ and 6.39×10⁻⁹, respectively, in the meta-analysis. The other seven had not been reported as shared loci and thus were tested in additional CelD (3,149 cases and 4,714 controls) and CD (1,835 cases and 1,669 controls) cohorts. Two of these loci, TAGAP and PUS10, showed significant evidence of replication (Bonferroni corrected p-values <0.0071) in the combined CelD and CD replication cohorts and were firmly established as shared risk loci of genome-wide significance, with overall combined p-values of 1.55×10⁻¹⁰ and 1.38×10⁻¹¹ respectively. Through a meta-analysis of GWAS data from CD and CelD, we have identified four shared risk loci: PTPN2, IL18RAP, TAGAP, and PUS10. The combined analysis of the two datasets provided the power, lacking in the individual GWAS for single diseases, to detect shared loci with a relatively small effect.     

Protective role of TIRAP functional variant against development of coronary artery disease

Coronary artery disease (CAD) is the leading cause of sudden death worldwide. Inflammation is proved to be an important player in development of the CAD. Inflammation is directly regulated by the Toll like receptors (TLRs). Susceptibility of CAD is influenced by genetic variations within TLRs and the proteins involved in its signaling cascade. The TIRAP/MAL {TIR domain containing adaptor protein / MyD88 (myeloid differentiation primary response gene 88) adaptor-like} exhibits maximum genetic variations of all adaptor proteins involved in TLR signaling cascade. Susceptibility to a number of diseases can be influenced due to presence of S180L single nucleotide polymorphism (SNP) of TIRAP/MAL. This study was conducted to investigate the functional role of this well characterized S180L polymorphism on susceptibility to CAD among Pakistani patients. A total of 146 Pakistani CAD patients and 147 controls were genotyped by Amplification Refractory Mutation System-Polymerase Chain Reaction (ARMS-PCR) and the data was analyzed by using 2-tailed Chi square (x²) test. The p value ≤ 0.05 was considered to be significant.Significantly high frequency of homozygous L180L genotype was observed among healthy subjects as compared to the CAD patients [24 (16%) vs 7 (5%); x² 11.85; p = 0.003]. Moreover, the allele frequency of the minor allele; 180L was observed to be significantly higher among controls than the CAD patients, having same direction of association [156 (53%) vs 131 (45%); OR (95% CI) = 0.7198 (0.520–0.996); p < 0.05).Our results indicate that protective effect of L180L; a coding variant of TIRAP/MAL against CAD is discernible.     

Overlapping genetic susceptibility variants between three autoimmune disorders: rheumatoid arthritis, type 1 diabetes and coeliac disease

Introduction

Genome wide association studies, replicated by numerous well powered validation studies, have revealed a large number of loci likely to play a role in susceptibility to many multifactorial diseases. It is now well established that some of these loci are shared between diseases with similar aetiology. For example, a number of autoimmune diseases have been associated with variants in the PTPN22, TNFAIP3 and CTLA4 genes. Here we have attempted to define overlapping genetic variants between rheumatoid arthritis (RA), type 1 diabetes (T1D) and coeliac disease (CeD).

Methods

We selected eight SNPs previously identified as being associated with CeD and six T1D-associated SNPs for validation in a sample of 3,962 RA patients and 3,531 controls. Genotyping was performed using the Sequenom MassArray platform and comparison of genotype and allele frequencies between cases and controls was undertaken. A trend test P-value < 0.004 was regarded as significant.

Results

We found statistically significant evidence for association of the TAGAP locus with RA (P = 5.0 × 10^-4). A marker at one other locus, C1QTNF6, previously associated with T1D, showed nominal association with RA in the current study but did not remain statistically significant at the corrected threshold.

Conclusions

In exploring the overlap between T1D, CeD and RA, there is strong evidence that variation within the TAGAP gene is associated with all three autoimmune diseases. Interestingly a number of loci appear to be specific to one of the three diseases currently studied suggesting that they may play a role in determining the particular autoimmune phenotype at presentation.     

PTPN22 genetic variation: evidence for multiple variants associated with rheumatoid arthritis

The minor allele of the R620W missense single-nucleotide polymorphism (SNP) (rs2476601) in the hematopoietic-specific protein tyrosine phosphatase gene, PTPN22, has been associated with multiple autoimmune diseases, including rheumatoid arthritis (RA). These genetic data, combined with biochemical evidence that this SNP affects PTPN22 function, suggest that this phosphatase is a key regulator of autoimmunity. To determine whether other genetic variants in PTPN22 contribute to the development of RA, we sequenced the coding regions of this gene in 48 white North American patients with RA and identified 15 previously unreported SNPs, including 2 coding SNPs in the catalytic domain. We then genotyped 37 SNPs in or near PTPN22 in 475 patients with RA and 475 individually matched controls (sample set 1) and selected a subset of markers for replication in an additional 661 patients with RA and 1,322 individually matched controls (sample set 2). Analyses of these results predict 10 common (frequency >1%) PTPN22 haplotypes in white North Americans. The sole haplotype found to carry the previously identified W620 risk allele was strongly associated with disease in both sample sets, whereas another haplotype, identical at all other SNPs but carrying the R620 allele, showed no association. R620W, however, does not fully explain the association between PTPN22 and RA, since significant differences between cases and controls persisted in both sample sets after the haplotype data were stratified by R620W. Additional analyses identified two SNPs on a single common haplotype that are associated with RA independent of R620W, suggesting that R620W and at least one additional variant in the PTPN22 gene region influence RA susceptibility.     

Systematic Analysis of the Gene Expression in the Livers of Nonalcoholic Steatohepatitis: Implications on Potential Biomarkers and Molecular Pathological Mechanism

Non-alcoholic steatohepatitis (NASH) is a severe form of non-alcoholic fatty liver disease (NAFLD). The molecular pathological mechanism of NASH is poorly understood. Recently, high throughput data such as microarray data together with bioinformatics methods have become a powerful way to identify biomarkers and to investigate pathogenesis of diseases. Taking advantage of well characterized microarray datasets of NASH livers, we performed a systematic analysis of potential biomarkers and possible pathological mechanism of NASH from a bioinformatics perspective.CodeLink Human Whole Genome Bioarrays were analyzed to find differentially expressed genes (DEGs) between controls and NASH patients. Four methods were used to identify DEGs and the intersection of DEGs identified by these methods was subsequently used for both biomarker prediction and molecular pathological mechanism analysis. For biomarker prediction, rank aggregation was used to rank DEGs identified by all these methods according to their significance of different expression. Alcohol dehydrogenase 4 (ADH4) exhibited the highest rank suggesting the most significant differential expression between normal and disease condition. Together with the previous report demonstrating the association between ADH4 and the pathogenesis of NASH, our data suggest that ADH4 could be a potential biomarker for NASH. For molecular pathological mechanism analysis, two clusters of highly correlated annotation terms and genes in these terms were identified based on the intersection of DEGs. Then, pathways enriched with these genes were identified to construct the network. Using this network, both for the first time, amino acid catabolism is implicated to play a pivotal role and urea cycle is implicated to be involved in the development of NASH.The results of our study identified potential biomarkers and suggested possible molecular pathological mechanism of NASH. These findings provide a comprehensive and systematic understanding of the pathogenesis of NASH and may facilitate the diagnosis, prevention and treatment of NASH.