Impaired NFKBIE gene function decreases cellular uptake of methotrexate by down-regulating SLC19A1 expression in a human rheumatoid arthritis cell line

Assessment of genetic polymorphisms within nuclear factor-κB signaling pathway genes in rheumatoid arthritis: Evidence for replication and genetic interaction

OBJECTIVE

This study was performed to replicate the associations of genetic polymorphisms within nuclear factor-κB (NF-κB) signaling pathway genes with rheumatoid arthritis (RA), and to further examine genetic interactions in a Chinese population.

METHODS

A total of eleven single-nucleotide polymorphisms (SNPs) were genotyped in 594 RA patients and 604 healthy controls.

RESULTS

Genetic association analysis revealed that NFKBIE rs2233434, TNIP1 rs10036748 and BLK rs13277113 were significantly associated with RA, cyclic citrullinated peptide (CCP)-positive RA and rheumatoid factor (RF)-positive RA, and TNFAIP3 rs2230926 was significantly associated with CCP-positive RA. Significant additive interaction was observed between NFKB1 rs28362491 and IKBKE rs12142086 (RERI = 0.76, 95% CI 0.13-1.38; AP = 0.57, 95% CI 0.11-1.03), NFKBIE rs2233434 and BLK rs13277113 (RERI = 1.41, 95% CI 0.88-1.94; AP = 0.85, 95% CI 0.50-1.20), NFKBIL rs2071592 and TNIP1 rs10036748 (RERI = 0.59, 95% CI 0.17-1.02; AP = 0.46, 95% CI 0.05-0.87), UBE2L3 rs5754217 and TNFSF4 rs2205960 (RERI = 0.50, 95% CI 0.16-0.84; AP = 0.57, 95% CI 0.09-1.05). Significant multiplicative interaction was detected between BLK rs13277113 and UBE2L3 rs5754217 (p = 0.02), BLK rs13277113 and TNFSF4 rs2205960 (p = 0.03).

CONCLUSIONS

Our results lent further support to the role of NF-κB signaling pathway in the pathogenesis of RA from a genetic perspective.

Identification of Candidate Risk Factor Genes for Human Idelalisib Toxicity Using a Collaborative Cross Approach

Idelalisib is a phosphatidylinositol 3-kinase inhibitor highly selective for the delta isoform that has shown good efficacy in treating chronic lymphocytic leukemia and follicular lymphoma. In clinical trials, however, idelalisib was associated with rare, but potentially serious liver and lung toxicities. In this study, we used the Collaborative Cross (CC) mouse population to identify genetic factors associated with the drug response that may inform risk management strategies for idelalisib in humans. Eight male mice (4 matched pairs) from 50 CC lines were treated once daily for 14 days by oral gavage with either vehicle or idelalisib at a dose selected to achieve clinically relevant peak plasma concentrations (150 mg/kg/day). The drug was well tolerated across all CC lines, and there were no observations of overt liver injury. Differences across CC lines were seen in drug concentration in plasma samples collected at the approximate Tmax on study Days 1, 7, and 14. There were also small but statistically significant treatment-induced alterations in plasma total bile acids and microRNA-122, and these may indicate early hepatocellular stress required for immune-mediated hepatotoxicity in humans. Idelalisib treatment further induced significant elevations in the total cell count of terminal bronchoalveolar lavage fluid, which may be analogous to pneumonitis observed in the clinic. Genetic mapping identified loci associated with interim plasma idelalisib concentration and the other 3 treatment-related endpoints. Thirteen priority candidate quantitative trait genes identified in CC mice may now guide interrogation of risk factors for adverse drug responses associated with idelalisib in humans.

Methylation analysis of the SLC19A1 promoter region in Chinese children with acute lymphoblastic leukaemia

WHAT IS KNOWN AND OBJECTIVE

Reduced folate carrier 1 (RFC1), which is encoded by the human solute carrier family 19 member 1 (SLC19A1) gene, plays an essential role in the cellular uptake of methotrexate (MTX). RFC1 expression is regulated by genetic variations and epigenetic modifications. The aim of the present study was to investigate the methylation status of the SLC19A1 promoter in peripheral blood and its association with MTX levels and toxicities in children with acute lymphoblastic leukaemia (ALL).

METHODS

Serum MTX concentrations were measured using a fluorescence polarization immunoassay. Methylation quantification for SLC19A1 promoter region #17 was performed by Sequenom MassARRAY in 52 paediatric ALL patients.

RESULTS AND DISCUSSION

Overall, the investigated region of the SLC19A1 promoter was in a hypermethylated state. No significant associations were detected between the methylation levels of six CpG units in the SLC19A1 promoter region #17 and clinical parameters of patients with ALL, including sex, age, immunotype and risk stratification. The methylation level of CpG_10 showed a significant positive correlation with MTX 24 hours after the initiation of infusion. No significant differences in the methylation levels of six CpG units were observed between patients with and without MTX toxicities. Due to the small sample size of this study, there was a high chance of false-positive results. A large-scale study would be required to confirm these preliminary results.

WHAT IS NEW AND CONCLUSION

Our preliminary results suggested the hypermethylated status of the SLC19A1 promoter in children with ALL. The methylation levels of the SLC19A1 promoter might affect MTX exposure. These findings have implications for the mechanisms underlying the variability of MTX responses in childhood ALL.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

Exploration of the immune-related signature and immune infiltration analysis for breast ductal and lobular carcinoma

Background

In this study, we aimed to explore the tumour associated immune signature of breast cancer (BC) and conduct integrative analyses with immune infiltrates in BC.

Methods

We downloaded the transcriptome profiling and clinical data of BC from The Cancer Genome Atlas (TCGA) database. The list of immune-related signatures was from the Innate database. The limma package was utilized to conduct the normalization, and we screened the differential immune signatures in BC. A univariate Cox regression model and the LASSO method were used to find the hub prognostic immune genes. The TAIG risk model was calculated based on the multivariate Cox regression results, and a receiver operating characteristic (ROC) curve was generated to assess the predictive power of TAIG. Moreover, we also conducted a correlation analysis between TAIG and the clinical characteristics. Additionally, we utilized the METABRIC cohort as the validation data set. The TIMER database is a comprehensive resource for performing systematic analyses of immune infiltrates across various malignancies. We evaluated the associations of immune signatures with several immune cells based on TIMER. Furthermore, we used the CIBERSORT algorithm to determine the fractions of immune cells in each sample and compared the differential distributions of immune infiltrates between two TAIG groups using the Wilcoxon rank-sum test.

Results

A total of 1,178 samples were obtained from the TCGA-BRCA database, but only 1,045 breast tumour samples were matched with complete transcriptome expression data. Meanwhile, we collected a total of 1,094 BC patients from the METABRIC cohort. We found a list of 1,399 differential immune signatures associated with survival, and functional analysis revealed that these genes participated in cytokine-cytokine receptor interactions, Th1 and Th2 cell differentiation and the JAK-STAT signalling pathway. The TAIG risk model was established from the multivariate Cox analysis, and we observed that high TAIG levels correlated with poor survival outcomes based on Kaplan-Meier analysis. The Kruskal-Wallis test suggested that high TAIG levels correlated with high AJCC-TNM stages and advanced pathological stages (P<0.01). We validated the well robustness of TAIG in METABRIC cohort and 5-year AUC reached up to 0.829. Moreover, we further uncovered the associations of hub immune signatures with immune cells and calculated the immune cell fractions in specific tumour samples based on gene signature expression. Last, we used the Wilcoxon rank-sum test to compare the differential immune density in the two groups and found that several immune cells had a significantly lower infiltrating density in the high TAIG groups, including CD8⁺ T cells (P=0.031), memory resting CD4⁺ T cells (P=0.026), M0 macrophages (P=0.023), and M2 macrophages (P=0.048).

Conclusions

In summary, we explored the immune signature of BC and constructed a TAIG risk model to predict prognosis. Moreover, we integrated the identified immune signature with tumour-infiltrating immune cells and found adverse associations between the TAIG levels and immune cell infiltrating density.

Methotrexate pharmacogenetics in the treatment of rheumatoid arthritis

For many decades, methotrexate (MXT) has remained the drug of choice in the treatment of rheumatoid arthritis (RA). Unfortunately, a considerable number of patients do not achieve an appropriate therapeutic response. Pharmacogenetics studies do not give usable results regarding differences in MTX response among RA patients. The mechanism of MTX action in RA is not completely understood. We present and discuss data regarding the molecular basis of folate and adenosine pathways, the most obvious MTX targets, to explain possible causes of therapy failure. The molecular basis of the disease could also have an impact on therapy outcomes and in this review we explore this. Finally, we make a short review of available pharmacogenetics study results.

Membrane-Spanning Protein Genetic Polymorphisms Related to Methotrexate Therapeutic Outcomes in a Chinese Rheumatoid Arthritis Population

Methotrexate (MTX) is a first-line disease-modifying antirheumatic drug for rheumatoid arthritis (RA), but individual variation in treatment response remains unexplained. The differences in drug efficacy and adverse drug reactions may be caused by genetic variations. We investigated the effects of single-nucleotide polymorphisms (SNPs) in 2 genes encoding membrane-spanning proteins, namely, reduced folate carrier-1 RFC-1/SLC19A1 (G>A [rs7499], A>G [rs2838956] and 180G>A [rs1051266]) and adenosine triphosphate-binding cassette B1 (rs1045642). Tagged SNPs were genotyped in 162 patients with RA in China. Then, we analyzed the relationships between these SNPs and therapeutic outcomes related to MTX in Chinese RA patients. No significant associations were found between the RFC-1/SLC19A1 (G>A [rs7499] and A>G [rs2838956]) and adenosine triphosphate-binding cassette B1 (rs1045642) gene polymorphisms and the response to MTX in RA patients. However, MTX-related toxicity was associated with one SNP, RFC-1 rs1051266 AA vs GG (odds ratio, 6.523; 95% confidence interval, 1.596-26.565; P = .009). SLC19A1 A>G rs2838956 showed a trend toward a significant association (odds ratio, 0.377; 95% confidence interval, 0.124-1.143; P = .085) with toxicity. Our results suggest that the RFC-1 80G>A (rs1051266) SNP exerts a potentially protective effect against the risk of adverse drug reactions in Chinese RA patients treated with MTX. Further studies are required to validate these findings.

SKA1 induces de novo MTX-resistance in osteosarcoma through inhibiting FPGS transcription

De novo methotrexate (MTX)-resistance, whose underlying mechanism remains largely unknown, usually leads to very poor prognosis in patients with osteosarcoma (OS). In this study, we established the de novo MTX-resistant OS cell line SF-86 and identified the candidate gene spindle and kinetochore associated complex subunit 1 (SKA1) as potentially related to de novo MTX-resistance. Analysis of a cohort of 95 OS patients demonstrated that SKA1 overexpression significantly correlated with de novo MTX-resistance and poor 5-year survival. Mechanistically, SKA1 overexpression lead to a downregulation of folylpoly-γ-glutamate synthetase (FPGS), a key enzyme that converts MTX into its active form, MTX-PG. We further demonstrated that SKA1 interacts with DNA-directed RNA polymerase II subunit RPB3. ChIP analysis revealed that RPB3 binds the promoter region of the FPGS gene and triggers FPGS transcription upon MTX treatment in SW1353, a MTX-sensitive OS cell line lacking endogenous SKA1 expression. On the contrary, this process is blocked in SF-86 cells due to the formation of an inhibitory SKA1-RPB3 complex. Furthermore, downregulation of SKA1 levels restores MTX sensitivity in SF-86. Collectively, our study has established the de novo MTX-resistant cell line SF-86 and identified SKA1 as a novel regulator of FPGS, playing a key role in the development of de novo MTX-resistance in OS.