Altered mRNA Splicing, Chondrocyte Gene Expression and Abnormal Skeletal Development due to SF3B4 Mutations in Rodriguez Acrofacial Dysostosis

Pharmacological Inhibition of the Spliceosome SF3b Complex by Pladienolide-B Elicits Craniofacial Developmental Defects in Mouse and Zebrafish

BACKGROUND

Mutations in genes encoding spliceosome components result in craniofacial structural defects in humans, referred to as spliceosomopathies. The SF3b complex is a crucial unit of the spliceosome, but model organisms generated through genetic modification of the complex do not perfectly mimic the phenotype of spliceosomopathies. Since the phenotypes are suggested to be determined by the extent of spliceosome dysfunction, an alternative experimental system that can seamlessly control SF3b function is needed.

METHODS

To establish another experimental system for model organisms elucidating relationship between spliceosome function and human diseases, we administered Pladienolide-B (PB), a SF3b complex inhibitor, to mouse and zebrafish embryos and assessed resulting phenotypes.

RESULTS

PB-treated mouse embryos exhibited neural tube defect and exencephaly, accompanied by apoptosis and reduced cell proliferation in the neural tube, but normal structure in the midface and jaw. PB administration to heterozygous knockout mice of Sf3b4, a gene coding for a SF3b component, influenced the formation of cranial neural crest cells (CNCCs). Despite challenges in continuous PB administration and a high death rate in mice, PB was stably administered to zebrafish embryos, resulting in prolonged survival. Brain, cranial nerve, retina, midface, and jaw development were affected, mimicking spliceosomopathy phenotypes. Additionally, alterations in cell proliferation, cell death, and migration of CNCCs were detected.

CONCLUSIONS

We demonstrated that zebrafish treated with PB exhibited phenotypes similar to those observed in human spliceosomopathies. This experimental system may serve as a valuable research tool for understanding spliceosome function and human diseases.

ADAM interact with large protein complexes to regulate Histone modification, gene expression and splicing

Cranial neural crest (CNC) cells are key stem cells that contribute to most of the facial structures in vertebrates. ADAM ( A D isintegrin A nd M etalloprotease) proteins are essential for the induction and migration of the CNC. We have shown that Adam13 associates with the transcription factor Arid3a to regulate gene expression. Here we show that Adam13 modulates Histone modifications in the CNC. We show that Arid3a binding to the tfap2α promoter depends on the presence of Adam13. This association promotes the expression of one tfap2α variant expressed in the CNC that uniquely activates the expression of gene critical for CNC migration. We show that both Adam13 and human ADAM9 associate with proteins involved in histone modification and RNA splicing, a function critically affected by the loss of Adam13. We propose that ADAMs may act as extracellular sensors to modulate chromatin availability, leading to changes in gene expression and splicing.

Etiology of craniofacial and cardiac malformations in a mouse model of SF3B4-related syndromes

Pathogenic variants in SF3B4, a component of the U2 snRNP complex important for branchpoint sequence recognition and splicing, are responsible for the acrofacial disorders Nager and Rodriguez Syndrome, also known as SF3B4-related syndromes. Patients exhibit malformations in the head, face, limbs, vertebrae as well as the heart. To uncover the etiology of craniofacial malformations found in SF3B4-related syndromes, mutant mouse lines with homozygous deletion of Sf3b4 in neural crest cells (NCC) were generated. Like in human patients, these embryos had craniofacial and cardiac malformations with variable expressivity and penetrance. The severity and survival of Sf3b4 NCC mutants was modified by the level of Sf3b4 in neighboring non-NCC. RNA sequencing analysis of heads of embryos prior to morphological abnormalities revealed significant changes in expression of genes forming the NCC regulatory network, as well as an increase in exon skipping. Additionally, several key histone modifiers involved in craniofacial and cardiac development showed increased exon skipping. Increased exon skipping was also associated with use of a more proximal branch point, as well as an enrichment in thymidine bases in the 50 bp around the branch points. We propose that decrease in Sf3b4 causes changes in the expression and splicing of transcripts required for proper craniofacial and cardiac development, leading to abnormalities.

Human stem cell model of neural crest cell differentiation reveals a requirement of SF3B4 in survival, maintenance, and differentiation

In vitro modeling is a powerful approach to investigate the pathomechanisms driving human congenital conditions. Here we use human embryonic stem cells (hESCs) to model Nager and Rodriguez syndromes, two craniofacial conditions characterized by hypoplastic neural crest-derived craniofacial bones, caused by pathogenic variants of SF3B4, a core component of the spliceosome. We observed that siRNA-mediated knockdown of SF3B4 interferes with the production of hESC-derived neural crest cells, as seen by a marked reduction in neural crest gene expression. This phenotype is associated with an increase in neural crest cell apoptosis and premature neuronal differentiation. Altogether these results point at a role of SF3B4 in neural crest cell survival, maintenance, and differentiation. We propose that the dysregulation of these processes may contribute to Nager/Rodriguez syndrome associated craniofacial defects.

The role of telomerase reverse transcriptase in the mitochondrial protective functions of Angiotensin-(1-7) in diabetic CD34+ cells

Angiotensin (Ang)-(1-7) stimulates vasoprotective functions of diabetic (DB) CD34+ hematopoietic stem/progenitor cells partly by decreasing reactive oxygen species (ROS), increasing nitric oxide (NO) levels and decreasing TGFβ1 secretion. Telomerase reverse transcriptase (TERT) translocates to mitochondria and regulates ROS generation. Alternative splicing of TERT results in variants α-, β- and α-β-TERT, which may oppose functions of full-length (FL) TERT. This study tested if the protective functions of Ang-(1-7) or TGFβ1-silencing are mediated by mitoTERT and that diabetes decreases FL-TERT expression by inducing splicing. CD34+ cells were isolated from the peripheral blood mononuclear cells of nondiabetic (ND, n = 68) or DB (n = 74) subjects. NO and mitoROS levels were evaluated by flow cytometry. TERT splice variants and mitoDNA-lesions were characterized by qPCR. TRAP assay was used for telomerase activity. Decoy peptide was used to block mitochondrial translocation (mitoXTERT). TERT inhibitor or mitoXTERT prevented the effects of Ang-(1-7) on NO or mitoROS levels in DB-CD34+ cells. FL-TERT expression and telomerase activity were lower and mitoDNA-lesions were higher in DB cells compared to ND and were reversed by Ang-(1-7) or TGFβ1-silencing. The prevalence of TERT splice variants, with predominant β-TERT expression, was higher and the expression of FL-TERT was lower in DB cells (n = 25) compared to ND (n = 30). Ang-(1-7) or TGFβ1-silencing decreased TERT-splicing and increased FL-TERT. Blocking of β-splicing increased FL-TERT and protected mitoDNA in DB-cells. The findings suggest that diabetes induces TERT-splicing in CD34+ cells and that β-TERT splice variant largely contributes to the mitoDNA oxidative damage.

JANUS, a spliceosome-associated protein, promotes miRNA biogenesis in Arabidopsis

MicroRNAs (miRNAs) are important regulators of genes expression. Their levels are precisely controlled through modulating the activity of the microprocesser complex (MC). Here, we report that JANUS, a homology of the conserved U2 snRNP assembly factor in yeast and human, is required for miRNA accumulation. JANUS associates with MC components Dicer-like 1 (DCL1) and SERRATE (SE) and directly binds the stem-loop of pri-miRNAs. In a hypomorphic janus mutant, the activity of DCL1, the numbers of MC, and the interaction of primary miRNA transcript (pri-miRNAs) with MC are reduced. These data suggest that JANUS promotes the assembly and activity of MC through its interaction with MC and/or pri-miRNAs. In addition, JANUS modulates the transcription of some pri-miRNAs as it binds the promoter of pri-miRNAs and facilitates Pol II occupancy of at their promoters. Moreover, global splicing defects are detected in janus. Taken together, our study reveals a novel role of a conserved splicing factor in miRNA biogenesis.

Defects of the spliceosomal gene SNRPB affect osteo- and chondro-differentiation

Although gene splicing occurs throughout the body, the phenotype of spliceosomal defects is largely limited to specific tissues. Cerebro-costo-mandibular syndrome (CCMS) is one such spliceosomal disease, which presents as congenital skeletal dysmorphism and is caused by mutations of SNRPB gene encoding Small Nuclear Ribonucleoprotein Polypeptides B/B' (SmB/B'). This study employed in vitro cell cultures to monitor osteo- and chondro-differentiation and examined the role of SmB/B' in the differentiation process. We found that low levels of SmB/B' by knockdown or mutations of SNRPB led to suppressed osteodifferentiation in Saos-2 osteoprogenitor-like cells, which was accompanied by affected splicing of Dlx5. On the other hand, low SmB/B' led to promoted chondrogenesis in HEPM mesenchymal stem cells. Consistent with other reports, osteogenesis was promoted by the Wnt/β-catenin pathway activator and suppressed by Wnt and BMP blockers, whereas chondrogenesis was promoted by Wnt inhibitors. Suppressed osteogenic markers by SNRPB knockdown were partly rescued by Wnt/β-catenin pathway activation. Reporter analysis revealed that suppression of SNRPB results in attenuated Wnt pathway and/or enhanced BMP pathway activities. SNRPB knockdown altered splicing of TCF7L2 which impacts Wnt/β-catenin pathway activities. This work helps unravel the mechanism underlying CCMS whereby reduced expression of spliceosomal proteins causes skeletal phenotypes.

The role of alternative pre-mRNA splicing in cancer progression

Alternative pre-mRNA splicing is a critical mechanism that generates multiple mRNA from a single gene, thereby increasing the diversity of the proteome. Recent research has highlighted the significance of specific splicing isoforms in cellular processes, particularly in regulating cell numbers. In this review, we examine the current understanding of the role of alternative splicing in controlling cancer cell growth and discuss specific splicing factors and isoforms and their molecular mechanisms in cancer progression. These isoforms have been found to intricately control signaling pathways crucial for cell cycle progression, proliferation, and apoptosis. Furthermore, studies have elucidated the characteristics and functional importance of splicing factors that influence cell numbers. Abnormal expression of oncogenic splicing isoforms and splicing factors, as well as disruptions in splicing caused by genetic mutations, have been implicated in the development and progression of tumors. Collectively, these findings provide valuable insights into the complex interplay between alternative splicing and cell proliferation, thereby suggesting the potential of alternative splicing as a therapeutic target for cancer.

De novo PHF5A variants are associated with craniofacial abnormalities, developmental delay, and hypospadias

PURPOSE

The SF3B splicing complex is composed of SF3B1-6 and PHF5A. We report a developmental disorder caused by de novo variants in PHF5A.

METHODS

Clinical, genomic, and functional studies using subject-derived fibroblasts and a heterologous cellular system were performed.

RESULTS

We studied 9 subjects with congenital malformations, including preauricular tags and hypospadias, growth abnormalities, and developmental delay who had de novo heterozygous PHF5A variants, including 4 loss-of-function (LOF), 3 missense, 1 splice, and 1 start-loss variant. In subject-derived fibroblasts with PHF5A LOF variants, wild-type and variant PHF5A mRNAs had a 1:1 ratio, and PHF5A mRNA levels were normal. Transcriptome sequencing revealed alternative promoter use and downregulated genes involved in cell-cycle regulation. Subject and control fibroblasts had similar amounts of PHF5A with the predicted wild-type molecular weight and of SF3B1-3 and SF3B6. SF3B complex formation was unaffected in 2 subject cell lines.

CONCLUSION

Our data suggest the existence of feedback mechanisms in fibroblasts with PHF5A LOF variants to maintain normal levels of SF3B components. These compensatory mechanisms in subject fibroblasts with PHF5A or SF3B4 LOF variants suggest disturbed autoregulation of mutated splicing factor genes in specific cell types, that is, neural crest cells, during embryonic development rather than haploinsufficiency as pathomechanism.

New pathogenic variant in DLX5: New clues for a clinical spectrum from split-hand-foot malformation to fibular aplasia, tibial campomelia and oligosyndactyly

Introduction: FATCO (Fibular Aplasia, Tibial Campomelia and Oligosyndactyly) is a very infrequent skeletal dysplasia classified within the limb hypoplasia-reduction defects group whose genetic cause has not yet been identified. The advent of next-generation sequencing is enabling the diagnosis of diseases with no previously known genetic cause. Methods: We performed a thorough autopsy on a fetus whose pregnancy was legally terminated due to severe malformations detected by ultrasound. A trio exome was run to identify the genetic cause and risk of recurrence. Previous literature of similar cases was systematically searched. Results: Anatomopathological analyses revealed complete fibular aplasia, shortened and campomelic tibia, absent ankle joint, club right foot and a split foot malformation, leading to the diagnosis of FATCO. Exome sequencing showed that the female fetus carried a de novo nonsense variant in DLX5. The literature search permitted the collection of information on 43 patients with FATCO, the majority of whom were males diagnosed postnatally. In most cases, lower limbs were affected exclusively, but in 39.5% of cases the upper limbs were also affected. Conclusion: The pathologies associated with DLX5 variants encompass a wide spectrum of manifestations ranging from abnormalities exclusively in the hands and feet to long bones such as the tibia and fibula.

History and highlights of the teratological collection in the Narrenturm, Vienna (Austria)

The collection of the Narrenturm in Vienna houses and maintains more than 50,000 objects including approximately 1200 teratological specimens; making it one of the biggest collections of specimens from human origin in Europe. The existence of this magnificent collection-representing an important resource for dysmorphology research, mostly awaiting contemporary diagnoses-is not widely known in the scientific community. Here, we show that the Narrenturm harbors a wealth of specimens with (exceptionally) rare congenital anomalies. These museums can be seen as physical repositories of human malformation, covering hundreds of years of dedicated collecting and preserving, thereby creating unique settings that can be used to expand our knowledge of developmental conditions that have to be preserved for future generations of scientists.

The emerging significance of splicing in vertebrate development

Splicing is a crucial regulatory node of gene expression that has been leveraged to expand the proteome from a limited number of genes. Indeed, the vast increase in intron number that accompanied vertebrate emergence might have aided the evolution of developmental and organismal complexity. Here, we review how animal models for core spliceosome components have provided insights into the role of splicing in vertebrate development, with a specific focus on neuronal, neural crest and skeletal development. To this end, we also discuss relevant spliceosomopathies, which are developmental disorders linked to mutations in spliceosome subunits. Finally, we discuss potential mechanisms that could underlie the tissue-specific phenotypes often observed upon spliceosome inhibition and identify gaps in our knowledge that, we hope, will inspire further research.

SF3B4 Frameshift Variants Represented a More Severe Clinical Manifestation in Nager Syndrome

Nager syndrome (NS) is a rare disease marked with craniofacial and preaxial limb anomalies. In this report, we summarized the current evidence to determine a possible genotype-phenotype association among NS individuals. Twenty-four articles comprising of 84 NS (including 9 patients with a severe form of NS [Rodriguez syndrome]) patients were examined, of which 76% were caused by variants in SF3B4 (OMIM *605593, Splicing Factor 3B, Subunit 4). Within the SF3B4 gene, variants located in exon 3 commonly occurred (20%) from a total identified variant, while hotspot location was identified in exon 1 (12%), and primarily occurred as frameshift variants (64%). Thirty-five distinct pathogenic variants within SF3B4 gene were identified with two common sites, c.1A > G and c.1060dupC in exons 1 and 5, respectively. Although no significant genotype-phenotype association was found, it is notable that patients with frameshift SF3B4 variants and predicted to lead to nonsense-mediated RNA decay (NMD) of the transcripts tended to have a more severe clinical manifestation. Additionally, patients harboring variants in exons 2 and 3 displayed a higher proportion of cardiac malformations. Taken together, this article summarizes the pathogenic variants observed in SF3B4 and provides a possible genotype-phenotype relationship in this disease.

Emerging roles of spliceosome in cancer and immunity

Precursor messenger RNA (pre-mRNA) splicing is catalyzed by an intricate ribonucleoprotein complex called the spliceosome. Although the spliceosome is considered to be general cell “housekeeping” machinery, mutations in core components of the spliceosome frequently correlate with cell- or tissue-specific phenotypes and diseases. In this review, we expound the links between spliceosome mutations, aberrant splicing, and human cancers. Remarkably, spliceosome-targeted therapies (STTs) have become efficient anti-cancer strategies for cancer patients with splicing defects. We also highlight the links between spliceosome and immune signaling. Recent studies have shown that some spliceosome gene mutations can result in immune dysregulation and notable phenotypes due to mis-splicing of immune-related genes. Furthermore, several core spliceosome components harbor splicing-independent immune functions within the cell, expanding the functional repertoire of these diverse proteins.

Molecular mechanisms of hearing loss in Nager syndrome

Nager syndrome is a rare human developmental disorder characterized by hypoplastic neural crest-derived craniofacial bones and limb defects. Mutations in SF3B4 gene, which encodes a component of the spliceosome, are a major cause for Nager. A review of the literature indicates that 45% of confirmed cases are also affected by conductive, sensorineural or mixed hearing loss. Conductive hearing loss is due to defective middle ear ossicles, which are neural crest derived, while sensorineural hearing loss typically results from defective inner ear or vestibulocochlear nerve, which are both derived from the otic placode. Animal model of Nager syndrome indicates that upon Sf3b4 knockdown cranial neural crest progenitors are depleted, which may account for the conductive hearing loss in these patients. To determine whether Sf3b4 plays a role in otic placode formation we analyzed the impact of Sf3b4 knockdown on otic development. Sf3b4-depleted Xenopus embryos exhibited reduced expression of several pan-placodal genes six1, dmrta1 and foxi4.1. We confirmed the dependence of placode genes expression on Sf3b4 function in animal cap explants expressing noggin, a BMP antagonist critical to induce placode fate in the ectoderm. Later in development, Sf3b4 morphant embryos had reduced expression of pax8, tbx2, otx2, bmp4 and wnt3a at the otic vesicle stage, and altered otic vesicle development. We propose that in addition to the neural crest, Sf3b4 is required for otic development, which may account for sensorineural hearing loss in Nager syndrome.

Alterations of RNA splicing patterns in esophagus squamous cell carcinoma

Alternative splicing (AS) is an important biological process for regulating the expression of various isoforms from a single gene and thus to promote proteome diversity. In this study, RNA-seq data from 15 pairs of matched esophageal squamous cell carcinoma (ESCC) and normal tissue samples as well as two cell lines were analyzed. AS events with significant differences were identified between ESCC and matched normal tissues, which were re-annotated to find protein coding genes or non-coding RNAs. A total of 45,439 AS events were found. Of these, 6019 (13.25%) significant differentially AS events were identified. Exon skipping (SE) events occupied the largest proportion of abnormal splicing events. Fifteen differential splicing events with the same trends of ΔΨ values in ESCC tissues, as well in the two cell lines were found. Four pathways and 20 biological processes related to pro-metastasis cell junction and migration were significantly enriched for the differentially spliced genes. The upregulated splicing factor SF3B4, which regulates 92 gene splicing events, could be a potential prognostic factor of ESCC. Differentially spliced genes, including HNRNPC, VCL, ZNF207, KIAA1217, TPM1 and CALD1 are shown with a sashimi plot. These results suggest that cell junction- and migration-related biological processes are influenced by AS abnormalities, and aberrant splicing events can be affected by splicing factor expression changes. The involved splicing factor SF3B4 was found to be a survival-related gene in ESCC and is presumed to regulate AS in multiple cancers. In summary, we identified significant differentially expressed AS events which may be related to the development of ESCC.

An Alternative Splicing Program for Mouse Craniofacial Development

Alternative splicing acts as a fundamental mechanism to increase the number of functional transcripts that can be derived from the genome - and its appropriate regulation is required to direct normal development, differentiation, and physiology, in many species. Recent studies have highlighted that mutation of splicing factors, resulting in the disruption of alternative splicing, can have profound consequences for mammalian craniofacial development. However, there has been no systematic analysis of the dynamics of differential splicing during the critical period of face formation with respect to age, tissue layer, or prominence. Here we used deep RNA sequencing to profile transcripts expressed in the developing mouse face for both ectodermal and mesenchymal tissues from the three facial prominences at critical ages for facial development, embryonic days 10.5, 11.5, and 12.5. We also derived separate expression data from the nasal pit relating to the differentiation of the olfactory epithelium for a total of 60 independent datasets. Analysis of these datasets reveals the differential expression of multiple genes, but we find a similar number of genes are regulated only via differential splicing, indicating that alternative splicing is a major source of transcript diversity during facial development. Notably, splicing changes between tissue layers and over time are more prevalent than between prominences, with exon skipping the most common event. We next examined how the variation in splicing correlated with the expression of RNA binding proteins across the various datasets. Further, we assessed how binding sites for splicing regulatory molecules mapped with respect to intron exon boundaries. Overall these studies help define an alternative splicing regulatory program that has important consequences for facial development.

Modelling the developmental spliceosomal craniofacial disorder Burn-McKeown syndrome using induced pluripotent stem cells

The craniofacial developmental disorder Burn-McKeown Syndrome (BMKS) is caused by biallelic variants in the pre-messenger RNA splicing factor gene TXNL4A/DIB1. The majority of affected individuals with BMKS have a 34 base pair deletion in the promoter region of one allele of TXNL4A combined with a loss-of-function variant on the other allele, resulting in reduced TXNL4A expression. However, it is unclear how reduced expression of this ubiquitously expressed spliceosome protein results in craniofacial defects during development. Here we reprogrammed peripheral mononuclear blood cells from a BMKS patient and her unaffected mother into induced pluripotent stem cells (iPSCs) and differentiated the iPSCs into induced neural crest cells (iNCCs), the key cell type required for correct craniofacial development. BMKS patient-derived iPSCs proliferated more slowly than both mother- and unrelated control-derived iPSCs, and RNA-Seq analysis revealed significant differences in gene expression and alternative splicing. Patient iPSCs displayed defective differentiation into iNCCs compared to maternal and unrelated control iPSCs, in particular a delay in undergoing an epithelial-to-mesenchymal transition (EMT). RNA-Seq analysis of differentiated iNCCs revealed widespread gene expression changes and mis-splicing in genes relevant to craniofacial and embryonic development that highlight a dampened response to WNT signalling, the key pathway activated during iNCC differentiation. Furthermore, we identified the mis-splicing of TCF7L2 exon 4, a key gene in the WNT pathway, as a potential cause of the downregulated WNT response in patient cells. Additionally, mis-spliced genes shared common sequence properties such as length, branch point to 3’ splice site (BPS-3’SS) distance and splice site strengths, suggesting that splicing of particular subsets of genes is particularly sensitive to changes in TXNL4A expression. Together, these data provide the first insight into how reduced TXNL4A expression in BMKS patients might compromise splicing and NCC function, resulting in defective craniofacial development in the embryo.

Cancer-Associated Substitutions in RNA Recognition Motifs of PUF60 and U2AF65 Reveal Residues Required for Correct Folding and 3' Splice-Site Selection

U2AF65 (U2AF2) and PUF60 (PUF60) are splicing factors important for recruitment of the U2 small nuclear ribonucleoprotein to lariat branch points and selection of 3′ splice sites (3′ss). Both proteins preferentially bind uridine-rich sequences upstream of 3′ss via their RNA recognition motifs (RRMs). Here, we examined 36 RRM substitutions reported in cancer patients to identify variants that alter 3′ss selection, RNA binding and protein properties. Employing PUF60- and U2AF65-dependent 3′ss previously identified by RNA-seq of depleted cells, we found that 43% (10/23) and 15% (2/13) of independent RRM mutations in U2AF65 and PUF60, respectively, conferred splicing defects. At least three RRM mutations increased skipping of internal U2AF2 (~9%, 2/23) or PUF60 (~8%, 1/13) exons, indicating that cancer-associated RRM mutations can have both cis- and trans-acting effects on splicing. We also report residues required for correct folding/stability of each protein and map functional RRM substitutions on to existing high-resolution structures of U2AF65 and PUF60. These results identify new RRM residues critical for 3′ss selection and provide relatively simple tools to detect clonal RRM mutations that enhance the mRNA isoform diversity.

Spliceosomopathies and neurocristopathies: Two sides of the same coin?

Mutations in core components of the spliceosome are responsible for a group of syndromes collectively known as spliceosomopathies. Patients exhibit microcephaly, micrognathia, malar hypoplasia, external ear anomalies, eye anomalies, psychomotor delay, intellectual disability, limb, and heart defects. Craniofacial malformations in these patients are predominantly found in neural crest cells-derived structures of the face and head. Mutations in eight genes SNRPB, RNU4ATAC, SF3B4, PUF60, EFTUD2, TXNL4, EIF4A3, and CWC27 are associated with craniofacial spliceosomopathies. In this review, we provide a brief description of the normal development of the head and the face and an overview of mutations identified in genes associated with craniofacial spliceosomopathies. We also describe a model to explain how and when these mutations are most likely to impact neural crest cells. We speculate that mutations in a subset of core splicing factors lead to disrupted splicing in neural crest cells because these cells have increased sensitivity to inefficient splicing. Hence, disruption in splicing likely activates a cellular stress response that includes increased skipping of regulatory exons in genes such as MDM2 and MDM4, key regulators of P53. This would result in P53-associated death of neural crest cells and consequently craniofacial malformations associated with spliceosomopathies.

The SF3b complex: splicing and beyond

The SF3b complex is an intrinsic component of the functional U2 small nuclear ribonucleoprotein (snRNP). As U2 snRNP enters nuclear pre-mRNA splicing, SF3b plays key roles in recognizing the branch point sequence (BPS) and facilitating spliceosome assembly and activation. Since the discovery of SF3b, substantial progress has been made in elucidating its molecular mechanism during splicing. In addition, numerous recent studies indicate that SF3b and its components are engaged in various molecular and cellular events that are beyond the canonical role in splicing. This review summarizes the current knowledge on the SF3b complex and highlights its multiple roles in splicing and beyond.

SF3b4: A Versatile Player in Eukaryotic Cells

Spliceosomes are large protein-RNA complexes regulating pre-mRNA processing in eukaryotes. SF3b4 encodes a core subunit of the U2-type spliceosome, loss- or gain-of-function of which often associates with abnormal cell growth, leading to tumorigenesis. Homologs of SF3b4 in other phyla are also essential. In this review, we summarize recent findings on the function of SF3b4. Importantly, we highlight the versatile roles of SF3b4, not only as a component for pre-mRNA splicing, but also as a regulator for transcription, translation, and cell signaling. Recent studies of SF3b4 homologs in different species across evolution will facilitate a better understanding of human diseases caused by the malfunction of SF3b4, such as Nager syndrome (NS) and cancer, in the future.

Heterozygous mutation of the splicing factor Sf3b4 affects development of the axial skeleton and forebrain in mouse

BACKGROUND

Splicing factor 3B subunit 4 (SF3B4) is a causative gene of an acrofacial dysostosis, Nager syndrome. Although in vitro analyses of SF3B4 have proposed multiple noncanonical functions unrelated to splicing, less information is available based on in vivo studies using model animals.

RESULTS

We performed expression and functional analyses of Sf3b4 in mice. The mouse Sf3b4 transcripts were found from two-cell stage, and were ubiquitously present during embryogenesis with high expression levels in several tissues such as forming craniofacial bones and brain. In contrast, expression of a pseudogene-like sequence of mouse Sf3b4 (Sf3b4_ps) found by in silico survey was not detected up to embryonic day 10. We generated a Sf3b4 knockout mouse using CRISPR-Cas9 system. The homozygous mutant mouse of Sf3b4 was embryonic lethal. The heterozygous mutant of Sf3b4 mouse (Sf3b4^+/- ) exhibited smaller body size compared to the wild-type from postnatal to adult period, as well as homeotic posteriorization of the vertebral morphology and flattened calvaria. The flattened calvaria appears to be attributable to mild microcephaly due to a lower cell proliferation rate in the forebrain.

CONCLUSIONS

Our study suggests that Sf3b4 controls anterior-posterior patterning of the axial skeleton and guarantees cell proliferation for forebrain development in mice.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.

Allele-specific binding of RNA-binding proteins reveals functional genetic variants in the RNA

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs.

Developmental processes regulate craniofacial variation in disease and evolution

Variation in development mediates phenotypic differences observed in evolution and disease. Although the mechanisms underlying phenotypic variation are still largely unknown, recent research suggests that variation in developmental processes may play a key role. Developmental processes mediate genotype-phenotype relationships and consequently play an important role regulating phenotypes. In this review, we provide an example of how shared and interacting developmental processes may explain convergence of phenotypes in spliceosomopathies and ribosomopathies. These data also suggest a shared pathway to disease treatment. We then discuss three major mechanisms that contribute to variation in developmental processes: genetic background (gene-gene interactions), gene-environment interactions, and developmental stochasticity. Finally, we comment on evolutionary alterations to developmental processes, and the evolution of disease buffering mechanisms.

Genes uniquely expressed in human growth plate chondrocytes uncover a distinct regulatory network

Background

Chondrogenesis is the earliest stage of skeletal development and is a highly dynamic process, integrating the activities and functions of transcription factors, cell signaling molecules and extracellular matrix proteins. The molecular mechanisms underlying chondrogenesis have been extensively studied and multiple key regulators of this process have been identified. However, a genome-wide overview of the gene regulatory network in chondrogenesis has not been achieved.

Results

In this study, employing RNA sequencing, we identified 332 protein coding genes and 34 long non-coding RNA (lncRNA) genes that are highly selectively expressed in human fetal growth plate chondrocytes. Among the protein coding genes, 32 genes were associated with 62 distinct human skeletal disorders and 153 genes were associated with skeletal defects in knockout mice, confirming their essential roles in skeletal formation. These gene products formed a comprehensive physical interaction network and participated in multiple cellular processes regulating skeletal development. The data also revealed 34 transcription factors and 11,334 distal enhancers that were uniquely active in chondrocytes, functioning as transcriptional regulators for the cartilage-selective genes.

Conclusions

Our findings revealed a complex gene regulatory network controlling skeletal development whereby transcription factors, enhancers and lncRNAs participate in chondrogenesis by transcriptional regulation of key genes. Additionally, the cartilage-selective genes represent candidate genes for unsolved human skeletal disorders.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4378-y) contains supplementary material, which is available to authorized users.

snRNP proteins in health and disease

Split gene architecture of most human genes requires removal of intervening sequences by mRNA splicing that occurs on large multiprotein complexes called spliceosomes. Mutations compromising several spliceosomal components have been recorded in degenerative syndromes and haematological neoplasia, thereby highlighting the importance of accurate splicing execution in homeostasis of assorted adult tissues. Moreover, insufficient splicing underlies defective development of craniofacial skeleton and upper extremities. This review summarizes recent advances in the understanding of splicing factor function deduced from cryo-EM structures. We combine these data with the characterization of splicing factors implicated in hereditary or somatic disorders, with a focus on potential functional consequences the mutations may elicit in spliceosome assembly and/or performance. Given aberrant splicing or perturbations in splicing efficiency substantially underpin disease pathogenesis, profound understanding of the mis-splicing principles may open new therapeutic vistas. In three major sections dedicated to retinal dystrophies, hereditary acrofacial syndromes, and haematological malignancies, we delineate the noticeable variety of conditions associated with dysfunctional splicing and accentuate recurrent patterns in splicing defects.

MED resulting from recessively inherited mutations in the gene encoding calcium-activated nucleotidase CANT1

Multiple Epiphyseal Dysplasia (MED) is a relatively mild skeletal dysplasia characterized by mild short stature, joint pain, and early-onset osteoarthropathy. Dominantly inherited mutations in COMP, MATN3, COL9A1, COL9A2, and COL9A3, and recessively inherited mutations in SLC26A2, account for the molecular basis of disease in about 80-85% of the cases. In two families with recurrent MED of an unknown molecular basis, we used exome sequencing and candidate gene analysis to identify homozygosity for recessively inherited missense mutations in CANT1, which encodes calcium-activated nucleotidase 1. The MED phenotype is thus allelic to the more severe Desbuquois dysplasia phenotype and the results identify CANT1 as a second locus for recessively inherited MED.

Genome-wide association studies identified multiple genetic loci for body size at four growth stages in Chinese Holstein cattle

The growth and maturity of cattle body size affect not only feed efficiency, but also productivity and longevity. Dissecting the genetic architecture of body size is critical for cattle breeding to improve both efficiency and productivity. The volume and weight of body size are indicated by several measurements. Among them, Heart Girth (HG) and Hip Height (HH) are the most important traits. They are widely used as predictors of body weight (BW). Few association studies have been conducted for HG and HH in cattle focusing on single growth stage. In this study, we extended the Genome-wide association studies to a full spectrum of four growth stages (6-, 12-, 18-, and 24-months after birth) in Chinese Holstein heifers. The whole genomic single nucleotide polymorphisms (SNPs) were obtained from the Illumina BovineSNP50 v2 BeadChip genotyped on 3,325 individuals. Estimated breeding values (EBVs) were derived for both HG and HH at the four different ages and analyzed separately for GWAS by using the Fixed and random model Circuitous Probability Unification (FarmCPU) method. In total, 27 SNPs were identified to be significantly associated with HG and HH at different growth stages. We found 66 candidate genes located nearby the associated SNPs, including nine genes that were known as highly related to development and skeletal and muscular growth. In addition, biological function analysis was performed by Ingenuity Pathway Analysis and an interaction network related to development was obtained, which contained 16 genes out of the 66 candidates. The set of putative genes provided valuable resources and can help elucidate the genomic architecture and mechanisms underlying growth traits in dairy cattle.