Homeobox NKX2-3 promotes marginal-zone lymphomagenesis by activating B-cell receptor signalling and shaping lymphocyte dynamics

Quantitative Analysis of NKX2-3 Expression in Human Colon: An Immunohistochemical Study

In mice, Nkx2-3 homeodomain transcription factor defines the vascular specification of secondary and tertiary lymphoid tissues of the intestines. In human studies, polymorphisms in NKX2-3 have been identified as a susceptibility factor in inflammatory bowel diseases, whereas in mice, its absence is associated with protection against experimental colitis and enhanced intestinal epithelial proliferation. Here, we investigated the expression of NKX2-3 in normal, polyp, and adenocarcinoma human colon samples using immunohistochemistry and quantitative morphometry, correlating its expression with endothelial and mesenchymal stromal markers. Our results revealed that the expression of NKX2-3 is regionally confined to the lamina propria and lamina muscularis mucosae, and its production is restricted mostly to endothelial cells and smooth muscle cells with variable co-expression of CD34, alpha smooth muscle antigen (αSMA), and vascular adhesion protein-1 (VAP-1). The frequency of NKX2-3-positive cells and intensity of expression correlated inversely with aging. Furthermore, in most colorectal carcinoma samples, we observed a significant reduction of NKX2-3 expression. These findings indicate that the NKX2-3 transcription factor is produced by both endothelial and non-endothelial tissue constituents in the colon, and its expression changes during aging and in colorectal malignancies. (J Histochem Cytochem XX: XXX-XXX, XXXX).

Transcription factor Nkx2-3 maintains the self-renewal of hematopoietic stem cells by regulating mitophagy

Hematopoietic stem cells (HSCs) reside at the top of the hematopoietic hierarchy, exhibiting a unique capacity to self-renew and differentiate into all blood cells throughout the lifetime. However, how to prevent HSC exhaustion during long-term hematopoietic output is not fully understood. Here, we show that the homeobox transcription factor Nkx2-3 is required for HSC self-renewal by preserving metabolic fitness. We found that Nkx2-3 is preferentially expressed in HSCs with excessive regenerative potential. Mice with conditional deletion of Nkx2-3 displayed a reduced HSC pool and long-term repopulating capacity as well as increased sensitivity to irradiation and 5-flurouracil treatment due to impaired HSC quiescence. In contrast, overexpression of Nkx2-3 improved HSC function both in vitro and in vivo. Furthermore, mechanistic studies revealed that Nkx2-3 can directly control the transcription of the critical mitophagy regulator ULK1, which is essential for sustaining metabolic homeostasis in HSCs by clearing activated mitochondria. More importantly, a similar regulatory role of NKX2-3 was observed in human cord blood-derived HSCs. In conclusion, our data demonstrate an important role of the Nkx2-3/ULK1/mitophagy axis in regulating the self-renewal of HSCs, therefore providing a promising strategy to improve the function of HSCs in the clinic.

Tumor-suppressive role of the musculoaponeurotic fibrosarcoma gene in colorectal cancer

Somatic cell reprogramming using the microRNAs miR-200c, miR-302s, and miR-369s leads to increased expression of cyclin-dependent kinase inhibitors in human colorectal cancer (CRC) cells and suppressed tumor growth. Here, we investigated whether these microRNAs inhibit colorectal tumorigenesis in CPC;Apc mice, which are prone to colon and rectal polyps. Repeated administration of microRNAs inhibited polyp formation. Microarray analysis indicated that c-MAF, which reportedly shows oncogene-like behavior in multiple myeloma and T cell lymphoma, decreased in tumor samples but increased in microRNA-treated normal mucosa. Immunohistochemistry identified downregulation of c-MAF as an early tumorigenesis event in CRC, with low c-MAF expression associated with poor prognosis. Of note, c-MAF expression and p53 protein levels were inversely correlated in CRC samples. c-MAF knockout led to enhanced tumor formation in azoxymethane/dextran sodium sulfate-treated mice, with activation of cancer-promoting genes. c-MAF may play a tumor-suppressive role in CRC development.

Transcription factor TFII-I fine tunes innate properties of B lymphocytes

The ubiquitously expressed transcription factor TFII-I is a multifunctional protein with pleiotropic roles in gene regulation. TFII-I associated polymorphisms are implicated in Sjögren's syndrome and Lupus in humans and, germline deletion of the Gtf2i gene in mice leads to embryonic lethality. Here we report a unique role for TFII-I in homeostasis of innate properties of B lymphocytes. Loss of Gtf2i in murine B lineage cells leads to an alteration in transcriptome, chromatin landscape and associated transcription factor binding sites, which exhibits myeloid-like features and coincides with enhanced sensitivity to LPS induced gene expression. TFII-I deficient B cells also show increased switching to IgG3, a phenotype associated with inflammation. These results demonstrate a role for TFII-I in maintaining immune homeostasis and provide clues for GTF2I polymorphisms associated with B cell dominated autoimmune diseases in humans.

Preneoplastic somatic mutations including MYD88L265P in lymphoplasmacytic lymphoma

Normal cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B lymphocytes before lymphoma remains uncertain. We sequenced multiple stages of the B lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of the high prevalence of mutated MYD88. We observed similar accumulation of random mutations in B lineages from both cohorts and unexpectedly found MYD88^L265P in normal precursor and mature B lymphocytes from patients with lymphoma. We uncovered genetic and transcriptional pathways driving malignant transformation and leveraged these to model lymphoplasmacytic lymphoma in mice, based on mutated MYD88 in B cell precursors and BCL2 overexpression. Thus, MYD88^L265P is a preneoplastic event, which challenges the current understanding of lymphomagenesis and may have implications for early detection of B cell lymphomas.

Transcriptomic analysis of spleen B cell revealed the molecular basis of bursopentin on B cell differentiation

The bursa of Fabricius, the acknowledged humoral immune organ unique to birds, plays a vital role in B cell development. Bursopentin (BP5) derived from the bursa is reported to induce the development and formation of B cells. However, the mechanism of BP5 on B cell differentiation is still unclear. In this paper, total B lymphocytes from mice immunized with H9N2 subtype AIV vaccine were stimulated with BP5. The results show that BP5 at the experimental dosages promoted B cell differentiation, including the total B cells, activated B cells, differentiated B cells, mature B cells and plasma cells. Then, the in vivo immune experiment proved that the percentages of activated and differentiated B cells from mice immunized with AIV vaccine and 0.25 mg/mL BP5 were increased. To investigate the molecular mechanism of BP5 on B cell differentiation, the gene expression profiles of B cells purified from the spleen cells of mice immunized with AIV vaccine and BP5 were detected following RNA sequencing technology. The results show that BP5 at 0.05 and 0.25 mg/mL induced the enrichment of various biological functions, and stimulated five common significant enrichment pathways in B cells from the immunized mice. Additionally, 120 and 59 differentially expressed genes (DEG) represented transcriptional factors in B cells following 0.05 and 0.25 mg/mL BP5 immunization, respectively. In summary, these results suggest that BP5 regulates various gene expression involved in regulation of B cell development, which provides the knowledge required for additional studies on B cell differentiation in response to bursal-derived peptides and also provides an important experimental basis for improving vaccine immunity.

Single-cell profiling of infiltrating B cells and tertiary lymphoid structures in the TME of gastric adenocarcinomas

The occurrence and development of gastric adenocarcinoma (gADC) is closely related to the interaction between tumor cells and immune cells in the tumor microenvironment (TME). Our objective was to characterize the repertoire of immune cells in the TME of gADC. To analyze the transcriptomic, immune, and spatial information of TME in gADC, we constructed single-cell RNA sequencing, 10 × Genomics V(D)J analysis, multiple immunofluorescence techniques, and OSCmap analysis of 49,765 single cells in seven samples from four gADC patients. Our integrative analysis of B cells demonstrated that a large number of mucosal associated lymphoid tissue (MALT)-B cells were detected in the gADC tissues, which have mature tertiary lymphatic structures (mTLSs), and almost no MALT-B cells in peripheral blood sample. Moreover, MALT-B cells are a class of IgA⁺ plasma cells, which are characterized with high expression of complement pathway activation-related genes. Next, natural killer T (NKT) cells mainly exist in gADC tissues accompanied by mTLSs. This study also classified monocytes/macrophages and epithelial cells into benign and malignant types. Interestingly, CSOmap (q < .05) and multiple immunofluorescence (p < .05) results indicated more types of immune cells can be enriched in tissues with mTLSs than normal TLSs, and the density of mTLSs were higher than normal TLSs. Our findings provide novel insights for the signature of immune cells and tumor cells in the TME of gADC with TLSs and highlight the potential importance of IgA-mediated humoral immunity in gADC patients with TLSs.

An aberrant DNA methylation signature for predicting the prognosis of head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a common malignancy worldwide with a poor prognosis. DNA methylation is an epigenetic modification that plays a critical role in the etiology and pathogenesis of HNSCC. The current study aimed to develop a predictive methylation signature based on bioinformatics analysis to improve the prognosis and optimize therapeutic outcome in HNSCC. Clinical information and methylation sequencing data of patients with HNSCC were downloaded from The Cancer Genome Atlas database. The R package was used to identify differentially methylated genes (DMGs) between HNSCC and adjacent normal tissues. We identified 22 DMGs associated with 246 differentially methylated sites. Patients with HNSCC were classified into training and test groups. Cox regression analysis was used to build a risk score formula based on the five methylation sites (cg26428455, cg13754259, cg17421709, cg19229344, and cg11668749) in the training group. The Kaplan–Meier survival curves showed that the overall survival (OS) rates were significantly different between the high‐ and low‐risk groups sorted by the signature in the training group (median: 1.38 vs. 1.57 years, log‐rank test, p < 0.001). The predictive power was then validated in the test group (median: 1.34 vs. 1.75 years, log‐rank test, p < 0.001). The area under the receiver operating characteristic curve (area under the curve) based on the signature for predicting the 5‐year survival rates, was 0.7 in the training and 0.73 in test groups, respectively. The results of multivariate Cox regression analysis showed that the riskscore (RS) signature based on the five methylation sites was an independent prognostic tool for OS prediction in patients. In addition, a predictive nomogram model that incorporated the RS signature and patient clinical information was developed. The innovative methylation signature‐based model developed in our study represents a robust prognostic tool for guiding clinical therapy and predicting the OS in patients with HNSCC.

Novel Autophagy-Related Gene Signature Investigation for Patients With Oral Squamous Cell Carcinoma

The correlation between autophagy defects and oral squamous cell carcinoma (OSCC) has been previously studied, but only based on a limited number of autophagy-related genes in cell lines or animal models. The aim of the present study was to analyze differentially expressed autophagy-related genes through The Cancer Genome Atlas (TCGA) database to explore enriched pathways and potential biological function. Based on TCGA database, a signature composed of four autophagy-related genes (CDKN2A, NKX2-3, NRG3, and FADD) was established by using multivariate Cox regression models and two Gene Expression Omnibus datasets were applied for external validation. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to study the function of autophagy-related genes and their pathways. The most significant GO and KEGG pathways were enriched in several key pathways that were related to the progression of autophagy and OSCC. Furthermore, a prognostic risk score was constructed based on the four genes; patients were then divided into two groups (i.e., high risk and low risk) in terms of the median of risk score. Prognosis of the two groups and results showed that patients at the low-risk group had a much better prognosis than those at the high-risk group, regardless of whether they were in the training datasets or validation datasets. Multivariate Cox regression results indicated that the risk score of the autophagy-related gene signatures could greatly predict the prognosis of patients after controlling for several clinical covariates. The findings of the present study revealed that autophagy-related gene signatures play an important role in OSCC and are potential prognostic biomarkers and therapeutic targets.

Histone acetylation dynamics modulates chromatin conformation and allele-specific interactions at oncogenic loci

In cancer cells, enhancer hijacking mediated by chromosomal alterations and/or increased deposition of acetylated histone H3 lysine 27 (H3K27ac) can support oncogene expression. However, how the chromatin conformation of enhancer-promoter interactions is affected by these events is unclear. In the present study, by comparing chromatin structure and H3K27ac levels in normal and lymphoma B cells, we show that enhancer-promoter-interacting regions assume different conformations according to the local abundance of H3K27ac. Genetic or pharmacological depletion of H3K27ac decreases the frequency and the spreading of these interactions, altering oncogene expression. Moreover, enhancer hijacking mediated by chromosomal translocations influences the epigenetic status of the regions flanking the breakpoint, prompting the formation of distinct intrachromosomal interactions in the two homologous chromosomes. These interactions are accompanied by allele-specific gene expression changes. Overall, our work indicates that H3K27ac dynamics modulates interaction frequency between regulatory regions and can lead to allele-specific chromatin configurations to sustain oncogene expression.

NKL-Code in Normal and Aberrant Hematopoiesis

Simple Summary

Gene codes represent expression patterns of closely related genes in particular tissues, organs or body parts. The NKL-code describes the activity of NKL homeobox genes in the hematopoietic system. NKL homeobox genes encode transcription factors controlling basic developmental processes. Therefore, aberrations of this code may contribute to deregulated hematopoiesis including leukemia and lymphoma. Normal and abnormal activities of NKL homeobox genes are described and mechanisms of (de)regulation, function, and diseases exemplified.

Abstract

We have recently described physiological expression patterns of NKL homeobox genes in early hematopoiesis and in subsequent lymphopoiesis and myelopoiesis, including terminally differentiated blood cells. We thereby systematized differential expression patterns of eleven such genes which form the so-called NKL-code. Due to the developmental impact of NKL homeobox genes, these data suggest a key role for their activity in normal hematopoietic differentiation processes. On the other hand, the aberrant overexpression of NKL-code-members or the ectopical activation of non-code members have been frequently reported in lymphoid and myeloid leukemia/lymphoma, revealing the oncogenic potential of these genes in the hematopoietic compartment. Here, I provide an overview of the NKL-code in normal hematopoiesis and instance mechanisms of deregulation and oncogenic functions of selected NKL genes in hematologic cancers. As well as published clinical studies, our conclusions are based on experimental work using hematopoietic cell lines which represent useful models to characterize the role of NKL homeobox genes in specific tumor types.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Identification of novel, clonally stable, somatic mutations targeting transcription factors PAX5 and NKX2-3, the epigenetic regulator LRIF1, and BRAF in a case of atypical B-cell chronic lymphocytic leukemia harboring a t(14;18)(q32;q21)

Diagnosis of B-cell chronic lymphocytic leukemia (B-CLL) is usually straightforward, involving clinical, immunophenotypic (Matutes score), and (immuno)genetic analyses (to refine patient prognosis for treatment). CLL cases with atypical presentation (e.g., Matutes ≤ 3) are also encountered, and for these diseases, biology and prognostic impact are less clear. Here we report the genomic characterization of a case of atypical B-CLL in a 70-yr-old male patient; B-CLL cells showed a Matutes score of 3, chromosomal translocation t(14;18)(q32;q21) (BCL2/IGH), mutated IGHV, deletion 17p, and mutations in BCL2, NOTCH1 (subclonal), and TP53 (subclonal). Quite strikingly, a novel PAX5 mutation that was predicted to be loss of function was also seen. Exome sequencing identified, in addition, a potentially actionable BRAF mutation, together with novel somatic mutations affecting the homeobox transcription factor NKX2-3, known to control B-lymphocyte development and homing, and the epigenetic regulator LRIF1, which is implicated in chromatin compaction and gene silencing. Neither NKX2-3 nor LRIF1 mutations, predicted to be loss of function, have previously been reported in B-CLL. Sequencing confirmed the presence of these mutations together with BCL2, NOTCH1, and BRAF mutations, with the t(14;18)(q32;q21) translocation, in the initial diagnostic sample obtained 12 yr prior. This is suggestive of a role for these novel mutations in B-CLL initiation and stable clonal evolution, including upon treatment withdrawal. This case extends the spectrum of atypical B-CLL with t(14;18)(q32;q21) and highlights the value of more global precision genomics for patient follow-up and treatment in these patients.

Genetic and clinical basis for two distinct subtypes of primary Sjögren's syndrome

OBJECTIVES

Clinical presentation of primary Sjögren's syndrome (pSS) varies considerably. A shortage of evidence-based objective markers hinders efficient drug development and most clinical trials have failed to reach primary endpoints.

METHODS

We performed a multicentre study to identify patient subgroups based on clinical, immunological and genetic features. Targeted DNA sequencing of 1853 autoimmune-related loci was performed. After quality control, 918 patients with pSS, 1264 controls and 107 045 single nucleotide variants remained for analysis. Replication was performed in 177 patients with pSS and 7672 controls.

RESULTS

We found strong signals of association with pSS in the HLA region. Principal component analysis of clinical data distinguished two patient subgroups defined by the presence of SSA/SSB antibodies. We observed an unprecedented high risk of pSS for an association in the HLA-DQA1 locus of odds ratio 6.10 (95% CI: 4.93, 7.54, P=2.2×10-62) in the SSA/SSB-positive subgroup, while absent in the antibody negative group. Three independent signals within the MHC were observed. The two most significant variants in MHC class I and II respectively, identified patients with a higher risk of hypergammaglobulinaemia, leukopenia, anaemia, purpura, major salivary gland swelling and lymphadenopathy. Replication confirmed the association with both MHC class I and II signals confined to SSA/SSB antibody positive pSS.

CONCLUSION

Two subgroups of patients with pSS with distinct clinical manifestations can be defined by the presence or absence of SSA/SSB antibodies and genetic markers in the HLA locus. These subgroups should be considered in clinical follow-up, drug development and trial outcomes, for the benefit of both subgroups.

Identification of Immune-Related Prognostic Biomarkers Associated with HPV-Positive Head and Neck Squamous Cell Carcinoma

Background

As a type of malignant tumor, head and neck squamous cell carcinoma (HNSCC) seriously threatens human health. This study is aimed at constructing a new, reliable prognostic model.

Method

The gene expression profile data of HNSCC patients were downloaded from the Gene Expression Omnibus and The Cancer Genome Atlas databases. The immune-related differentially expressed genes (IRDEGs) related to HNSCC were identified. We then used Cox regression analysis and least absolute shrinkage and selection operator (LASSO) analysis to explore IRDEGs related to the HNSCC prognosis and to construct and validate a risk scoring model and used ESTIMATE to evaluate tumor immune infiltration in HNSCC patients. Finally, we validated IGSF5 expression and function in HNSCC cells.

Results

A total of 1,195 IRDEGs were found from the GSE65858 dataset. Thirty-one of the 1,195 IRDEGs were associated with the prognosis of HNSCC. Nine key IRDEGs were further selected using the LASSO method, and a risk scoring model was established for predicting the survival of HNSCC patients. According to the risk scoring model, the prognosis of patients in the high-risk group was worse than that of the low-risk group; the high-risk group had significantly higher immune scores than the low-risk group; and between the high- and low-risk samples, there were significant differences in the proportion of 10 types of cells, including naive cells, plasma cells, and resting CD4⁺ memory T cells. IGSF5 has low expression in HNSCC, and overexpression of IGSF5 significantly impaired HNSCC cell proliferation.

Conclusion

This prognostic risk assessment model can help systematically evaluate the survival prognosis of HNSCC patients and provides a new research direction for the improvement of the survival prognosis of HNSCC patients in clinical practice.

Significance of TP53 mutation in treatment and prognosis in head and neck squamous cell carcinoma

Background: TP53 is ranked as the most common mutated gene in head and neck squamous cell carcinoma (HNSCC). Results: The status of TP53 mutation was investigated on International Cancer Genome Consortium and The Cancer Genome Atlas database and TP53-related differentially expressed genes were screened out from transcriptome data and mutation information. A TP53-related prognostic gene signature (TIMP4, ONECUT2, CGNL1, DMRTA2 and NKX2.3) was constructed based on Cox regression analysis and LASSO algorithm. Univariate and multivariate analyses were carried out to identify promising prognosticators for HNSCC. Conclusion: Our findings provide a well-rounded landscape of TP53 mutation in HNSCC and pave the groundwork for developing innovative and effective cancer treatment methods for HNSCC.

Development of a Prognostic Signature Based on Autophagy-related Genes for Head and Neck Squamous Cell Carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a malignant tumor with relative low survival rate. Increasingly evidences have emphasized the importance of autophagy in cancer initiation, progression, and the responses to cancer treatment.

AIM OF THE STUDY

This study aimed to investigate the potential biological and prognostic significance of autophagy-related genes (ARGs) in HNSCC patients.

METHODS

We collected a list of ARGs from Human Autophagy Database and obtained expression profiles and clinical information of HNSCC samples from the Cancer Genome Atlas (TCGA) portal. Differential expression analysis and functional enrichment analysis were performed by R software. The prognostic value of differentially expressed ARGs was detected by Cox regression analysis and prognosis-related ARGs were subjected to LASSO regression analysis. Univariate and multivariate Cox regression analysis were applied to identify promising independent prognosticators for HNSCC.

RESULTS

A total of 35 differentially expressed ARGs were screened out and functional enrichment analysis results indicated these genes were mainly associated with autophagy-related biological processes and pathways. Seven prognosis-related ARGs (ITGA3, CDKN2A, FADD, NKX2-3, BAK1, CXCR4, and HSPB8) were selected to construct a risk signature, which proved to be effective in predicting the survival rate of HNSCC patients. Moreover, univariate analysis showed risk score, tumor stage, T stage, and N stage were negatively correlated with patient overall survival and the multivariate Cox regression analysis results indicated risk score, age, and N stage was significantly associated with patient prognosis.

CONCLUSIONS

Our findings may provide novel evidences for the diagnosis and prognosis evaluation for HNSCC.

Targeting pathogenic mechanisms in marginal zone lymphoma: from concepts and beyond

Marginal zone lymphoma (MZL) represents a group of three distinct though overlapping lymphoid malignancies that includes extranodal, nodal and splenic marginal lymphoma. MZL patients usually present an indolent clinical course, although the disease remains largely incurable, save early stage disease that might be irradiated. Therapeutic advances have been limited due to the small patient population, and have largely been adapted from other indolent lymphomas. Here, we discuss the numerous targets and pathways which may offer the prospect of directly inhibiting the mechanisms identified promoting and sustaining marginal zone lymphomagenesis. In particular, we focus on the agents that may have at least a theoretical application in the disease. Various dysregulated pathways converge to produce an overarching stimulation of nuclear factor κB (NF-κB) and the MYD88-IRAK4 axis, which can be thus leveraged or targeting B-cell receptor signaling through BTK inhibitors (such as ibrutinib, zanubrutinib, acalabrutinib) and PI3K inhibitors (such as idelalisib, copanlisib, duvelisib umbralisib) or via more novel agents in development such as MALT1 inhibitors, SMAC mimetics, NIK inhibitors, IRAK4 or MYD88 inhibitors. NOTCH signaling is also crucial for marginal zone cells, but no clinical data are available with NOTCH inhibitors such as the γ-secretase inhibitor PF-03084014 or the NICD inhibitor CB-103. The hypermethylation phenotype, the overexpression of the PRC2-complex or the presence of TET2 mutations reported in MZL subsets make epigenetic agents (demethylating agents, EZH2 inhibitors, HDAC inhibitors) also potential therapeutic tools for MZL patients.

Ameliorated Autoimmune Arthritis and Impaired B Cell Receptor-Mediated Ca2+ Influx in Nkx2-3 Knock-out Mice

B cells play a crucial role in the pathogenesis of rheumatoid arthritis. In Nkx2-3-deficient mice (Nkx2-3^−/−) the spleen’s histological structure is fundamentally changed; therefore, B cell homeostasis is seriously disturbed. Based on this, we were curious, whether autoimmune arthritis could be induced in Nkx2-3^−/− mice and how B cell activation and function were affected. We induced arthritis with immunization of recombinant human proteoglycan aggrecan G1 domain in Nkx2-3^−/− and control BALB/c mice. We followed the clinical picture, characterized the radiological changes, the immune response, and intracellular Ca²⁺ signaling of B cells. Incidence of the autoimmune arthritis was lower, and the disease severity was milder in Nkx2-3^−/− mice than in control BALB/c mice. The radiological changes were in line with the clinical picture. In Nkx2-3^−/− mice, we measured decreased antigen-induced proliferation and cytokine production in spleen cell cultures; in the sera, we found less anti-CCP-IgG2a, IL-17 and IFNγ, but more IL-1β, IL-4 and IL-6. B cells isolated from the lymph nodes of Nkx2-3^−/− mice showed decreased intracellular Ca²⁺ signaling compared to those isolated from BALB/c mice. Our findings show that the transcription factor Nkx2-3 might regulate the development of autoimmune arthritis most likely through modifying B cell activation.

The identification of autophagy-related genes in the prognosis of oral squamous cell carcinoma

OBJECTIVES

The aim of this study was to identify prognostic autophagy-related genes and assess the ability of these genes to predict clinical outcomes in oral squamous cell carcinoma (OSCC).

SUBJECTS AND METHODS

The details of the human autophagy-related genes were obtained from the Human Autophagy Database. The Cancer Genome Atlas database was used to obtain the gene expression profiles and clinical data of patients. Prediction of biological functions of differentially expressed genes was performed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Prognosis-related genes were identified by Cox regression analysis, and the coefficient was applied to construct a prognostic risk score model. The median of the risk score was applied to distinguish between high- and low-risk groups. The Gene Expression Omnibus database, qRT-PCR and immunohistochemistry were used to validate the expression of key genes.

RESULTS

KEGG analyses revealed that differentially expressed genes were mainly enriched in autophagy-related pathways and virus infection. BAK1, BID, NKX2-3 and SPHK1 were identified. The risk score model showed that the high-risk score had poorer overall survival (Kaplan-Meier analysis, p = 1.79 × 10^-7 ). SPHK1 was upregulated in OSCC tissues and cells, and NXK2-3 was downregulated.

CONCLUSIONS

Autophagy-related gene expression profiles may be a potential biomarker for OSCC prognosis.

VLA-4 Expression and Activation in B Cell Malignancies: Functional and Clinical Aspects

Lineage commitment and differentiation of hematopoietic cells takes place in well-defined microenvironmental surroundings. Communication with other cell types is a vital prerequisite for the normal functions of the immune system, while disturbances in this communication support the development and progression of neoplastic disease. Integrins such as the integrin very late antigen-4 (VLA-4; CD49d/CD29) control the localization of healthy as well as malignant B cells within the tissue, and thus determine the patterns of organ infiltration. Malignant B cells retain some key characteristics of their normal counterparts, with B cell receptor (BCR) signaling and integrin-mediated adhesion being essential mediators of tumor cell homing, survival and proliferation. It is thus not surprising that targeting the BCR pathway using small molecule inhibitors has proved highly effective in the treatment of B cell malignancies. Attenuation of BCR-dependent lymphoma–microenvironment interactions was, in this regard, described as a main mechanism critically contributing to the efficacy of these agents. Here, we review the contribution of VLA-4 to normal B cell differentiation on the one hand, and to the pathophysiology of B cell malignancies on the other hand. We describe its impact as a prognostic marker, its interplay with BCR signaling and its predictive role for novel BCR-targeting therapies, in chronic lymphocytic leukemia and beyond.

Whole Exome and Transcriptome RNA-Sequencing Model for the Diagnosis of Prostate Cancer

In our previous study, we developed a genome-wide DNA methylation model for the diagnosis of prostate cancer, and we pointed out that a considerable average error is associated with the current method for the diagnosis of prostate cancer, which is predicated on pathological assessment of biopsied tissue. In this study, we utilized whole exome and transcriptome RNA-sequencing (RNA-seq) data that were derived from 468 tumor samples and 51 normal samples of prostatic tissue, and we analyzed over 20,000 genes per sample. We were able to develop a mathematical model that classified tumor tissue versus normal tissue with a high accuracy. The overall sensitivity was 97.01%, and the overall specificity was 94.12%. The input variables to the model were the mRNA expression values of the following nine genes: ANGPT1, MED21, AOX1, PLP2, HPN, HPN-AS1, EPHA10, NKX2-3, and LRFN1. The model was validated with unknown samples, with a 10-fold cross-validation, and a leave-one-out cross-validation. We present here a genomic model, based on a whole exome and transcriptome RNA-seq analysis of biopsied prostatic tissue, that could be utilized in the diagnosis of prostate cancer.

Deregulated NKL Homeobox Genes in B-Cell Lymphoma

Recently, we have described physiological expression patterns of NKL homeobox genes in early hematopoiesis and in subsequent lymphopoiesis. We identified nine genes which constitute the so-called NKL-code. Aberrant overexpression of code-members or ectopically activated non-code NKL homeobox genes are described in T-cell leukemia and in T- and B-cell lymphoma, highlighting their oncogenic role in lymphoid malignancies. Here, we introduce the NKL-code in normal hematopoiesis and focus on deregulated NKL homeobox genes in B-cell lymphoma, including HLX, MSX1 and NKX2-2 in Hodgkin lymphoma; HLX, NKX2-1 and NKX6-3 in diffuse large B-cell lymphoma; and NKX2-3 in splenic marginal zone lymphoma. Thus, the roles of various members of the NKL homeobox gene subclass are considered in normal and pathological hematopoiesis in detail.

Epstein-Barr virus (EBV) activates NKL homeobox gene HLX in DLBCL

NKL homeobox genes encode developmental transcription factors regulating basic processes in cell differentiation. According to their physiological expression pattern in early hematopoiesis and lymphopoiesis, particular members of this homeobox gene subclass constitute an NKL-code. B-cell specific NKL-code genes generate a regulatory network and their deregulation is implicated in B-cell lymphomagenesis. Epstein-Barr virus (EBV) infects B-cells and influences the activity of signalling pathways including JAK/STAT and several genes encoding developmental regulators. Therefore, EBV-infection impacts the pathogenesis and the outcome of B-cell malignancies including Hodgkin lymphoma and diffuse large B-cell lymphoma (DLBCL). Here, we isolated EBV-positive and EBV-negative subclones from the DLBCL derived cell line DOHH-2. These subclones served as models to investigate the role of EBV in deregulation of the B-cell specific NKL-code members HHEX, HLX, MSX1 and NKX6-3. We showed that the EBV-encoded factors LMP1 and LMP2A activated the expression of HLX via STAT3. HLX in turn repressed NKX6-3, SPIB and IL4R which normally mediate plasma cell differentiation. In addition, HLX repressed the pro-apoptotic factor BCL2L11/BIM and hence supported cell survival. Thus, EBV aberrantly activated HLX in DLBCL, thereby disturbing both B-cell differentiation and apoptosis. The results of our study appreciate the pathogenic role of EBV in NKL homeobox gene deregulation and B-cell malignancies.

PD-1/PD-L1 immune checkpoint and p53 loss facilitate tumor progression in activated B-cell diffuse large B-cell lymphomas

Refractory or relapsed diffuse large B-cell lymphoma (DLBCL) often associates with the activated B-cell-like (ABC) subtype and genetic alterations that drive constitutive NF-κB activation and impair B-cell terminal differentiation. Here, we show that DNA damage response by p53 is a central mechanism suppressing the pathogenic cooperation of IKK2ca-enforced canonical NF-κB and impaired differentiation resulting from Blimp1 loss in ABC-DLBCL lymphomagenesis. We provide evidences that the interplay between these genetic alterations and the tumor microenvironment select for additional molecular addictions that promote lymphoma progression, including aberrant coexpression of FOXP1 and the B-cell mutagenic enzyme activation-induced deaminase, and immune evasion through major histocompatibility complex class II downregulation, PD-L1 upregulation, and T-cell exhaustion. Consistently, PD-1 blockade cooperated with anti-CD20-mediated B-cell cytotoxicity, promoting extended T-cell reactivation and antitumor specificity that improved long-term overall survival in mice. Our data support a pathogenic cooperation among NF-κB-driven prosurvival, genetic instability, and immune evasion mechanisms in DLBCL and provide preclinical proof of concept for including PD-1/PD-L1 blockade in combinatorial immunotherapy for ABC-DLBCL.

IL-22-Independent Protection from Colitis in the Absence of Nkx2.3 Transcription Factor in Mice

The transcription factor Nkx2.3 regulates the vascular specification of Peyer patches in mice through determining endothelial addressin preference and may function as a susceptibility factor in inflammatory bowel diseases in humans. We wished to analyze the role of Nkx2.3 in colonic solitary intestinal lymphoid tissue composition and in colitis pathogenesis. We studied the colonic solitary intestinal lymphoid tissue of Nkx2.3-deficient mice with immunofluorescence and flow cytometry. Colitis was induced in mice using 2.5% dextran sodium sulfate, and severity was assessed with histology, flow cytometry, and quantitative PCR. We found that the lack of Nkx2.3 impairs maturation of isolated lymphoid follicles and attenuates dextran sodium sulfate-induced colitis independent of endothelial absence of mucosal addressin cell-adhesion molecule-1 (MAdCAM-1), which was also coupled with enhanced colonic epithelial regeneration. Although we observed increased numbers of group 3 innate lymphoid cells and Th17 cells and enhanced transcription of IL-22, Ab-mediated neutralization of IL-22 did not abolish the protection from colitis in Nkx2.3-deficient mice. Nkx2.3-/- hematopoietic cells could not rescue wild-type mice from colitis. Using LacZ-Nkx2.3 reporter mice, we found that Nkx2.3 expression was restricted to VAP-1+ myofibroblast-like pericryptal cells. These results hint at a previously unknown stromal role of Nkx2.3 as driver of colitis and indicate that Nkx2.3+ stromal cells play a role in epithelial cell homeostasis.

NKL homeobox gene NKX2-2 is aberrantly expressed in Hodgkin lymphoma

NKL homeobox genes encode basic transcriptional regulators of cell and tissue differentiation. Recently, we described a hematopoietic NKL-code comprising nine specific NKL homeobox genes expressed in normal hematopoietic stem cells, lymphoid progenitors and during lymphopoiesis, highlighting their physiological role in the development of T-, B- and NK-cells. Here, we identified aberrant expression of the non-hematopoietic neural NKL homeobox gene NKX2-2 in about 12% of both, classical Hodgkin lymphoma (HL) and nodular lymphocyte predominant (NLP) HL patients. The NKX2-2 expressing NLPHL-derived cell line DEV served as a model by analysing chromosomal configurations and expression profiling data to reveal activating mechanisms and downstream targets of this developmental regulator. While excluding chromosomal rearrangements at the locus of NKX2-2 we identified t(3;14)(p21;q32) resulting in overexpression of the IL17 receptor gene IL17RB via juxtaposition with the IGH-locus. SiRNA-mediated knockdown experiments demonstrated that IL17RB activated NKX2-2 transcription. Overexpression of IL17RB-cofactor DAZAP2 via chromosomal gain of 12q13 and deletion of its proteasomal inhibitor SMURF2 at 17q24 supported expression of NKX2-2. IL17RB activated transcription factors FLI1 and FOXG1 which in turn mediated NKX2-2 expression. In addition, overexpressed chromatin-modulator AUTS2 contributed to NKX2-2 activation as well. Downstream analyses indicated that NKX2-2 inhibits transcription of lymphoid NKL homeobox gene MSX1 and activates expression of basic helix-loop-helix factor NEUROD1 which may disturb B-cell differentiation processes via reported interaction with TCF3/E2A. Taken together, our data reveal ectopic activation of a neural gene network in HL placing NKX2-2 at its hub, highlighting a novel oncogenic impact of NKL homeobox genes in B-cell malignancies.

NKL homeobox gene activities in B-cell development and lymphomas

Homeobox genes encode transcription factors which regulate basic processes in development and cell differentiation. Several members of the NKL subclass are deregulated in T-cell progenitors and support leukemogenesis. We have recently described particular expression patterns of nine NKL homeobox genes in early hematopoiesis and T-cell development. Here, we screened NKL homeobox gene activities in normal B-cell development and extended the NKL-code to include this lymphoid lineage. Analysis of public expression profiling datasets revealed that HHEX and NKX6-3 were the only members differentially active in naïve B-cells, germinal center B-cells, plasma cells and memory B-cells. Subsequent examination of different types of B-cell malignancies showed both aberrant overexpression of NKL-code members and ectopic activation of subclass members physiologically silent in lymphopoiesis including BARX2, DLX1, EMX2, NKX2-1, NKX2-2 and NKX3-2. Based on these findings we performed detailed studies of the B-cell specific NKL homeobox gene NKX6-3 which showed enhanced activity in patient subsets of follicular lymphoma, mantle cell lymphoma and diffuse large B-cell lymphoma (DLBCL), and in three DLBCL cell lines to serve as in vitro models. While excluding genomic and chromosomal rearrangements at the locus of NKX6-3 (8p11) promoter studies demonstrated that B-cell factors MYB and PAX5 activated NKX6-3 transcription. Furthermore, aberrant BMP7/SMAD1-signalling and deregulated expression of chromatin complex components AUTS2 and PCGF5 promoted NKX6-3 activation. Finally, NKL homeobox genes HHEX, HLX, MSX1 and NKX6-3 were expressed in B-cell progenitors and generated a regulatory gene network in cell lines which we propose may provide physiological support for NKL-code formation in early B-cell development. Together, we identified an NKL-code in B-cell development whose violation may deregulate differentiation and promote malignant transformation.

The transcription factor NKX2-3 mediates p21 expression and ectodysplasin-A signaling in the enamel knot for cusp formation in tooth development

Tooth morphogenesis is initiated by reciprocal interactions between the ectoderm and neural crest-derived mesenchyme. During tooth development, tooth cusps are regulated by precise control of proliferation of cell clusters, termed enamel knots, that are present among dental epithelial cells. The interaction of ectodysplasin-A (EDA) with its receptor, EDAR, plays a critical role in cusp formation by these enamel knots, and mutations of these genes is a cause of ectodermal dysplasia. It has also been reported that deficiency in Nkx2-3, encoding a member of the NK2 homeobox family of transcription factors, leads to cusp absence in affected teeth. However, the molecular role of NKX2-3 in tooth morphogenesis is not clearly understood. Using gene microarray analysis in mouse embryos, we found that Nkx2-3 is highly expressed during tooth development and increased during the tooth morphogenesis, especially during cusp formation. We also demonstrate that NKX2-3 is a target molecule of EDA and critical for expression of the cell cycle regulator p21 in the enamel knot. Moreover, NKX2-3 activated the bone morphogenetic protein (BMP) signaling pathway by up-regulating expression levels of Bmp2 and Bmpr2 in dental epithelium and decreased the expression of the dental epithelial stem cell marker SRY box 2 (SOX2). Together, our results indicate that EDA/NKX2-3 signaling is essential for enamel knot formation during tooth morphogenesis in mice.

Richter transformation driven by Epstein-Barr virus reactivation during therapy-related immunosuppression in chronic lymphocytic leukaemia

The increased risk of Richter transformation (RT) in patients with chronic lymphocytic leukaemia (CLL) due to Epstein-Barr virus (EBV) reactivation during immunosuppressive therapy with fludarabine other targeted agents remains controversial. Among 31 RT cases classified as diffuse large B-cell lymphoma (DLBCL), seven (23%) showed EBV expression. In contrast to EBV^- tumours, EBV⁺ DLBCLs derived predominantly from IGVH-hypermutated CLL, and they also showed CLL-unrelated IGVH sequences more frequently. Intriguingly, despite having different cellular origins, clonally related and unrelated EBV⁺ DLBCLs shared a previous history of immunosuppressive chemo-immunotherapy, a non-germinal centre DLBCL phenotype, EBV latency programme type II or III, and very short survival. These data suggested that EBV reactivation during therapy-related immunosuppression can transform either CLL cells or non-tumoural B lymphocytes into EBV⁺ DLBCL. To investigate this hypothesis, xenogeneic transplantation of blood cells from 31 patients with CLL and monoclonal B-cell lymphocytosis (MBL) was performed in Rag2^-/- IL2γc^-/- mice. Remarkably, the recipients' impaired immunosurveillance favoured the spontaneous outgrowth of EBV⁺ B-cell clones from 95% of CLL and 64% of MBL patients samples, but not from healthy donors. Eventually, these cells generated monoclonal tumours (mostly CLL-unrelated but also CLL-related), recapitulating the principal features of EBV⁺ DLBCL in patients. Accordingly, clonally related and unrelated EBV⁺ DLBCL xenografts showed indistinguishable cellular, virological and molecular features, and synergistically responded to combined inhibition of EBV replication with ganciclovir and B-cell receptor signalling with ibrutinib in vivo. Our study underscores the risk of RT driven by EBV in CLL patients receiving immunosuppressive therapies, and provides the scientific rationale for testing ganciclovir and ibrutinib in EBV⁺ DLBCL. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Aberrant expression of NKL homeobox gene HLX in Hodgkin lymphoma

NKL homeobox genes are basic regulators of cell and tissue differentiation, many acting as oncogenes in T-cell leukemia. Recently, we described an hematopoietic NKL-code comprising six particular NKL homeobox genes expressed in hematopoietic stem cells and lymphoid progenitors, unmasking their physiological roles in the development of these cell types. Hodgkin lymphoma (HL) is a B-cell malignancy showing aberrant activity of several developmental genes resulting in disturbed B-cell differentiation. To examine potential concordances in abnormal lymphoid differentiation of T- and B-cell malignancies we analyzed the expression of the hematopoietic NKL-code associated genes in HL, comprising HHEX, HLX, MSX1, NKX2-3, NKX3-1 and NKX6-3. Our approach revealed aberrant HLX activity in 8 % of classical HL patients and additionally in HL cell line L-540. Accordingly, to identify upstream regulators and downstream target genes of HLX we used L-540 cells as a model and performed chromosome and genome analyses, comparative expression profiling and functional assays via knockdown and overexpression experiments therein. These investigations excluded chromosomal rearrangements of the HLX locus at 1q41 and demonstrated that STAT3 operated directly as transcriptional activator of the HLX gene. Moreover, subcellular analyses showed highly enriched STAT3 protein in the nucleus of L-540 cells which underwent cytoplasmic translocation by repressing deacetylation. Finally, HLX inhibited transcription of B-cell differentiation factors MSX1, BCL11A and SPIB and of pro-apoptotic factor BCL2L11/BIM, thereby suppressing Etoposide-induced cell death. Collectively, we propose that aberrantly expressed NKL homeobox gene HLX is part of a pathological gene network in HL, driving deregulated B-cell differentiation and survival.

Nkx2-3-A Slippery Slope From Development Through Inflammation Toward Hematopoietic Malignancies

The development of peripheral lymphoid tissues from the mesoderm is the result of a complex convergence combining lymphohematopoietic differentiation with the local specification of nonhematopoietic mesenchymal components. Although the various transcriptional regulators with fate-determining effects in diversifying the mobile leukocyte subsets have been thoroughly studied and identified, the tissue-specific determinants promoting the regional differentiation of resident mesenchyme are less understood. Of these factors, various members of the NK-class Nkx paralogues have emerged as key regulators for the organogenesis of spleen and mucosal lymphoid tissues, and recent data have also indicated their involvement in various pathological events, including gut inflammation and hematopoietic malignancies. Here, we summarize available data on the roles of Nkx2-3 in lymphoid tissue development and discuss its possible value as a developmental marker and disease-associated pathogenic trait.

Distinguishing three subtypes of hematopoietic cells based on gene expression profiles using a support vector machine

Hematopoiesis is a complicated process involving a series of biological sub-processes that lead to the formation of various blood components. A widely accepted model of early hematopoiesis proceeds from long-term hematopoietic stem cells (LT-HSCs) to multipotent progenitors (MPPs) and then to lineage-committed progenitors. However, the molecular mechanisms of early hematopoiesis have not been fully characterized. In this study, we applied a computational strategy to identify the gene expression signatures distinguishing three types of closely related hematopoietic cells collected in recent studies: (1) hematopoietic stem cell/multipotent progenitor cells; (2) LT-HSCs; and (3) hematopoietic progenitor cells. Each cell in these cell types was represented by its gene expression profile among a total number of 20,475 genes. The expression features were analyzed by a Monte-Carlo Feature Selection (MCFS) method, resulting in a feature list. Then, the incremental feature selection (IFS) and a support vector machine (SVM) optimized with a sequential minimum optimization (SMO) algorithm were employed to access the optimal classifier with the highest Matthews correlation coefficient (MCC) value of 0.889, in which 6698 features were used to represent cells. In addition, through an updated program of MCFS method, seventeen decision rules can be obtained, which can classify the three cell types with an overall accuracy of 0.812. Using a literature review, both the rules and the top features used for building the optimal classifier were confirmed to be commonly used or potential biological markers for distinguishing the three cell types of HSPCs. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

Targeting the anion exchanger 2 with specific peptides as a new therapeutic approach in B lymphoid neoplasms

Regulatory T (Treg) cells can weaken antitumor immune responses, and inhibition of their function appears to be a promising therapeutic approach in cancer patients. Mice with targeted deletion of the gene encoding the Cl^-/HCO₃^- anion exchanger AE2 (also termed SLC4A2), a membrane-bound carrier involved in intracellular pH regulation, showed a progressive decrease in the number of Treg cells. We therefore challenged AE2 as a potential target for tumor therapy, and generated linear peptides designed to bind the third extracellular loop of AE2, which is crucial for its exchange activity. Peptide p17AE2 exhibited optimal interaction ability and indeed promoted apoptosis in mouse and human Treg cells, while activating effector T-cell function. Interestingly, this linear peptide also induced apoptosis in different types of human leukemia, lymphoma and multiple myeloma cell lines and primary malignant samples, while it showed only moderate effects on normal B lymphocytes. Finally, a macrocyclic AE2 targeting peptide exhibiting increased stability in vivo was effective in mice xenografted with B-cell lymphoma. These data suggest that targeting the anion exchanger AE2 with specific peptides may represent an effective therapeutic approach in B-cell malignancies.

TCRα rearrangements identify a subgroup of NKL-deregulated adult T-ALLs associated with favorable outcome

T-cell acute lymphoblastic leukemia (T-ALL) results from leukemic transformation of T-cell precursors arrested at specific differentiation stages, including an 'early-cortical' thymic maturation arrest characterized by expression of cytoplasmic TCRβ but no surface T-cell receptor (TCR) and frequent ectopic expression of the TLX1/3 NK-like homeotic proteins (NKL). We designed a TCRα VJC PCR to identify clonal TCRα rearrangements in 32% of 127 T-ALLs, including 0/52 immature/TCRγδ lineage cases and 41/75 (55%) TCRαβ lineage cases. Amongst the latter, TCRα rearrangements were not identified in 30/54 (56%) of IMβ/pre-αβ early-cortical T-ALLs, of which the majority (21/30) expressed TLX1/3. We reasoned that the remaining T-ALLs might express other NKL proteins, so compared transcript levels of 46 NKL in T-ALL and normal thymic subpopulations. Ectopic overexpression of 10 NKL genes, of which six are unreported in T-ALL (NKX2-3, BARHL1, BARX2, EMX2, LBX2 and MSX2), was detectable in 17/104 (16%) T-ALLs. Virtually all NKL overexpressing T-ALLs were TCRα unrearranged and ectopic NKL transcript expression strongly repressed Eα activity, suggesting that ectopic NKL expression is the major determinant in early-cortical thymic T-ALL maturation arrest. This immunogenetic T-ALL subtype, defined by TCRβ VDJ but no TCRα VJ rearrangement, is associated with a favorable outcome in GRAALL-treated adult T-ALLs.

NKL homeobox gene activities in hematopoietic stem cells, T-cell development and T-cell leukemia

T-cell acute lymphoblastic leukemia (T-ALL) cells represent developmentally arrested T-cell progenitors, subsets of which aberrantly express homeobox genes of the NKL subclass, including TLX1, TLX3, NKX2-1, NKX2-5, NKX3-1 and MSX1. Here, we analyzed the transcriptional landscape of all 48 members of the NKL homeobox gene subclass in CD34+ hematopoietic stem and progenitor cells (HSPCs) and during lymphopoiesis, identifying activities of nine particular genes. Four of these were expressed in HSPCs (HHEX, HLX1, NKX2-3 and NKX3-1) and three in common lymphoid progenitors (HHEX, HLX1 and MSX1). Interestingly, our data indicated downregulation of NKL homeobox gene transcripts in late progenitors and mature T-cells, a phenomenon which might explain the oncogenic impact of this group of genes in T-ALL. Using MSX1-expressing T-ALL cell lines as models, we showed that HHEX activates while HLX1, NKX2-3 and NKX3-1 repress MSX1 transcription, demonstrating the mutual regulation and differential activities of these homeobox genes. Analysis of a public T-ALL expression profiling data set comprising 117 patient samples identified 20 aberrantly activated members of the NKL subclass, extending the number of known NKL homeobox oncogene candidates. While 7/20 genes were also active during hematopoiesis, the remaining 13 showed ectopic expression. Finally, comparative analyses of T-ALL patient and cell line profiling data of NKL-positive and NKL-negative samples indicated absence of shared target genes but instead highlighted deregulation of apoptosis as common oncogenic effect. Taken together, we present a comprehensive survey of NKL homeobox genes in early hematopoiesis, T-cell development and T-ALL, showing that these genes generate an NKL-code for the diverse stages of lymphoid development which might be fundamental for regular differentiation.