Differentiation

Background: In humans, haploinsufficiency of TWIST1 leads to Saethre-Chotzen Syndrome, characterized by limb abnormalities and craniofacial defects. In mouse, Twist1 promotes limb bud outgrowth via maintenance of the apical ectodermal ridge (AER), in addition to its craniofacial roles. In zebrafish, two twist1 paralogs, twist1a and twist1b, are known to influence cranial neural crest and cranial suture development but have not yet been implicated in limb formation. Results: We find that twist1a and twist1b act redundantly during zebrafish pectoral fin development. While twist1b mutants have normal fins, twist1a mutants have small fins, and twist1a;twist1b double mutants lack fins. These defects appear to originate with impaired establishment of the pectoral fin field, followed by failure of AER formation, likely due to reduced fgf10a expression in the fin mesenchyme. Anterior-posterior patterning of the fin bud is relatively normal in twist1a mutants and twist1a;twist1b double mutants, albeit with an anterior extension of shha expression in some embryos. Conclusions: Together, our results show that twist1a and twist1b play crucial and overlapping roles in zebrafish pectoral fin development, echoing roles played by Twist1 in the mouse limb bud and highlighting an early impact of twist1 genes on the initial formation of the zebrafish forelimb field.

Inherited retinal disorders are significant contributors of blindness worldwide. Mutations in Peripherin-2 (PRPH2), a highly conserved vertebrate tetraspanin membrane protein responsible for formation and maintenance of OS morphology, have been shown to cause diverse types of inherited photoreceptor cell (PRC) disorders including but not limited to Leber congenital amaurosis, cone-rod dystrophy, and retinitis pigmentosa. In this study we used a cone-rich diurnal zebrafish model to characterize the loss of PRPH2 function. Of the four PRPH2 zebrafish orthologs only prph2a and prph2b were found to be expressed in PRCs. CRISPR-mediated single mutants of prph2a and prph2b did not yield striking rod or cone phenotypes. Double prph2a/2b mutants exhibited early loss of all cone cells, preceded by cone outer segment disorganization in the form of whorls akin to the phenotypes observed in PRPH2+/- mice. Surprisingly rod photoreceptor cells were not affected and in fact exhibited a striking lengthening of rod OSs with normal disc formation. Overgrowth of rod OSs proceeded up to 1 year, but no degeneration was observed. To determine how rod OS can persist without prhp2a/b we targeted rom1a and rom1b using CRISPR. Injection of rom1a/b crRNA resulted in complete loss of both rod and cone OSs in the prph2a/b double mutants. Surprisingly, inhibition of rom1a/b alone resulted in the loss of rod but not cone OSs. These findings suggest that unlike in mammals, zebrafish rom1a/b is essential for rod OS formation while prph2a/b is essential for cone OSs.

NANOG is a conserved regulator of epiblast identity and is essential for hypoblast development and embryo survival after the blastocyst stage in mice. However, its role during post-blastocyst development remains poorly understood in non-rodent mammals. Here, we investigated the function of NANOG during ovine embryogenesis using CRISPR/Cas9-mediated gene ablation. NANOG ablation did not affect blastocyst formation or the initial development of the epiblast. However, NANOG-deficient blastocysts contained significantly fewer SOX17+ hypoblast cells by day (D) 8, followed by reduced hypoblast migration during early post-hatching stages in D12 in vitro embryos. These defects were transient, as hypoblast development was comparable between NANOG-deficient and wild-type conceptuses during subsequent in vivo conceptus elongation from embryonic day (E) 11 onwards. In contrast, although extraembryonic membranes development and conceptus elongation proceeded normally, NANOG-deficient conceptuses lacked SOX2+ epiblast cells and failed to form an embryonic disc (ED) from E11 onwards. These findings demonstrate that NANOG is dispensable for hypoblast development but is required for epiblast maintenance and ED formation during conceptus elongation in sheep.

Hox genes play central roles in specifying positional identities along the vertebrate anterior-posterior axis. In mice, genetic analyses have demonstrated that Hox genes distributed among the four Hox clusters contribute to vertebral patterning, with extensive functional redundancy among paralogous genes. Our previous genetic analysis in zebrafish identified important roles for HoxB- and HoxC-related genes in specifying anterior vertebral identities, whereas the contributions of HoxA- and HoxD-related genes remained unresolved. Here, we examined adult zebrafish carrying extensive combinations of hoxaa, hoxab, and hoxda cluster deletions and generated five-gene homozygous mutants carrying frameshift mutations in hoxa3a, hoxa4a, hoxa5a, hoxd3a, and hoxd4a. X-ray micro-CT analysis revealed no obvious alterations in anterior vertebral morphology in either the compound cluster mutants or the five-gene mutants. These results indicate that HoxA/D genes make only a limited detectable contribution to anterior vertebral patterning in zebrafish. Together with our previous findings, they suggest that vertebral patterning functions are distributed unevenly among zebrafish Hox clusters, with a predominant contribution from HoxB/C-related genes.

Retinal development is orchestrated by a network of transcription factors that guide multipotent retinal progenitor cells (RPCs) to fate-committed lineages, ultimately producing seven major retinal cell classes. Among these, retinal ganglion cells (RGCs) serve as the sole output neurons of the retina, relaying visual and non-visual information to the brain. RGC specification requires a cascade of transcriptional regulators, including the basic helix-loop-helix (bHLH) factor Atoh7, which confers competence to RPCs, and Pou4f2 and Isl1, which drive terminal differentiation and subtype diversification. Previous studies demonstrate that the SoxC group transcription factors (Sox4, Sox11, and Sox12) are also involved in RGC genesis, but their precise integration into the Atoh7-driven regulatory hierarchy remains undefined. To address this question, we used the retina-specific Vsx2-Cre line to generate Sox4/Sox11 double conditional knockout (dcKO) and Sox4/Sox11/Atoh7 triple knockout (tKO) mice. Immunohistochemistry revealed profound lineage disruption and reduced progenitor proliferation and survival in both dcKO and tKO retinas; not only RGC genesis but also that of horizontal and amacrine (H&As) cells were severely compromised, whereas photoreceptor cells (PHCs) production increased. These results indicate that Sox4 and Sox11 are involved in the coordinated generation of the different early retinal lineages. Like the Atoh7-null retina, RGC precursors still formed in the SoxC dcKO retina, but their genesis was almost completely abolished in the tKO retina. Our findings indicate that the SoxC factors act in parallel with Atoh7 as a major upstream regulatory input to initiate RGC fate. Bulk RNA-seq revealed the SoxC-dependent transcriptional programs and signaling pathways and confirmed the lineage changes demonstrated by marker analysis. CUT&Tag analysis identified the genome-wide binding sites and thereby the target genes of Sox11, further illuminating the mechanisms underlying functions of the SoxC factors in multiple retinal cell states/types during development.

The MAF family transcription factor Nrl is a key regulator of vertebrate rod photoreceptor formation/differentiation. The expression of Nrl in nascent photoreceptors has been proposed to act as a master regulator to initiate a rod program and to repress a default cone photoreceptor sister fate. The evidence for rod-specific expression of Nrl has been most strongly supported by mouse transgenic models that use a transcriptional promoter element for Nrl. However, this element has not been rigorously evaluated to have rod lineage-specific activity or as a validated proxy for endogenous Nrl expression, despite its wide use under this assumption. Here we identify that the Nrl promoter is not exclusively active in rod photoreceptors but is also transiently active within a high proportion of the differentiating cone photoreceptor population. Thus, the Nrl promoter cannot be used to distinguish rods from cones at developmental timepoints. Furthermore, multiple lines of evidence using immunofluorescence and reanalysis of published single cell RNA-seq/ATAC-seq datasets support the presence of endogenous Nrl gene expression in early differentiating cone photoreceptors. Taken together, these results suggest that proposed gene regulatory networks that position Nrl as a master regulator of rod versus cone fate are inadequate as currently constructed. This impacts our understanding of visual system evolution and has implications for the development of therapeutic strategies of cell replacement for human blindness.

Background: Viable pediatric human growth plate (GP) tissue is rarely available for translational research, limiting direct investigation of human longitudinal bone growth and pediatric growth disorders. In this proof-of-concept study, we aimed to determine whether it is feasible to establish a clinically integrated ex vivo human GP model using tissue obtained during routine percutaneous epiphysiodesis (PE) procedures in adolescents treated for extreme tall stature or leg length difference due to trauma. Methods: GP tissue and cells were collected during PE and processed using protocols adapted from established methods of human osteoarthritic cartilage processing within the RAAK study. Feasibility was assessed by evaluating tissue collection, cell isolation, contamination rate, monolayer expansion, and generation of three-dimensional cartilage pellets. Proliferation of GP-derived chondrocytes was compared with osteoarthritis-derived articular chondrocytes, and histological assessment was performed to evaluate cartilage-like matrix formation. Results: Across consecutive surgical procedures, viable GP tissue could be obtained reproducibly, with only few samples failing to yield cells and no relevant contamination issues. Isolated GP chondrocytes expanded successfully in two-dimensional culture and showed a strong early proliferative response compared with RAAK-derived chondrocytes. In addition, GP-derived cells formed three-dimensional organoids and histology confirmed cartilage-like matrix deposition supporting their capacity to generate neo-cartilage tissue in vitro. Conclusion: This study demonstrates feasibility to obtain, culture, and functionally assess viable human GP tissue from routine PE surgery. As such, the Leiden ex vivo human GP model provides a unique platform to study local mechanisms of endochondral bone growth, link genetic determinants of height to functional GP biology, and support future therapeutic research in pediatric growth disorders.

Quantifying skeletal mineralization phenotypes in larval fish is complicated by the natural curvature of the notochord and by sample-to-sample variability in orientation, staining and imaging. Consequently, many studies rely on summary measures such as vertebral counts or total stain intensity. Here we present SCAMP (Spinal Calcification & Mineralization Profiler), an open-source, GUI-based Python tool that computationally straightens the curved notochord of Alizarin Red S-stained fish larvae and generates standardized mineralization profiles along the spinal axis. This approach reduces positional and shape variability, allowing direct, quantitative comparison of calcification patterns within and between experimental cohorts, without requiring programming expertise. We validate SCAMP using a zebrafish model of Pseudoxanthoma elasticum (abcc6aelu15/elu15), recovering genotype-specific differences in the intensity, extent and spatial distribution of ectopic calcification. Using SCAMP, we further show that inorganic pyrophosphate (PPi) supplementation of the medium suppresses ectopic notochord calcification, alters the anterior-posterior distribution of mineralized regions in homozygous mutants, and promotes mineralization at physiological vertebral sites. We also show that methylene blue, a routine antifungal additive in fish medium, reduces baseline calcification, with the most pronounced effects observed in heterozygous controls. SCAMP is freely available and has the potential to be adapted to other fish species used in skeletal and mineralization research.

MyoD plays a central role in determining the skeletal muscle lineage in vertebrate embryos. The core enhancer (CE) and distal regulatory region (DRR) are the only known MyoD enhancers, and together, they recapitulate all major aspects of MyoD expression in the embryo. However, knocking out each enhancer individually has only modest effects on MyoD expression. Here, we show that embryos lacking both enhancers maintain muscle-specific MyoD expression, indicating the existence of unknown MyoD regulatory elements. Precision run-on sequencing together with available ChIP-seq and DNase I hypersensitivity datasets identified three new candidate enhancer regions within 96 kb of MyoD 5 flanking sequences. Transgenic analysis revealed that DNA elements at -36 and -60 kb are active in all muscle-forming regions, each recapitulating aspects of endogenous MyoD expression. Enhancer activities in muscle regulatory factor-deficient mice suggest that they are components of the auto- and cross-regulatory circuitry that maintains MyoD expression. Analysis of Hi-C data showed that the entire -96 kb region constitutes a loop domain, within which multiple interactions between elements and with the MyoD gene were detected. A larger loop domain delimited by CTCF sites was also identified from -96 kb to +220 kb relative to the MyoD transcriptional start site. These data indicate that the newly identified enhancers are key components of a cis regulatory network that controls the activation and maintenance of MyoD expression in the embryo. One sentence summaryCis regulation of MyoD transcription during development

Background: Primary embryonic dorsal root ganglion (DRG) explants offer an accessible and physiologically relevant organotypic in vitro assay platform that preserves native three-dimensional cellular interactions and supports robust neurite outgrowth, enabling studies of sensory neuron development and neurite outgrowth, as well as controlled assessment of neurotoxicity, and pharmacological compound effects. Methods: Here we provide a comprehensive step-by-step method to micro-dissect and culture embryonic chicken DRGs and quantitatively measure neurite outgrowth. The protocol covers embryo preparation, vertebral column exposure, bilateral DRG microdissection, explant culturing, longitudinal phase-contrast imaging, and neurite image analysis. Results: DRG explants maintained their overall ganglion morphology and remained attached to the culture surface throughout the five-day culture period. Untreated control explants exhibited limited neurite outgrowth, whereas nerve growth factor (NGF) treated explants developed progressively longer and denser radial neurite networks. Neurite outgrowth was readily visualized by live phase contrast microscopy and quantified using an Image-based neurite tracing method, demonstrating the responsiveness of the explant model to NGF and its suitability for quantitative assessment of neurite outgrowth. Conclusion: This simple, low-cost protocol provides an in vitro assay platform for neurotoxicity, neuroregeneration, drug screening, and biomaterial evaluation.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a complex multisystem disorder. Its hallmark feature, post-exertional malaise (PEM), is a delayed and prolonged worsening of symptoms following even minor exertion. As PEM commonly includes muscle pain, weakness, exercise intolerance and impaired recovery, skeletal muscle dysfunction may contribute to its pathophysiology, but the underlying molecular processes remain poorly understood. We generated induced pluripotent stem cell-derived skeletal muscle myotubes from individuals with ME/CFS (n = 6) and healthy controls (n = 7) to investigate disease-associated transcriptomic alterations and identify potential therapeutic candidates. Myogenic progenitor cells from both groups differentiated into elongated, multinucleated myotubes expressing titin and myosin heavy chain, while qPCR profiling of myogenic differentiation markers further supported successful differentiation. Gene-level analysis identified three differentially expressed genes at an FDR-adjusted p < 0.05: ERVH48-1, CDX4 and SLC2A14. Isoform-switching analysis identified 16 genes with significant differential isoform usage, including CACNA1H, MIEF2, SMAD2 and CHD4. Gene set enrichment analysis identified broader pathway-level alterations involving RNA processing, protein synthesis, proteostasis, mitochondrial quality control, metabolism and muscle-related processes. Library of Integrated Network-based Cellular Signatures Phase 2 (LINCS2) analysis further prioritised 14 compounds associated with mitochondrial metabolism, glucose and lipid regulation, PPAR signalling, calcium and excitatory signalling, neuroimmune activity and cellular stress responses. Overall, ME/CFS iPSC-derived myotubes showed modest gene-level differences alongside broader isoform- and pathway-level alterations, supporting their use for investigating skeletal muscle biology and prioritising candidates for functional validation.

Myotonic dystrophy type 1 (DM1) is caused by a CTG expansion in the 3' untranslated region of the DMPK gene that results in the expression of expanded CUG RNA which trap splicing factors resulting in widespread alternative splicing dysregulation. A panel of these splicing events, or splicing index (SI), is currently used to measure splicing dysregulation in DM1 patients. Despite significant interest in splicing, little is known about its progression over time in DM1. The objective of this study was to evaluate the changes in SI over 3 years in a well-phenotyped cohort of DM1 participants and identify factors contributing to these changes. RNA was extracted from vastus lateralis (VL) biopsies of 20 DM1 participants taken 3 years apart and the SI score was measured using targeted RNA sequencing. Baseline SI scores correlated with percent predicted muscle strength of the ankle dorsiflexors, knee extensors and hand grip. On average, SI scores worsened over 3 years, but the progression varied markedly among participants. Adult/juvenile phenotypes and lower baseline SI scores were among the factors associated with higher progression of SI scores over time. These findings emphasize the importance of integrating clinical phenotype and baseline molecular severity when interpreting longitudinal SI change, particularly in the context of clinical trials.

Multiple hypotheses have been formulated to explain differences in tissue repair ability across vertebrates. One hypothesis posits that the accessibility of actively cycling stromal cells within uninjured tissue confers access to a proliferative population in response to tissue damage. This hypothesis further suggests that animals with an indeterminate growth mode possess an actively cycling cell population necessary for growth that can be readily accessed for tissue regeneration. Moreover, the absence of an actively cycling population in connective tissue provides a mechanism that restricts regeneration in animals with determinate growth whose cells are refractory to cell cycle progression and proliferation to produce new tissue for morphogenesis. Here, we explore this paradigm using an EdU-BrdU pulse chase strategy in four different vertebrate species: two with determinate (Acomys dimidiatus and Mus musculus) and two with indeterminate modes of growth (Danio rerio and Ambystoma mexicanum). We find that although indeterminate growers do possess a small population of actively cycling cells, this population does not contribute to regeneration. Moreover, we found that while Acomys does not possess a population of actively cycling stromal cells, cells re-enter the cell cycle de novo in these animals to contribute to regeneration. Furthermore, testing this hypothesis allowed us to ask whether tissue injury could stimulate cell cycle re-entry - a so-called primed state - in cells at distance from the injury site in these four species and we did not find evidence of such priming in stromal or epidermal tissue. HighlightsO_LICell cycle re-entry is a common response to injury in regenerative and non-regenerative vertebrates that is independent of actively cycling stromal cells in uninjured connective tissue C_LIO_LIActively cycling cells do not contribute to regenerative healing in spiny mice, axolotls or zebrafish. C_LIO_LIOur data do not support systemic cell cycle activation in response to injury. C_LI

In CNS endothelial cells (ECs), VEGF signaling promotes vascular permeability, and Wnt/beta-catenin signaling reduces vascular permeability by controlling the gene expression program for the blood-brain barrier. Here we show, using genetic mosaics, that, in mouse brain ECs in vivo, an increase in Wnt/beta-catenin signaling produces an increase in VEGFR1 levels and a decrease in VEGFR2 levels, and a decrease in Wnt/beta-catenin signaling produces a decrease in VEGFR1 levels and an increase in VEGFR2 levels. As VEGFR1 functions as a decoy receptor to reduce VEGF signaling through VEGFR2, these data imply that Wnt/beta-catenin signaling acts at the receptor level to reduce VEGF signaling. In HEK/293T cells, VEGF signaling is suppressed by Wnt/beta-catenin signaling in a dose-dependent manner, with nearly complete suppression observed at levels of Wnt/beta-catenin signaling that produce little or no change in VEGF receptor levels. These data reveal two mechanisms by which VEGF signaling is regulated by Wnt/beta-catenin signaling.

Scoliosis is a three-dimensional deformity of the vertebral column that commonly emerges during growth, yet the biological mechanisms that initiate and drive curve progression remain poorly understood. Serotonin (5-hydroxytryptamine [5-HT]) has long been linked to scoliosis through its relationship with pineal melatonin signaling and through genetic associations with human disease, but its role in spinal morphogenesis has remained unclear. Here, using zebrafish, we show that serotonergic signaling is required for normal spine morphology. Transient exposure to 5-HT during development induces later spinal abnormalities, while loss of Tph2, the principal enzyme required for neuronal 5-HT synthesis, causes progressive three-dimensional spinal curvature during juvenile growth, mimicking aspects of adolescent idiopathic scoliosis (AIS). tph2 mutants undergo initially normal vertebral patterning while targeted ablation of the pineal gland does not induce scoliosis, arguing against a pineal or early skeletal origin for the 5-HT-associated spinal curves. The Reissner fiber, which has been linked to spinal curvature in zebrafish, also formed normally in tph2 mutants. Instead, mutants showed a marked reduction in Fev-expressing intraspinal serotonergic neurons together with increased locomotor activity. Mutants exhibited longer swim bouts and shorter periods of inactivity between bouts. Importantly, these motor defects were evident shortly before the onset of overt spinal curvature. Together, these findings identify neuronal serotonin as a critical regulator of spinal stability and support a model in which altered serotonergic control of motor function changes the mechanical environment experienced by the growing spine, increasing its susceptibility to progressive curvature.

Objective: Kabuki syndrome type 2 (KS2) is a rare X-linked disorder caused by pathogenic variants in KDM6A and is frequently associated with hearing impairment. Because auditory phenotypes include sensorineural, conductive, and mixed hearing loss, we examined whether Kdm6a deficiency directly impairs cochlear morphology or auditory function in mice. Methods: We generated Nestin-Cre; Kdm6a conditional mutants and a germline-derived global Kdm6a deletion model using Sox2-Cre-mediated germline recombination. Auditory function and cochlear morphology were assessed by auditory brainstem responses (ABRs), distortion product otoacoustic emissions (DPOAEs), and histological analyses. Results: Whole-body hemizygous or homozygous Kdm6a deletion resulted in embryonic lethality. Analyses therefore focused on Nestin-Cre; Kdm6aflox/Y males and females carrying a germline-derived heterozygous Kdm6a deletion. In Nestin-Cre mutants, ABR thresholds, wave I amplitudes and latencies, and DPOAE measures at 8, 16, and 32 kHz were comparable with littermate controls at 4 and 8 weeks. In heterozygous deletion females, ABR thresholds, wave I amplitudes and latencies at 4-64 kHz and DPOAE measures at 8, 16, and 32 kHz were comparable with male and female wild-type controls. No overt cochlear pathology was observed in either model. Conclusion: Kdm6a deficiency in the viable mouse models examined did not produce overt cochlear pathology or baseline auditory dysfunction. These findings suggest that partial Kdm6a loss alone is insufficient to produce a robust primary cochlear phenotype in young mice and highlight the potential contributions of non-cochlear, developmental, central, and later-onset mechanisms to hearing impairment in KS2. Clinically, these findings support careful evaluation for potentially treatable middle-ear disease.

The neural crest and neurogenic placodes, which arise from the neural plate border, give rise to morphological characteristics that distinguish vertebrates from invertebrates (Gans & Northcutt, 1983). Recent studies have suggested that embryos of ascidians, a group of tunicates that are the closest invertebrate relatives of vertebrates, possess cells that share an evolutionary origin with vertebrate neural crest cells (Abitua et al., 2012; Fatieieva et al., 2025; Ishida & Satou, 2024; Stolfi et al., 2015; Todorov et al., 2024; Waki et al., 2015) and neurogenic placode cells (Abitua et al., 2015; Ikeda et al., 2013; Liu et al., 2023; Liu & Satou, 2019; Manni et al., 2004; Mazet et al., 2005; Papadogiannis et al., 2022; Wagner & Levine, 2012). To dissect the neural plate border of ascidian embryos at the molecular level, and to gain deeper insights into evolutionary origins of the neural crest and neurogenic placodes, we comprehensively analyzed expression patterns of transcription factor genes at single-cell resolution. We demonstrated that the ascidian neural plate border consists of three domains with distinct gene expression profiles: the anterior, inner lateral, and outer lateral domains. A cross-species comparison of transcriptomes from ascidians and zebrafish suggests that the anterior domain is homologous to zebrafish neurogenic placodes, the inner lateral domain to the neural crest and tail bud and the outer lateral domain to the median fin fold ectoderm. We propose that a tripartite neural plate border was present in the last common ancestor of vertebrates and tunicates, providing a blueprint for evolutionary emergence of vertebrate morphological novelties.

Pathogenic variants in POLR1C, which encodes a shared subunit of RNA Polymerases (Pols) I and III, cause Treacher Collins syndrome (TCS) and POLR3-related leukodystrophy. While Pol I and Tp53-dependent mechanisms have been implicated in the pathogenesis of TCS, the basis of hypomyelination in POLR1C-associated POLR3-related leukodystrophy remains incompletely understood. Here, we show that polr1c mutant zebrafish exhibit reduced myelination in addition to previously described craniofacial anomalies. Oligodendrocyte precursor cells exhibit increased activation of the Tp53 pathway; however, these cells do not undergo apoptosis. Consistent with this finding, tp53 inhibition reduces cell death in polr1c mutants but fails to restore myelination, indicating that myelination deficits are not driven by Tp53-dependent progenitor loss in this model. polr1c mutants also exhibit reduced rRNA transcription by Pol I and reduced expression of some Pol III-transcribed tRNAs. Altogether, these data indicate distinct tissue-specific responses to polr1c deficiency and suggest persistent impairment of rRNA transcription contributes to deficient myelin development. These findings expand the developmental consequences of polr1c loss and advance our understanding of the molecular basis of POLR1C-associated diseases.

Osteoarthritis (OA) disproportionately affects women, and estrogen has been implicated in cartilage and joint homeostasis, yet its effects on meniscal fibrochondrocytes (MFCs), and whether those effects differ by sex, remain poorly defined. We reanalyzed a publicly available RNA-sequencing dataset (Gene Expression Omnibus, GSE199087) comprising human MFCs from a male and a female donor treated with 17{beta}-estradiol (E2) or vehicle. Differential expression, Gene Set Enrichment Analysis, Enrichr, and iLINCS were integrated to identify the transcriptional programs modulated by E2 in each sex. In female MFCs, E2 upregulated pathways governing DNA replication, cell-cycle progression, and genomic maintenance (e.g., MCM8, BRCA1, RAD51), while downregulating inflammatory and extracellular matrix-degrading genes, including MMP1, IL12A, CXCL12, and MYD88. In male MFCs, E2 instead upregulated chromatin remodeling and developmental signaling programs, led by SRCAP, ERCC6, and NOTCH1, and downregulated antigen presentation and mitochondrial genes. These divergent responses indicate that E2 engages distinct, sex-specific transcriptional programs in MFCs, with the female profile favoring proliferation and matrix preservation. Although derived from a limited sample and therefore hypothesis-generating, these findings nominate candidate mechanisms underlying sex differences in meniscal biology and OA susceptibility, and warrant validation in larger, sex-balanced cohorts. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/740612v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@16026fforg.highwire.dtl.DTLVardef@1a048acorg.highwire.dtl.DTLVardef@3246corg.highwire.dtl.DTLVardef@4c598b_HPS_FORMAT_FIGEXP M_FIG C_FIG

Mammalian skeletal muscle research predominantly focuses on locomotor muscles, and feeding related muscles remain less extensively characterized despite their role in mastication, mandibular stabilization, and swallowing. In this study, we investigated the transcriptomic specialization of three functionally and developmentally unique rat muscles: the anterior digastric (AD), a jaw opening muscle; the superficial masseter (SM), a jaw closing muscle; and the Sternohyoid (SH), a non-mandibular muscle involved in swallowing. Differential gene expression and weighted gene co-expression network analysis were used to characterize the transcription level features associated with their distinct roles. Our results indicated that all three muscles predominantly expressed fast-twitch contractile isoforms. However, the AD showed lower overall expression of several contractile gene families, including myosin heavy chain, myosin light chain, and tropomyosin isoforms, while exhibiting elevated expression of slow/oxidative myosin isoforms like Myh7 and Myh2. Network analysis revealed that modules correlated with AD are strongly enriched for fatty acid catabolism, mitochondrial energy production, and vascular/extracellular matrix remodeling. Additionally, AD and SM shared a distinct gene set compared to SH, highlighting their common developmental origin from the first branchial arch. Our findings show that the rat feeding related muscles possess unique transcriptomic profiles shaped by their contractile functions, developmental origins, and metabolic functions.