Genome-wide association studies and gene expression profiles of rheumatoid arthritis: An analysis

Genetic polymorphisms of ACE I/D, IL-1β G > A and IL-4 VNTR among Egyptian subjects with rheumatoid arthritis

BACKGROUND

Various reports examined the contribution of ACE I/D, IL-1β G > A and IL-4 VNTR with the susceptibility to RA but with conflicting findings. The goal of this study is to assess the impact of these three variants with the susceptibility, clinical and biochemical markers in addition to different composite indices of RA.

SUBJECTS AND METHODS

This case-control survey enclosed 120 RA Egyptian patients who were emulated with 150 healthy controls from the same territory. Their genomic DNA was genotyped utilising the PCR approach.

RESULTS

RA patients with ACE allele (D) and IL-4 VNTR allele (2 R) were expressively higher emulated with healthy controls (p < .05). Nevertheless, RA patients with IL-1β (A allele) failed to achieve an apparent significant emulated with healthy controls (p > .05).

CONCLUSIONS

This work specifies a noteworthy association for ACE I/D and IL-4 VNTR but not IL-1β G > A polymorphisms with RA susceptibility among Egyptian subjects.

Functional Analysis of Autoantibody Signatures in Rheumatoid Arthritis

For the identification of antigenic protein biomarkers for rheumatoid arthritis (RA), we conducted IgG profiling on high density protein microarrays. Plasma IgG of 96 human samples (healthy controls, osteoarthritis, seropositive and seronegative RA, n = 24 each) and time-series plasma of a pristane-induced arthritis (PIA) rat model (n = 24 total) were probed on AIT’s 16k protein microarray. To investigate the analogy of underlying disease pathways, differential reactivity analysis was conducted. A total of n = 602 differentially reactive antigens (DIRAGs) at a significance cutoff of p < 0.05 were identified between seropositive and seronegative RA for the human samples. Correlation with the clinical disease activity index revealed an inverse correlation of antibodies against self-proteins found in pathways relevant for antigen presentation and immune regulation. The PIA model showed n = 1291 significant DIRAGs within acute disease. Significant DIRAGs for (I) seropositive, (II) seronegative and (III) PIA were subjected to the Reactome pathway browser which also revealed pathways relevant for antigen presentation and immune regulation; of these, seven overlapping pathways had high significance. We therefore conclude that the PIA model reflects the biological similarities of the disease pathogenesis. Our data show that protein array analysis can elucidate biological differences and pathways relevant in disease as well be a useful additional layer of omics information.

A 9 mRNAs-based diagnostic signature for rheumatoid arthritis by integrating bioinformatic analysis and machine-learning

Background

Rheumatoid arthritis (RA) is an autoimmune rheumatic disease that carries a substantial burden for both patients and society. Early diagnosis of RA is essential to prevent disease progression and select an optimal therapeutic strategy. However, RA diagnosis is challenging, partly due to a lack of reliable biomarkers. Here, we aimed to explore the diagnostic signature and establish a predictive model of RA.

Methods

The mRNA expression profiling data of GSE17755, containing blood samples of 112 RA patients and 53 healthy control patients, were obtained from the Gene Expression Omnibus (GEO) database, followed by differential expression, GO (Gene Ontology), and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis. A PPI network was constructed to select candidate hub genes, then logistic regression and random forest models were established based on the identified genes.

Results

Significantly, we identified 52 differentially expressed genes (DEGs), including 16 upregulated genes and 36 downregulated genes in RA samples compared with control samples. GO and KEGG analysis showed that several immune-related cellular processes were particularly enriched. We identified nine hub genes in the PPI network, including CFL1, COTL1, ACTG1, PFN1, LCP1, LCK, HLA-E, FYN, and HLA-DRA. The logistic regression and random forest models based on the nine identified genes reliably distinguished the RA samples from the healthy samples with substantially high AUC.

Conclusion

The diagnostic logistic regression and random forest models based on nine hub genes reliably predicted the occurrence of RA. Our findings could provide new insights into RA diagnostics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-020-02180-w.

LILRB1 Intron 1 Has a Polymorphic Regulatory Region That Enhances Transcription in NK Cells and Recruits YY1

LILRB1 is a highly polymorphic receptor expressed by subsets of innate and adaptive immune cells associated with viral and autoimmune diseases and targeted by pathogens for immune evasion. LILRB1 expression on human NK cells is variegated, and the frequency of LILRB1⁺ cells differs among people. However, little is known about the processes and factors mediating LILRB1 transcription in NK cells. LILRB1 gene expression in lymphoid and myeloid cells arises from two distinct promoters that are separated by the first exon and intron. In this study, we identified a polymorphic 3-kb region within LILRB1 intron 1 that is epigenetically marked as an active enhancer in human lymphoid cells and not monocytes. This region possesses multiple YY1 sites, and complexes of the promoter/enhancer combination were isolated using anti-YY1 in chromatin immunoprecipitation-loop. CRISPR-mediated deletion of the 3-kb region lowers LILRB1 expression in human NKL cells. Together, these results indicate the enhancer in intron 1 binds YY1 and suggest YY1 provides a scaffold function enabling enhancer function in regulating LILRB1 gene transcription in human NK cells.

HLA-C: An Accomplice in Rheumatic Diseases

Human leukocyte antigen c (HLA-C) is a polymorphic membrane protein encoded by the HLA-C gene in the class I major histocompatibility complex. HLA-C plays an essential role in protection against cancer and viruses but has also been implicated in allograft rejection, preeclampsia, and autoimmune disease. This review summarizes reports and proposed mechanisms for the accessory role of HLA-C in rheumatic diseases. Historically, contributions of HLA-C to rheumatic diseases were eclipsed by the stronger association with HLA-DRB1 alleles containing the "shared epitope" with rheumatoid arthritis. Larger genetic association studies and more powerful analytical approaches have revealed independent associations of HLA-C with rheumatic disease-associated phenotypes, including development of anticitrullinated peptide antibodies. HLA-C functions by presenting antigens to T cells and by binding activatory and inhibitory receptors on natural killer (NK) cells, but the exact mechanisms by which the HLA-C locus contributes to autoimmunity are largely undefined. Studies have suggested that HLA-C and NK cell receptor polymorphisms may predict responsiveness to pharmacotherapy. Understanding the mechanisms of the role of HLA-C in rheumatic disease could uncover therapeutic targets or guide precision pharmacologic treatments.

Characterization of the full-length human Grb7 protein and a phosphorylation representative mutant

It is well known the dimerization state of receptor tyrosine kinases (RTKs), in conjunction with binding partners such as the growth factor receptor bound protein 7 (Grb7) protein, plays an important role in cell signaling regulation. Previously, we proposed, downstream of RTKs, that the phosphorylation state of Grb7SH2 domain tyrosine residues could control Grb7 dimerization, and dimerization may be an important regulatory step in Grb7 binding to RTKs. In this manner, additional dimerization-dependent regulation could occur downstream of the membrane-bound kinase in RTK-mediated signaling pathways. Extrapolation to the full-length (FL) Grb7 protein, and the ability to test this hypothesis further, has been hampered by the availability of large quantities of pure and stable FL protein. Here, we report the biophysical characterization of the FL Grb7 protein and also a mutant representing a tyrosine-phosphorylated Grb7 protein form. Through size exclusion chromatography and analytical ultracentrifugation, we show the phosphorylated-tyrosine-mimic Y492E-FL-Grb7 protein (Y492E-FL-Grb7) is essentially monomeric at expected physiological concentrations. It has been shown previously the wild-type FL Grb7(WT-FLGrb7) protein is dimeric with a dissociation constant (Kd) of approximately 11μM. Our studies here measure a FL protein dimerization K_d of WT-FL-Grb7 within one order of magnitude at approximately 1μM. The approximate size and shape of the WT-FL-Grb7 in comparison the tyrosine-phosphorylation mimic Y492E-FL-Grb7 protein was determined by dynamic light scattering methods. In vitro phosphorylation of the Grb7SH2 domain indicates only one of the available tyrosine residues is phosphorylated, suggesting the same phosphorylation pattern could be relevant in the FL protein. The biophysical characterization studies in total are interpreted with a view towards understanding the functionally active Grb7 protein conformation.

Maternal pre-pregnancy obesity, offspring cord blood DNA methylation, and offspring cardiometabolic health in early childhood: an epigenome-wide association study

Pre-pregnancy obesity is an established risk factor for adverse sex-specific cardiometabolic health in offspring. Epigenetic alterations, such as in DNA methylation (DNAm), are a hypothesized link; however, sex-specific epigenomic targets remain unclear. Leveraging data from the Newborn Epigenetics Study (NEST) cohort, linear regression models were used to identify CpG sites in cord blood leukocytes associated with pre-pregnancy obesity in 187 mother-female and 173 mother-male offsprings. DNAm in cord blood was measured using the Illumina HumanMethylation450k BeadChip. Replication analysis was conducted among the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Associations between pre-pregnancy obesity-associated CpG sites and offspring BMI z-score (BMIz) and blood pressure (BP) percentiles at 4-5-years of age were also examined. Maternal pre-pregnacy obesity was associated with 876 CpGs in female and 293 CpGs in male offspring (false discovery rate <5%). Among female offspring, 57 CpG sites, including the top 18, mapped to the TAPBP gene (range of effect estimates: -0.83% decrease to 4.02% increase in methylation). CpG methylation differences in the TAPBP gene were also observed among males (range of effect estimates: -0.30% decrease to 2.59% increase in methylation). While technically validated, none of the TAPBP CpG sites were replicated in ALSPAC. In NEST, methylation differences at CpG sites of the TAPBP gene were associated with BMI z-score (cg23922433 and cg17621507) and systolic BP percentile (cg06230948) in female and systolic (cg06230948) and diastolic (cg03780271) BP percentile in male offspring. Together, these findings suggest sex-specific effects, which, if causal, may explain observed sex-specific effects of maternal obesity.

Genetic variants associated with patent ductus arteriosus in extremely preterm infants

OBJECTIVE

Patent ductus arteriosus (PDA) is a commonly observed condition in preterm infants. Prior studies have suggested a role for genetics in determining spontaneous ductal closure. Using samples from a large neonatal cohort we tested the hypothesis that common genetic variations are associated with PDA in extremely preterm infants.

STUDY DESIGN

Preterm infants (n = 1013) enrolled at NICHD Neonatal Research Network sites were phenotyped for PDA. DNA was genotyped for 1634 single nucleotide polymorphisms (SNPs) from candidate genes. Analyses were adjusted for ancestral eigenvalues and significant epidemiologic variables.

RESULTS

SNPs in several genes were associated with the clinical diagnosis of PDA and with surgical ligation in extremely preterm neonates diagnosed with PDA (p < 0.01). None of the associations were significant after correction for multiple comparisons.

CONCLUSION

We identified several common genetic variants associated with PDA. These findings may inform further studies on genetic risk factors for PDA in preterm infants.

Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis

Objectives

The aim of this study was to identify key pathological genes in osteoarthritis (OA).

Methods

We searched and downloaded mRNA expression data from the Gene Expression Omnibus database to identify differentially expressed genes (DEGs) of joint synovial tissues from OA and normal individuals. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analyses were used to assess the function of identified DEGs. The protein-protein interaction (PPI) network and transcriptional factors (TFs) regulatory network were used to further explore the function of identified DEGs. The quantitative real-time polymerase chain reaction (qRT-PCR) was applied to validate the result of bioinformatics analysis. Electronic validation was performed to verify the expression of selected DEGs. The diagnosis value of identified DEGs was accessed by receiver operating characteristic (ROC) analysis.

Results

A total of 1085 DEGs were identified. KEGG pathway analysis displayed that Wnt was a significantly enriched signalling pathway. Some hub genes with high interactions such as USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, and ZNF423 were identified in the PPI and TFs network. The results of qRT-PCR showed that GADD45B, ADAMTS1, and TFAM were down-regulated in joint synovial tissues of OA, which was consistent with the bioinformatics analysis. The expression levels of USP46, CPVL, FOSL2, and PTGS1 in electronic validation were compatible with the bio-informatics result. CPVL and TFAM had a potential diagnostic value for OA based on the ROC analysis.

Conclusion

The deregulated genes including USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, ZNF423, ADAMTS1, and TFAM might be involved in the pathology of OA.

Cite this article: X. Zhang, Y. Bu, B. Zhu, Q. Zhao, Z. Lv, B. Li, J. Liu. Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis. Bone Joint Res 2018;7:298–307. DOI: 10.1302/2046-3758.74.BJR-2017-0245.R1.

Grb7 protein RA domain oligomerization

The growth factor receptor bound protein 7 (Grb7) is an adaptor protein that is often coamplified with the erythroblastosis oncogene B 2 receptor in 20% to 30% of breast cancer patients. Grb7 overexpression has been linked to increased cell migration and cancer metastasis. The ras associating and pleckstrin homology domain region of Grb7 has been reported to interact with various other downstream signaling proteins such as four and half Lin11, Isl-1, Mec-3 (LIM) domains isoform 2 and filamin α. These interactions are believed to play a role in regulating Grb7-mediated cell migration function. The full-length Grb7 protein has been shown to dimerize, and the oligomeric state of the Grb7SH2 domain has been extensively studied; however, the oligomerization state of the ras associating and pleckstrin homology domains, and the importance of this oligomerization in Grb7 function, is yet to be fully known. In this study, we characterize the oligomeric state of the Grb7RA domain using size exclusion chromatography, nuclear magnetic resonance, nuclear relaxation studies, glutaraldehyde cross linking, and dynamic light scattering. We report the Grb7RA domain can exist in transient multimeric forms and, based upon modeling results, postulate the potential role of Grb7RA domain oligomerization in Grb7 function.

IL-32 promoter SNP rs4786370 predisposes to modified lipoprotein profiles in patients with rheumatoid arthritis

Patients with rheumatoid arthritis (RA) are at higher risk of developing cardiovascular diseases (CVD). Interleukin (IL)-32 has previously been shown to be involved in the pathogenesis of RA and might be linked to the development of atherosclerosis. However, the exact mechanism linking IL-32 to CVD still needs to be elucidated. The influence of a functional genetic variant of IL-32 on lipid profiles and CVD risk was therefore studied in whole blood from individuals from the NBS cohort and RA patients from 2 independent cohorts. Lipid profiles were matched to the specific IL-32 genotypes. Allelic distribution was similar in all three groups. Interestingly, significantly higher levels of high density lipoprotein cholesterol (HDLc) were observed in individuals from the NBS cohort and RA patients from the Nijmegen cohort homozygous for the C allele (p = 0.0141 and p = 0.0314 respectively). In contrast, the CC-genotype was associated with elevated low density lipoprotein cholesterol (LDLc) and total cholesterol (TC) in individuals at higher risk for CVD (plaque positive) (p = 0.0396; p = 0.0363 respectively). Our study shows a functional effect of a promoter single-nucleotide polymorphism (SNP) in IL32 on lipid profiles in RA patients and individuals, suggesting a possible protective role of this SNP against CVD.