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

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.

•

Osteoblasts express EMT/MET signatures during zebrafish fin regeneration

•

De/re-differentiating osteoblasts express cdh11, an EMT marker

•

A subset of osteoprogenitors express the MET marker esrp1

•

Our scRNA-seq data can be explored online

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

Simple Summary

Colorectal cancer (CRC) ranks third for incidence and second for number of deaths among cancer types worldwide. Poor patient survival due to inadequate response to currently available treatment regimens points to the urgent requirement for personalized therapy in CRC patients. Our aim was to provide mechanistic insights into the pro-tumorigenic role of the RNA-binding protein ESRP1, which is highly expressed in a subset of CRC patients. We show that, in CRC cells, ESRP1 binds to and has the same trend in expression as RAC1b, a well-known tumor promoter. Thus, RAC1b may be a potential therapeutic target in ESRP1-overexpressing CRC.

Abstract

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.

A conditional multi-trait sequence GWAS discovers pleiotropic candidate genes and variants for sheep wool, skin wrinkle and breech cover traits

Background

Imputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to “fly strike” (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle.

Results

A stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed.

Conclusions

The wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00651-0.

Analysis of immune subtypes across the epithelial-mesenchymal plasticity spectrum

Epithelial-mesenchymal plasticity plays a critical role in many solid tumor types as a mediator of metastatic dissemination and treatment resistance. In addition, there is also a growing appreciation that the epithelial/mesenchymal status of a tumor plays a role in immune evasion and immune suppression. A deeper understanding of the immunological features of different tumor types has been facilitated by the availability of large gene expression datasets and the development of methods to deconvolute bulk RNA-Seq data. These resources have generated powerful new ways of characterizing tumors, including classification of immune subtypes based on differential expression of immunological genes. In the present work, we combine scoring algorithms to quantify epithelial-mesenchymal plasticity with immune subtype analysis to understand the relationship between epithelial plasticity and immune subtype across cancers. We find heterogeneity of epithelial-mesenchymal transition (EMT) status both within and between cancer types, with greater heterogeneity in the expression of EMT-related factors than of MET-related factors. We also find that specific immune subtypes have associated EMT scores and differential expression of immune checkpoint markers.

Alternative splicing associated with cancer stemness in kidney renal clear cell carcinoma

Backgroud

Cancer stemness is associated with metastases in kidney renal clear cell carcinoma (KIRC) and negatively correlates with immune infiltrates. Recent stemness evaluation methods based on the absolute expression have been proposed to reveal the relationship between stemness and cancer. However, we found that existing methods do not perform well in assessing the stemness of KIRC patients, and they overlooked the impact of alternative splicing. Alternative splicing not only progresses during the differentiation of stem cells, but also changes during the acquisition of the stemness features of cancer stem cells. There is an urgent need for a new method to predict KIRC-specific stemness more accurately, so as to provide help in selecting treatment options.

Methods

The corresponding RNA-Seq data were obtained from the The Cancer Genome Atlas (TCGA) data portal. We also downloaded stem cell RNA sequence data from the Progenitor Cell Biology Consortium (PCBC) Synapse Portal. Independent validation sets with large sample size and common clinic pathological characteristics were obtained from the Gene Expression Omnibus (GEO) database. we constructed a KIRC-specific stemness prediction model using an algorithm called one-class logistic regression based on the expression and alternative splicing data to predict stemness indices of KIRC patients, and the model was externally validated. We identify stemness-associated alternative splicing events (SASEs) by analyzing different alternative splicing event between high- and low- stemness groups. Univariate Cox and multivariable logistic regression analysisw as carried out to detect the prognosis-related SASEs respectively. The area under curve (AUC) of receiver operating characteristic (ROC) was performed to evaluate the predictive values of our model.

Results

Here, we constructed a KIRC-specific stemness prediction model with an AUC of 0.968,and to provide a user-friendly interface of our model for KIRC stemness analysis, we have developed KIRC Stemness Calculator and Visualization (KSCV), hosted on the Shiny server, can most easily be accessed via web browser and the url https://jiang-lab.shinyapps.io/kscv/. When applied to 605 KIRC patients, our stemness indices had a higher correlation with the gender, smoking history and metastasis of the patients than the previous stemness indices, and revealed intratumor heterogeneity at the stemness level. We identified 77 novel SASEs by dividing patients into high- and low- stemness groups with significantly different outcome and they had significant correlations with expression of 17 experimentally validated splicing factors. Both univariate and multivariate survival analysis demonstrated that SASEs closely correlated with the overall survival of patients.

Conclusions

Basing on the stemness indices, we found that not only immune infiltration but also alternative splicing events showed significant different at the stemness level. More importantly, we highlight the critical role of these differential alternative splicing events in poor prognosis, and we believe in the potential for their further translation into targets for immunotherapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08470-8.

Formative pluripotent stem cells show features of epiblast cells poised for gastrulation

The pluripotency of mammalian early and late epiblast could be recapitulated by naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), respectively. However, these two states of pluripotency may not be sufficient to reflect the full complexity and developmental potency of the epiblast during mammalian early development. Here we report the establishment of self-renewing formative pluripotent stem cells (fPSCs) which manifest features of epiblast cells poised for gastrulation. fPSCs can be established from different mouse ESCs, pre-/early-gastrula epiblasts and induced PSCs. Similar to pre-/early-gastrula epiblasts, fPSCs show the transcriptomic features of formative pluripotency, which are distinct from naïve ESCs and primed EpiSCs. fPSCs show the unique epigenetic states of E6.5 epiblast, including the super-bivalency of a large set of developmental genes. Just like epiblast cells immediately before gastrulation, fPSCs can efficiently differentiate into three germ layers and primordial germ cells (PGCs) in vitro. Thus, fPSCs highlight the feasibility of using PSCs to explore the development of mammalian epiblast.

A cell-to-patient machine learning transfer approach uncovers novel basal-like breast cancer prognostic markers amongst alternative splice variants

BACKGROUND

Breast cancer is amongst the 10 first causes of death in women worldwide. Around 20% of patients are misdiagnosed leading to early metastasis, resistance to treatment and relapse. Many clinical and gene expression profiles have been successfully used to classify breast tumours into 5 major types with different prognosis and sensitivity to specific treatments. Unfortunately, these profiles have failed to subclassify breast tumours into more subtypes to improve diagnostics and survival rate. Alternative splicing is emerging as a new source of highly specific biomarkers to classify tumours in different grades. Taking advantage of extensive public transcriptomics datasets in breast cancer cell lines (CCLE) and breast cancer tumours (TCGA), we have addressed the capacity of alternative splice variants to subclassify highly aggressive breast cancers.

RESULTS

Transcriptomics analysis of alternative splicing events between luminal, basal A and basal B breast cancer cell lines identified a unique splicing signature for a subtype of tumours, the basal B, whose classification is not in use in the clinic yet. Basal B cell lines, in contrast with luminal and basal A, are highly metastatic and express epithelial-to-mesenchymal (EMT) markers, which are hallmarks of cell invasion and resistance to drugs. By developing a semi-supervised machine learning approach, we transferred the molecular knowledge gained from these cell lines into patients to subclassify basal-like triple negative tumours into basal A- and basal B-like categories. Changes in splicing of 25 alternative exons, intimately related to EMT and cell invasion such as ENAH, CD44 and CTNND1, were sufficient to identify the basal-like patients with the worst prognosis. Moreover, patients expressing this basal B-specific splicing signature also expressed newly identified biomarkers of metastasis-initiating cells, like CD36, supporting a more invasive phenotype for this basal B-like breast cancer subtype.

CONCLUSIONS

Using a novel machine learning approach, we have identified an EMT-related splicing signature capable of subclassifying the most aggressive type of breast cancer, which are basal-like triple negative tumours. This proof-of-concept demonstrates that the biological knowledge acquired from cell lines can be transferred to patients data for further clinical investigation. More studies, particularly in 3D culture and organoids, will increase the accuracy of this transfer of knowledge, which will open new perspectives into the development of novel therapeutic strategies and the further identification of specific biomarkers for drug resistance and cancer relapse.

Silencing of SRRM4 suppresses microexon inclusion and promotes tumor growth across cancers

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3–27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Using data from The Cancer Genome Atlas, this study shows that the splicing factor SRRM4 and its program of differentiation-promoting microexons are downregulated across tumor types with remarkable consistency, providing tumors with a proliferative advantage.