On the NF-Y regulome as in ENCODE (2019)

The role(s) of NF-Y in development and differentiation

NF-Y is a conserved sequence-specific trimeric Transcription Factor -TF- binding to the CCAAT element. We review here the role(s) in development, from pre-implantation embryo to terminally differentiated tissues, by rationalizing and commenting on genetic, genomic, epigenetic and biochemical studies. This effort brings to light the impact of NF-YA isoforms on stemness and differentiation, as well as binding to distal vs promoter proximal sites and connections with selected TFs.

Genomic binding of NF-Y in mouse and human cells

NF-Y is a Transcription Factor that regulates transcription through binding to the CCAAT-box. To understand its strategy, we analyzed 16 ChIP-seq datasets from human and mouse cells. Shared loci, mostly located in promoters of expressed genes of cell cycle, metabolism and gene expression pathways, are associated with histone marks of active chromatin and specific modules of TFs. Other peaks are in enhancers and Transposable Elements -TE- of retroviral origin in human and mouse. We evaluated the relationship with USF1, a common synergistic partner in promoters and MLT1 TEs, upon NF-YB inactivation: USF1 binding decreases in promoters, modestly in MLT1, suggesting a pioneering role of NF-Y in formers, not in the latters. These data define a common set of NF-Y functional targets across different mammalian cell types, suggesting a pioneering role in promoters with respect to TEs.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development

STUDY QUESTION

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development?

SUMMARY ANSWER

Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce.

STUDY DESIGN, SIZE, DURATION

This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome.

MAIN RESULTS AND THE ROLE OF CHANCE

A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation.

LARGE SCALE DATA

All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501.

LIMITATIONS, REASONS FOR CAUTION

Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells.

WIDER IMPLICATIONS OF THE FINDINGS

The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

NF-YA isoforms with alternative splicing of exon-5 in Aves

NF-YA, the regulatory subunit of the trimeric CCAAT-binding transcription factor NF-Y, is present in vertebrates in two major alternative spliced isoforms: NF-YAl and NF-YAs, differing for the presence of exon-3. NF-YAx, a third isoform without exon-3/-5, was reported only in human neuronal cells and tumors. These events affect the Trans-Activation Domain. We provide here evidence for the expression of NF-YAx and for the existence of a new isoform, NF-YAg, skipping only exon-5. These isoforms are abundant in Aves, but not in reptiles, and are the prevalent transcripts in the initial phases of embryo development in chicken. Finally, we analyzed NF-YAg and NF-YAx amino acid sequence using AlphaFold: absence of exon-5 denotes a global reduction of β-stranded elements, while removal of the disordered exon-3 sequence has limited effects on TAD architecture. These data identify an expanded program of NF-YA isoforms within the TAD in Aves, implying a role during early development.

The NF-Y splicing signature controls hybrid EMT and ECM-related pathways to promote aggressiveness of colon cancer

Aberrant splicing events are associated with colorectal cancer (CRC) and provide new opportunities for tumor diagnosis and treatment. The expression of the splice variants of NF-YA, the DNA binding subunit of the transcription factor NF-Y, is deregulated in multiple cancer types compared to healthy tissues. NF-YAs and NF-YAl isoforms differ in the transactivation domain, which may result in distinct transcriptional programs. In this study, we demonstrated that the NF-YAl transcript is higher in aggressive mesenchymal CRCs and predicts shorter patients' survival. In 2D and 3D conditions, CRC cells overexpressing NF-YAl (NF-YAl^high) exhibit reduced cell proliferation, rapid single cell amoeboid-like migration, and form irregular spheroids with poor cell-to-cell adhesion. Compared to NF-YAs^high, NF-YAl^high cells show changes in the transcription of genes involved in epithelial-mesenchymal transition, extracellular matrix and cell adhesion. NF-YAl and NF-YAs bind similarly to the promoter of the E-cadherin gene, but oppositely regulate its transcription. The increased metastatic potential of NF-YAl^high cells in vivo was confirmed in zebrafish xenografts. These results suggest that the NF-YAl splice variant could be a new CRC prognostic factor and that splice-switching strategies may reduce metastatic CRC progression.

Phylogeny of NF-YA trans-activation splicing isoforms in vertebrate evolution

NF-Y is a trimeric pioneer Transcription Factor (TF) whose target sequence -the CCAAT box- is present in ~25% of mammalian promoters. We reconstruct the phylogenetic history of the regulatory NF-YA subunit in vertebrates. We find that in addition to the remarkable conservation of the subunits-interaction and DNA-binding parts, the Transcriptional Activation Domain (TAD) is also conserved (>90% identity among bony vertebrates). We infer the phylogeny of the alternatively spliced exon-3 and partial splicing events of exon-7 -7N and 7C- revealing independent clade-specific losses of these regions. These isoforms shape the TAD. Absence of exon-3 in basal deuterostomes, cartilaginous fishes and hagfish, but not in lampreys, suggests that the "short" isoform is primordial, with emergence of exon-3 in chordates. Exon 7N was present in the vertebrate common ancestor, while 7C is a molecular innovation of teleost fishes. RNA-seq analysis in several species confirms expression of all these isoforms. We identify 3 blocks of amino acids in the TAD shared across deuterostomes, yet structural predictions and sequence analyses suggest an evolutionary drive for maintenance of an Intrinsically Disordered Region -IDR- within the TAD. Overall, these data help reconstruct the logic for alternative splicing of this essential eukaryotic TF.

The USR domain of USF1 mediates NF-Y interactions and cooperative DNA binding

The trimeric CCAAT-binding NF-Y is a "pioneer" Transcription Factor -TF- known to cooperate with neighboring TFs to regulate gene expression. Genome-wide analyses detected a precise stereo-alignment -10/12 bp- of CCAAT with E-box elements and corresponding colocalization of NF-Y with basic-Helix-Loop-Helix (bHLH) TFs. We dissected here NF-Y interactions with USF1 and MAX. USF1, but not MAX, cooperates in DNA binding with NF-Y. NF-Y and USF1 synergize to activate target promoters. Reconstruction of complexes by structural means shows independent DNA binding of MAX, whereas USF1 has extended contacts with NF-Y, involving the USR, a USF-specific amino acid sequence stretch required for trans-activation. The USR is an intrinsically disordered domain and adopts different conformations based on E-box-CCAAT distances. Deletion of the USR abolishes cooperative DNA binding with NF-Y. Our data indicate that the functionality of certain unstructured domains involves adapting to small variation in stereo-alignments of the multimeric TFs sites.

Genome-wide screening and identification of nuclear Factor-Y family genes and exploration their function on regulating abiotic and biotic stress in potato (Solanum tuberosum L.)

The Nuclear Factor-Y (NF-Y) transcription factor (TF), which includes three distinct subunits (NF-YA, NF-YB and NF-YC), is known to manipulate various aspects of plant growth, development, and stress responses. Although the NF-Y gene family was well studied in many species, little is known about their functions in potato. In this study, a total of 37 potato NF-Y genes were identified, including 11 StNF-YAs, 20 StNF-YBs, and 6 StNF-YCs. The genetic features of these StNF-Y genes were investigated by comparing their evolutionary relationship, intron/exon organization and motif distribution pattern. Multiple alignments showed that all StNF-Y proteins possessed clearly conserved core regions that were flanked by non-conserved sequences. Gene duplication analysis indicated that nine StNF-Y genes were subjected to tandem duplication and eight StNF-Ys arose from segmental duplication events. Synteny analysis suggested that most StNF-Y genes (33 of 37) were orthologous to potato's close relative tomato (Solanum lycopersicum L.). Tissue-specific expression of the StNF-Y genes suggested their potential roles in controlling potato growth and development. The role of StNF-Ys in regulating potato responses to abiotic stress (ABA, drought and salinity) was also confirmed: twelve StNF-Y genes were up-regulated and another two were down-regulated under different abiotic treatments. In addition, genes responded differently to pathogen challenges, suggesting that StNF-Y genes may play distinct roles under certain biotic stress. In summary, insights into the evolution of NF-Y family members and their functions in potato development and stress responses are provided.