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

Molecular cloning of PRD-like homeobox genes expressed in bovine oocytes and early IVF embryos

BACKGROUND

Embryonic genome activation (EGA) is a critical step in early embryonic development, as it marks the transition from relying on maternal factors to the initiation of transcription from embryo's own genome. The factors associated with EGA are not well understood and need further investigation. PRD-like (PRDL) homeodomain transcription factors (TFs) are considered to play crucial roles in this early event during development but these TFs have evolved differently, even within mammalian lineages. Different numbers of PRDL TFs have been predicted in bovine (Bos taurus); however, their divergent evolution requires species-specific confirmation and functional investigations.

RESULTS

In this study, we conducted molecular cloning of mRNAs for the PRDL TFs ARGFX, DUXA, LEUTX, NOBOX, TPRX1, TPRX2, and TPRX3 in bovine oocytes or in vitro fertilized (IVF) preimplantation embryos. Our results confirmed the expression of PRDL TF genes in early bovine development at the cDNA level and uncovered their structures. For each investigated PRDL TF gene, we isolated at least one homeodomain-encoding cDNA fragment, indicative of DNA binding and thus potential role in transcriptional regulation in developing bovine embryos. Additionally, our cDNA cloning approach allowed us to reveal breed-related differences in bovine, as evidenced by the identification of a high number of single nucleotide variants (SNVs) across the PRDL class homeobox genes. Subsequently, we observed the prediction of the 9aa transactivation domain (9aaTAD) motif in the putative protein sequence of TPRX3 leading us to conduct functional analysis of this gene. We demonstrated that the TPRX3 overexpression in bovine fibroblast induces not only protein-coding genes but also short noncoding RNAs involved in splicing and RNA editing. We supported this finding by identifying a shared set of genes between our and published bovine early embryo development datasets.

CONCLUSIONS

Providing full-length cDNA evidence for previously predicted homeobox genes that belong to PRDL class improves the annotation of the bovine genome. Updating the annotation with seven developmentally-important genes will enhance the accuracy of RNAseq analysis with datasets derived from bovine preimplantation embryos. In addition, the absence of TPRX3 in humans highlights the species-specific and TF-specific regulation of biological processes during early embryo development.

Rise and SINE: roles of transcription factors and retrotransposons in zygotic genome activation

In sexually reproducing organisms, life begins with the fusion of transcriptionally silent gametes, the oocyte and sperm. Although initiation of transcription in the embryo, known as zygotic genome activation (ZGA), is universally required for development, the transcription factors regulating this process are poorly conserved. In this Perspective, we discuss recent insights into the mechanisms of ZGA in totipotent mammalian embryos, namely ZGA regulation by several transcription factors, including by orphan nuclear receptors (OrphNRs) such as the pioneer transcription factor NR5A2, and by factors of the DUX, TPRX and OBOX families. We performed a meta-analysis and compiled a list of pan-ZGA genes, and found that most of these genes are indeed targets of the above transcription factors. Remarkably, more than a third of these ZGA genes appear to be regulated both by OrphNRs such as NR5A2 and by OBOX proteins, whose motifs co-occur in SINE B1 retrotransposable elements, which are enriched near ZGA genes. We propose that ZGA in mice is activated by recruitment of multiple transcription factors to SINE B1 elements that function as enhancers, and discuss a potential relevance of this mechanism to Alu retrotransposable elements in human ZGA.

Kick-starting the zygotic genome: licensors, specifiers, and beyond

Zygotic genome activation (ZGA), the first transcription event following fertilization, kickstarts the embryonic program that takes over the control of early development from the maternal products. How ZGA occurs, especially in mammals, is poorly understood due to the limited amount of research materials. With the rapid development of single-cell and low-input technologies, remarkable progress made in the past decade has unveiled dramatic transitions of the epigenomes, transcriptomes, proteomes, and metabolomes associated with ZGA. Moreover, functional investigations are yielding insights into the key regulators of ZGA, among which two major classes of players are emerging: licensors and specifiers. Licensors would control the permission of transcription and its timing during ZGA. Accumulating evidence suggests that such licensors of ZGA include regulators of the transcription apparatus and nuclear gatekeepers. Specifiers would instruct the activation of specific genes during ZGA. These specifiers include key transcription factors present at this stage, often facilitated by epigenetic regulators. Based on data primarily from mammals but also results from other species, we discuss in this review how recent research sheds light on the molecular regulation of ZGA and its executors, including the licensors and specifiers.

LEUTX regulates porcine embryonic genome activation in somatic cell nuclear transfer embryos

Emerging evidence highlights the regulatory role of paired-like (PRD-like) homeobox transcription factors (TFs) in embryonic genome activation (EGA). However, the majority of PRD-like genes are lost in rodents, thus prompting an investigation into PRD-like TFs in other mammals. Here, we showed that PRD-like TFs were transiently expressed during EGA in human, monkey, and porcine fertilized embryos, yet they exhibited inadequate expression in their cloned embryos. This study, using pig as the research model, identified LEUTX as a key PRD-like activator of porcine EGA through genomic profiling and found that LEUTX overexpression restored EGA failure and improved preimplantation development and cloning efficiency in porcine cloned embryos. Mechanistically, LEUTX opened EGA-related genomic regions and established histone acetylation via recruiting acetyltransferases p300 and KAT2A. These findings reveal the regulatory mechanism of LEUTX to govern EGA in pigs, which may provide valuable insights into the study of early embryo development for other non-rodent mammals.

Interaction network of human early embryonic transcription factors

Embryonic genome activation (EGA) occurs during preimplantation development and is characterized by the initiation of de novo transcription from the embryonic genome. Despite its importance, the regulation of EGA and the transcription factors involved in this process are poorly understood. Paired-like homeobox (PRDL) family proteins are implicated as potential transcriptional regulators of EGA, yet the PRDL-mediated gene regulatory networks remain uncharacterized. To investigate the function of PRDL proteins, we are identifying the molecular interactions and the functions of a subset family of the Eutherian Totipotent Cell Homeobox (ETCHbox) proteins, seven PRDL family proteins and six other transcription factors (TFs), all suggested to participate in transcriptional regulation during preimplantation. Using mass spectrometry-based interactomics methods, AP-MS and proximity-dependent biotin labeling, and chromatin immunoprecipitation sequencing we derive the comprehensive regulatory networks of these preimplantation TFs. By these interactomics tools we identify more than a thousand high-confidence interactions for the 21 studied bait proteins with more than 300 interacting proteins. We also establish that TPRX2, currently assigned as pseudogene, is a transcriptional activator.

OBOX regulates mouse zygotic genome activation and early development

Zygotic genome activation (ZGA) activates the quiescent genome to enable the maternal-to-zygotic transition1,2. However, the identity of transcription factors that underlie mammalian ZGA in vivo remains elusive. Here we show that OBOX, a PRD-like homeobox domain transcription factor family (OBOX1-OBOX8)3-5, are key regulators of mouse ZGA. Mice deficient for maternally transcribed Obox1/2/5/7 and zygotically expressed Obox3/4 had a two-cell to four-cell arrest, accompanied by impaired ZGA. The Obox knockout defects could be rescued by restoring either maternal and zygotic OBOX, which suggests that maternal and zygotic OBOX redundantly support embryonic development. Chromatin-binding analysis showed that Obox knockout preferentially affected OBOX-binding targets. Mechanistically, OBOX facilitated the 'preconfiguration' of RNA polymerase II, as the polymerase relocated from the initial one-cell binding targets to ZGA gene promoters and distal enhancers. Impaired polymerase II preconfiguration in Obox mutants was accompanied by defective ZGA and chromatin accessibility transition, as well as aberrant activation of one-cell polymerase II targets. Finally, ectopic expression of OBOX activated ZGA genes and MERVL repeats in mouse embryonic stem cells. These data thus demonstrate that OBOX regulates mouse ZGA and early embryogenesis.

The NFkB activation domain is 14-amino-acid-long variant of the 9aaTAD

The nine-amino-acid transactivation domains (9aaTAD) was identified in numerous transcription factors including Gal4, p53, E2A, MLL, c-Myc, N-Myc, and also in SP, KLF, and SOX families. Most of the 9aaTAD domains interact with the KIX domain of transcription mediators MED15 and CBP to activate transcription. The NFkB activation domain occupied the same position on the KIX domain as the 9aaTADs of MLL, E2A, and p53. Binding of the KIX domain is established by the two-point interaction involving 9aaTAD positions p3-4 and p6-7. The NFkB primary binding region (positions p3-4) is almost identical with MLL and E2A, but secondary NFkB binding region differs by the position and engages the distal NFkB region p10-11. Thus, the NFkB activation domain is five amino acids longer than the other 9aaTADs. The NFkB activation domain includes an additional region, which we called the Omichinski Insert extending activation domain length to 14 amino acids. By deletion, we demonstrated that Omichinski Insert is an entirely non-essential part of NFkB activation domain. In summary, we recognized the NFkB activation domain as prolonged 9aaTAD conserved in evolution from humans to amphibians.

Pediatric-type high-grade neuroepithelial tumors with CIC gene fusion share a common DNA methylation signature

Pediatric neoplasms in the central nervous system (CNS) show extensive clinical and molecular heterogeneity and are fundamentally different from those occurring in adults. Molecular genetic testing contributes to accurate diagnosis and enables an optimal clinical management of affected children. Here, we investigated a rare, molecularly distinct type of pediatric high-grade neuroepithelial tumor (n = 18), that was identified through unsupervised visualization of genome-wide DNA methylation array data, together with copy number profiling, targeted next-generation DNA sequencing, and RNA transcriptome sequencing. DNA and/or RNA sequencing revealed recurrent fusions involving the capicua transcriptional repressor (CIC) gene in 10/10 tumor samples analyzed, with the most common fusion being CIC::LEUTX (n = 9). In addition, a CIC::NUTM1 fusion was detected in one of the tumors. Apart from the detected fusion events, no additional oncogenic alteration was identified in these tumors. The histopathological review demonstrated a morphologically heterogeneous group of high-grade neuroepithelial tumors with positive immunostaining for markers of glial differentiation in combination with weak and focal expression of synaptophysin, CD56 and CD99. All tumors were located in the supratentorial compartment, occurred during childhood (median age 8.5 years) and typically showed early relapses. In summary, we expand the spectrum of pediatric-type tumors of the CNS by reporting a previously uncharacterized group of rare high-grade neuroepithelial tumors that share a common DNA methylation signature and recurrent gene fusions involving the transcriptional repressor CIC. Downstream functional consequences of the fusion protein CIC::LEUTX and potential therapeutic implications need to be further investigated.

Comprehensive characterization of the embryonic factor LEUTX

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

Translatome and transcriptome co-profiling reveals a role of TPRXs in human zygotic genome activation

Translational regulation plays a critical role during the oocyte-to-embryo transition (OET) and zygotic genome activation (ZGA). Here, we integrated ultra-low-input Ribo-seq with mRNA-seq to co-profile the translatome and transcriptome in human oocytes and early embryos. Comparison with mouse counterparts identified widespread differentially translated genes functioning in epigenetic reprogramming, transposon defense, and small RNA biogenesis, in part driven by species-specific regulatory elements in 3' untranslated regions. Moreover, PRD-like homeobox transcription factors, including TPRXL, TPRX1, and TPRX2, are highly translated around ZGA. TPRX1/2/L knockdown leads to defective ZGA and preimplantation development. Ectopically expressed TPRXs bind and activate key ZGA genes in human embryonic stem cells. These data reveal the conservation and divergence of translation landscapes during OET and identify critical regulators of human ZGA.

PRD-class homeobox genes in bovine early embryos: function, evolution and overlapping roles

Eutherian Totipotent Cell Homeobox (ETCHbox) genes are mammalian-specific PRD-class homeobox genes with conserved expression in the preimplantation embryo but fast-evolving and highly divergent sequences. Here, we exploit an ectopic expression approach to examine the role of bovine ETCHbox genes and show that ARGFX and LEUTX homeodomain proteins upregulate genes normally expressed in the blastocyst; the identities of the regulated genes suggest that, in vivo, the ETCHbox genes play a role in coordinating the physical formation of the blastocyst structure. Both genes also downregulate genes expressed earlier during development and genes associated with an undifferentiated cell state, possibly via the JAK/STAT pathway. We find evidence that bovine ARGFX and LEUTX have overlapping functions, in contrast to their antagonistic roles in humans. Finally, we characterize a mutant bovine ARGFX allele which eliminates the homeodomain and show that homozygous mutants are viable. These data support the hypothesis of functional overlap between ETCHbox genes within a species, roles for ETCHbox genes in blastocyst formation and the change of their functions over evolutionary time.

DUX4 is a multifunctional factor priming human embryonic genome activation

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Dynamic Molecular Evolution of Mammalian Homeobox Genes: Duplication, Loss, Divergence and Gene Conversion Sculpt PRD Class Repertoires

The majority of homeobox genes are highly conserved across animals, but the eutherian-specific ETCHbox genes, embryonically expressed and highly divergent duplicates of CRX, are a notable exception. Here we compare the ETCHbox genes of 34 mammalian species, uncovering dynamic patterns of gene loss and tandem duplication, including the presence of a large tandem array of LEUTX loci in the genome of the European rabbit (Oryctolagus cuniculus). Despite extensive gene gain and loss, all sampled species possess at least two ETCHbox genes, suggesting their collective role is indispensable. We find evidence for positive selection and show that TPRX1 and TPRX2 have been the subject of repeated gene conversion across the Boreoeutheria, homogenising their sequences and preventing divergence, especially in the homeobox region. Together, these results are consistent with a model where mammalian ETCHbox genes are dynamic in evolution due to functional overlap, yet have collective indispensable roles.

A novel BRD4-LEUTX fusion in a pediatric sarcoma with epithelioid morphology and diffuse S100 expression

Malignant epithelioid soft tissue tumors encompass a wide spectrum of lesions. Among them, Epithelioid Malignant Peripheral Nerve Sheath Tumors (MPNST) constitute a distinct subgroup, accounting for <5% of all MPNST. Epithelioid MPNST are infrequently associated with neurofibromatosis type 1, occasionally arise in a schwannoma and show diffuse S100 and CD34 expression, often combined with INI-1 loss. However, the molecular mechanisms underlying the tumorigenesis of epithelioid MPNST remain largely unknown. We describe a case of a 10-year-old girl with an epithelioid malignancy of the orbit. The tumor proved positive for S100, CD34 and SOX10, and, although INI-1 expression was maintained, the overall features suggested the possibility of an epithelioid MPNST, arising in an unusual location. NGS analysis revealed a novel in-frame BRD4-LEUTX fusion gene. LEUTX plays an important role in embryonal genome activation and its expression is mostly suppressed postnatally. We were able to detect increased levels of LEUTX transcript in the tumor, indicating that BRD4-LEUTX fusion leads to LEUTX re-activation. To our knowledge, this fusion has never been reported previously. Whether the current case represents an example of epithelioid MPNST or a distinct tumor entity remains to be determined.

Relationship of DUX4 and target gene expression in FSHD myocytes

Facioscapulohumeral dystrophy (FSHD) is associated with the upregulation of the DUX4 transcription factor and its target genes. However, low-frequency DUX4 upregulation in patient myocytes is difficult to detect and examining the relationship and dynamics of DUX4 and target gene expression has been challenging. Using RNAScope in situ hybridization with highly specific probes, we detect the endogenous DUX4 and target gene transcripts in situ in patient skeletal myotubes during 13-day differentiation in vitro. We found that the endogenous DUX4 transcripts primarily localize as foci in one or two nuclei as compared with the accumulation of the recombinant DUX4 transcripts in the cytoplasm. We also found the continuous increase of DUX4 and target gene-positive myotubes after Day 3, arguing against its expected immediate cytotoxicity. Interestingly, DUX4 and target gene expression become discordant later in differentiation with the increase of DUX4-positive/target gene-negative as well as DUX4-negative/target gene-positive myotubes. Depletion of DUX4-activated transcription factors, DUXA and LEUTX, specifically repressed a DUX4-target gene, KDM4E, later in differentiation, suggesting that after the initial activation by DUX4, target genes themselves contribute to the maintenance of downstream gene expression. Together, the study provides important new insights into the dynamics of the DUX4 transcriptional network in FSHD patient myocytes.