An Alternative Splicing Program for Mouse Craniofacial Development

Splicing Defects and Cell Death Cause SF3B2-Linked Craniofacial Microsomia

Craniofacial microsomia (CFM) is a genetically and phenotypically heterogeneous disorder characterized by hypoplasia of facial tissue that is often asymmetric. Affected tissues typically include the ears (external and internal), mandible, and maxilla, but various extracranial anomalies have also been reported. Loss-of-function variants in the SF3B2 gene have recently been reported in 8 cases of CFM, representing one of the more common genetic causes identified to date. To better define the full phenotypic spectrum associated with variants in SF3B2, we report novel loss-of-function variants in SF3B2 in 5 new families with CFM. Furthermore, to determine the mechanism by which SF3B2 loss-of-function perturbs craniofacial development, we established sf3b2-null mutant zebrafish, which exhibited severe deficiencies in craniofacial cartilage and bone progenitors due to elevated apoptosis and reduced proliferation of cranial neural crest cells. In addition, we generated a heterozygous truncating variant of SF3B2 in human induced pluripotent stem cells using CRISPR/Cas9 gene editing. Differentiation of these cells into neural crest cells was accompanied by increased cell death and reduced proliferation. RNA sequencing of sf3b2 mutant zebrafish revealed widespread disruption of mRNA splicing, including mdm2, a key regulator of Tp53-mediated apoptosis. Genetic inhibition of tp53 in sf3b2 mutants demonstrated that tp53 inhibition reduces early cell death but does not improve proliferation or craniofacial cartilage development. Therefore, our functional studies indicate that widespread mRNA splicing disruption, in addition to Tp53-dependent cell death, contributes to the craniofacial features observed in SF3B2-related CFM.

Gene expression patterns of the developing human face at single cell resolution reveal cell type contributions to normal facial variation and disease risk

Craniofacial development gives rise to the complex structures of the face and involves the interplay of diverse cell types. Despite its importance, our understanding of human-specific craniofacial developmental mechanisms and their genetic underpinnings remains limited. Here, we present a comprehensive single-nucleus RNA sequencing (snRNA-seq) atlas of human craniofacial development from craniofacial tissues of 24 embryos that span six key time points during the embryonic period (4-8 post-conception weeks). This resource resolves the transcriptional dynamics of seven major cell types and uncovers distinct major cell types, including muscle progenitors and cranial neural crest cells (CNCCs), as well as dozens of subtypes of ectoderm and mesenchyme. Comparative analyses reveal substantial conservation of major cell types, alongside human biased differences in gene expression programs. CNCCs, which play a crucial role in craniofacial morphogenesis, exhibit the lowest marker gene conservation, underscoring their evolutionary plasticity. Spatial transcriptomics further localizes cell populations, providing a detailed view of their developmental roles and anatomical context. We also link these developmental processes to genetic variation, identifying cell type-specific enrichments for common variants associated with facial morphology and rare variants linked to orofacial clefts. Intriguingly, Neanderthal-introgressed sequences are enriched near genes with biased expression in cartilage and specialized ectodermal subtypes, suggesting their contribution to modern human craniofacial features. This atlas offers unprecedented insights into the cellular and genetic mechanisms shaping the human face, highlighting conserved and distinctly human aspects of craniofacial biology. Our findings illuminate the developmental origins of craniofacial disorders, the genetic basis of facial variation, and the evolutionary legacy of ancient hominins. This work provides a foundational resource for exploring craniofacial biology, with implications for developmental genetics, evolutionary biology, and clinical research into congenital anomalies.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

The Tumor Suppressor Adenomatous Polyposis Coli (apc) Is Required for Neural Crest-Dependent Craniofacial Development in Zebrafish

Neural crest (NC) is a unique vertebrate cell type arising from the border of the neural plate and epidermis that gives rise to diverse tissues along the entire body axis. Roberto Mayor and colleagues have made major contributions to our understanding of NC induction, delamination, and migration. We report that a truncating mutation of the classical tumor suppressor Adenomatous Polyposis Coli (apc) disrupts craniofacial development in zebrafish larvae, with a marked reduction in the cranial neural crest (CNC) cells that contribute to mandibular and hyoid pharyngeal arches. While the mechanism is not yet clear, the altered expression of signaling molecules that guide CNC migration could underlie this phenotype. For example, apcmcr/mcr larvae express substantially higher levels of complement c3, which Mayor and colleagues showed impairs CNC cell migration when overexpressed. However, we also observe reduction in stroma-derived factor 1 (sdf1/cxcl12), which is required for CNC migration into the head. Consistent with our previous work showing that APC directly enhances the activity of glycogen synthase kinase 3 (GSK-3) and, independently, that GSK-3 phosphorylates multiple core mRNA splicing factors, we identify 340 mRNA splicing variations in apc mutant zebrafish, including a splice variant that deletes a conserved domain in semaphorin 3f (sema3f), an axonal guidance molecule and a known regulator of CNC migration. Here, we discuss potential roles for apc in CNC development in the context of some of the seminal findings of Mayor and colleagues.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Systematic analysis of alternative splicing in time course data using Spycone

MOTIVATION

During disease progression or organism development, alternative splicing may lead to isoform switches that demonstrate similar temporal patterns and reflect the alternative splicing co-regulation of such genes. Tools for dynamic process analysis usually neglect alternative splicing.

RESULTS

Here, we propose Spycone, a splicing-aware framework for time course data analysis. Spycone exploits a novel IS detection algorithm and offers downstream analysis such as network and gene set enrichment. We demonstrate the performance of Spycone using simulated and real-world data of SARS-CoV-2 infection.

AVAILABILITY AND IMPLEMENTATION

The Spycone package is available as a PyPI package. The source code of Spycone is available under the GPLv3 license at https://github.com/yollct/spycone and the documentation at https://spycone.readthedocs.io/en/latest/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptome Analysis Explored the Differential Genes’ Expression During the Development of the Stropharia rugosoannulata Fruiting Body

Stropharia rugosoannulata (S. rugosoannulata) is a fungus with great edible and nutritional values; however, the development mechanism of its fruiting body has not been studied. Thus, this study aimed to analyze the differentially expressed genes (DEGs) in four stages; primordia stage (Sra1), young mushroom stage (Sra2), picking stage (Sra3), and opening umbrella stage (Sra4). Therefore, total RNA was extracted for further RNA-sequencing analysis. In three pairwise comparison groups (PCGs), Sra1 vs. Sra2, Sra2 vs. Sra3, and Sra3 vs. Sra4, a total of 3,112 DEGs were identified among the three PCGs. A GO analysis of the DEGs showed that there were 21 terms significantly enriched in Sra1 vs. Sra2 PCG. There was no significantly enriched GO term in the other two PCGs. Furthermore, KEGG pathway analysis showed that these DEGs were mainly enriched in glucose and amino acid metabolisms. Moreover we found that intron retention (IR) and the alternative 3′ splice site (A3SS) accounted for more than 80%. The development of the S. rugosoannulata fruiting body mainly involved glucose and amino acid metabolisms. IR and A3SS were the two main types of ASE, which played an important role in the development and maturation of the S. rugosoannulata fruiting body.

AP-2α and AP-2β cooperatively function in the craniofacial surface ectoderm to regulate chromatin and gene expression dynamics during facial development

The facial surface ectoderm is essential for normal development of the underlying cranial neural crest cell populations, providing signals that direct appropriate growth, patterning, and morphogenesis. Despite the importance of the ectoderm as a signaling center, the molecular cues and genetic programs implemented within this tissue are understudied. Here, we show that removal of two members of the AP-2 transcription factor family, AP-2α and AP-2ß, within the early embryonic ectoderm of the mouse leads to major alterations in the craniofacial complex. Significantly, there are clefts in both the upper face and mandible, accompanied by fusion of the upper and lower jaws in the hinge region. Comparison of ATAC-seq and RNA-seq analyses between controls and mutants revealed significant changes in chromatin accessibility and gene expression centered on multiple AP-2 binding motifs associated with enhancer elements within these ectodermal lineages. In particular, loss of these AP-2 proteins affects both skin differentiation as well as multiple signaling pathways, most notably the WNT pathway. We also determined that the mutant clefting phenotypes that correlated with reduced WNT signaling could be rescued by Wnt1 ligand overexpression in the ectoderm. Collectively, these findings highlight a conserved ancestral function for AP-2 transcription factors in ectodermal development and signaling, and provide a framework from which to understand the gene regulatory network operating within this tissue that directs vertebrate craniofacial development.

The Role of RNA-Binding Proteins in Vertebrate Neural Crest and Craniofacial Development

Cranial neural crest (NC) cells delaminate from the neural folds in the forebrain to the hindbrain during mammalian embryogenesis and migrate into the frontonasal prominence and pharyngeal arches. These cells generate the bone and cartilage of the frontonasal skeleton, among other diverse derivatives. RNA-binding proteins (RBPs) have emerged as critical regulators of NC and craniofacial development in mammals. Conventional RBPs bind to specific sequence and/or structural motifs in a target RNA via one or more RNA-binding domains to regulate multiple aspects of RNA metabolism and ultimately affect gene expression. In this review, we discuss the roles of RBPs other than core spliceosome components during human and mouse NC and craniofacial development. Where applicable, we review data on these same RBPs from additional vertebrate species, including chicken, Xenopus and zebrafish models. Knockdown or ablation of several RBPs discussed here results in altered expression of transcripts encoding components of developmental signaling pathways, as well as reduced cell proliferation and/or increased cell death, indicating that these are common mechanisms contributing to the observed phenotypes. The study of these proteins offers a relatively untapped opportunity to provide significant insight into the mechanisms underlying gene expression regulation during craniofacial morphogenesis.

Srsf3 mediates alternative RNA splicing downstream of PDGFRα signaling in the facial mesenchyme

Signaling through the platelet-derived growth factor receptor alpha (PDGFRα) is crucial for mammalian craniofacial development, although the mechanisms by which the activity of downstream intracellular effectors is regulated to mediate gene expression changes have not been defined. We find that the RNA-binding protein Srsf3 is phosphorylated at Akt consensus sites downstream of PI3K-mediated PDGFRα signaling in mouse palatal mesenchyme cells, leading to its nuclear translocation. We further demonstrate that ablation of Srsf3 in the mouse neural crest lineage leads to facial clefting due to defective cranial neural crest cell proliferation and survival. Finally, we show that Srsf3 regulates the alternative RNA splicing of transcripts encoding protein kinases in the mouse facial process mesenchyme to regulate PDGFRα-dependent intracellular signaling. Collectively, our findings reveal that alternative RNA splicing is an important mechanism of gene expression regulation downstream of PI3K/Akt-mediated PDGFRα signaling in the facial mesenchyme and identify Srsf3 as a critical regulator of craniofacial development.

FaceBase 3: analytical tools and FAIR resources for craniofacial and dental research

The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a 'big data' resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements.