The splicing regulators Esrp1 and Esrp2 direct an epithelial splicing program essential for mammalian development

Gene Duplication and Alternative Splicing as Evolutionary Drivers of Proteome Specialization

Animals comprise hundreds of cell types, each with specialized biological functions. However, many genes expressed in each cell type belong to widely conserved gene families with ancestrally ubiquitous expression. This raises a paradox: how have these genes evolved to shape cell type-specific traits without compromising their ancestral function in all other cells? This can be achieved through gene duplication and the origin of regulated, alternatively spliced exons, which generate new related proteins in the form of paralogous genes and alternative isoforms, respectively. Here, we explore how such new related proteins can contribute to the evolution of specific cell types while preserving broader ancestral roles. Specifically, we separately classify possible expression and functional fates for new related proteins and discuss their interplays and evolutionary likelihood. Our primary hypothesis is that expression specialization, mostly coupled with functional specialization, is the predominant fate for both paralogous genes and alternative isoforms throughout animal evolution.

Spatial gene regulatory networks driving cell state transitions during human liver disease

Liver fibrosis is a major cause of death worldwide. As a progressive step in chronic liver disease, fibrosis is almost always diagnosed too late with limited treatment options. Here, we uncover the spatial transcriptional landscape driving human liver fibrosis using single nuclei RNA and Assay for Transposase-Accessible Chromatin (ATAC) sequencing to deconvolute multi-cell spatial transcriptomic profiling in human liver cirrhosis. Through multi-modal data integration, we define molecular signatures driving cell state transitions in liver disease and define an impaired cellular response and directional trajectory between hepatocytes and cholangiocytes associated with disease remodelling. We identify pro-fibrogenic signatures in non-parenchymal cell subpopulations co-localised within the fibrotic niche and localise transitional cell states at the scar interface. This combined approach provides a spatial atlas of gene regulation and defines molecular signatures associated with liver disease for targeted therapeutics or as early diagnostic markers of progressive liver disease.

Genetic requirement for Esrp1/2 in vertebrate pituitary morphogenesis

The pituitary gland produces several hormones that regulate growth, metabolism, stress response, reproduction, and homeostasis. Congenital hypopituitarism is a deficiency in one or more pituitary hormones and encompasses a spectrum of clinical conditions. The pituitary has a complex embryonic origin with the oral ectoderm contributing the anterior lobe, and the neural ectoderm generating the posterior lobe. Pituitary abnormalities and growth deficiencies are associated with cleft palate however the developmental genetic connection between pituitary and orofacial cleft malformations remains to be determined. The epithelial RNA splicing regulators Esrp1 and Esrp2 are required for orofacial development in zebrafish, mice, and humans, and loss of function of these genes results in a cleft palate. Here we present a detailed developmental analysis of the genetic requirement for Esrp1/2 in pituitary morphogenesis in mouse and zebrafish. Further, we describe a patient with cleft palate and hypopituitarism that harbors a nucleotide variant in the RNA binding domain of ESRP2. The discovery of this key function for Esrp1/2 in pituitary formation has significant fundamental and clinical implications for understanding congenital hypopituitarism and craniofacial anomalies.

Esrp1 and Esrp2 regulate the stability of tmc1/2a mRNAs in zebrafish sensory hair cells

RNA-binding proteins (RBPs) are important for post-transcriptional RNA processing, including pre-mRNA alternative splicing, mRNA stability, and translation. Several RBPs have been shown to play pivotal roles in the inner ear, whose dysfunction leads to auditory and/or balance impairments. Epithelial splicing-regulatory protein 1 (ESRP1) regulates alternative splicing and mRNA stability, and mutations in ESRP1 gene have been associated with sensorineural hearing loss in humans. In Esrp1 knockout mouse embryos, alternative splicing of its target genes such as Fgfr2 is impaired, which eventually result in cochlear development deficits. However, Esrp1 knockout mice die soon after birth because of complications from cleft-lip and palate defects, impeding further investigations at later postnatal ages. In the present study, we explored the role of ESRP1 in hearing using zebrafish as a model. We showed that esrp1 and its paralog esrp2 are expressed in the inner ear and certain anterior lateral line (ALL) neuromasts. Furthermore, our data suggested that Esrp1 and Esrp2 are required for the mechano-electrical transduction (MET) function of hair cells. RNA sequencing results indicated a significant decrease in the levels of several mRNAs in esrp1/2 double knockout larvae. Among the dysregulated genes are tmc1 and tmc2a, which encode essential subunits of the MET complex. Further investigations demonstrated that Esrp1/2 could directly bind to tmc1 and tmc2a mRNAs and affect their stability. Taken together, we showed here that Esrp1 and Esrp2 regulate the MET function of zebrafish sensory hair cells by modulating the stability of tmc1 and tmc2a mRNAs.Significance statement ESRP1 is an important RNA-binding protein, whose malfunction has been associated with hearing loss in humans. Esrp1 knockout affects alternative splicing of its target mRNAs such as Fgfr2, eventually leading to cochlear development deficits in mice. However, Esrp1 knockout mice die soon after birth, precluding further investigations at later postnatal ages. In this study, we explored the role of ESRP1 in hearing using zebrafish as a model. Our results demonstrated that esrp1 and its paralog esrp2 are expressed in the zebrafish inner ear, and that esrp1/esrp2 double knockout compromised the mechano-electrical transduction (MET) function of hair cells. Additionally, we successfully identified tmc1 and tmc2a mRNAs as the targets of Esrp1/2, which encode essential subunits of the MET complex.

ESRP2-microRNA-122 axis promotes the postnatal onset of liver polyploidization and maturation

Hepatocyte polyploidy and maturity are critical to acquiring specialized liver functions. Multiple intracellular and extracellular factors influence ploidy, but how they cooperate temporally to steer liver polyploidization and maturation or how post-transcriptional mechanisms integrate into these paradigms is unknown. Here, we identified an important regulatory hierarchy in which postnatal activation of epithelial splicing regulatory protein 2 (ESRP2) stimulates processing of liver-specific microRNA (miR-122) to facilitate polyploidization, maturation, and functional competence of hepatocytes. By determining transcriptome-wide protein-RNA interactions in vivo and integrating them with single-cell and bulk hepatocyte RNA-seq data sets, we delineated an ESRP2-driven RNA processing program that drives sequential replacement of fetal-to-adult transcript isoforms. Specifically, ESRP2 binds the primary miR-122 host gene transcript to promote its processing/biogenesis. Combining constitutive and inducible ESRP2 gain- and loss-of-function mouse models with miR-122 rescue experiments, we demonstrated that timed activation of ESRP2 augments the miR-122-driven program of cytokinesis failure, ensuring the proper onset and extent of hepatocyte polyploidization.

Genome-wide landscape of genetic diversity, runs of homozygosity, and runs of heterozygosity in five Alpine and Mediterranean goat breeds

BACKGROUND

Goat breeds in the Alpine area and Mediterranean basin exhibit a unique genetic heritage shaped by centuries of selection and adaptability to harsh environments. Understanding their adaptive traits can aid breeding programs target enhanced resilience and productivity, especially as we are facing important climate and agriculture challenges. To this aim the genomic architecture of 480 goats belonging to five breeds (i.e., Saanen [SAA], Camosciata delle Alpi [CAM], Murciano-Granadina [MUR], Maltese [MAL], Sarda [SAR]) reared in the Sardinia Island were genotyped and their genomic architecture evaluated to find molecular basis of adaptive traits. Inbreeding, runs of homozygosity (ROH) and runs of heterozygosity (ROHet) were identified. Finally, candidate genes in the ROH and ROHet regions were explored through a pathway analysis to assess their molecular role.

RESULTS

In total, we detected 10,341 ROH in the SAA genome, 11,063 ROH in the CAM genome, 12,250 ROH in the MUR genome, 8,939 ROH in the MAL genome, and 18,441 ROH in the SAR genome. Moreover, we identified 4,087 ROHet for SAA, 3,360 for CAM, 2,927 for MUR, 3,701 for MAL, and 3,576 for SAR, with SAR having the highest heterozygosity coefficient. Interestingly, when computing the inbreeding coefficient using homozygous segment (FROH), SAA showed the lowest value while MAL the highest one, suggesting the need to improve selecting strategies to preserve genetic diversity within the population. Among the most significant candidate genes, we identified several ones linked to different physiological functions, such as milk production (e.g., DGAT1, B4GALT1), immunity (GABARAP, GPS2) and adaptation to environment (e.g., GJA3, GJB2 and GJB6).

CONCLUSIONS

This study highlighted the genetic diversity within and among five goat breeds. The high levels of ROH identified in some breeds might indicate high levels of inbreeding and a lack in genetic variation, which might negatively impact the animal population. Conversely, high levels of ROHet might indicate regions of the genetic diversity, beneficial for breed health and resilience. Therefore, these findings could aid breeding programs in managing inbreeding and preserving genetic diversity.

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.

PDGFRα signaling regulates Srsf3 transcript binding to affect PI3K signaling and endosomal trafficking

Signaling through the platelet-derived growth factor receptor alpha (PDGFRα) plays a critical role in craniofacial development. Phosphatidylinositol 3-kinase (PI3K)/Akt is the primary effector of PDGFRα signaling during mouse skeletal development. We previously demonstrated that Akt phosphorylates the RNA-binding protein serine/arginine-rich splicing factor 3 (Srsf3) downstream of PI3K-mediated PDGFRα signaling in mouse embryonic palatal mesenchyme (MEPM) cells, leading to its nuclear translocation. We further showed that ablation of Srsf3 in the murine neural crest lineage results in severe midline facial clefting and widespread alternative RNA splicing (AS) changes. Here, we demonstrated via enhanced UV-crosslinking and immunoprecipitation of MEPM cells that PDGF-AA stimulation leads to preferential binding of Srsf3 to exons and loss of binding to canonical Srsf3 CA-rich motifs. Through the analysis of complementary RNA-seq data, we showed that Srsf3 activity results in the preferential inclusion of exons with increased GC content and lower intron to exon length ratio. We found that Srsf3 activity downstream of PDGFRα signaling leads to retention of the receptor in early endosomes and increases in downstream PI3K-mediated Akt signaling. Taken together, our findings reveal that growth factor-mediated phosphorylation of an RNA-binding protein underlies gene expression regulation necessary for mammalian craniofacial development.

Down-regulation of ESRP2 inhibits breast cancer cell proliferation via inhibiting cyclinD1

Epithelial splicing regulatory protein 2 (ESRP2),an important alternative splicing protein of mRNA, is reported to have a dual role in tumors, which can promote or inhibit the occurrence and development of tumors. However, the function and mechanism of ESRP2 in breast cancer (BC) remain unclear. The distribution of ESRP2 expression in breast cancer and the correlation between ESRP2 expression and the overall survival rate were detected by The Cancer Genome Atlas (TCGA) database. Gene Ontology(GO)analysis, containing biological process, cellular components, and molecular function, was utilized to evaluate the potential mechanism of ESRP2 in breast cancer. The ESRP2 expression in breast cancer cell lines was detected by real-time quantitative PCR analysis (RT-qPCR) and western blotting. Cell clone was performed to examine the proliferation of ESRP2 knockdown in MCF-7 cells. The cell cycle was measured by flow cytometry assays. The role of ESRP2 knockdown in synergistic effect with chemotherapeutic agents was also determined by MTT assay. Bioinformatics analysis demonstrated that the ESRP2 gene was elevated in breast cancer cells and its overexpression was strongly correlated with shorter overall survival. GO analysis revealed that ESRP2 expression was related to cell proliferation. ESRP2 mRNA and protein expression were elevated in breast cancer cell lines, compared to the normal human breast cell line MCF-10 A. Dwon-regulation of ESRP2 inhibited cell proliferation and promoted the sensitivity of chemotherapy drug, Cisplatin(DDP) and Paclitaxel (TAXOL), in MCF-7 cells.Additionally, ESRP2 knockdown obstructed the cell cycle at the G1 phase and caused a decrease in cyclinD1 protein expression. These findings reveal that ESRP2 is highly expressed in breast cancer and is correlated with poor prognosis in breast cancer patients. ESRP2 knockdown can inhibit MCF-7 cell proliferation by arresting the cell cycle at the G1 phase and promoting the sensitivity of chemotherapy drugs (DDP and TAXOL)in MCF-7 cells. ESRP2 may be required for the regulation of breast cancer progression, as well as a critical target for the clinical treatment of breast cancer.

Calcium deposition in chicken eggshells: role of host genetics and gut microbiota

Eggshell is predominantly composed of calcium carbonate, making up about 95% of its composition. Eggshell quality is closely related to the amount of calcium deposition in the shell, which requires chickens to maintain a robust state of calcium metabolism. In this study, we introduced a novel parameter, Total Eggshell Weight (TESW), which measures the total weight of eggshells produced by chickens over a period of 10 consecutive d, providing valuable information on the intensity of calcium metabolism in chickens. Genome-wide association study (GWAS) was conducted to explore the genetic determinants of eggshell calcification in a population of 570 Rhode Island Red laying hens at 90 wk of age. This study revealed a significant association between a specific SNP (rs14249431) and TESW. Additionally, using random forest modeling and 2-tailed testing, we identified 3 genera, Lactobacillus in the jejunum, Lactobacillus, and Fournierella in the cecum, that exhibited a significant association with TESW. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis of claudin-1 and occludin genes in individuals with low TESW and high abundance of jejunal Lactobacillus confirmed that the inhibitory effect of jejunal Lactobacillus on calcium uptake was achieved through the up-regulation of tight junctions in intestinal epithelial cells. Notably, both host and microbial factors influence TESW, displaying a mutually influential relationship between them. The microbiome-wide Genome-Wide Association Study (mb-GWAS) identified significant associations between these 3 genera and specific genomic variants, such as rs316115020 and rs316420452 on chromosome 5, rs313198529 on chromosome 11, linked to Lactobacillus in the cecum. Moreover, rs312552529 on chromosome 1 exhibited potential association with Fournierella in the cecum. This study highlights the influence of host genetics and gut microbiota on calcium deposition in eggshells during the late laying phase, providing a foundational reference for studying calcium metabolism in hens.

PDGFRα signaling regulates Srsf3 transcript binding to affect PI3K signaling and endosomal trafficking

Signaling through the platelet-derived growth factor receptor alpha (PDGFRα) plays a critical role in craniofacial development, as mutations in PDGFRA are associated with cleft lip/palate in humans and Pdgfra mutant mouse models display varying degrees of facial clefting. Phosphatidylinositol 3-kinase (PI3K)/Akt is the primary effector of PDGFRα signaling during skeletal development in the mouse. We previously demonstrated that Akt phosphorylates the RNA-binding protein serine/arginine-rich splicing factor 3 (Srsf3) downstream of PI3K-mediated PDGFRα signaling in mouse embryonic palatal mesenchyme (MEPM) cells, leading to its nuclear translocation. We further showed that ablation of Srsf3 in the murine neural crest lineage results in severe midline facial clefting, due to defects in proliferation and survival of cranial neural crest cells, and widespread alternative RNA splicing (AS) changes. Here, we sought to determine the molecular mechanisms by which Srsf3 activity is regulated downstream of PDGFRα signaling to control AS of transcripts necessary for craniofacial development. We demonstrated via enhanced UV-crosslinking and immunoprecipitation (eCLIP) of MEPM cells that PDGF-AA stimulation leads to preferential binding of Srsf3 to exons and loss of binding to canonical Srsf3 CA-rich motifs. Through the analysis of complementary RNA-seq data, we showed that Srsf3 activity results in the preferential inclusion of exons with increased GC content and lower intron to exon length ratio. Moreover, we found that the subset of transcripts that are bound by Srsf3 and undergo AS upon PDGFRα signaling commonly encode regulators of PI3K signaling and early endosomal trafficking. Functional validation studies further confirmed that Srsf3 activity downstream of PDGFRα signaling leads to retention of the receptor in early endosomes and increases in downstream PI3K-mediated Akt signaling. Taken together, our findings reveal that growth factor-mediated phosphorylation of an RNA-binding protein underlies gene expression regulation necessary for mammalian craniofacial development.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

ESRP2-microRNA-122 axis directs the postnatal onset of liver polyploidization and maturation

Hepatocyte polyploidy and maturity are critical to acquiring specialized liver functions. Multiple intra- and extracellular factors influence ploidy, but how they cooperate temporally to steer liver polyploidization and maturation or how post-transcriptional mechanisms integrate into these paradigms is unknown. Here, we identified an important regulatory hierarchy in which postnatal activation of Epithelial-Splicing-Regulatory-Protein-2 (ESRP2) stimulates biogenesis of liver-specific microRNA (miR-122), thereby facilitating polyploidization, maturation, and functional competence of hepatocytes. By determining transcriptome-wide protein-RNA interactions in vivo and integrating them with single-cell and bulk hepatocyte RNA-seq datasets, we delineate an ESRP2-driven RNA processing program that drives sequential replacement of fetal-to-adult transcript isoforms. Specifically, ESRP2 binds the primary miR-122 host gene transcript to promote its processing/biogenesis. Combining constitutive and inducible ESRP2 gain- and loss-of-function mice models with miR-122 rescue experiments, we demonstrate that timed activation of ESRP2 augments miR-122-driven program of cytokinesis failure, ensuring proper onset and extent of hepatocyte polyploidization.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

ACLY alternative splicing correlates with cancer phenotypes

ATP-citrate lyase (ACLY) links carbohydrate and lipid metabolism and provides nucleocytosolic acetyl-CoA for protein acetylation. ACLY has two major splice isoforms: the full-length canonical "long" isoform and an uncharacterized "short" isoform in which exon 14 is spliced out. Exon 14 encodes 10 amino acids within an intrinsically disordered region and includes at least one dynamically phosphorylated residue. Both isoforms are expressed in healthy tissues to varying degrees. Analysis of human transcriptomic data revealed that the percent spliced in (PSI) of exon 14 is increased in several cancers and correlated with poorer overall survival in a pan-cancer analysis, though not in individual tumor types. This prompted us to explore potential biochemical and functional differences between ACLY isoforms. Here, we show that there are no discernible differences in enzymatic activity or stability between isoforms or phosphomutants of ACLY in vitro. Similarly, both isoforms and phosphomutants were able to rescue ACLY functions, including fatty acid synthesis and bulk histone acetylation, when re-expressed in Acly knockout cells. Deletion of Acly exon 14 in mice did not overtly impact development or metabolic physiology nor did it attenuate tumor burden in a genetic model of intestinal cancer. Notably, expression of epithelial splicing regulatory protein 1 (ESRP1) is highly correlated with ACLY PSI. We report that ACLY splicing is regulated by ESRP1. In turn, both ESRP1 expression and ACLY PSI are correlated with specific immune signatures in tumors. Despite these intriguing patterns of ACLY splicing in healthy and cancer tissues, functional differences between the isoforms remain elusive.

Alternative splicing in EMT and TGF-β signaling during cancer progression

Epithelial to mesenchymal transition (EMT) is a physiological process during development where epithelial cells transform to acquire mesenchymal characteristics, which allows them to migrate and colonize secondary tissues. Many cellular signaling pathways and master transcriptional factors exert a myriad of controls to fine tune this vital process to meet various developmental and physiological needs. Adding to the complexity of this network are post-transcriptional and post-translational regulations. Among them, alternative splicing has been shown to play important roles to drive EMT-associated phenotypic changes, including actin cytoskeleton remodeling, cell-cell junction changes, cell motility and invasiveness. In advanced cancers, transforming growth factor-β (TGF-β) is a major inducer of EMT and is associated with tumor cell metastasis, cancer stem cell self-renewal, and drug resistance. This review aims to provide an overview of recent discoveries regarding alternative splicing events and the involvement of splicing factors in the EMT and TGF-β signaling. It will emphasize the importance of various splicing factors involved in EMT and explore their regulatory mechanisms.

Characterization of Alternative Splicing in High-Risk Wilms' Tumors

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.

Alternative mRNA splicing events and regulators in epidermal differentiation

Alternative splicing (AS) of messenger RNAs occurs in ∼95% of multi-exon human genes and generates diverse RNA and protein isoforms. We investigated AS events associated with human epidermal differentiation, a process crucial for skin function. We identified 6,413 AS events, primarily involving cassette exons. We also predicted 34 RNA-binding proteins (RBPs) regulating epidermal AS, including 19 previously undescribed candidate regulators. From these results, we identified FUS as an RBP that regulates the balance between keratinocyte proliferation and differentiation. Additionally, we characterized the function of a cassette exon AS event in MAP3K7, which encodes a kinase involved in cell signaling. We found that a switch from the short to long isoform of MAP3K7, triggered during differentiation, enforces the demarcation between proliferating basal progenitors and overlying differentiated strata. Our findings indicate that AS occurs extensively in the human epidermis and has critical roles in skin homeostasis.

ESRP1 controls biogenesis and function of a large abundant multiexon circRNA

While the majority of circRNAs are formed from infrequent back-splicing of exons from protein coding genes, some can be produced at quite high level and in a regulated manner. We describe the regulation, biogenesis and function of circDOCK1(2-27), a large, abundant circular RNA that is highly regulated during epithelial-mesenchymal transition (EMT) and whose formation depends on the epithelial splicing regulator ESRP1. CircDOCK1(2-27) synthesis in epithelial cells represses cell motility both by diverting transcripts from DOCK1 mRNA production to circRNA formation and by direct inhibition of migration by the circRNA. HITS-CLIP analysis and CRISPR-mediated deletions indicate ESRP1 controls circDOCK1(2-27) biosynthesis by binding a GGU-containing repeat region in intron 1 and detaining its splicing until Pol II completes its 157 kb journey to exon 27. Proximity-dependent biotinylation (BioID) assay suggests ESRP1 may modify the RNP landscape of intron 1 in a way that disfavours communication of exon 1 with exon 2, rather than physically bridging exon 2 to exon 27. The X-ray crystal structure of RNA-bound ESRP1 qRRM2 domain reveals it binds to GGU motifs, with the guanines embedded in clamp-like aromatic pockets in the protein.

Proteomic analysis of pulmonary arteries and lung tissues from dogs affected with pulmonary hypertension secondary to degenerative mitral valve disease

In dogs with degenerative mitral valve disease (DMVD), pulmonary hypertension (PH) is a common complication characterized by abnormally elevated pulmonary arterial pressure (PAP). Pulmonary arterial remodeling is the histopathological changes of pulmonary artery that has been recognized in PH. The underlying mechanisms that cause this arterial remodeling are poorly understood. This study aimed to perform shotgun proteomics to investigate changes in protein expression in pulmonary arteries and lung tissues of DMVD dogs with PH compared to normal control dogs and DMVD dogs without PH. Tissue samples were collected from the carcasses of 22 small-sized breed dogs and divided into three groups: control (n = 7), DMVD (n = 7) and DMVD+PH groups (n = 8). Differentially expressed proteins were identified, and top three upregulated and downregulated proteins in the pulmonary arteries of DMVD dogs with PH including SIK family kinase 3 (SIK3), Collagen type I alpha 1 chain (COL1A1), Transforming growth factor alpha (TGF-α), Apoptosis associated tyrosine kinase (AATYK), Hepatocyte growth factor activator (HGFA) and Tyrosine-protein phosphatase non-receptor type 13 (PTPN13) were chosen. Results showed that some of the identified proteins may play a role in the pathogenesis of pulmonary arterial remodeling. This study concluded shotgun proteomics has potential as a tool for exploring candidate proteins associated with the pathogenesis of PH secondary to DMVD in dogs.

Spatial and Temporal Relationship between Epithelial-Mesenchymal Transition (EMT) and Stem Cells in Cancer

BACKGROUND

Epithelial-mesenchymal transition (EMT) is often linked with carcinogenesis. However, EMT is also important for embryo development and only reactivates in cancer. Connecting how EMT occurs during embryonic development and in cancer could help us further understand the root mechanisms of cancer diseases.

CONTENT

There are key regulatory elements that contribute to EMT and the induction and maintenance of stem cell properties during embryogenesis, tissue regeneration, and carcinogenesis. Here, we explore the implications of EMT in the different stages of embryogenesis and tissue development. We especially highlight the necessity of EMT in the mesodermal formation and in neural crest cells. Through EMT, these cells gain epithelial-mesenchymal plasticity (EMP). With this transition, crucial morphological changes occur to progress through the metastatic cascade as well as tissue regeneration after an injury. Stem-like cells, including cancer stem cells, are generated from EMT and during this process upregulate factors necessary for stem cell maintenance. Hence, it is important to understand the key regulators allowing stem cell awakening in cancer, which increases plasticity and promotes treatment resistance, to develop strategies targeting this cell population and improve patient outcomes.

SUMMARY

EMT involves multifaceted regulation to allow the fluidity needed to facilitate adaptation. This regulatory mechanism, plasticity, involves many cooperating transcription factors. Additionally, posttranslational modifications, such as splicing, activate the correct isoforms for either epithelial or mesenchymal specificity. Moreover, epigenetic regulation also occurs, such as acetylation and methylation. Downstream signaling ultimately results in the EMT which promotes tissue generation/regeneration and cancer progression.

Loss of ESRP2 Activates TAK1-MAPK Signaling through the Fetal RNA-Splicing Program to Promote Hepatocellular Carcinoma Progression

Tumors usually display fetal-like characteristics, and many oncofetal proteins have been identified. However, fetal-like reprogramming of RNA splicing in hepatocellular carcinoma (HCC) is poorly understood. Here, it is demonstrated that the expression of epithelial splicing regulatory protein 2 (ESRP2), an RNA splicing factor, is suppressed in fetal hepatocytes and HCC, in parallel with tumor progression. By combining RNA-Seq with splicing analysis, it is identified that ESRP2 controls the fetal-to-adult switch of multiple splice isoforms in HCC. Functionally, ESRP2 suppressed cell proliferation and migration by specifically switching the alternative splicing (AS) of the TAK1 gene and restraining the expression of the fetal and oncogenic isoform, TAK1_ΔE12. Notably, aberrant TAK1 splicing led to the activation of p38MAPK signaling and predicted poor prognosis in HCC patients. Further investigation revealed that TAK1_ΔE12 protein interacted closely with TAB3 and formed liquid condensation in HCC cells, resulting in p38MAPK activation, enhanced cell migration, and accelerated tumorigenesis. Loss of ESRP2 sensitized HCC cells to TAK1 kinase inhibitor (TAK1i), promoting pyroptotic cell death and CD8+ T cell infiltration. Combining TAK1i with immune checkpoint therapy achieved potent tumor regression in mice. Overall, the findings reveal a previously unexplored onco-fetal reprogramming of RNA splicing and provide novel therapeutic avenues for HCC.

Epithelial specific splicing regulator proteins as emerging oncogenes in aggressive prostate cancer

Prostate cancer progression is connected to the activity of conventional oncogenes and tumour suppressors and driven by circulating steroid hormones. A key issue has been how to identify and care for aggressively developing prostate tumours. Here we discuss how expression of the splicing regulators ESRP1 and ESRP2, and how their role as "masterminds" of epithelial splicing patterns, have been identified as markers of aggressively proliferating prostate primary tumours. We suggest that the origin of prostate cancer within epithelial cells, and the subsequent association of ESRP1 and ESRP2 expression with more aggressive disease progression, identify ESRP1 and ESRP2 as lineage survival oncogenes. To move this field on in the future it will be important to identify the gene expression targets controlled by ESRP1/2 that regulate prostate cancer proliferation. Potential future therapies could be designed to target ESRP1 and ESRP2 protein activity or their regulated splice isoforms in aggressive prostate tumours. Design of these therapies is potentially complicated by the risk of producing a more mesenchymal splicing environment that might promote tumour metastasis.

Differential exon usage of developmental genes is associated with deregulated epigenetic marks

Alternative exon usage is known to affect a large portion of genes in mammalian genomes. Importantly, different splice isoforms sometimes possess distinctly different protein functions. Here, we analyzed data from the Human Epigenome Atlas for 11 different human adult tissues and for 8 cultured cells that mimic early developmental stages. We found a significant enrichment of cases where differential usage of exons in various developmental stages of human cells and tissues is associated with differential epigenetic modifications in the flanking regions of individual exons. Many of the genes that were differentially regulated at the exon level and showed deregulated histone marks at the respective exon flanks are functionally associated with development and metabolism.

Therapeutic Targeting of RNA Splicing in Cancer

RNA splicing is a key regulatory step in the proper control of gene expression. It is a highly dynamic process orchestrated by the spliceosome, a macro-molecular machinery that consists of protein and RNA components. The dysregulation of RNA splicing has been observed in many human pathologies ranging from neurodegenerative diseases to cancer. The recent identification of recurrent mutations in the core components of the spliceosome in hematologic malignancies has advanced our knowledge of how splicing alterations contribute to disease pathogenesis. This review article will discuss our current understanding of how aberrant RNA splicing regulation drives tumor initiation and progression. We will also review current therapeutic modalities and highlight emerging technologies designed to target RNA splicing for cancer treatment.

The discovery, function, and regulation of epithelial splicing regulatory proteins (ESRP) 1 and 2

Alternative splicing is a broad and evolutionarily conserved mechanism to diversify gene expression and functionality. The process relies on RNA binding proteins (RBPs) to recognize and bind target sequences in pre-mRNAs, which allows for the inclusion or skipping of various alternative exons. One recently discovered family of RBPs is the epithelial splicing regulatory proteins (ESRP) 1 and 2. Here, we discuss the structure and physiological function of the ESRPs in a variety of contexts. We emphasize the current understanding of their splicing activities, using the classic example of fibroblast growth factor receptor 2 mutually exclusive splicing. We also describe the mechanistic roles of ESRPs in coordinating the splicing and functional output of key signaling pathways that support the maintenance of, or shift between, epithelial and mesenchymal cell states. In particular, we highlight their functions in the development of mammalian limbs, the inner ear, and craniofacial structure while discussing the genetic and biochemical evidence that showcases their conserved roles in tissue regeneration, disease, and cancer pathogenesis.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

Transcriptomic analysis of the upper lip and primary palate development in mice

Background: Normal fusion of the upper lip and primary palate is a complex process involving a series of characteristic and orderly regulated cellular events. Cleft lip with or without palate (CL/P), one of the most common congenital malformations, may be induced by abnormalities in any of these events. However, less is known about the precise regulatory process in the fusion of the upper lip and primary palate. Methods: Lambdoidal junction tissues of mice from embryonic days 10.5, 11.5, and 12.5- three key fusion stages-were acquired for RNA sequencing. Results: Gene expression profiles in distinct fusion stages of mice were identified. Some of the differentially expressed genes (DEGs) have been reported to affect upper lip and primary palate development. However, other DEGs, such as Krt5, Pax1, Ambn, Hey2, and Tnmd, have not previously been investigated. Gene set enrichment analysis (GSEA) of these DEGs revealed the sequential intensification of Wnt, PI3K-Akt, MAPK, Hippo, and TGF-beta signaling pathways and identified relatively highly expressed genes including Tnn, Wnt3a, and Wnt16. We also observed substantial alternative splicing events during the fusion of the upper lip and primary palate and identified potentially important genes including Gtpbp8, Armcx1, Tle3, and Numa1. Protein-protein interaction (PPI) network analysis identified a series of hub genes, including Col1a2, Fos, Bmp2, Shh, Col1a1, Wnt3a, Anxa1, Gem, etc. Conclusion: Overall, the results of this study provided a comprehensive analysis of the development of the upper lip and primary palate. Our work provides insight into future studies of normal upper lip and primary palate development and the etiology of CL/P.

Downregulated ESRP1/2 promotes lung metastasis of bladder carcinoma through altering FGFR2 splicing and macrophage polarization

Introduction

Lung metastasis occurs in parts of the bladder carcinoma (BC) patients but represents the highest severity and a poor outcome of the disease. The molecular mechanism underlying lung metastasis of BC is not fully understood. Fibroblast growth factor receptor 2 (FGFR2) signaling plays a substantial role in the BC cell growth and invasion. In this study, we assessed the regulation of the alternative splicing of FGFR2 by epithelial splicing regulatory proteins (ESRPs) in lung metastasis of BC.

Methods

Gene profile of BC in comparison with adjacent non-tumor bladder tissue was obtained from GEO public database to analyze the levels of differentiated genes and pathways. Moreover, the association of ESRP1 or ESRP2 with lung metastasis of BC was analyzed on our own clinic samples. The effects of altered expression of ESRP1 or ESRP2 on alternative splicing of FGFR2 IIIb and IIIc, which represents epithelial and mesenchymal-like splicing, were analyzed on BC cell lines T24 and RT4. The in vivo effects of ESRP1 or ESRP2 on lung metastasis of BC were assessed in mice subcutaneously grafted with ESRP1/2-modified BC labeled with fluorescent and luciferase reporters.

Results

We detected significant reduction of ESRP1 and ESRP2 in BC in public database of BC specimens. Moreover, analysis on our own specimens also showed strong downregulation of ESRP1 or ESRP2 in BC, and the latter was more pronounced in cases with lung metastasis. In vitro, altered levels of ESRP1 or ESRP2 caused a switch of FGFR2 splicing between FGFR2-IIIb and FGFR2-IIIc, resulting in changes in tumor cell growth and metastatic potential. In vivo, re-expression of ESRP1 or ESRP2 in BC cells not only inhibited the growth of the xenografted tumor formation in nude mice, but also reduced the occurrence of lung metastasis, partially through altering polarization of tumor-associated macrophages.

Conclusion

Our data thus suggest that reduction in ESRP1 or ESRP2 promotes lung metastasis of BC through altering FGFR2 splicing and macrophage polarization.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.

Characterization of alternative mRNA splicing in cultured cell populations representing progressive stages of human fetal kidney development

Nephrons are the functional units of the kidney. During kidney development, cells from the cap mesenchyme-a transient kidney-specific progenitor state-undergo a mesenchymal to epithelial transition (MET) and subsequently differentiate into the various epithelial cell types that create the tubular structures of the nephron. Faults in this transition can lead to a pediatric malignancy of the kidney called Wilms' tumor that mimics normal kidney development. While human kidney development has been characterized at the gene expression level, a comprehensive characterization of alternative splicing is lacking. Therefore, in this study, we performed RNA sequencing on cell populations representing early, intermediate, and late developmental stages of the human fetal kidney, as well as three blastemal-predominant Wilms' tumor patient-derived xenografts. Using this newly generated RNAseq data, we identified a set of transcripts that are alternatively spliced between the different developmental stages. Moreover, we found that cells from the earliest developmental stage have a mesenchymal splice-isoform profile that is similar to that of blastemal-predominant Wilms' tumor xenografts. RNA binding motif enrichment analysis suggests that the mRNA binding proteins ESRP1, ESRP2, RBFOX2, and QKI regulate alternative mRNA splicing during human kidney development. These findings illuminate new molecular mechanisms involved in human kidney development and pediatric kidney cancer.

The global Protein-RNA interaction map of ESRP1 defines a post-transcriptional program that is essential for epithelial cell function

The epithelial splicing regulatory proteins, ESRP1 and ESRP2, are essential for mammalian development through the regulation of a global program of alternative splicing of genes involved in the maintenance of epithelial cell function. To further inform our understanding of the molecular functions of ESRP1, we performed enhanced crosslinking immunoprecipitation coupled with high-throughput sequencing (eCLIP) in epithelial cells of mouse epidermis. The genome-wide binding sites of ESRP1 were integrated with RNA-Seq analysis of alterations in splicing and total gene expression that result from epidermal ablation of Esrp1 and Esrp2. These studies demonstrated that ESRP1 functions in splicing regulation occur primarily through direct binding in a position-dependent manner to promote either exon inclusion or skipping. In addition, we also identified widespread binding of ESRP1 in 3' and 5' untranslated regions (UTRs) of genes involved in epithelial cell function, suggesting that its post-transcriptional functions extend beyond splicing regulation.

Emerging Roles of RNA-Binding Proteins in Inner Ear Hair Cell Development and Regeneration

RNA-binding proteins (RBPs) regulate gene expression at the post-transcriptional level. They play major roles in the tissue- and stage-specific expression of protein isoforms as well as in the maintenance of protein homeostasis. The inner ear is a bi-functional organ, with the cochlea and the vestibular system required for hearing and for maintaining balance, respectively. It is relatively well documented that transcription factors and signaling pathways are critically involved in the formation of inner ear structures and in the development of hair cells. Accumulating evidence highlights emerging functions of RBPs in the post-transcriptional regulation of inner ear development and hair cell function. Importantly, mutations of splicing factors of the RBP family and defective alternative splicing, which result in inappropriate expression of protein isoforms, lead to deafness in both animal models and humans. Because RBPs are critical regulators of cell proliferation and differentiation, they present the potential to promote hair cell regeneration following noise- or ototoxin-induced damage through mitotic and non-mitotic mechanisms. Therefore, deciphering RBP-regulated events during inner ear development and hair cell regeneration can help define therapeutic strategies for treatment of hearing loss. In this review, we outline our evolving understanding of the implications of RBPs in hair cell formation and hearing disease with the aim of promoting future research in this field.

Antigens Expressed by Breast Cancer Cells Undergoing EMT Stimulate Cytotoxic CD8+ T Cell Immunity

Antigenic differences formed by alterations in gene expression and alternative splicing are predicted in breast cancer cells undergoing epithelial to mesenchymal transition (EMT) and the reverse plasticity known as MET. How these antigenic differences impact immune interactions and the degree to which they can be exploited to enhance immune responses against mesenchymal cells is not fully understood. We utilized a master microRNA regulator of EMT to alter mesenchymal-like EO771 mammary carcinoma cells to a more epithelial phenotype. A computational approach was used to identify neoantigens derived from the resultant differentially expressed somatic variants (SNV) and alternative splicing events (neojunctions). Using whole cell vaccines and peptide-based vaccines, we find superior cytotoxicity against the more-epithelial cells and explore the potential of neojunction-derived antigens to elicit T cell responses through experiments designed to validate the computationally predicted neoantigens. Overall, results identify EMT-associated splicing factors common to both mouse and human breast cancer cells as well as immunogenic SNV- and neojunction-derived neoantigens in mammary carcinoma cells.

Revisiting the embryogenesis of lip and palate development

Clefts of the lip and palate (CLP), the major causes of congenital facial malformation globally, result from failure of fusion of the facial processes during embryogenesis. With a prevalence of 1 in 500-2500 live births, CLP causes major morbidity throughout life as a result of problems with facial appearance, feeding, speaking, obstructive apnoea, hearing and social adjustment and requires complex, multi-disciplinary care at considerable cost to healthcare systems worldwide. Long-term outcomes for affected individuals include increased mortality compared with their unaffected siblings. The frequent occurrence and major healthcare burden imposed by CLP highlight the importance of dissecting the molecular mechanisms driving facial development. Identification of the genetic mutations underlying syndromic forms of CLP, where CLP occurs in association with non-cleft clinical features, allied to developmental studies using appropriate animal models is central to our understanding of the molecular events underlying development of the lip and palate and, ultimately, how these are disturbed in CLP.

The X-linked splicing regulator MBNL3 has been co-opted to restrict placental growth in eutherians

Understanding the regulatory interactions that control gene expression during the development of novel tissues is a key goal of evolutionary developmental biology. Here, we show that Mbnl3 has undergone a striking process of evolutionary specialization in eutherian mammals resulting in the emergence of a novel placental function for the gene. Mbnl3 belongs to a family of RNA-binding proteins whose members regulate multiple aspects of RNA metabolism. We find that, in eutherians, while both Mbnl3 and its paralog Mbnl2 are strongly expressed in placenta, Mbnl3 expression has been lost from nonplacental tissues in association with the evolution of a novel promoter. Moreover, Mbnl3 has undergone accelerated protein sequence evolution leading to changes in its RNA-binding specificities and cellular localization. While Mbnl2 and Mbnl3 share partially redundant roles in regulating alternative splicing, polyadenylation site usage and, in turn, placenta maturation, Mbnl3 has also acquired novel biological functions. Specifically, Mbnl3 knockout (M3KO) alone results in increased placental growth associated with higher Myc expression. Furthermore, Mbnl3 loss increases fetal resource allocation during limiting conditions, suggesting that location of Mbnl3 on the X chromosome has led to its role in limiting placental growth, favoring the maternal side of the parental genetic conflict.

Application of zebrafish in the study of craniomaxillofacial developmental anomalies

Craniomaxillofacial developmental anomalies are one of the most prevalent congenital defects worldwide and could result from any disruption of normal development processes, which is generally influenced by interactions between genes and the environment. Currently, with the advances in genetic screening strategies, an increasing number of novel variants and their roles in orofacial diseases have been explored. Zebrafish is recognized as a powerful animal model, and its homologous genes and similar oral structure and development process provide an ideal platform for studying the contributions of genetic and environmental factors to human craniofacial malformations. Here, we reviewed zebrafish models for the study of craniomaxillofacial developmental anomalies, such as human nonsyndromic cleft lip with or without an affected palate and jaw and tooth developmental anomalies. Due to its potential for gene expression and regulation research, zebrafish may provide new perspectives for understanding craniomaxillofacial diseaseand its treatment.

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.

Mouse models in palate development and orofacial cleft research: Understanding the crucial role and regulation of epithelial integrity in facial and palate morphogenesis

Cleft lip and cleft palate are common birth defects resulting from genetic and/or environmental perturbations of facial development in utero. Facial morphogenesis commences during early embryogenesis, with cranial neural crest cells interacting with the surface ectoderm to form initially partly separate facial primordia consisting of the medial and lateral nasal prominences, and paired maxillary and mandibular processes. As these facial primordia grow around the primitive oral cavity and merge toward the ventral midline, the surface ectoderm undergoes a critical differentiation step to form an outer layer of flattened and tightly connected periderm cells with a non-stick apical surface that prevents epithelial adhesion. Formation of the upper lip and palate requires spatiotemporally regulated inter-epithelial adhesions and subsequent dissolution of the intervening epithelial seam between the maxillary and medial/lateral nasal processes and between the palatal shelves. Proper regulation of epithelial integrity plays a paramount role during human facial development, as mutations in genes encoding epithelial adhesion molecules and their regulators have been associated with syndromic and non-syndromic orofacial clefts. In this chapter, we summarize mouse genetic studies that have been instrumental in unraveling the mechanisms regulating epithelial integrity and periderm differentiation during facial and palate development. Since proper epithelial integrity also plays crucial roles in wound healing and cancer, understanding the mechanisms regulating epithelial integrity during facial development have direct implications for improvement in clinical care of craniofacial patients.

Single-cell resolution of MET- and EMT-like programs in osteoblasts during zebrafish fin regeneration

Zebrafish regenerate fin rays following amputation through epimorphic regeneration, a process that has been proposed to involve the epithelial-to-mesenchymal transition (EMT). We performed single-cell RNA sequencing (scRNA-seq) to elucidate osteoblastic transcriptional programs during zebrafish caudal fin regeneration. We show that osteoprogenitors are enriched with components associated with EMT and its reverse, mesenchymal-to-epithelial transition (MET), and provide evidence that the EMT markers cdh11 and twist2 are co-expressed in dedifferentiating cells at the amputation stump at 1 dpa, and in differentiating osteoblastic cells in the regenerate, the latter of which are enriched in EMT signatures. We also show that esrp1, a regulator of alternative splicing in epithelial cells that is associated with MET, is expressed in a subset of osteoprogenitors during outgrowth. This study provides a single cell resource for the study of osteoblastic cells during zebrafish fin regeneration, and supports the contribution of MET- and EMT-associated components to this process.

The epithelial splicing regulator ESRP2 is epigenetically repressed by DNA hypermethylation in Wilms tumour and acts as a tumour suppressor

Wilms tumour (WT), an embryonal kidney cancer, has been extensively characterised for genetic and epigenetic alterations, but a proportion of WTs still lack identifiable abnormalities. To uncover DNA methylation changes critical for WT pathogenesis, we compared the epigenome of foetal kidney with two WT cell lines, filtering our results to remove common cancer-associated epigenetic changes and to enrich for genes involved in early kidney development. This identified four hypermethylated genes, of which ESRP2 (epithelial splicing regulatory protein 2) was the most promising for further study. ESRP2 was commonly repressed by DNA methylation in WT, and this occurred early in WT development (in nephrogenic rests). ESRP2 expression was reactivated by DNA methyltransferase inhibition in WT cell lines. When ESRP2 was overexpressed in WT cell lines, it inhibited cellular proliferation in vitro, and in vivo it suppressed tumour growth of orthotopic xenografts in nude mice. RNA-seq of the ESRP2-expressing WT cell lines identified several novel splicing targets. We propose a model in which epigenetic inactivation of ESRP2 disrupts the mesenchymal to epithelial transition in early kidney development to generate WT.

Orofacial Cleft and Mandibular Prognathism-Human Genetics and Animal Models

Many complex molecular interactions are involved in the process of craniofacial development. Consequently, the network is sensitive to genetic mutations that may result in congenital malformations of varying severity. The most common birth anomalies within the head and neck are orofacial clefts (OFCs) and prognathism. Orofacial clefts are disorders with a range of phenotypes such as the cleft of the lip with or without cleft palate and isolated form of cleft palate with unilateral and bilateral variations. They may occur as an isolated abnormality (nonsyndromic-NSCLP) or coexist with syndromic disorders. Another cause of malformations, prognathism or skeletal class III malocclusion, is characterized by the disproportionate overgrowth of the mandible with or without the hypoplasia of maxilla. Both syndromes may be caused by the presence of environmental factors, but the majority of them are hereditary. Several mutations are linked to those phenotypes. In this review, we summarize the current knowledge regarding the genetics of those phenotypes and describe genotype-phenotype correlations. We then present the animal models used to study these defects.

Alternative RNA splicing in tumour heterogeneity, plasticity and therapy

Alternative splicing is a process by which a single gene is able to encode multiple different protein isoforms. It is regulated by the inclusion or exclusion of introns and exons that are joined in different patterns prior to protein translation, thus enabling transcriptomic and proteomic diversity. It is now widely accepted that alternative splicing is dysregulated across nearly all cancer types. This widespread dysregulation means that nearly all cellular processes are affected - these include processes synonymous with the hallmarks of cancer - evasion of apoptosis, tissue invasion and metastasis, altered cellular metabolism, genome instability and drug resistance. Emerging evidence indicates that the dysregulation of alternative splicing also promotes a permissive environment for increased tumour heterogeneity and cellular plasticity. These are fundamental regulators of a patient's response to therapy. In this Review, we introduce the mechanisms of alternative splicing and the role of aberrant splicing in cancer, with particular focus on newfound evidence of alternative splicing promoting tumour heterogeneity, cellular plasticity and altered metabolism. We discuss recent in vivo models generated to study alternative splicing and the importance of these for understanding complex tumourigenic processes. Finally, we review the effects of alternative splicing on immune evasion, cell death and genome instability, and how targeting these might enhance therapeutic efficacy.

Regulation of Alternative Splicing during Epithelial-Mesenchymal Transition

Alternative splicing is an essential mechanism of gene regulation, giving rise to remarkable protein diversity in higher eukaryotes. Epithelial-mesenchymal transition (EMT) is a developmental process that plays an essential role in metazoan embryogenesis. Recent studies have revealed that alternative splicing serves as a fundamental layer of regulation that governs cells to undergo EMT. In this review, we summarize recent findings on the functional impact of alternative splicing in EMT and EMT-associated activities. We then discuss the regulatory mechanisms that control alternative splicing changes during EMT.

Single-cell characterization of CRISPR-modified transcript isoforms with nanopore sequencing

We developed a single-cell approach to detect CRISPR-modified mRNA transcript structures. This method assesses how genetic variants at splicing sites and splicing factors contribute to alternative mRNA isoforms. We determine how alternative splicing is regulated by editing target exon-intron segments or splicing factors by CRISPR-Cas9 and their consequences on transcriptome profile. Our method combines long-read sequencing to characterize the transcript structure and short-read sequencing to match the single-cell gene expression profiles and gRNA sequence and therefore provides targeted genomic edits and transcript isoform structure detection at single-cell resolution.

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.

The RNA-Binding Protein ESRP1 Modulates the Expression of RAC1b in Colorectal Cancer Cells

RNA binding proteins are well recognized as critical regulators of tumorigenic processes through their capacity to modulate RNA biogenesis, including alternative splicing, RNA stability and mRNA translation. The RNA binding protein Epithelial Splicing Regulatory Protein 1 (ESRP1) can act as a tumor suppressor or promoter in a cell type- and disease context-dependent manner. We have previously shown that elevated expression of ESRP1 in colorectal cancer cells can drive tumor progression. To gain further insights into the pro-tumorigenic mechanism of action of ESRP1, we performed cDNA microarray analysis on two colorectal cells lines modulated for ESRP1 expression. Intriguingly, RAC1b was highly expressed, both at mRNA and protein levels, in ESRP1-overexpressing cells, while the opposite trend was observed in ESRP1-silenced CRC cells. Moreover, RAC1 and RAC1b mRNA co-immunoprecipitate with ESRP1 protein. Silencing of RAC1b expression significantly reduced the number of soft agar colonies formed by ESRP1-overexpressing cells, suggesting that ESRP1 acted, at least partially, through RAC1b in its tumor-promoting activities in CRC cells. Thus, our data provide molecular cues on targetable candidates in CRC cases with high ESRP1 expression.

Genome-wide association studies of preweaning growth and in vivo carcass composition traits in Esme sheep

Sheep are considered as a major contributor of global food security. Moreover, sheep preweaning growth traits as well as in vivo carcass composition traits such as ultrasonic measurements of Longissimus dorsi muscle depth (UMD) and back-fat thickness (UFD) are crucially important indicators of meat yield and hot carcass composition. Despite their relative importance for productivity and profitability of a sheep production system, detected QTL for these traits are quite scarce. Therefore, we implemented GWAS for these traits using animal mixed model-based association approach provided by GenABEL in Esme sheep. Three genome-wide and 14 individual chromosome-wide associated SNPs were discovered. As a result, ESRP1, LOC105613082, ZNF641, DUSP5, TEAD1, SMOX, PTPRT, RALYL, POM121C, PHIP, LOC101106051, ZIM3, PEG3, TRPC7, FBXL4, LOC105610397, LOC105616489 and DNAAF2 were suggested as candidates. Some of the discovered genes and involved pathways were already annotated to contribute growth and development in various species including human, mice and cattle. All in all, the results of this study are expected to strongly contribute to shed a light on the underlying molecular mechanisms behind growth and carcass composition traits, with potential implications on studies aiming faster genetic improvement, targeted low-resolution SNP panel designs and genome-editing studies.

Epithelial Splicing Regulatory Protein 1 Is Overexpressed in Breast Cancer and Predicts Poor Prognosis for Breast Cancer Patients

Background

Epithelial splicing regulatory proteins (ESRPs), including ESRP1 and ESRP2, are important proteins for alternative splicing of mRNAs and are reported to promote or inhibit the progression of some tumors. However, the effects of ESRPs in breast cancer are still unknown.

Material/Methods

In this study, we detected the transcriptional level and alterations of ESRP1 in patients with breast cancer based on the Oncomine, Gene Expression Profiling Interactive Analysis, Gene Expression-Based Outcome for Breast Cancer Online, and cBioPortal databases. Using immunohistochemistry and quantitative polymerase chain reaction, the expression pattern of ESRP1 in breast cancer was analyzed. Analysis of the clinicopathological characteristics and function of ESRP1 in breast cancer were actualized through the University of Alabama Cancer database and Database for Annotation, Visualization and Integrated Discovery. Using the Kaplan-Meier plotter, the prognostic values of ESRP1 in patients with breast cancer were analyzed. The Encyclopedia of RNA Interactomes database was used to predict miRNAs that regulated ESRP1.

Results

We found that ESRP1 was significantly overexpressed in patients with breast cancer, compared with patients without breast cancer, and had statistically significant clinicopathological characteristics. Kaplan-Meier plotter analysis indicated that the elevated expression of ESRP1 was associated with poor prognosis in patients with breast cancer. Furthermore, hsa-miR-181c-5p was identified to be potentially involved in the regulation of ESRP1.

Conclusions

These results suggest that ESRP1 is a valuable target for the precise treatment of breast cancer and a potential biomarker for the prognosis of patients with breast cancer.

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.