Identification of key genes implicated in the suppressive function of human FOXP3+CD25+CD4+ regulatory T cells through the analysis of time-series data

Abstract

Human forkhead box P3 (FOXP3)+ cluster of differentiation (CD)25+CD4+ regulatory T cells (Tregs) are a type of T cell that express CD4, CD25 and FOXP3, which are critical for maintaining immune homeostasis. The present study aimed to determine the mechanisms underlying Treg function. The GSE11292 dataset was downloaded from the Gene Expression Omnibus, which included data from Treg cells at 19 time points (0-360 min) with an equal interval of 20 min, and corresponding repeated samples. However, data for Treg cells at time point 120 min were missing. Using the Mfuzz package, the key genes were identified by clustering analysis. Subsequently, regulatory networks and protein-protein interaction (PPI) networks were constructed and merged into integrated networks using Cytoscape software. Using Database for Annotation, Visualization and Integrated Discover software, enrichment analyses were performed for the genes involved in the PPI networks. Cluster 1 (including 292 genes), cluster 2 (including 111 genes), cluster 3 (including 194 genes) and cluster 4 (including 103 genes) were obtained from the clustering analysis. GAPDH (degree, 40) in cluster 1, Janus kinase 2 (JAK2) (degree, 10) and signal transducer and activator of transcription 5A (STAT5A) (degree, 9) in cluster 3, and tumor necrosis factor (TNF) (degree, 26) and interleukin 2 (IL2) (degree, 22) in cluster 4 had higher degrees in the PPI networks. In addition, it was indicated that several genes may have a role in Treg function by targeting other genes [e.g. microRNA (miR)-146b-3p→TNF, miR-146b-5p→TNF, miR-142-5p→TNF and tripartite motif containing 28 (TRIM28)→GAPDH]. Enrichment analyses indicated that IL2 and TNF were enriched in the immune response and T cell receptor signaling pathway. In conclusion, GAPDH targeted by TRIM28, TNF targeted by miR-146b-3p, miR-146b-5p and miR-142-5p, in addition to JAK2, IL2, and STAT5A may serve important roles in Treg function.