An NF-κB Transcription-Factor-Dependent Lineage-Specific Transcriptional Program Promotes Regulatory T Cell Identity and Function.

Both conventional T (Tconv) cells and regulatory T (Treg) cells are activated through ligation of the T cell receptor (TCR) complex, leading to the induction of the transcription factor NF-κB. In Tconv cells, NF-κB regulates expression of genes essential for T cell activation, proliferation, and function. However the role of NF-κB in Treg function remains unclear. We conditionally deleted canonical NF-κB members p65 and c-Rel in developing and mature Treg cells and found they have unique but partially redundant roles. c-Rel was critical for thymic Treg development while p65 was essential for mature Treg identity and maintenance of immune tolerance. Transcriptome and NF-κB p65 binding analyses demonstrated a lineage specific, NF-κB-dependent transcriptional program, enabled by enhanced chromatin accessibility. These dual roles of canonical NF-κB in Tconv and Treg cells highlight the functional plasticity of the NF-κB signaling pathway and underscores the need for more selective strategies to therapeutically target NF-κB.

A pan-cancer analysis implicates human NKIRAS1 as a tumor-suppressor gene.

The NF-κB family of transcription factors and the Ras family of small GTPases are important mediators of proproliferative signaling that drives tumorigenesis and carcinogenesis. The κB-Ras proteins were previously shown to inhibit both NF-κB and Ras activation through independent mechanisms, implicating them as tumor suppressors with potentially broad relevance to human cancers. In this study, we have used two mouse models to establish the relevance of the κB-Ras proteins for tumorigenesis. Additionally, we have utilized a pan-cancer bioinformatics analysis to explore the role of the κB-Ras proteins in human cancers. Surprisingly, we find that the genes encoding κB-Ras 1 (NKIRAS1) and κB-Ras 2 (NKIRAS2) are rarely down-regulated in tumor samples with oncogenic Ras mutations. Reduced expression of human NKIRAS1 alone is associated with worse prognosis in at least four cancer types and linked to a network of genes implicated in tumorigenesis. Our findings provide direct evidence that loss of NKIRAS1 in human tumors that do not carry oncogenic RAS mutations is associated with worse clinical outcomes.

FCRL3 expression is upregulated and closely correlates with TIGIT expression in regulatory T cells of patients with systemic lupus erythematosus.

Using data from single-cell RNA sequencing and flow cytometry, we initially examined the expression of FCRL3, finding it to be elevated and positively associated with TIGIT expression in the regulatory T cells of patients with systemic lupus erythematosus. This also suggests that the co-expression of FCRL3 and TIGIT warrants further attention.

G Protein-Coupled Receptor 56 Characterizes CTLs and Reflects the Progression of Lung Cancer Patients.

CTLs play important roles in host immune responses to tumors. CD4 CTLs are characterized by their ability to secrete cytotoxic effector molecules, such as granzyme B and perforin, and kill target cells in a MHC class II-restricted manner. However, the cell surface markers of CD4 CTLs remain unknown, which hinders their separation and research on their function. In this study, we performed a bioinformatics analysis and experimental validation that revealed that G protein-coupled receptor 56 (GPR56) is a cell surface marker that can be used to characterize CD4 CTLs. We found that GPR56 and granzyme B were coexpressed in extremely high levels in human peripheral blood T cells, and that anti-GPR56 stimulation significantly upregulated the expression of granzyme B in both CD4+GPR56+ and CD8+GPR56+ T cells. These findings suggest that GPR56 expression and the GPR56 signaling pathway could contribute directly to the toxic function of either CD4+ or CD8+ T cells. We also used GPR56 as a biomarker to investigate the clinical significance of CD4 CTLs. GPR56+ T cell levels were increased in patients with lung cancer, and GPR56 expression was significantly correlated with lung cancer progression. A further analysis revealed an increase in exhausted cell states in lung cancer patients because of upregulation of programmed cell death protein 1 expression in GPR56+ T cells. The findings of this study suggest that GPR56 characterizes the cytotoxic states of either CD4+ or CD8+ T cells.

Identification of rat Vstm1 with conservative anti-inflammatory effect between rat and human homologs.

Human VSTM1 (also known as SIRL1) is an inhibitory immune checkpoint receptor involved in leukocyte activation. Identification of the homologous genes in other species, such as mice and rats, will undoubtedly contribute to functional studies and clinical applications. Here, we successfully cloned the Vstm1 gene in rats, as supported by high-throughput sequencing data. However, Vstm1 is degenerated to a pseudogene in the mouse genome. Rat Vstm1 mRNA contains a complete open reading frame (ORF) of 630 nucleotides encoding 209 amino acids. Rat Vstm1 is highly expressed in bone marrow, especially in granulocytes. The expression levels of Vstm1 gradually increase with the development of granulocytes in bone marrow but are downregulated in response to inflammatory stimuli. Rat VSTM1 does not have an immunoreceptor tyrosine-based inhibitory motif (ITIM), however, it shows a conservative function of inflammatory inhibition with human VSTM1, and both are anti-correlated with many inflammatory cytokines, such as IL-1α and TNF-α. In bone marrow-derived macrophages (BMDMs), either rat or human VSTM1 suppressed the secretion of inflammatory cytokines in response to LPS stimulation. Further analysis in lung cancer microenvironment revealed that VSTM1 is mainly expressed in myeloid cells, anti-correlated with inflammatory cytokines and associated with tumor development and metastasis.

VNN2-expressing circulating monocytes exhibit unique functional characteristics and are decreased in patients with primary Sjögren's syndrome.

OBJECTIVE: This study aims to elucidate the significance of VNN2 expression in peripheral blood monocytes and its clinical relevance in primary Sjögren's syndrome (pSS). METHODS: We investigated VNN2 expression by analyzing single-cell RNA sequencing (scRNA-seq) data from peripheral blood mononuclear cells. Flow cytometry was used to detect and compare VNN2 expression in total monocytes, classical monocytes (cMo), intermediate monocytes (iMo) and non-classical monocytes (ncMo). Additionally, we examined the expression of HLA, ICAM1, CD62L, ITGAM, S100A8, S100A9, CCR2, CCR6, CX3CR1 and CXCR3 in VNN2+ and VNN2- cells. We analyzed the correlation between VNN2 expression and clinical indicators and assessed the clinical utility of VNN2+ monocytes in pSS diagnosis using receiver operating characteristic curves. RESULTS: We observed high VNN2 expression in monocytes, with significantly higher levels in CD14++ monocytes compared to ncMo. VNN2+ monocytes exhibited decreased expression of HLA and CD62L and increased expression of ICAM1, ITGAM, S100A8, S100A9, CCR2, CCR6, CX3CR1 and CXCR3 compared to VNN2- monocytes. Although scRNA-seq data showed that VNN2 mRNA was upregulated, cell surface expression of VNN2 was decreased in monocytes from pSS patients compared to healthy controls. The reduced levels of VNN2+ monocyte subpopulations in pSS patients were negatively correlated with anti-ribosome antibody levels and positively correlated with complement 4 levels. Detection of VNN2 expression in monocytes can aid in the auxiliary diagnosis of pSS. CONCLUSION: Monocytes expressing cell surface VNN2 are significantly reduced in pSS patients. This suggests a potential role for VNN2 in pSS development and its potential use as a diagnostic marker for pSS.

Evaluating imputation methods for single-cell RNA-seq data.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) enables the high-throughput profiling of gene expression at the single-cell level. However, overwhelming dropouts within data may obscure meaningful biological signals. Various imputation methods have recently been developed to address this problem. Therefore, it is important to perform a systematic evaluation of different imputation algorithms. RESULTS: In this study, we evaluated 11 of the most recent imputation methods on 12 real biological datasets from immunological studies and 4 simulated datasets. The performance of these methods was compared, based on numerical recovery, cell clustering and marker gene analysis. Most of the methods brought some benefits on numerical recovery. To some extent, the performance of imputation methods varied among protocols. In the cell clustering analysis, no method performed consistently well across all datasets. Some methods performed poorly on real datasets but excellent on simulated datasets. Surprisingly and importantly, some methods had a negative effect on cell clustering. In marker gene analysis, some methods identified potentially novel cell subsets. However, not all of the marker genes were successfully imputed in gene expression, suggesting that imputation challenges remain. CONCLUSIONS: In summary, different imputation methods showed different effects on different datasets, suggesting that imputation may have dataset specificity. Our study reveals the benefits and limitations of various imputation methods and provides a data-driven guidance for scRNA-seq data analysis.

Omic horizon expression: a database of gene expression based on RNA sequencing data.

BACKGROUND: Gene expression profiles have important significance for gene expression characteristics and further functional studies. More attention has been given to the expression databases in humans and mice, but less attention has been given to rats, while rat models also play an irreplaceable role in biomedical experiments. RESULTS: To depict the rat gene expression profiles in mRNA expression levels, we analyzed over 2,700 RNA sequencing (RNA-Seq) samples from 48 tissues, 40 primary cell types and 25 cell lines; and then mapped them to the latest version of the rat genome reference, mRatBN7.2. Based on these datasets and reanalysis, we constructed a new database, the Omic Horizon Expression Database ( http://immudb.bjmu.edu.cn/expression.html ), which allows expressional profile query of over 25,000 rat genes based on non-redundant gene symbols. The database supports requests using gene symbols (or alias), Ensemble and Entrez gene IDs. Gene expression profiles can be queried in three categories: tissues, primary cells and cell lines. Application examples including expression profiling and comparison, as well as identification of novel rat genes, were illustrated to show the utility of the database. CONCLUSIONS: As an omic resource, the Omic Horizon Expression Database provides horizons of gene expression profiles across various tissues and cells, which greatly facilitates the identification of rat genes as well as functional clues.

PDK1 Is Required for Maintenance of CD4+ Foxp3+ Regulatory T Cell Function.

Regulatory T (Treg) cells have an essential role in maintaining immune homeostasis, in part by suppressing effector T cell functions. Phosphoinositide-dependent kinase 1 (PDK1) is a pleiotropic kinase that acts as a key effector downstream of PI3K in many cell types. In T cells, PDK1 has been shown to be critical for activation of NF-κB and AKT signaling upon TCR ligation and is therefore essential for effector T cell activation, proliferation, and cytokine production. Using Treg cell-specific conditional deletion, we now demonstrate that PDK1 is also essential for Treg cell suppressive activity in vivo. Ablation of Pdk1 specifically in Treg cells led to systemic, lethal, scurfy-like inflammation in mice. Genome-wide analysis confirmed that PDK1 is essential for the regulation of key Treg cell signature gene expression and, further, suggested that PDK1 acts primarily to control Treg cell gene expression through regulation of the canonical NF-κB pathway. Consistent with these results, the scurfy-like phenotype of mice lacking PDK1 in Treg cells was rescued by enforced activation of NF-κB downstream of PDK1. Therefore, PDK1-mediated activation of the NF-κB signaling pathway is essential for regulation of Treg cell signature gene expression and suppressor function.

FAM19A5/TAFA5, a novel neurokine, plays a crucial role in depressive-like and spatial memory-related behaviors in mice.

FAM19A5/TAFA5 is a member of the family with sequence similarity 19 with unknown function in emotional and cognitive regulation. Here, we reported that FAM19A5 was highly expressed in the embryonic and postnatal mouse brain, especially in the hippocampus. Behaviorally, genetic deletion of Fam19a5 resulted in increased depressive-like behaviors and impaired hippocampus-dependent spatial memory. These behavioral alterations were associated with the decreased expression of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors and N-methyl-D-aspartic acid receptors, as well as significantly reduced glutamate release and neuronal activity in the hippocampus. Subsequently, these changes led to the decreased density of dendritic spines. In recent years, the roles of chronic stress participating in the development of depression have become increasingly clear, but the mechanism remains to be elucidated. We found that the levels of FAM19A5 in plasma and hippocampus of chronic stress-treated mice were significantly decreased whereas overexpression of human FAM19A5 selectively in the hippocampus could attenuate chronic stress-induced depressive-like behaviors. Taken together, our results revealed for the first time that FAM19A5 plays a key role in the regulation of depression and spatial cognition in the hippocampus. Furthermore, our study provided a new mechanism for chronic stress-induced depression, and also provided a potential biomarker for the diagnosis and a new strategy for the treatment of depression.

CMTM7 plays key roles in TLR-induced plasma cell differentiation and p38 activation in murine B-1 B cells.

Terminal differentiation of B cells into antibody-secreting cells is the foundation of humoral immune response. B-1 cells, which are different from B-2 cells, preferentially differentiate into plasma cells. CMTM7 is a MARVEL-domain-containing membrane protein predominantly expressed in B cells that plays an important role in B-1a cell development. The present study assessed CMTM7 function in response to antigen stimulation. Following immunization with T cell-dependent and T cell-independent antigens, Cmtm7-deficient mice exhibited decreased IgM but normal IgG responses in vivo. In vitro stimulation with LPSs induced Cmtm7-/- B-1 cell activation, whereas proliferation was marginally reduced. Notably, Cmtm7 deficiency markedly suppressed plasma cell differentiation in response to TLR agonists, accompanied by a decrease in IgM and IL-10 production. At the molecular level, loss of Cmtm7 repressed the downregulation of Pax5 and the upregulation of Xbp1, Irf4, and Prdm1. Furthermore, p38 phosphorylation was inhibited in Cmtm7-/- B-1 cells. Experiments using a p38 inhibitor revealed that p38 activation was essential for the terminal differentiation of B-1 cells, suggesting that Cmtm7 contributes to B-1 cell differentiation by maintaining p38 activation. Overall, the data reveal the crucial functions of CMTM7 in TLR-induced terminal differentiation and p38 activation in B-1 cells.

Essential role for Cmtm7 in cell-surface phenotype, BCR signaling, survival and Igµ repertoire of splenic B-1a cells.

B-1a cells represent a distinct B cell population with unique phenotype, self-renewing capacity and restricted Igµ repertoire. They primarily locate in body cavity and also exist in spleen. The different subpopulations of B-1a cells are heavily affected by local environment. Our previous studies revealed that MARVEL-domain-containing membrane protein, CMTM7, was involved in B-1a cell development. Here, we focused its influence on peritoneal and splenic B-1a cells. Unlike peritoneal B-1a cells, we found that splenic Cmtm7-/- B-1a cells expressed higher level of CD5, CD80 and CD86 compared with WT counterparts. They also exhibited an enhanced tonic BCR signals in steady state. Though the cell viability was unaffected in vitro, Cmtm7 knockout markedly promoted splenic B-1a cell apoptosis in situ, which was likely associated with down-regulation of Il-5rα. With regard to Igµ repertoire, peritoneal and splenic Cmtm7-/- B-1a cells exhibit similar changes exemplified by the loss of VH11 and gain of VH12, whereas an increase in VH1 usage and skewed J segments from JH1 to JH2 and JH4 families could only be detected within splenic Cmtm7-/- B-1a cells. Overall, these data indicate that Cmtm7 functions differently in peritoneal and splenic B-1a cells and plays a more important role in splenic cells.

Cmtm7 knockout inhibits B-1a cell development at the transitional (TrB-1a) stage.

Innate-like B-1a cells are an important cell population for production of natural IgM and interleukin-10 (IL-10), and act as the first line against pathogens. We determined that CMTM7 is essential for B-1a cell development. Following Cmtm7 (CKLF-like MARVEL transmembrane domain-containing 7) knockout, B-1a cell numbers decreased markedly in all investigated tissues. Using a bone marrow and fetal liver adoptive transfer model and conditional knockout mice, we showed that the reduction of B-1a cells resulted from B-cell-intrinsic defects. Because of B-1a cell loss, Cmtm7-deficient mice produced less IgM and IL-10, and were more susceptible to microbial sepsis. Self-renewal and homeostasis of mature B-1a cells in Cmtm7-/- mice were not impaired, suggesting the effect of Cmtm7 on B-1a cell development. Further investigations demonstrated that the function of Cmtm7 in B-1a cell development occurred at the specific transitional B-1a (TrB-1a) stage. Cmtm7 deficiency resulted in a slow proliferation and high cell death rate of TrB-1a cells. Thus, Cmtm7 controls B-1a cell development at the transitional stage.

Novel Adipokine, FAM19A5, Inhibits Neointima Formation After Injury Through Sphingosine-1-Phosphate Receptor 2.

BACKGROUND: Obesity plays crucial roles in the development of cardiovascular diseases. However, the mechanisms that link obesity and cardiovascular diseases remain elusive. Compelling evidence indicates that adipokines play an important role in obesity-related cardiovascular diseases. Here, we found a new adipokine-named family with sequence similarity 19, member A5 (FAM19A5), a protein with unknown function that was predicted to be distantly related to the CC-chemokine family. We aimed to test whether adipose-derived FAM19A5 regulates vascular pathology on injury. METHODS: DNA cloning, protein expression, purification, and N-terminal sequencing were applied to characterize FAM19A5. Adenovirus infection and siRNA transfection were performed to regulate FAM19A5 expression. Balloon and wire injury were performed in vivo on the rat carotid arteries and mouse femoral arteries, respectively. Bioinformatics analysis, radioactive ligand-receptor binding assays, receptor internalization, and calcium mobilization assays were used to identify the functional receptor for FAM19A5. RESULTS: We first characterized FAM19A5 as a secreted protein, and the first 43 N-terminal amino acids were the signal peptides. Both FAM19A5 mRNA and protein were abundantly expressed in the adipose tissue but were downregulated in obese mice. Overexpression of FAM19A5 markedly inhibited vascular smooth muscle cell proliferation and migration and neointima formation in the carotid arteries of balloon-injured rats. Accordingly, FAM19A5 silencing in adipocytes significantly promoted vascular smooth muscle cell activation. Adipose-specific FAM19A5 transgenic mice showed greater attenuation of neointima formation compared with wild-type littermates fed with or without Western-style diet. We further revealed that sphingosine-1-phosphate receptor 2 was the functional receptor for FAM19A5, with a dissociation constant (Kd) of 0.634 nmol/L. Inhibition of sphingosine-1-phosphate receptor 2 or its downstream G12/13-RhoA signaling circumvented the suppressive effects of FAM19A5 on vascular smooth muscle cell proliferation and migration. CONCLUSIONS: We revealed that a novel adipokine, FAM19A5, was capable of inhibiting postinjury neointima formation via sphingosine-1-phosphate receptor 2-G12/13-RhoA signaling. Downregulation of FAM19A5 during obesity may trigger cardiometabolic diseases.

FAM19A1 is a new ligand for GPR1 that modulates neural stem-cell proliferation and differentiation.

FAM19A1 is a member of the family with sequence similarity 19 with unknown function. FAM19A1 mRNA expression is restricted to the CNS. Here, we report that FAM19A1 is a classic secretory protein, and expression levels correlate with brain development, increasing from embryonic d 12.5, peaking between postnatal d (P)1 and P7 and decreasing at wk 8. The adult hippocampus is a region of FAM19A1 high expression. Recombinant FAM19A1 suppressed the proliferation and self-renewal of neural stem cells (NSCs) and altered the lineage progression of NSCs with promoted neuron differentiation and suppressed astrocyte differentiation. Although GPCR 1 (GPR1) has been reported to be expressed in the CNS, its functions in the brain remain unclear. We identified GPR1 to be a functional receptor for FAM19A1. FAM19A1 interacted with GPR1 via the N-terminal domain (GPR1-ND), and its NSC modulatory functions required the Rho-associated protein kinase (ROCK) /ERK1/2 and ROCK/signal transducer and activator of transcription 3 signaling pathways. GPR1-ND that selectively bound to FAM19A1 neutralized the effects of FAM19A1 on NSC functions. Taken together, our results show, for the first time to our knowledge, that FAM19A1 is a novel regulatory factor of the proliferation and differentiation of NSCs, and identified a novel mechanism by which GPCR mediates the effects of FAM19A1 on NSC functions that may be important for brain development and neurogenesis. Additional exploration of the functions of FAM19A1 and GPR1 in the CNS may broaden the range of therapeutic options available for major brain disorders.-Zheng, C., Chen, D., Zhang, Y., Bai, Y., Huang, S., Zheng, D., Liang, W., She, S., Peng, X., Wang, P., Mo, X., Song, Q., Lv, P., Huang, J., Ye, R. D., Wang, Y. FAM19A1 is a new ligand for GPR1 that modulates neural stem-cell proliferation and differentiation.

Virtual Sorting Has a Distinctive Advantage in Identification of Anticorrelated Genes and Further Negative Regulators of Immune Cell Subpopulations.

Immune cells are highly plastic in both gene expression and cell phenotype. We have established a method of gene expressional plasticity and virtual sorting to evaluate immune cell subpopulations and their characteristic genes in human CD4+ T cells. In this study, we continued to investigate the informatics mechanism on the effectiveness of virtual sorting. We found that virtual sorting had an overall positive correlation to the Pearson correlation in the identification of positively correlated genes. However, owing to nonlinear biological anticorrelation, virtual sorting showed a distinctive advantage for anticorrelated genes, suggesting an important role in the identification of negative regulators. In addition, based on virtual sorting results, we identified two basic gene sets among highly plastic genes, i.e., highly plastic cell cycle-associated molecules and highly plastic immune and defense response-associated molecules. Genes within each set tended to be positively connected, but genes between two sets were often anticorrelated. Further analysis revealed preferential transcription factor binding motifs existed between highly plastic cell cycle-associated molecules and highly plastic immune and defense response-associated molecules. Our results strongly suggested predetermined regulation, which was called an immune cell internal phenotype, should exist and could be mined by virtual sorting analysis. This provided efficient functional clues to study immune cell phenotypes and their regulation. Moreover, the current substantial virtual sorting results in both CD4+ T cells and B cells provide a useful resource for big-data-driven experimental studies and knowledge discoveries.

CMTM4 inhibits cell proliferation and migration via AKT, ERK1/2, and STAT3 pathway in colorectal cancer.

CMTM4 (CKLF-like MARVEL transmembrane domain containing 4), a potential tumor suppressor gene, is involved in several types of malignancies. It has been reported to be downregulated and exhibit anti-tumorigenic activities by regulating cell growth and cell cycle in clear cell renal cell carcinoma. It has also been identified as a tumor suppressor in hepatocellular carcinoma (HCC), and its negative expression is a risk factor for poor prognosis of HCC patients. In the present study, an integrated bioinformatics analysis based on The Cancer Genome Atlas (TCGA) database showed that CMTM4 was frequently reduced in colorectal cancer (CRC) and high expression of CMTM4 was associated with increased overall survival rates. Based on these findings, we adopted gain-of-function and lost-of-function strategies using SW480 and HT29 CRC cell lines which have relatively low and high endogenous CMTM4 levels, respectively. We observed impeded cell proliferation and migration upon overexpression of CMTM4 in SW480 cells, and the opposite effects were observed upon knockdown of CMTM4 in HT-29 cells. Cell signaling pathways essential for CRC progression were then examined, and the phosphorylation levels of AKT, ERK1/2, and STAT3 were found to be decreased by CMTM4 overexpression in SW480 cells and elevated by CMTM4 silencing in HT29 cells. Their inhibitors were used to validate that the three signaling pathways contributed to the inhibitory effects of CMTM4 on CRC cells. Taken together, our results suggest that CMTM4 plays a tumor suppressive role in CRC.

Identification of FAM3D as a new endogenous chemotaxis agonist for the formyl peptide receptors.

The family with sequence similarity 3 (FAM3) gene family is a cytokine-like gene family with four members FAM3A, FAM3B, FAM3C and FAM3D. In this study, we found that FAM3D strongly chemoattracted human peripheral blood neutrophils and monocytes. To identify the FAM3D receptor, we used chemotaxis, receptor internalization, Ca(2+) flux and radioligand-binding assays in FAM3D-stimulated HEK293 cells that transiently expressed formyl peptide receptor (FPR)1 or FPR2 to show that FAM3D was a high affinity ligand of these receptors, both of which were highly expressed on the surface of neutrophils, and monocytes and macrophages. After being injected into the mouse peritoneal cavity, FAM3D chemoattracted CD11b+ Ly6G+ neutrophils in a short time. In response to FAM3D stimulation, phosphorylated ERK1/2 and phosphorylated p38 MAPK family proteins were upregulated in the mouse neutrophils, and this increase was inhibited upon treatment with an inhibitor of FPR1 or FPR2. FAM3D has been reported to be constitutively expressed in the gastrointestinal tract. We found that FAM3D expression increased significantly during colitis induced by dextran sulfate sodium. Taken together, we propose that FAM3D plays a role in gastrointestinal homeostasis and inflammation through its receptors FPR1 and FPR2.

Electronic Sorting of Immune Cell Subpopulations Based on Highly Plastic Genes.

Immune cells are highly heterogeneous and plastic with regard to gene expression and cell phenotype. In this study, we categorized genes into those with low and high gene plasticity, and those categories revealed different functions and applications. We proposed that highly plastic genes could be suited for the labeling of immune cell subpopulations; thus, novel immune cell subpopulations could be identified by gene plasticity analysis. For this purpose, we systematically analyzed highly plastic genes in human and mouse immune cells. In total, 1,379 human and 883 mouse genes were identified as being extremely plastic. We also expanded our previous immunoinformatic method, electronic sorting, which surveys big data to perform virtual analysis. This approach used correlation analysis and took dosage changes into account, which allowed us to identify the differentially expressed genes. A test with human CD4(+) T cells supported the method's feasibility, effectiveness, and predictability. For example, with the use of human nonregulatory T cells, we found that FOXP3(hi)CD4(+) T cells were highly expressive of certain known molecules, such as CD25 and CTLA4, and that this process of investigation did not require isolating or inducing these immune cells in vitro. Therefore, the sorting process helped us to discover the potential signature genes or marker molecules and to conduct functional evaluations for immune cell subpopulations. Finally, in human CD4(+) T cells, 747 potential immune cell subpopulations and their candidate signature genes were identified, which provides a useful resource for big data-driven knowledge discoveries.

Cytokine-like 1 Chemoattracts Monocytes/Macrophages via CCR2.

Cytokine-like 1 (CYTL1) is a novel potential cytokine that was first identified in CD34(+) cells derived from bone marrow and cord blood, and it was also found using our immunogenomics strategy. The immunobiological functions of CYTL1 remain largely unknown, and its potential receptor(s) has not been identified. A previous proposed hypothesis suggested that CYTL1 had structural similarities with CCL2 and that CCR2 was a potential receptor of CYTL1. In this study, we verify that CYTL1 possesses chemotactic activity and demonstrate that its functional receptor is CCR2B using a series of experiments performed in HEK293 cells expressing CCR2B or CCR2B-EGFP, including chemotaxis, receptor internalization, and radioactive binding assays. CYTL1 chemoattracts human monocytes but not PBLs, and its chemotactic activity toward monocytes is dependent on the CCR2B-ERK pathway. Furthermore, both human and mouse recombinant CYTL1 protein have chemotactic effects on macrophages from wild-type mice but not from Ccr2(-/-) mice. Furthermore, the chemotactic activity of CYTL1 is sensitive to pertussis toxin. All of the above data confirm that CCR2B is a functional receptor of CYTL1.