RNA Binding Properties of SOX Family Members.

SOX proteins are a family of transcription factors (TFs) that play critical functions in sex determination, neurogenesis, and chondrocyte differentiation, as well as cardiac, vascular, and lymphatic development. There are 20 SOX family members in humans, each sharing a 79-residue L-shaped high mobility group (HMG)-box domain that is responsible for DNA binding. SOX2 was recently shown to interact with long non-coding RNA and large-intergenic non-coding RNA to regulate embryonic stem cell and neuronal differentiation. The RNA binding region was shown to reside within the HMG-box domain; however, the structural details of this binding remain unclear. Here, we show that all SOX family members, except group H, interact with RNA. Our mutational experiments demonstrate that the disordered C-terminal region of the HMG-box domain plays an important role in RNA binding. Further, by determining a high-resolution structure of the HMG-box domain of the group H family member SOX30, we show that despite differences in RNA binding ability, SOX30 shares a very similar secondary structure with other SOX protein HMG-box domains. Together, our study provides insight into the interaction of SOX TFs with RNA.

Identification and expression analysis of Sox family genes in echinoderms.

The Sox gene family, a collection of transcription factors widely distributed throughout the animal kingdom, plays a crucial role in numerous developmental processes. Echinoderms occupy a pivotal position in many research fields, such as neuroscience, sex determination and differentiation, and embryonic development. However, to date, no comprehensive study has been conducted to characterize and analyze Sox genes in echinoderms. In the present study, the evolution and expression of Sox family genes across 11 echinoderms were analyzed using bioinformatics methods. The results revealed a total of 70 Sox genes, with counts ranging from 5 to 8 across different echinoderms. Phylogenetic analysis revealed that the identified Sox genes could be categorized into seven distinct classes: the SoxB1 class, SoxB2 class, SoxC class, SoxD class, SoxE class, SoxF class and SoxH class. Notably, the SoxB1, SoxB2, and SoxF genes were ubiquitously present in all the echinoderms studied, which suggests that these genes may be conserved in echinoderms. The spatiotemporal expression patterns observed for Sox genes in the three echinoderms indicated that various Sox members perform distinct functional roles. Notably, SoxB1 is likely involved in echinoderm ovary development, while SoxH may play a crucial role in testis development in starfish and sea cucumber. In general, the present investigation provides a molecular foundation for exploring the Sox gene in echinoderms, providing a valuable resource for future phylogenetic and genomic studies.

HMG-B transcription factors of unicellular opisthokonts and their relatedness to the Sox-Tcf/Lef-Mata proteins of Metazoa and fungi.

A phylogenetic analysis of transcription factors of the Sox-Tcf/Lef-Mata (STM) family of the HMG-B superfamily was carried out in order to clarify the evolutionary roots of the Wnt signaling pathway in unicellular organisms. The data set for analysis included protein sequences of metazoans, fungi, unicellular opisthokonts, apusomonads and amoebozoans. The topology of the phylogenetic tree suggests that STM-related proteins arose in the common ancestor of Opisthokonta and Amoebozoa, two of amoebozoan STM proteins are sister-related to opisthokont ones and the three known lineages of STM transcription factors (STM family in narrow sence) are found in Opisthokonta only. Of these, the holozoan Sox protein branch is the result of either the first or second branching, that originated in the common ancestor of Opisthokonta. The lineage containing Tcf/Lef proteins (holozoan) and the lineage containing Mata proteins (holomycotan) are sister. They derived either at the time of the Holozoa and Holomycota divergence or originate from two paralogs of the common ancestor of Opisthokonta, which arose after the separation of the Sox lineage. Interaction with Armadillo-like proteins may be an original feature of the STM protein family and existed in the unicellular ancestors of multicellular animals; a connection is possible between the presence of Mata-related proteins in Aphelidium protococcorum and specific genome feature of this species.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum).

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

SOX2 and SOX9 Expression in Developing Postnatal Opossum (Monodelphis domestica) Cortex.

(1) Background: Central nervous system (CNS) development is characterized by dynamic changes in cell proliferation and differentiation. Key regulators of these transitions are the transcription factors such as SOX2 and SOX9. SOX2 is involved in the maintenance of progenitor cell state and neural stem cell multipotency, while SOX9, expressed in neurogenic niches, plays an important role in neuron/glia switch with predominant expression in astrocytes in the adult brain. (2) Methods: To validate SOX2 and SOX9 expression patterns in developing opossum (Monodelphis domestica) cortex, we used immunohistochemistry (IHC) and the isotropic fractionator method on fixed cortical tissue from comparable postnatal ages, as well as dissociated primary neuronal cultures. (3) Results: Neurons positive for both neuronal (TUJ1 or NeuN) and stem cell (SOX2) markers were identified, and their presence was confirmed with all methods and postnatal age groups (P4-6, P6-18, and P30) analyzed. SOX9 showed exclusive staining in non-neuronal cells, and it was coexpressed with SOX2. (4) Conclusions: The persistence of SOX2 expression in developing cortical neurons of M. domestica during the first postnatal month implies the functional role of SOX2 during neuronal differentiation and maturation, which was not previously reported in opossums.

SOX family transcription factors as therapeutic targets in wound healing: A comprehensive review.

The SOX family plays a vital role in determining the fate of cells and has garnered attention in the fields of cancer research and regenerative medicine. It also shows promise in the study of wound healing, as it actively participates in the healing processes of various tissues such as skin, fractures, tendons, and the cornea. However, our understanding of the mechanisms behind the SOX family's involvement in wound healing is limited compared to its role in cancer. Gaining insight into its role, distribution, interaction with other factors, and modifications in traumatized tissues could provide valuable new knowledge about wound healing. Based on current research, SOX2, SOX7, and SOX9 are the most promising members of the SOX family for future interventions in wound healing. SOX2 and SOX9 promote the renewal of cells, while SOX7 enhances the microvascular environment. The SOX family holds significant potential for advancing wound healing research. This article provides a comprehensive review of the latest research advancements and therapeutic tools related to the SOX family in wound healing, as well as the potential benefits and challenges of targeting the SOX family for wound treatment.

[The Mechanism of miR-1294 Targeting SOX15 to Regulate Wnt/ß-catenin Signaling Pathway and Promote the Proliferation of Acute Lymphoblastic Leukemia Cells in Children].

OBJECTIVE: To explore the effect of abnormal miRNA expression on the proliferation of pediatric acute lymphoblastic leukemia (ALL) cells and its related mechanism. METHODS: 15 children with ALL and 15 healthy subjects were collected from the Second Affiliated Hospital of Hainan Medical University from July 2018 to March 2021. MiRNA sequencing was performed on their bone marrow cells, and validated using qRT-PCR. MiR-1294 and miR-1294-inhibitory molecule (miR-1294-inhibitor) were transfected into Nalm-6 cells, and the proliferation of Nalm-6 cells was detected by CCK-8 and colony formation assays. Western blot and ELISA were used to detect apoptosis of Nalm-6 cells. Biological prediction of miR-1294 was performed to find the target gene, which was verified by luciferase reporter assay. Si-SOX15 was transfected into Nalm-6 cells, Western blot was used to detect the expression of Wnt signaling pathway-related proteins and to verify the effect of si-SOX15 on the proliferation and apoptosis of Nalm-6 cells. RESULTS: Compared with healthy subjects, 22 miRNAs were significantly upregulated in bone marrow cells of ALL patients, of which miR-1294 was the most significantly upregulated. In addition, the expression level of SOX15 gene was significantly reduced in bone marrow cells of ALL patients. Compared with the NC group, the miR-1294 group showed increased protein expression levels of Wnt3a and ß-catenin, faster cell proliferation, and more colony-forming units, while caspase-3 protein expression level and cell apoptosis were reduced. Compared with the NC group, the miR-1294-inhibitor group showed reduced protein expression levels of Wnt3a and ß-catenin, slower cell proliferation, and fewer colony-forming units, while caspase-3 protein expression level was increased and apoptosis rate was elevated. miR-1294 had a complementary base-pair with the 3'UTR region of SOX15 , and miR-1294 directly targeted SOX15 . The expression of miR-1294 was negatively correlated with SOX15 in ALL cells. Compared with the si-NC group, the si-SOX15 group showed increased protein expression levels of Wnt3a and ß-catenin, accelerated cell proliferation, and decreased caspase-3 protein expression level and cell apoptosis rate. CONCLUSION: MiR-1294 can target and inhibit SOX15 expression, thus activating the Wnt/ß-Catenin signaling pathway to promote the proliferation of ALL cells, inhibit cell apoptosis, and ultimately affect the disease progression.

Single-cell transcriptome profiling reveals distinct expression patterns among genes in the mouse incisor dental pulp.

SOX transcription factors play key roles in cell differentiation and cell fate determination during development. Using single-cell RNA-sequencing data, we examined the expression profiles of Sox genes in the mouse incisor dental pulp. Our analysis showed that Sox4, Sox5, Sox9, Sox11, and Sox12 are mainly expressed in mesenchymal stem/stromal cells (MSCs) representing osteogenic cells at different stages of differentiation. We found that in several MSCs, Sox genes co-expressed with regulatory genes such as Sp7, Satb2, Msx1, Snai2, Dlx1, Twist2, and Tfap2a. In addition, Sox family genes colocalized with Runx2 and Lef1, which are highly enriched in MSCs undergoing osteoblast differentiation. A protein interaction network analysis uncovered that CREBBP, CEBPB, TLE1, TWIST1, and members of the HDAC and SMAD families are interacting partners of RUNX2 and LEF1 during skeletal development. Collectively, the distinct expression patterns of the SOX transcription factors suggest that they play essential regulatory roles in directing lineage-specific gene expression during differentiation of MSCs.

Zebrafish sox2 Is Required for the Swim Bladder Inflation by Controlling the Swim-Up Behavior.

The swim bladder functions to maintain the fish balance at a certain position under water. Although the motoneuron-dependent swim-up behavior is important for swim bladder inflation, the underlying molecular mechanism remains largely unknown. We generated a sox2 KO zebrafish using TALEN and found that the posterior chamber of the swim bladder was uninflated. The tail flick and the swim-up behavior were absent in the mutant zebrafish embryos and the behavior could not be accomplished. As the tail flick behavior is absent, the mutant larvae therefore cannot reach the water surface to gulp air, ultimately leading to the uninflation of the swim bladder. To understand the mechanism underlying the swim-up defects, we crossed the sox2 null allele in the background of Tg(huc:eGFP) and Tg(hb9:GFP). The deficiency of sox2 in zebrafish resulted in abnormal motoneuron axons in the regions of trunk, tail, and swim bladder. To identify the downstream target gene of sox2 to control the motor neuron development, we performed RNA sequencing on the transcriber of mutant embryos versus wild type embryos and found that the axon guidance pathway was abnormal in the mutant embryos. RT-PCR demonstrated that the expression of sema3bl, ntn1b, and robo2 were decreased in the mutants.

Genetic background of idiopathic neurodevelopmental delay patients with significant brain deviation volume.

BACKGROUND: Significant brain volume deviation is an essential phenotype in children with neurodevelopmental delay (NDD), but its genetic basis has not been fully characterized. This study attempted to analyze the genetic factors associated with significant whole-brain deviation volume (WBDV). METHODS: We established a reference curve based on 4222 subjects ranging in age from the first postnatal day to 18 years. We recruited only NDD patients without acquired etiologies or positive genetic results. Cranial magnetic resonance imaging (MRI) and clinical exome sequencing (2742 genes) data were acquired. A genetic burden test was performed, and the results were compared between patients with and without significant WBDV. Literature review analyses and BrainSpan analysis based on the human brain developmental transcriptome were performed to detect the potential role of genetic risk factors in human brain development. RESULTS: We recruited a total of 253 NDD patients. Among them, 26 had significantly decreased WBDV (<-2 standard deviations [SDs]), and 14 had significantly increased WBDV (>+2 SDs). NDD patients with significant WBDV had higher rates of motor development delay (49.8% [106/213] vs . 75.0% [30/40], P  = 0.003) than patients without significant WBDV. Genetic burden analyses found 30 genes with an increased allele frequency of rare variants in patients with significant WBDV. Analyses of the literature further demonstrated that these genes were not randomly identified: burden genes were more related to the brain development than background genes ( P  = 1.656e -9 ). In seven human brain regions related to motor development, we observed burden genes had higher expression before 37-week gestational age than postnatal stages. Functional analyses found that burden genes were enriched in embryonic brain development, with positive regulation of synaptic growth at the neuromuscular junction, positive regulation of deoxyribonucleic acid templated transcription, and response to hormone, and these genes were shown to be expressed in neural progenitors. Based on single cell sequencing analyses, we found TUBB2B gene had elevated expression levels in neural progenitor cells, interneuron, and excitatory neuron and SOX15 had high expression in interneuron and excitatory neuron. CONCLUSION: Idiopathic NDD patients with significant brain volume changes detected by MRI had an increased prevalence of motor development delay, which could be explained by the genetic differences characterized herein.

Evolutionary Landscape of SOX Genes to Inform Genotype-to-Phenotype Relationships.

The SOX transcription factor family is pivotal in controlling aspects of development. To identify genotype-phenotype relationships of SOX proteins, we performed a non-biased study of SOX using 1890 open-reading frame and 6667 amino acid sequences in combination with structural dynamics to interpret 3999 gnomAD, 485 ClinVar, 1174 Geno2MP, and 4313 COSMIC human variants. We identified, within the HMG (High Mobility Group)- box, twenty-seven amino acids with changes in multiple SOX proteins annotated to clinical pathologies. These sites were screened through Geno2MP medical phenotypes, revealing novel SOX15 R104G associated with musculature abnormality and SOX8 R159G with intellectual disability. Within gnomAD, SOX18 E137K (rs201931544), found within the HMG box of ~0.8% of Latinx individuals, is associated with seizures and neurological complications, potentially through blood-brain barrier alterations. A total of 56 highly conserved variants were found at sites outside the HMG-box, including several within the SOX2 HMG-box-flanking region with neurological associations, several in the SOX9 dimerization region associated with Campomelic Dysplasia, SOX14 K88R (rs199932938) flanking the HMG box associated with cardiovascular complications within European populations, and SOX7 A379V (rs143587868) within an SOXF conserved far C-terminal domain heterozygous in 0.716% of African individuals with associated eye phenotypes. This SOX data compilation builds a robust genotype-to-phenotype association for a gene family through more robust ortholog data integration.

The Role of SOX Transcription Factors in Ageing and Age-Related Diseases.

The quest for eternal youth and immortality is as old as humankind. Ageing is an inevitable physiological process accompanied by many functional declines that are driving factors for age-related diseases. Stem cell exhaustion is one of the major hallmarks of ageing. The SOX transcription factors play well-known roles in self-renewal and differentiation of both embryonic and adult stem cells. As a consequence of ageing, the repertoire of adult stem cells present in various organs steadily declines, and their dysfunction/death could lead to reduced regenerative potential and development of age-related diseases. Thus, restoring the function of aged stem cells, inducing their regenerative potential, and slowing down the ageing process are critical for improving the health span and, consequently, the lifespan of humans. Reprograming factors, including SOX family members, emerge as crucial players in rejuvenation. This review focuses on the roles of SOX transcription factors in stem cell exhaustion and age-related diseases, including neurodegenerative diseases, visual deterioration, chronic obstructive pulmonary disease, osteoporosis, and age-related cancers. A better understanding of the molecular mechanisms of ageing and the roles of SOX transcription factors in this process could open new avenues for developing novel strategies that will delay ageing and prevent age-related diseases.

Nanog organizes transcription bodies.

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells.1-3 How proteins of the transcriptional machinery come together to form such bodies, however, is unclear. Here, we take advantage of two large, isolated, and long-lived transcription bodies that reproducibly form during early zebrafish embryogenesis to characterize the dynamics of transcription body formation. Once formed, these transcription bodies are enriched for initiating and elongating RNA polymerase II, as well as the transcription factors Nanog and Sox19b. Analyzing the events leading up to transcription, we find that Nanog and Sox19b cluster prior to transcription. The clustering of transcription factors is sequential; Nanog clusters first, and this is required for the clustering of Sox19b and the initiation of transcription. Mutant analysis revealed that both the DNA-binding domain as well as one of the two intrinsically disordered regions of Nanog are required to organize the two bodies of transcriptional activity. Taken together, our data suggest that the clustering of transcription factors dictates the formation of transcription bodies.

SOX transcription factors direct TCF-independent WNT/ß-catenin responsive transcription to govern cell fate in human pluripotent stem cells.

WNT/ß-catenin signaling controls gene expression across biological contexts from development and stem cell homeostasis to diseases including cancer. How ß-catenin is recruited to distinct enhancers to activate context-specific transcription is unclear, given that most WNT/ß-catenin-responsive transcription is thought to be mediated by TCF/LEF transcription factors (TFs). With time-resolved multi-omic analyses, we show that SOX TFs can direct lineage-specific WNT-responsive transcription during the differentiation of human pluripotent stem cells (hPSCs) into definitive endoderm and neuromesodermal progenitors. We demonstrate that SOX17 and SOX2 are required to recruit ß-catenin to lineage-specific WNT-responsive enhancers, many of which are not occupied by TCFs. At TCF-independent enhancers, SOX TFs establish a permissive chromatin landscape and recruit a WNT-enhanceosome complex to activate SOX/ß-catenin-dependent transcription. Given that SOX TFs and the WNT pathway are critical for specification of most cell types, these results have broad mechanistic implications for the specificity of WNT responses across developmental and disease contexts.

SOX15 transcriptionally increases the function of AOC1 to modulate ferroptosis and progression in prostate cancer.

Amine oxidase copper-containing 1 (AOC1) is considered an oncogene in many types of tumors. Nevertheless, there have been no investigations of AOC1 and its regulatory mechanism in prostate cancer. Here, we reveal a novel action of AOC1 and a tumor suppressor mechanism in prostate cancer. AOC1 is downregulated in prostate cancer. Abatement of AOC1 in prostate cancer tissue is positively correlated with the tumor size, lymph node metastasis, and Gleason score for prostate cancer. Conversely, high expression of AOC1 is significantly associated with reduced proliferation and migration in prostate cancer both in vitro and in vivo. We show that the anticancer effect of AOC1 is mediated by its action on spermidine which leads to the activation of reactive oxygen species and ferroptosis. AOC1 expression in prostate cancer is positively regulated by the transcription factor SOX15. Therefore, SOX15 can transcriptionally promote AOC1 expression and strengthen this effect. Targeting AOC1 and SOX15 may be promising for the treatment of prostate cancer.

Identification of novel genes by targeted exome sequencing in Retinoblastoma.

BACKGROUND: Retinoblastoma (RB) is initiated by mutation in both alleles of RB1 gene. However, few cases may occur even in the absence of RB1 mutation suggesting the role of genes other than RB1. METHODOLOGY: The current study was planned to utilize targeted exome sequencing in Indian RB patients affected with unilateral non-familial RB. 75 unilateral RB patients below 5 years of age were enrolled. Genomic DNA was extracted from blood and tumor tissue. From peripheral blood DNA, all coding and exon/intron regions were amplified using PCR and direct sequencing. Cases which did not harbor pathogenic variants in peripheral blood DNA were further screened for mutations in their tumor tissue DNA using targeted exome sequencing. Three pathogenicity prediction tools (Mutation Taster, SIFT, and PolyPhen-2) were used to determine the pathogenicity of non-synonymous variations. An in-house bioinformatics pipeline was devised for the mutation screening by targeted exome sequencing. Protein modeling studies were also done to predict the effect of the mutations on the protein structure and function. RESULTS: Using the mentioned approach, we found two novel variants (g.69673_69674insT and g.48373314C>A) in RB1 gene in peripheral blood DNA. We also found novel variants in eight genes (RB1, ACAD11, GPR151, KCNA1, OTOR, SOX30, ARL11, and MYCT1) that may be associated with RB pathogenesis. CONCLUSION: The present study expands our current knowledge regarding the genomic landscape of RB and also highlights the importance of NGS technologies to detect genes and novel variants that may play an important role in cancer initiation, progression, and prognosis.

Down Regulation of Catsper1 Expression by Calmodulin Inhibitor (Calmidazolium): Possible Implications for Fertility.

The CatSper channel localizes exclusively in the flagella of sperm cells. The Catsper1 protein, together with three pore units, is essential for the CatSper Channel formation, which produces flagellum hyperactivation and confers sperm fertility. Catsper1 expression is dependent on Sox transcription factors, which can recognize in vitro at least three Sox binding sites on the promoter. Sox transcription factors have calmodulin-binding domains for nuclear importation. Calmodulin (CaM) is affected by the specific inhibitor calmidazolium (CMZ), which prevents the nuclear transport of Sox factors. In this work, we assess the regulation of the Catsper1 promoter in vivo by Sox factors in the murine testis and evaluate the effects of the inhibitor calmidazolium on the expression of the Casper genes, and the motility and fertility of the sperm. Catsper1 promoter has significant transcriptional activity in vivo; on the contrary, three Sox site mutants in the Catsper1 promoter reduced transcriptional activity in the testis. CaM inhibition affects Sox factor nuclear transport and has notable implications in the expression and production of Catsper1, as well as in the motility and fertility capability of sperm. The molecular mechanism described here might conform to the basis of a male contraceptive strategy acting at the transcriptional level by affecting the production of the CatSper channel, a fundamental piece of male fertility.

Pluripotency factors determine gene expression repertoire at zygotic genome activation.

Awakening of zygotic transcription in animal embryos relies on maternal pioneer transcription factors. The interplay of global and specific functions of these proteins remains poorly understood. Here, we analyze chromatin accessibility and time-resolved transcription in single and double mutant zebrafish embryos lacking pluripotency factors Pou5f3 and Sox19b. We show that two factors modify chromatin in a largely independent manner. We distinguish four types of direct enhancers by differential requirements for Pou5f3 or Sox19b. We demonstrate that changes in chromatin accessibility of enhancers underlie the changes in zygotic expression repertoire in the double mutants. Pou5f3 or Sox19b promote chromatin accessibility of enhancers linked to the genes involved in gastrulation and ventral fate specification. The genes regulating mesendodermal and dorsal fates are primed for activation independently of Pou5f3 and Sox19b. Strikingly, simultaneous loss of Pou5f3 and Sox19b leads to premature expression of genes, involved in regulation of organogenesis and differentiation.

The landscape of pioneer factor activity reveals the mechanisms of chromatin reprogramming and genome activation.

Upon fertilization, embryos undergo chromatin reprogramming and genome activation; however, the mechanisms that regulate these processes are poorly understood. Here, we generated a triple mutant for Nanog, Pou5f3, and Sox19b (NPS) in zebrafish and found that NPS pioneer chromatin opening at >50% of active enhancers. NPS regulate acetylation across core histones at enhancers and promoters, and their function in gene activation can be bypassed by recruiting histone acetyltransferase to individual genes. NPS pioneer chromatin opening individually, redundantly, or additively depending on sequence context, and we show that high nucleosome occupancy facilitates NPS pioneering activity. Nucleosome position varies based on the input of different transcription factors (TFs), providing a flexible platform to modulate pioneering activity. Altogether, our results illuminate the sequence of events during genome activation and offer a conceptual framework to understand how pioneer factors interpret the genome and integrate different TF inputs across cell types and developmental transitions.

Antibody-free profiling of transcription factor occupancy during early embryogenesis by FitCUT&RUN.

Key transcription factors (TFs) play critical roles in zygotic genome activation (ZGA) during early embryogenesis, whereas genome-wide occupancies of only a few factors have been profiled during ZGA due to the limitation of cell numbers or the lack of high-quality antibodies. Here, we present FitCUT&RUN, a modified CUT&RUN method, in which an Fc fragment of immunoglobulin G is used for tagging, to profile TF occupancy in an antibody-free manner and demonstrate its reliability and robustness using as few as 5000 K562 cells. We applied FitCUT&RUN to zebrafish undergoing embryogenesis to generate reliable occupancy profiles of three known activators of zebrafish ZGA: Nanog, Pou5f3, and Sox19b. By profiling the time-series occupancy of Nanog during zebrafish ZGA, we observed a clear trend toward a gradual increase in Nanog occupancy and found that Nanog occupancy prior to the major phase of ZGA is important for the activation of some early transcribed genes.