Survival of mature T cells in the periphery is intrinsically dependent on GIMAP1 in mice

Investigation and verification of GIMAP6 as a robust biomarker for prognosis and tumor immunity in lung adenocarcinoma

BACKGROUND AND AIM

According to previous reports, GTPase of immunity-associated protein 6 (GIMAP6) is essential for autophagy. However, it is unclear how GIMAP6 affects the development and tumor immunity of lung adenocarcinoma (LUAD).

METHODS

In the present study, the role of GIMAP6 in vivo and in vitro was examined using reverse transcription-quantitative PCR, western blotting, and Cell Counting Kit-8, colony formation and Transwell assays. Datasets from The Cancer Genome Atlas and Genotype-Tissue Expression databases were thoroughly analyzed using R software. A nomogram was created using GIMAP6 and prognostic characteristics. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and Gene Set Enrichment Analysis were applied to explore the potential mechanism of GIMAP6 in lung cancer. The link between GIMAP6 and the immunological landscape was studied using single-cell RNA sequencing datasets from Tumor Immune Estimation Resource (TIMER) 2.0 and Tumor Immune Single-cell Hub.

RESULTS

Patients with high GIMAP6 expression had improved overall and disease-specific survival compared with those patients with low GIMAP6 expression. According to the receiver operating characteristic and calibration curve, the nomogram based on T stage, N stage and GIMAP6 had predictive value for prognosis. According to functional enrichment analysis, GIMAP6 was primarily involved in T-cell receptor signaling pathway, chemokine signaling pathway, cytokine and cytokine receptor interaction. GIMAP6 was shown to be favorably linked with the infiltration of immune cells and immune-related molecules, including cytotoxic T-lymphocyte associated protein 4, programmed death-ligand 1, and T cell immunoreceptor with Ig and ITIM domains, by single-cell sequencing and TIMER2.0 analysis. The role of GIMAP6 in lung cancer cell proliferation, invasion, migration and immunity was experimentally verified.

CONCLUSION

These findings confirmed that GIMAP6 was an effective prognostic molecule that was involved in the regulation of the immune microenvironment of LUAD, and may become a predictor for the efficacy of immunotherapy.

Radiogenomic Features of GIMAP Family Genes in Clear Cell Renal Cell Carcinoma: An Observational Study on CT Images

GTPases of immunity-associated proteins (GIMAP) genes include seven functional genes and a pseudogene. Most of the GIMAPs have a role in the maintenance and development of lymphocytes. GIMAPs could inhibit the development of tumors by increasing the amount and antitumor activity of infiltrating immunocytes. Knowledge of key factors that affect the tumor immune microenvironment for predicting the efficacy of immunotherapy and establishing new targets in ccRCC is of great importance. A computed tomography (CT)-based radiogenomic approach was used to detect the imaging phenotypic features of GIMAP family gene expression in ccRCC. In this retrospective study we enrolled 193 ccRCC patients divided into two groups: ccRCC patients with GIMAP expression (n = 52) and ccRCC patients without GIMAP expression (n = 141). Several imaging features were evaluated on preoperative CT scan. A statistically significant correlation was found with absence of endophytic growth pattern (p = 0.049), tumor infiltration (p = 0.005), advanced age (p = 0.018), and high Fuhrman grade (p = 0.024). This study demonstrates CT imaging features of GIMAP expression in ccRCC. These results could allow the collection of data on GIMAP expression through a CT-approach and could be used for the development of a targeted therapy.

Comparative Characterization of Immune Response in Sheep with Caseous Lymphadenitis through Analysis of the Whole Blood Transcriptome

Caseous lymphadenitis (CL) is a chronic contagious disease that affects small ruminants and is characterized by the formation of pyogranulomas in lymph nodes and other organs. However, the pathogenesis of this disease and the response of the host genome to infection are not yet fully understood. This study aimed to investigate the whole blood transcriptome and evaluate differential gene expression during the later stages of CL in naturally infected ewes. The study included diseased, serologically positive (EP), exposed, serologically negative (EN) ewes from the same infected flock and healthy ewes (CN) from a different flock. RNA sequencing was performed using the Illumina NextSeq system, and differential gene expression was estimated using DESeq2 and Edge R approaches. The analysis identified 191 annotated differentially expressed genes (DEGs) in the EP group (102 upregulated and 89 downregulated) and 256 DEGs in the EN group (106 upregulated and 150 downregulated) compared to the CN group. Numerous immunoregulatory interactions between lymphoid and nonlymphoid cells were influenced in both EP and EN ewes. Immune DEGs were preferentially assigned to antigen presentation through the MHC complex, T lymphocyte-mediated immunity, and extracellular matrix interactions. Furthermore, the EP group showed altered regulation of cytokine and chemokine signaling and activation and recombination of B-cell receptors. Conversely, NF-kappa B signaling, apoptosis, and stress response were the main processes influenced in the EN group. In addition, statistically significant enrichment of the essential immune pathways of binding and uptake of ligands by scavenger receptors in EP and p53 signaling in the EN group was found. In conclusion, this study provides new insights into the disease course and host-pathogen interaction in naturally CL-infected sheep by investigating the blood transcriptome.

Computational gene expression analysis reveals distinct molecular subgroups of T-cell prolymphocytic leukemia

T-cell prolymphocytic leukemia (T-PLL) is a rare blood cancer with poor prognosis. Overexpression of the proto-oncogene TCL1A and missense mutations of the tumor suppressor ATM are putative main drivers of T-PLL development, but so far only little is known about the existence of T-PLL gene expression subtypes. We performed an in-depth computational reanalysis of 68 gene expression profiles of one of the largest currently existing T-PLL patient cohorts. Hierarchical clustering combined with bootstrapping revealed three robust T-PLL gene expression subgroups. Additional comparative analyses revealed similarities and differences of these subgroups at the level of individual genes, signaling and metabolic pathways, and associated gene regulatory networks. Differences were mainly reflected at the transcriptomic level, whereas gene copy number profiles of the three subgroups were much more similar to each other, except for few characteristic differences like duplications of parts of the chromosomes 7, 8, 14, and 22. At the network level, most of the 41 predicted potential major regulators showed subgroup-specific expression levels that differed at least in comparison to one other subgroup. Functional annotations suggest that these regulators contribute to differences between the subgroups by altering processes like immune responses, angiogenesis, cellular respiration, cell proliferation, apoptosis, or migration. Most of these regulators are known from other cancers and several of them have been reported in relation to leukemia (e.g. AHSP, CXCL8, CXCR2, ELANE, FFAR2, G0S2, GIMAP2, IL1RN, LCN2, MBTD1, PPP1R15A). The existence of the three revealed T-PLL subgroups was further validated by a classification of T-PLL patients from two other smaller cohorts. Overall, our study contributes to an improved stratification of T-PLL and the observed subgroup-specific molecular characteristics could help to develop urgently needed targeted treatment strategies.

Large−Scale Profiling of Extracellular Vesicles Identified miR−625−5p as a Novel Biomarker of Immunotherapy Response in Advanced Non−Small−Cell Lung Cancer Patients

Simple Summary

Immune checkpoint inhibitors (ICIs) have revolutionized the treatment of advanced non-small-cell lung cancer (NSCLC) leading to substantial improvement in survival time and quality of life. Nevertheless, the clinical benefit of treatment is still limited to a minority of patients, reflecting the need to identify novel noninvasive biomarkers to improve patient selection. Currently available markers such as PD-L1 expression have important limitations. In this study, we focused on extracellular vesicles (EV)-associated miRNAs produced by cancer cells and their microenvironment that can be easily detected in blood. In particular, after a large-scale screening of 799 EV-miRNAs, we identified EV-miR-625-5p as a novel independent biomarker of response and survival in ICI-treated NSCLC patients, in particular in patients with PD-L1 expression ≥ 50%. EV-miR-625-5p integrated with PDL-1 test could allow the clinician to identify in advance patients that would benefit from ICIs.

Abstract

Immune checkpoint inhibitors (ICIs) are largely used in the treatment of patients with advanced non-small-cell lung cancer (NSCLC). Novel biomarkers that provide biological information that could be useful for clinical management are needed. In this respect, extracellular vesicles (EV)-associated microRNAs (miRNAs) that are the principal vehicle of intercellular communication may be important sources of biomarkers. We analyzed the levels of 799 EV-miRNAs in the pretreatment plasma of 88 advanced NSCLC patients who received anti-PD-1 therapy as single agent. After data normalization, we used a two-step approach to identify candidate biomarkers associated to both objective response (OR) by RECIST and longer overall survival (OS). Univariate and multivariate analyses including known clinicopathologic variables and new findings were performed. In our cohort, 24/88 (27.3%) patients showed OR by RECIST. Median OS in the whole cohort was 11.5 months. In total, 196 EV-miRNAs out 799 were selected as expressed above background. After multiplicity adjustment, abundance of EV-miR-625-5p was found to be correlated with PD-L1 expression and significantly associated to OR by RECIST (p = 0.0366) and OS (p = 0.0031). In multivariate analysis, PD-L1 staining and EV-miR-625-5p levels were constantly associated to OR and OS. Finally, we showed that EV-miR-625-5p levels could discriminate patients with longer survival, in particular in the class expressing PD-L1 ≥50%. EV-miRNAs represent a source of relevant biomarkers. EV-miR-625-5p is an independent biomarker of response and survival in ICI-treated NSCLC patients, in particular in patients with PD-L1 expression ≥50%.

STAT1 maintains naïve CD8+ T cell quiescence by suppressing the type I IFN-STAT4-mTORC1 signaling axis

Naïve CD8⁺ T cell quiescence is maintained at a steady state. Although this state of quiescence involves various cell-intrinsic and cell-extrinsic regulators, the mechanisms underlying this regulation remain incompletely understood. Here, we found that signal transducer and activator of transcription 1 (STAT1), a key transcription factor downstream of interferon receptor (IFNR) signaling, plays a cell-intrinsic role in maintaining naïve CD8⁺ T cell quiescence. STAT1-deficient mice showed enhanced proliferation of peripheral naïve CD8⁺ T cells, which resulted in an abnormal increase in the number of CD44^hi memory/activated phenotype cells and an enlargement of secondary lymphoid tissues. This phenomenon was not observed in IFNR-deficient mice but was paradoxically dependent on type I interferon and its alternative signaling pathway via the STAT4–RagD–lysosomal mTORC1 axis. Collectively, these findings underline the importance of STAT1 in regulating the homeostasis of peripheral naïve CD8⁺ T cells by suppressing their responsiveness to homeostatic cues at a steady state.

The GIMAP Family Proteins: An Incomplete Puzzle

Overview: Long-term survival of T lymphocytes in quiescent state is essential to maintain their cell numbers in secondary lymphoid organs and in peripheral circulation. In the BioBreeding diabetes-prone strain of rats (BB-DP), loss of functional GIMAP5 (GTPase of the immune associated nucleotide binding protein 5) results in profound peripheral T lymphopenia. This discovery heralded the identification of a new family of proteins initially called Immune-associated nucleotide binding protein (IAN) family. In this review we will use ‘GIMAP’ to refer to this family of proteins. Recent studies suggest that GIMAP proteins may interact with each other and also be involved in the movement of the cellular cargo along the cytoskeletal network. Here we will summarize the current knowledge on the characteristics and functions of GIMAP family of proteins.

Mining TCGA Data for Key Biomarkers Related to Immune Microenvironment in Endometrial cancer by Immune Score and Weighted Correlation Network Analysis

Background: Endometrial cancer (EC) is one of the most lethal gynecological cancers around the world. The aim of this study is to identify the potential immune microenvironment-related biomarkers associated with the prognosis for EC.

Methods: RNA-seq data and clinical information of EC patients were derived from The Cancer Genome Atlas (TCGA). The immune score of each EC sample was obtained by ESTIMATE algorithm. Weighted gene co-expression network analysis (WGCNA) was used to identify the interesting module and potential key genes concerning the immune score. The expression patterns of the key genes were then verified via the GEPIA database. Finally, CIBERSORT was applied to evaluate the relative abundances of 22 immune cell types in EC.

Results: Immune scores were significantly associated with tumor grade and histology of EC, and high immune scores may exert a protective influence on the survival outcome for EC. WGCNA indicated that the black module was significantly correlated with the immune score. Function analysis revealed it mainly involved in those terms related to immune regulation and inflammatory response. Moreover, 11 key genes (APOL3, C10orf54, CLEC2B, GIMAP1, GIMAP4, GIMAP6, GIMAP7, GIMAP8, GYPC, IFFO1, TAGAP) were identified from the black module, validated by the GEPIA database, and revealed strong correlations with infiltration levels of multiple immune cell types, as was the prognosis of EC.

Conclusion: In this study, 11 key genes showed abnormal expressions and strong correlations with immune infiltration in EC, most of which were significantly associated with the prognosis of EC. These findings made them promising therapeutic targets for the treatment of EC.

Systemic analyses of expression patterns and clinical features for GIMAPs family members in lung adenocarcinoma

GTPase of immunity-associated proteins (GIMAPs) are frequently prescribed as important components of immune regulation complexes, which were known to play key roles in lung adenocarcinoma. However, little is known about the function of distinct GIMAPs in lung adenocarcinoma. To address this issue, this study investigated the biological function and pathway of GIMAPs in lung adenocarcinoma using multiple public databases. Absent expression of GIMAPs was found in lung adenocarcinoma at mRNA and protein levels. While a purity-corrected value uncovered that all GIMAPs were positively associated with the immune infiltration of lung adenocarcinoma. Furthermore, the expressions of GIMAPs were considered to be negatively associated with clinical cancer stages, patient’s gender and pathological tumor grades in patients with lung adenocarcinoma. Besides, higher mRNA expression of GIMAPs was significantly associated with longer overall survival of patients with lung adenocarcinoma. Taken together, these results may enable GIMAPs family members as diagnostic and survival biomarker candidates or even potential therapeutic targets for patients with lung adenocarcinoma.

Heterogeneity of human bone marrow and blood natural killer cells defined by single-cell transcriptome

Natural killer (NK) cells are critical to both innate and adaptive immunity. However, the development and heterogeneity of human NK cells are yet to be fully defined. Using single-cell RNA-sequencing technology, here we identify distinct NK populations in human bone marrow and blood, including one population expressing higher levels of immediate early genes indicative of a homeostatic activation. Functionally matured NK cells with high expression of CX3CR1, HAVCR2 (TIM-3), and ZEB2 represents terminally differentiated status with the unique transcriptional profile. Transcriptomic and pseudotime analyses identify a transitional population between CD56^bright and CD56^dim NK cells. Finally, a donor with GATA2^T354M mutation exhibits reduced percentage of CD56^bright NK cells with altered transcriptome and elevated cell death. These data expand our understanding of the heterogeneity and development of human NK cells.

Natural killer (NK) cells are important innate immune cells with diverse functions. Here the authors use single-cell RNA-sequencing of purified human bone marrow and peripheral blood NK cells to define five populations of NK cells with distinct transcriptomic profile to further our understanding of NK development and heterogeneity.

GIMAP6 is required for T cell maintenance and efficient autophagy in mice

The GTPases of the immunity-associated proteins (GIMAP) GTPases are a family of proteins expressed strongly in the adaptive immune system. We have previously reported that in human cells one member of this family, GIMAP6, interacts with the ATG8 family member GABARAPL2, and is recruited to autophagosomes upon starvation, suggesting a role for GIMAP6 in the autophagic process. To study this possibility and the function of GIMAP6 in the immune system, we have established a mouse line in which the Gimap6 gene can be inactivated by Cre-mediated recombination. In mice bred to carry the CD2Cre transgene such that the Gimap6 gene was deleted within the T and B cell lineages there was a 50–70% reduction in peripheral CD4⁺ and CD8⁺ T cells. Analysis of splenocyte-derived proteins from these mice indicated increased levels of MAP1LC3B, particularly the lipidated LC3-II form, and S405-phosphorylation of SQSTM1. Electron microscopic measurements of Gimap6^-/- CD4⁺ T cells indicated an increased mitochondrial/cytoplasmic volume ratio and increased numbers of autophagosomes. These results are consistent with autophagic disruption in the cells. However, Gimap6^-/- T cells were largely normal in character, could be effectively activated in vitro and supported T cell-dependent antibody production. Treatment in vitro of CD4⁺ splenocytes from GIMAP6^fl/flERT2Cre mice with 4-hydroxytamoxifen resulted in the disappearance of GIMAP6 within five days. In parallel, increased phosphorylation of SQSTM1 and TBK1 was observed. These results indicate a requirement for GIMAP6 in the maintenance of a normal peripheral adaptive immune system and a significant role for the protein in normal autophagic processes. Moreover, as GIMAP6 is expressed in a cell-selective manner, this indicates the potential existence of a cell-restricted mode of autophagic regulation.

Functional and biochemical characterization of a T cell-associated anti-apoptotic protein, GIMAP6

GTPases of immunity-associated proteins (GIMAPs) are expressed in lymphocytes and regulate survival/death signaling and cell development within the immune system. We found that human GIMAP6 is expressed primarily in T cell lines. By sorting human peripheral blood mononuclear cells and performing quantitative RT-PCR, GIMAP6 was found to be expressed in CD3+ cells. In Jurkat cells that had been knocked down for GIMAP6, treatment with hydrogen peroxide, FasL, or okadaic acid significantly increased cell death/apoptosis. Exogenous expression of GMAP6 protected Huh-7 cells from apoptosis, suggesting that GIMAP6 is an anti-apoptotic protein. Furthermore, knockdown of GIMAP6 not only rendered Jurkat cells sensitive to apoptosis but also accelerated T cell activation under phorbol 12-myristate 13-acetate/ionomycin treatment conditions. Using this experimental system, we also observed a down-regulation of p65 phosphorylation (Ser-536) in GIMAP6 knockdown cells, indicating that GIMAP6 might display anti-apoptotic function through NF-κB activation. The conclusion from the study on cultured T cells was corroborated by the analysis of primary CD3+ T cells, showing that specific knockdown of GIMAP6 led to enhancement of phorbol 12-myristate 13-acetate/ionomycin-mediated activation signals. To characterize the biochemical properties of GIMAP6, we purified the recombinant GIMAP6 to homogeneity and revealed that GIMAP6 had ATPase as well as GTPase activity. We further demonstrated that the hydrolysis activity of GIMAP6 was not essential for its anti-apoptotic function in Huh-7 cells. Combining the expression data, biochemical properties, and cellular features, we conclude that GIMAP6 plays a role in modulating immune function and that it does this by controlling cell death and the activation of T cells.