V-J combinations of T-cell receptor predict responses to erythropoietin in end-stage renal disease patients

Immune-associated biomarkers identification for diagnosing carotid plaque progression with uremia through systematical bioinformatics and machine learning analysis

BACKGROUND

Uremia is one of the most challenging problems in medicine and an increasing public health issue worldwide. Patients with uremia suffer from accelerated atherosclerosis, and atherosclerosis progression may trigger plaque instability and clinical events. As a result, cardiovascular and cerebrovascular complications are more likely to occur. This study aimed to identify diagnostic biomarkers in uremic patients with unstable carotid plaques (USCPs).

METHODS

Four microarray datasets (GSE37171, GSE41571, GSE163154, and GSE28829) were downloaded from the NCBI Gene Expression Omnibus database. The Limma package was used to identify differentially expressed genes (DEGs) in uremia and USCP. Weighted gene co-expression network analysis (WGCNA) was used to determine the respective significant module genes associated with uremia and USCP. Moreover, a protein-protein interaction (PPI) network and three machine learning algorithms were applied to detect potential diagnostic genes. Subsequently, a nomogram and a receiver operating characteristic curve (ROC) were plotted to diagnose USCP with uremia. Finally, immune cell infiltrations were further analyzed.

RESULTS

Using the Limma package and WGCNA, the intersection of 2795 uremia-related DEGs and 1127 USCP-related DEGs yielded 99 uremia-related DEGs in USCP. 20 genes were selected as candidate hub genes via PPI network construction. Based on the intersection of genes from the three machine learning algorithms, three hub genes (FGR, LCP1, and C5AR1) were identified and used to establish a nomogram that displayed a high diagnostic performance (AUC: 0.989, 95% CI 0.971-1.000). Dysregulated immune cell infiltrations were observed in USCP, showing positive correlations with the three hub genes.

CONCLUSION

The current study systematically identified three candidate hub genes (FGR, LCP1, and C5AR1) and established a nomogram to assist in diagnosing USCP with uremia using various bioinformatic analyses and machine learning algorithms. Herein, the findings provide a foothold for future studies on potential diagnostic candidate genes for USCP in uremic patients. Additionally, immune cell infiltration analysis revealed that the dysregulated immune cell proportions were identified, and macrophages could have a critical role in USCP pathogenesis.

Divergent Characteristics of T-Cell Receptor Repertoire Between Essential Hypertension and Aldosterone-Producing Adenoma

Aldosterone-producing adenoma (APA) is a benign adrenal tumor that results in persistent hyperaldosteronism. As one major subtype of primary aldosteronism, APA leads to secondary hypertension that is associated with immune dysregulation. However, how the adaptive immune system, particularly the T-cell population, is altered in APA patients remains largely unknown. Here, we performed TCR sequencing to characterize the TCR repertoire between two age-matched groups of patients: one with APA and the other one with essential hypertension (EH). Strikingly, we found a significant reduction of TCR repertoire diversity in the APA group. Analyses on TCR clustering and antigen annotation further showed that the APA group possessed lower diversity in TCR clonotypes with non-common antigen-specific features, compared with the EH group. In addition, our results indicated that the strength of correlation between generation probabilities and frequencies of TCR clonotypes was significantly higher in the APA group than that in the EH group. Finally, we observed that clinical features, including plasma aldosterone level, aldosterone–renin ratio, and blood sodium level, were positively associated with the strength of correlation between generation and abundance of TCR clonotypes in the APA group. Our findings unveiled the correlation between T-cell immune repertoire and APA, suggesting a critical role of such adrenal adenoma in the T-cell immunity of patients with hypertension.

pmTR database: population matched (pm) germline allelic variants of T-cell receptor (TR) loci

The IMGT database profiles the TR germline alleles for all four TR loci (TRA, TRB, TRG and TRD), however, it does not comprise of the information regarding population specificity and allelic frequencies of these germline alleles. The specificity of allelic variants to different human populations can, however, be a rich source of information when studying the genetic basis of population-specific immune responses in disease and in vaccination. Therefore, we meticulously identified true germline alleles enriched with complete TR allele sequences and their frequencies across 26 different human populations, profiled by "1000 Genomes data". We identified 205 TRAV, 249 TRBV, 16 TRGV and 5 TRDV germline alleles supported by at least four haplotypes. The diversity of germline allelic variants in the TR loci is the highest in Africans, while the majority of the Non-African alleles are specific to the Asian populations, suggesting a diverse profile of TR germline alleles in different human populations. Interestingly, the alleles in the IMGT database are frequent and common across all five super-populations. We believe that this new set of germline TR sequences represents a valuable new resource which we have made available through the new population-matched TR (pmTR) database, accessible via https://pmtrig.lumc.nl/ .

The role of genetic polymorphisms in STIM1 and ORAI1 for erythropoietin resistance in patients with renal failure

ABSTRACT

Anemia is a common complication in patients with renal failure. While erythropoietin is commonly used to treat anemia, some patients exhibit a poor response to erythropoietin. Since store-operated calcium channel (SOC) signaling is one of the erythropoietin activated pathways, we aimed to investigate the association between the genetic polymorphisms of SOC signaling pathway and erythropoietin resistance in patients with renal failure.Four tagging single nucleotide polymorphisms in STIM1 and five in ORAI1 were selected in this study. Genotyping was performed with the TaqMan Allelic Discrimination assay and the association of individual tagging single nucleotide polymorphisms with erythropoietin resistance was analyzed by multivariable adjusted random intercepts model.194 patients were enrolled in this study. The mean age of participants is 68 years, and 56% were men. The mean erythropoietin resistance index was 9.04 ± 4.51 U/Kg/week/g/dL. We found that patients with the AA genotype of rs1561876 in STIM1, and the CC or CT genotypes of rs6486795 in ORAI1, were associated with increased risk of erythropoietin resistance. Functional annotation of expression quantitative trait loci revealed that the AA genotype of rs1561876 in STIM1 has a relatively lower expression of ribonucleotide reductase catalytic subunit M1 in skeletal muscle, while the CC genotype of rs6486795 in ORAI1 has a relatively higher expression of ORAI1 in the whole blood and thyroid.Overall, we demonstrate a significant association between erythropoietin resistance and genetic polymorphisms of STIM1 and ORAI1. Annotation prediction revealed the importance of SOC-mediated calcium signaling for erythropoietin resistance.

Replacement of TCR Dβ With Immunoglobulin DH DSP2.3 Imposes a Tyrosine-Enriched TCR Repertoire and Adversely Affects T Cell Development

Enrichment for tyrosine in immunoglobulin CDR-H3 is due in large part to natural selection of germline immunoglobulin DH sequence. We have previously shown that when DH sequence is modified to reduce the contribution of tyrosine codons, epitope recognition is altered and B cell development, antibody production, autoantibody production, and morbidity and mortality following pathogen challenge are adversely affected. TCRβ diversity (Dβ) gene segment sequences are even more highly conserved than DH, with trout Dβ1 identical to human and mouse Dβ1. We hypothesized that natural selection of Dβ sequence also shapes CDR-B3 diversity and influences T cell development and T cell function. To test this, we used a mouse strain that lacked Dβ2 and contained a novel Dβ1 allele (DβYTL) that replaces Dβ1 with an immunoglobulin DH, DSP2.3. Unlike Dβ1, wherein glycine predominates in all three reading frames (RFs), in DSP2.3 there is enrichment for tyrosine in RF1, threonine in RF2, and leucine in RF3. Mature T cells using DβYTL expressed TCRs enriched at particular CDR-B3 positions for tyrosine but depleted of leucine. Changing Dβ sequence altered thymocyte and peripheral T cell numbers and the T cell response to an ovalbumin immunodominant epitope. The differences in tyrosine content might explain, at least in part, why TCRs are more polyspecific and of lower affinity for their cognate antigens than their immunoglobulin counterparts.

Preimmune Control of the Variance of TCR CDR-B3: Insights Gained From Germline Replacement of a TCR Dβ Gene Segment With an Ig D H Gene Segment

We have previously shown that the sequence of the immunoglobulin diversity gene segment (D H ) helps dictate the structure and composition of complementarity determining region 3 of the immunoglobulin heavy chain (CDR-H3). In order to test the role of germline D sequence on the diversity of the preimmune TCRβ repertoire of T cells, we generated a mouse with a mutant TCRβ DJC locus wherein the Dβ2-Jβ2 gene segment cluster was deleted and the remaining diversity gene segment, Dβ1 (IMGT:TRDB1), was replaced with DSP2.3 (IMGT:IGHD2-02), a commonly used B cell immunoglobulin D H gene segment. Crystallographic studies have shown that the length and thus structure of TCR CDR-B3 places amino acids at the tip of CDR-B3 in a position to directly interact with peptide bound to an MHC molecule. The length distribution of complementarity determining region 3 of the T cell receptor beta chain (CDR-B3) has been proposed to be restricted largely by MHC-specific selection, disfavoring CDR-B3 that are too long or too short. Here we show that the mechanism of control of CDR-B3 length depends on the Dβ sequence, which in turn dictates exonucleolytic nibbling. By contrast, the extent of N addition and the variance of created CDR3 lengths are regulated by the cell of origin, the thymocyte. We found that the sequence of the D and control of N addition collaborate to bias the distribution of CDR-B3 lengths in the pre-immune TCR repertoire and to focus the diversity provided by N addition and the sequence of the D on that portion of CDR-B3 that is most likely to interact with the peptide that is bound to the presenting MHC.

Characterization of T-Cell Receptor Repertoire in Patients with Rheumatoid Arthritis Receiving Biologic Therapies

Rheumatoid arthritis (RA) is a systematic autoimmune disease, predominantly causing chronic polyarticular inflammation and joint injury of patients. For the treatment of RA, biologic disease-modifying antirheumatic drugs (bDMARDs) have been used to reduce inflammation and to interfere with disease progression through targeting and mediating the immune system. Although the therapeutic effects of bDMARDs in RA patients have been widely reported, whether these drugs also play important roles in T-cell repertoire status is still unclear. We therefore designed the study to identify the role of T-cell repertoire profiles in RA patients with different types of bDMARD treatments. A high-throughput sequencing approach was applied to profile the T-cell receptor beta chain (TCRB) repertoire of circulating T lymphocytes in eight patients given adalimumab (anti-TNF-α) with/without the following use of either rituximab (anti-CD20) or tocilizumab (anti-IL6R). We subsequently analyzed discrepancies in the clonal diversity and CDR3 length distribution as well as usages of the V and J genes of TCRB repertoire and interrogated the association between repertoire diversity and disease activities followed by the treatment of bDMARDs in these RA patients. All groups of patients showed well-controlled DAS28 scores (<2.6) after different treatment regimens of drugs and displayed no significant statistical differences in repertoire diversity, distribution of CDR3 lengths, and usage of V and J genes of TCRB. Nonetheless, a trend between overall TCRB repertoire diversity and disease activity scores in all bDMARD-treated RA patients was observed. Additionally, age was found to be associated with repertoire diversity in RA patients treated with bDMARDs. Through the profiling of the TCR repertoire in RA patients receiving different biologic medications, our study indicated an inverse tendency between TCR repertoire diversity and disease activity after biologic treatment in RA patients.

Association of endothelin genetic variants and hospitalized infection complications in end-stage renal disease (ESRD) patients

BACKGROUND

Infection is the second most common cause of mortality for patients with end-stage renal disease (ESRD), accompanying with immune dysfunction. Endothelin (EDN) is known to be related to inflammation; however, it is unknown whether genetic variants of the EDN gene family are associated with increased risk of hospitalized infection events.

METHODS

Nineteen tagging single-nucleotide polymorphisms (tSNPs) of the EDN gene family were selected for genotyping a cohort of 190 ESRD patients. Patient demographics were recorded, the subtypes of infection events were identified, and association analysis between the EDN genetic variants and hospitalized infection events was performed.

RESULTS

In this study, 106 patients were hospitalized for infection events. The leading events were pneumonia, bacteremia, and cellulitis. The minor allele of rs260741, rs197173, and rs926632 SNPs of EDN3 were found to be associated with reduced risk of hospitalized bacteremia events.

CONCLUSIONS

The minor allele of rs260741, rs197173, and rs926632 in EDN3 were associated with reduced risk of hospitalized bacteremia events in ESRD patients.

Physiology and pathophysiology of renal erythropoietin-producing cells

Anemia is a common complication and contributes to increased morbidity and mortality in chronic kidney disease (CKD) patients. Whereas there has been a significant improvement of understanding the underlying mechanism of erythropoiesis, the treatment of renal anemia is still restricted to erythropoietin (EPO)-stimulating agents. The purpose of this article is to review the physiology of erythropoiesis, functional role of EPO and underlying molecular and cellular basis that regulate EPO production. Regulation of EPO production is at mRNA level. When anemia or hypoxia occurs, transcriptional factor, hypoxia-inducible factor (HIF), binds to EPO 5' hypoxic response element and EPO gene transcription increases. The renal EPO is mainly produced by pericytes. In CKD, pericytes transdifferentiate to myofibroblasts, and subsequently the ability of EPO production decreases, leading to renal anemia. Recent experimental and clinical studies show the promising efficacy of prolyl hydroxylase inhibitors in renal anemia through increasing EPO production by stabilizing HIF. Recent advances on epigenetics create a new field to study EPO gene expression at chromatin level. We will discuss the role of demethylating agent on restoring EPO expression, providing a novel approach to the treatment of renal anemia.