Diversity of human and mouse homeobox gene expression in development and adult tissues

The Diagnostic and Prognostic Values of HOXA Gene Family in Kidney Clear Cell Renal Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) is one of the most common cancers with high mortality worldwide. As members of the homeobox (HOX) family, homeobox-A (HOXA) genes have been reported to play an increasingly important role in tumorigenesis and the progression of multiple cancers. However, limited studies have investigated the potential diagnostic and prognostic roles of HOXA genes in KIRC. In this research, we explored the expression pattern of the HOXA gene family in KIRC progression by differential analysis of expression profiles from The Cancer Genome Atlas (TCGA). By using univariate Cox analysis and lasso regression analysis, we comprehensively evaluated the prognostic value of HOXA genes and eventually identified a prognostic risk model consisting of five HOXA genes (HOXA2, HOXA3, HOXA7, HOXA11, and HOXA13). The risk model was further validated as a novel independent prognostic factor for KIRC patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses, and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, to explore the potential mechanism of tumorigenesis and clinical application of KIRC, we also developed the HOXA-based competing endogenous RNA (ceRNA) regulatory network and machine learning classification model. Valproic acid and tretinoin were predicted to be the most promising small molecules to adjuvant treatment of KIRC by mining the CMAP and DGIdb drug database. Subsequently, pathway and functional enrichment analyses provided us with new ways to search for a possible mechanism of action of drugs. Taken together, our study demonstrated the nonnegligible role of HOXA genes in KIRC and constructed an effective prognostic and diagnostic model, which offers novel insights into KIRC prognosis.

A genome-wide identification and analysis of the homeobox genes in the brown planthopper, Nilaparvata lugens (Hemiptera: Delphacidae)

The homeobox family is a large and diverse superclass of genes, many of which act as transcription factors that play important roles in tissue differentiation and embryogenesis in animals. The brown planthopper (BPH), Nilaparvata lugens, is the most destructive pest of rice in Asia, and high fecundity contributes significantly to its ecological success in natural and agricultural habits. Here, we identified 94 homeobox genes in BPH, which could be divided into 75 gene families and 9 classes. This number is comparable to the number of homeobox genes found in the honeybee Apis mellifera, but is slightly less than in Drosophila or the red flour beetle Tribolium castaneum. A spatio-temporal analysis indicated that most BPH homeobox genes were expressed in a development and tissue-specific manner, of which 21 genes were highly expressed in ovaries. RNA interference (RNAi)-mediated functional assay showed that 22 homeobox genes were important for nymph development and the nymph to adult transition, whereas 67 genes were dispensable during this process. Fecundity assay showed that knockdown of 13 ovary-biased genes (zfh1, schlank, abd-A, Lim3_2, Lmxb, Prop, ap_1, Not, lab, Hmx, vis, Pknox, and C15) led to the reproductive defect. This is the first comprehensive investigation into homeobox genes in a hemipteran insect and thus helps us to understand the functional significance of homeobox genes in insect reproduction.

Heritable pattern of oxidized DNA base repair coincides with pre-targeting of repair complexes to open chromatin

Human genome stability requires efficient repair of oxidized bases, which is initiated via damage recognition and excision by NEIL1 and other base excision repair (BER) pathway DNA glycosylases (DGs). However, the biological mechanisms underlying detection of damaged bases among the million-fold excess of undamaged bases remain enigmatic. Indeed, mutation rates vary greatly within individual genomes, and lesion recognition by purified DGs in the chromatin context is inefficient. Employing super-resolution microscopy and co-immunoprecipitation assays, we find that acetylated NEIL1 (AcNEIL1), but not its non-acetylated form, is predominantly localized in the nucleus in association with epigenetic marks of uncondensed chromatin. Furthermore, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) revealed non-random AcNEIL1 binding near transcription start sites of weakly transcribed genes and along highly transcribed chromatin domains. Bioinformatic analyses revealed a striking correspondence between AcNEIL1 occupancy along the genome and mutation rates, with AcNEIL1-occupied sites exhibiting fewer mutations compared to AcNEIL1-free domains, both in cancer genomes and in population variation. Intriguingly, from the evolutionarily conserved unstructured domain that targets NEIL1 to open chromatin, its damage surveillance of highly oxidation-susceptible sites to preserve essential gene function and to limit instability and cancer likely originated ∼500 million years ago during the buildup of free atmospheric oxygen.

HOXA5 Expression Is Elevated in Breast Cancer and Is Transcriptionally Regulated by Estradiol

HOXA5 is a homeobox-containing gene associated with the development of the lung, gastrointestinal tract, and vertebrae. Here, we investigate potential roles and the gene regulatory mechanism in HOXA5 in breast cancer cells. Our studies demonstrate that HOXA5 expression is elevated in breast cancer tissues and in estrogen receptor (ER)-positive breast cancer cells. HOXA5 expression is critical for breast cancer cell viability. Biochemical studies show that estradiol (E2) regulates HOXA5 gene expression in cultured breast cancer cells in vitro. HOXA5 expression is also upregulated in vivo in the mammary tissues of ovariectomized female rats. E2-induced HOXA5 expression is coordinated by ERs. Knockdown of either ERα or ERβ downregulated E2-induced HOXA5 expression. Additionally, ER co-regulators, including CBP/p300 (histone acetylases) and MLL-histone methylases (MLL2, MLL3), histone acetylation-, and H3K4 trimethylation levels are enriched at the HOXA5 promoter in present E2. In summary, our studies demonstrate that HOXA5 is overexpressed in breast cancer and is transcriptionally regulated via estradiol in breast cancer cells.

Complex impact of DNA methylation on transcriptional dysregulation across 22 human cancer types

Accumulating evidence has demonstrated that transcriptional regulation is affected by DNA methylation. Understanding the perturbation of DNA methylation-mediated regulation between transcriptional factors (TFs) and targets is crucial for human diseases. However, the global landscape of DNA methylation-mediated transcriptional dysregulation (DMTD) across cancers has not been portrayed. Here, we systematically identified DMTD by integrative analysis of transcriptome, methylome and regulatome across 22 human cancer types. Our results revealed that transcriptional regulation was affected by DNA methylation, involving hundreds of methylation-sensitive TFs (MethTFs). In addition, pan-cancer MethTFs, the regulatory activity of which is generally affected by DNA methylation across cancers, exhibit dominant functional characteristics and regulate several cancer hallmarks. Moreover, pan-cancer MethTFs were found to be affected by DNA methylation in a complex pattern. Finally, we investigated the cooperation among MethTFs and identified a network module that consisted of 43 MethTFs with prognostic potential. In summary, we systematically dissected the transcriptional dysregulation mediated by DNA methylation across cancer types, and our results provide a valuable resource for both epigenetic and transcriptional regulation communities.

NKL homeobox gene activities in normal and malignant myeloid cells

Recently, we have documented a hematopoietic NKL-code mapping physiological expression patterns of NKL homeobox genes in early hematopoiesis and in lymphopoiesis, which spotlights genes deregulated in lymphoid malignancies. Here, we enlarge this map to include normal NKL homeobox gene expressions in myelopoiesis by analyzing public expression profiling data and primary samples from developing and mature myeloid cells. We thus uncovered differential activities of six NKL homeobox genes, namely DLX2, HHEX, HLX, HMX1, NKX3-1 and VENTX. We further examined public expression profiling data of 251 acute myeloid leukemia (AML) and 183 myelodysplastic syndrome (MDS) patients, thereby identifying 24 deregulated genes. These results revealed frequent deregulation of NKL homeobox genes in myeloid malignancies. For detailed analysis we focused on NKL homeobox gene NANOG, which acts as a stem cell factor and is correspondingly expressed alone in hematopoietic progenitor cells. We detected aberrant expression of NANOG in a small subset of AML patients and in AML cell line NOMO-1, which served as a model. Karyotyping and genomic profiling discounted rearrangements of the NANOG locus at 12p13. But gene expression analyses of AML patients and AML cell lines after knockdown and overexpression of NANOG revealed regulators and target genes. Accordingly, NKL homeobox genes HHEX, DLX5 and DLX6, stem cell factors STAT3 and TET2, and the NOTCH-pathway were located upstream of NANOG while NKL homeobox genes HLX and VENTX, transcription factors KLF4 and MYB, and anti-apoptosis-factor MIR17HG represented target genes. In conclusion, we have extended the NKL-code to the myeloid lineage and thus identified several NKL homeobox genes deregulated in AML and MDS. These data indicate a common oncogenic role of NKL homeobox genes in both lymphoid and myeloid malignancies. For misexpressed NANOG we identified an aberrant regulatory network, which contributes to the understanding of the oncogenic activity of NKL homeobox genes.

How Weird is The Worm? Evolution of the Developmental Gene Toolkit in Caenorhabditis elegans

Comparative developmental biology and comparative genomics are the cornerstones of evolutionary developmental biology. Decades of fruitful research using nematodes have produced detailed accounts of the developmental and genomic variation in the nematode phylum. Evolutionary developmental biologists are now utilising these data as a tool with which to interrogate the evolutionary basis for the similarities and differences observed in Nematoda. Nematodes have often seemed atypical compared to the rest of the animal kingdom-from their totally lineage-dependent mode of embryogenesis to their abandonment of key toolkit genes usually deployed by bilaterians for proper development-worms are notorious rule breakers of the bilaterian handbook. However, exploring the nature of these deviations is providing answers to some of the biggest questions about the evolution of animal development. For example, why is the evolvability of each embryonic stage not the same? Why can evolution sometimes tolerate the loss of genes involved in key developmental events? Lastly, why does natural selection act to radically diverge toolkit genes in number and sequence in certain taxa? In answering these questions, insight is not only being provided about the evolution of nematodes, but of all metazoans.

The Role of the HOXA Gene Family in Acute Myeloid Leukemia

The HOXA gene family is associated with various cancer types. However, the role of HOXA genes in acute myeloid leukemia (AML) have not been comprehensively studied. We compared the transcriptional expression, survival data, and network analysis of HOXA-associated signaling pathways in patients with AML using the ONCOMINE, GEPIA, LinkedOmics, cBioPortal, and Metascape databases. We observed that HOXA2-10 mRNA expression levels were significantly upregulated in AML and that high HOXA1-10 expression was associated with poor AML patient prognosis. The HOXA genes were altered in ~18% of the AML samples, either in terms of amplification, deep deletion, or elevated mRNA expression. The following pathways were modulated by HOXA gene upregulation: GO:0048706: embryonic skeletal system development; R-HSA-5617472: activation of HOX genes in anterior hindbrain development during early embryogenesis; GO:0060216: definitive hemopoiesis; hsa05202: transcriptional mis-regulation in cancer; and GO:0045638: negative regulation of myeloid cell differentiation, and they were significantly regulated due to alterations affecting the HOXA genes. This study identified HOXA3-10 genes as potential AML therapeutic targets and prognostic markers.

Direct and Indirect Targeting of HOXA9 Transcription Factor in Acute Myeloid Leukemia

HOXA9 (Homeobox A9) is a homeotic transcription factor known for more than two decades to be associated with leukemia. The expression of HOXA9 homeoprotein is associated with anterior-posterior patterning during embryonic development, and its expression is then abolished in most adult cells, with the exception of hematopoietic progenitor cells. The oncogenic function of HOXA9 was first assessed in human acute myeloid leukemia (AML), particularly in the mixed-phenotype associated lineage leukemia (MPAL) subtype. HOXA9 expression in AML is associated with aggressiveness and a poor prognosis. Since then, HOXA9 has been involved in other hematopoietic malignancies and an increasing number of solid tumors. Despite this, HOXA9 was for a long time not targeted to treat cancer, mainly since, as a transcription factor, it belongs to a class of protein long considered to be an "undruggable" target; however, things have now evolved. The aim of the present review is to focus on the different aspects of HOXA9 targeting that could be achieved through multiple ways: (1) indirectly, through the inhibition of its expression, a strategy acting principally at the epigenetic level; or (2) directly, through the inhibition of its transcription factor function by acting at either the protein/protein interaction or the protein/DNA interaction interfaces.

Transcriptome analysis reveals differentially expressed genes between human primary bone marrow mesenchymal stem cells and human primary dermal fibroblasts

Stromal cells have been widely used in biomedical research and disease modeling studies in vitro. The most commonly used types of stromal cells are mesenchymal stem cells and fibroblasts. Their cellular phenotypes and differentiation capabilities are quite similar and there are no specific distinction criteria. In order to identify transcriptomic differences between these 2 cell types, we performed next-generation RNA sequencing. Using the global gene expression profile and pathway analysis, we showed that human primary bone marrow mesenchymal stem cells and human primary dermal fibroblasts have different molecular signatures. We also identified critical transcription factors that are differentially expressed between these cells. We then proposed that homeobox genes and some other sequence-specific transcription factors are not only responsible for transcriptional differences, but also discriminate bone marrow mesenchymal stem cells and dermal fibroblasts at the transcriptional level.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

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.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Novel and divergent genes in the evolution of placental mammals

Analysis of genome sequences within a phylogenetic context can give insight into the mode and tempo of gene and protein evolution, including inference of gene ages. This can reveal whether new genes arose on particular evolutionary lineages and were recruited for new functional roles. Here, we apply MCL clustering with all-versus-all reciprocal BLASTP to identify and phylogenetically date 'Homology Groups' among vertebrate proteins. Homology Groups include new genes and highly divergent duplicate genes. Focusing on the origin of the placental mammals within the Eutheria, we identify 357 novel Homology Groups that arose on the stem lineage of Placentalia, 87 of which are deduced to play core roles in mammalian biology as judged by extensive retention in evolution. We find the human homologues of novel eutherian genes are enriched for expression in preimplantation embryo, brain, and testes, and enriched for functions in keratinization, reproductive development, and the immune system.

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.

A sister of NANOG regulates genes expressed in pre-implantation human development

The NANOG homeobox gene plays a pivotal role in self-renewal and maintenance of pluripotency in human, mouse and other vertebrate embryonic stem cells, and in pluripotent cells of the blastocyst inner cell mass. There is a poorly studied and atypical homeobox locus close to the Nanog gene in some mammals which could conceivably be a cryptic paralogue of NANOG, even though the loci share only 20% homeodomain identity. Here we argue that this gene, NANOGNB (NANOG Neighbour), is an extremely divergent duplicate of NANOG that underwent radical sequence change in the mammalian lineage. Like NANOG, the NANOGNB gene is expressed in pre-implantation embryos of human and cow; unlike NANOG, NANOGNB expression is restricted to 8-cell and morula stages, preceding blastocyst formation. When expressed ectopically in adult cells, human NANOGNB elicits gene expression changes, including downregulation of a set of genes that have an expression pulse at the 8-cell stage of pre-implantation development. We conclude that gene duplication and massive sequence divergence in mammals generated a novel homeobox gene that acquired new developmental roles complementary to those of Nanog.

Characterization of the porcine Nanog 5'-flanking region

Objective

Nanog, a homeodomain protein, has been investigated in humans and mice using embryonic stem cells (ESCs). Because of the limited availability of ESCs, few studies have reported the function and role of Nanog in porcine ESCs. Therefore, in this study, we investigated the location of the porcine Nanog chromosome and its basal promoter activity, which might have potential applications in development of ESCs specific marker as well as understanding its operating systems in the porcine.

Methods

To characterize the porcine Nanog promoter, the 5′-flanking region of Nanog was isolated from cells of mini-pig ears. BLAST database search showed that there are two porcine Nanog genomic loci, chromosome 1 and 5, both of which contain an exon with a start codon. Deletion mutants from the 5′-flanking region of both loci were measured using the Dual-Luciferase Reporter Assay System, and a fluorescence marker, green fluorescence protein.

Results

Promoter activity was detected in the sequences of chromosome 5, but not in those of chromosome 1. We identified the sequences from –99 to +194 that possessed promoter activity and contained transcription factor binding sites from deletion fragment analysis. Among the transcription factor binding sites, a Sp1 was found to play a crucial role in basal promoter activity, and point mutation of this site abolished its activity, confirming its role in promoter activity. Furthermore, gel shift analysis and chromatin immunoprecipitation analysis confirmed that Sp1 transcription factor binds to the Sp1 binding site in the porcine Nanog promoter. Taken together, these results show that Sp1 transcription factor is an essential element for porcine Nanog basal activity the same as in human and mouse.

Conclusion

We showed that the porcine Nanog gene is located on porcine chromosome 5 and its basal transcriptional activity is controlled by Sp1 transcription factor.

Conserved roles of mouse DUX and human DUX4 in activating cleavage-stage genes and MERVL/HERVL retrotransposons

To better understand transcriptional regulation during human oogenesis and pre-implantation development, we defined stage-specific transcription, which revealed the cleavage stage as highly distinctive. Here, we present multiple lines of evidence that a eutherian-specific, multi-copy retrogene, DUX4, encodes a transcription factor which activates hundreds of endogenous genes (e.g. ZSCAN4, ZFP352, KDM4E) and retroviral elements (MERVL/HERVL-family) that defines the cleavage-specific transcriptional programs in mouse and human. Remarkably, mouse Dux expression is both necessary and sufficient to convert mouse embryonic stem cells into two-cell embryo-like (‘2C-like’) cells, measured here by the reactivation of ‘2C’ genes and repeat elements, the loss of POU5F1 protein and chromocenters, and by the conversion of the chromatin landscape (assessed by ATAC-seq) to a state strongly resembling mouse two-cell embryos. Taken together, we propose mouse DUX and human DUX4 as major drivers of the cleavage/‘2C’ state.