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

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.

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.

Obox4 promotes zygotic genome activation upon loss of Dux

Once fertilized, mouse zygotes rapidly proceed to zygotic genome activation (ZGA), during which long terminal repeats (LTRs) of murine endogenous retroviruses with leucine tRNA primer (MERVL) are activated by a conserved homeodomain-containing transcription factor, DUX. However, Dux-knockout embryos produce fertile mice, suggesting that ZGA is redundantly driven by an unknown factor(s). Here, we present multiple lines of evidence that the multicopy homeobox gene, Obox4, encodes a transcription factor that is highly expressed in mouse two-cell embryos and redundantly drives ZGA. Genome-wide profiling revealed that OBOX4 specifically binds and activates MERVL LTRs as well as a subset of murine endogenous retroviruses with lysine tRNA primer (MERVK) LTRs. Depletion of Obox4 is tolerated by embryogenesis, whereas concomitant Obox4/Dux depletion markedly compromises embryonic development. Our study identified OBOX4 as a transcription factor that provides genetic redundancy to preimplantation development.

Detection of newly synthesized RNA reveals transcriptional reprogramming during ZGA and a role of Obox3 in totipotency acquisition

Zygotic genome activation (ZGA) after fertilization enables the maternal-to-zygotic transition. However, the global view of ZGA, particularly at initiation, is incompletely understood. Here, we develop a method to capture and sequence newly synthesized RNA in early mouse embryos, providing a view of transcriptional reprogramming during ZGA. Our data demonstrate that major ZGA gene activation begins earlier than previously thought. Furthermore, we identify a set of genes activated during minor ZGA, the promoters of which show enrichment of the Obox factor motif, and find that Obox3 or Obox5 overexpression in mouse embryonic stem cells activates ZGA genes. Notably, the expression of Obox factors is severely impaired in somatic cell nuclear transfer (SCNT) embryos, and restoration of Obox3 expression corrects the ZGA profile and greatly improves SCNT embryo development. Hence, our study reveals dynamic transcriptional reprogramming during ZGA and underscores the crucial role of Obox3 in facilitating totipotency acquisition.

Role of maternal spiralian-specific homeobox gene SPILE-E in the specification of blastomeres along the animal-vegetal axis during the early cleavage stages of mollusks

Spiralians, one of the major clades of bilaterians, share a unique development known as spiralian development, characterized by the formation of tiers of cells called quartets, which exhibit different developmental potentials along the animal-vegetal axis. Recently, spiralian-specific TALE-type homeobox genes (SPILE) have been identified, some of which show zygotic and staggered expression patterns along the animal-vegetal axis and function in quartet specification in mollusks. However, it is unclear which maternal molecular components control the zygotic expression of these transcription factors. In this study, we focused on SPILE-E, a maternal transcription factor, and investigated its expression and function in mollusks. We found that the maternal and ubiquitous expression of SPILE-E in the cleavage stages is conserved in molluskan species, including limpets, mussels, and chitons. We knocked down SPILE-E in limpets and revealed that the expression of transcription factors specifically expressed in the first quartet (1q2 ; foxj1b) and second quartet (2q; SPILE-B) was abolished, whereas the macromere-quartet marker (SPILE-C) was ectopically expressed in 1q2 in SPILE-E morphants. Moreover, we showed that the expression of SPILE-A, which upregulates SPILE-B but represses SPILE-C expression, decreased in SPILE-E morphants. Consistent with changes in the expression pattern of the above transcription factors, SPILE-E-morphant larvae exhibited patchy or complete loss of expression of marker genes of ciliated cells and shell fields, possibly reflecting incomplete specification of 1q2 and 2q. Our results provide a molecular framework for quartet specification and highlight the importance of maternal lineage-specific transcription factors in the development and evolution of spiralians.

Human 8-cell embryo-like cells from pluripotent stem cells

The totipotent embryo initiates transcription during zygotic or embryonic genome activation (EGA, ZGA). ZGA occurs at the 8-cell stage in humans and its failure leads to developmental arrest. Understanding the molecular pathways underlying ZGA and totipotency is essential to comprehend human development. Recently, human 8-cell-like cells (8CLCs) have been discovered in vitro that resemble the 8-cell embryo. 8CLCs exist among naive pluripotent stem cells and can be induced genetically or chemically. Their ZGA-like transcriptome, transposable element activation, 8-cell embryo-specific protein expression, and developmental properties make them an exceptional model system to study early embryonic cell-state transitions and human totipotency programs in vitro.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

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.

Dynamic evolutionary history of spiralian-specific TALE homeobox genes in mollusks

Homeobox genes play essential roles in the early development of many animals. Although the repertoire of most homeobox genes, including three amino acid loop extension (TALE)-type homeobox genes, is conserved in animals, spiralian-TALE (SPILE) genes are a notable exception. In this study, SPILE genes were extracted from the genomic data of 22 mollusk species and classified into four clades (-A/C, -B, -D, and -E) to determine which SPILE genes exhibit dynamic repertoire changes. While SPILE-D and -E duplications were rarely observed, SPILE-B duplication was observed in the bivalve lineage and SPILE-A/C duplication was observed in multiple clades. Conversely, most or all SPILE genes were lost in cephalopods and in some gastropod lineages. SPILE gene expression patterns were also analyzed in multiple mollusk species using publicly available RNA-seq data. The majority of SPILE genes examined, particularly those in the A/C- and B-clades, were specifically expressed during early development, suggesting that most SPILE genes exert specific roles in early development. This comprehensive cataloging and characterization revealed a dynamic evolutionary history, including SPILE-A/C and -B gene duplications and the loss of SPILE genes in several lineages. Furthermore, this study provides a useful resource for studying the molecular mechanism of spiralian early development and the evolution of young and lineage-specific transcription factors.