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.

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.

[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.

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.

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.

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.

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.

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.

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.

Differential expression of members of SOX family of transcription factors in failing human hearts.

The Sry-related high-mobility-group box (SOX) gene family, with 20 known transcription factors in humans, plays an essential role during development and disease processes. Several SOX proteins (SOX4, 11, and 9) are required for normal heart morphogenesis. SOX9 was shown to contribute to cardiac fibrosis. However, differential expression of other SOXs and their roles in the failing human myocardium have not been explored. Here, we used the whole-transcriptome sequencing (RNA-seq), gene co-expression, and meta-analysis to examine whether any SOX factors might play a role in the failing human myocardium. RNA-seq analysis was performed for cardiac tissue samples from heart failure (HF) patients due to dilated cardiomyopathy (DCM), or hypertrophic cardiomyopathy (HCM) and healthy donors (NF). The RNA levels of 20 SOX genes from RNA-seq data were extracted and compared to the 3 groups. Four SOX genes whose RNA levels were significantly upregulated in DCM or HCM compared to NF. However, only SOX4 and SOX8 proteins were markedly increased in the HF groups. A moderate to strong correlation was observed between the RNA level of SOX4/8 and fibrotic genes among each individual. Gene co-expression network analysis identified genes associated and respond similarly to perturbations with SOX4 in cardiac tissues. Using a meta-analysis combining epigenetics and genome-wide association data, we reported several genomic variants associated with HF phenotype linked to SOX4 or SOX8. In summary, our results implicate that SOX4 and SOX8 have a role in cardiomyopathy, leading to HF in humans. The molecular mechanism associated with them in HF warrants further investigation.

Sox-positive cell population in the adult cerebellum increases upon tissue degeneration.

Adult neurogenesis is well-described in the subventricular and subgranular zones of the mammalian brain. Recent observations that resident glia express stem cell markers in some areas of the brain not traditionally associated with neurogenesis hint to a possible role in tissue repair. The Bergmann glia (BG) population in the cerebellum displays markers and in vitro features associated with neural stem cells (NSC), however the physiological relevance of this phenotypic overlap remains unclear in the absence of established in vivo evidence of tissue regeneration in the adult cerebellum. Here, this BG population was analysed in the adult cerebellum of different species and showed conservation of NSC-associated marker expression including Sox1, Sox2 and Sox9, in chick, primate and mouse cerebellum tissue. NSC-like cells isolated from adult mouse cerebellum showed slower growth when compared to lateral ventricle NSC, as well as differences upon differentiation. In a mouse model of cerebellar degeneration, progressive Purkinje cell loss was linked to cerebellar cortex disorganisation and a significant increase in Sox-positive cells compared to matching controls. These results show that this Sox-positive population responds to cerebellar tissue disruption, suggesting it may represent a mobilisable cellular resource for targeted strategies to promote tissue repair.

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.

Unraveling the genotype by environment interaction in a thermosensitive fish with a polygenic sex determination system.

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to high-temperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.