The extracellular matrix protein TGFBI promotes myofibril bundling and muscle fibre growth in the zebrafish embryo

Identification and immune landscape of sarcopenia-related molecular clusters in inflammatory bowel disease by machine learning and integrated bioinformatics

Sarcopenia, a prevalent comorbidity of inflammatory bowel disease (IBD), is characterized by diminished skeletal muscle mass and strength. Nevertheless, the underlying interconnected mechanisms remain elusive. This study identified distinct expression patterns of sarcopenia-associated genes (SRGs) across individuals with IBD and in samples of normal tissue. By analyzing SRG expression profiles, we effectively segregated 541 IBD samples into three distinct clusters, each marked by its unique immune landscape. To unravel the transcriptional disruptions underlying these clusters, the Weighted Gene Co-expression Network Analysis (WGCNA) algorithm was employed to spotlight key genes linked to each cluster. A diagnostic model based on four key genes (TIMP1, PLAU, PHLDA1, TGFBI) was established using Random Forest and LASSO (least absolute shrinkage and selection operator) algorithms, and validated with the GSE179285 dataset. Moreover, the GSE112366 dataset facilitated the exploration of gene expression dynamics within the ileum mucosa of UC patients pre- and post-Ustekinumab treatment. Additionally, insights into the intricate relationship between immune cells and these pivotal genes were gleaned from the single-cell RNA dataset GSE162335. In conclusion, our findings collectively underscored the pivotal role of sarcopenia-related genes in the pathogenesis of IBD. Their potential as robust biomarkers for future diagnostic and therapeutic strategies is particularly promising, opening avenues for a deeper understanding and improved management of these interconnected conditions.

Quantitative Image Analysis of Axonal Morphology in In Vivo Model

Quantifying axonal branching is crucial for understanding neural circuit function, developmental and regeneration processes and disease mechanisms. Factors that regulate patterns of axonal arborization and tune neuronal circuits are investigated for their implication in various disorders in brain connectivity. The lack of a reliable and user-friendly method makes the quantitative analysis of axon morphology difficult. Specifically, methods to visualize and quantify the complex axon arborization are challenging to implement and apply practically. Our study was aimed at developing a robust but simple method of quantification that used ImageJ 2D analysis and compared it with Imaris visualization and analysis of 3D images. We used zebrafish fluorescent transgenic lines to perform in vivo imaging of developing motor neuron axons that adequately reflected the complexity of axonal networks. Our new method, developed on ImageJ, is easy and fast, giving access to new information such as collateral distribution along the axonal shaft. This study describes step-by-step procedures that can be easily applied to a variety of organisms and in vitro systems. Our study provides a basis for further exploration of neural circuits to gain new insights into neuronal disorders and potential therapeutic interventions.

The myokine Fibcd1 is an endogenous determinant of myofiber size and mitigates cancer-induced myofiber atrophy

Decline in skeletal muscle cell size (myofiber atrophy) is a key feature of cancer-induced wasting (cachexia). In particular, atrophy of the diaphragm, the major muscle responsible for breathing, is an important determinant of cancer-associated mortality. However, therapeutic options are limited. Here, we have used Drosophila transgenic screening to identify muscle-secreted factors (myokines) that act as paracrine regulators of myofiber growth. Subsequent testing in mouse myotubes revealed that mouse Fibcd1 is an evolutionary-conserved myokine that preserves myofiber size via ERK signaling. Local administration of recombinant Fibcd1 (rFibcd1) ameliorates cachexia-induced myofiber atrophy in the diaphragm of mice bearing patient-derived melanoma xenografts and LLC carcinomas. Moreover, rFibcd1 impedes cachexia-associated transcriptional changes in the diaphragm. Fibcd1-induced signaling appears to be muscle selective because rFibcd1 increases ERK activity in myotubes but not in several cancer cell lines tested. We propose that rFibcd1 may help reinstate myofiber size in the diaphragm of patients with cancer cachexia.

Pannexin 1 Regulates Skeletal Muscle Regeneration by Promoting Bleb-Based Myoblast Migration and Fusion Through a Novel Lipid Based Signaling Mechanism

Adult skeletal muscle has robust regenerative capabilities due to the presence of a resident stem cell population called satellite cells. Muscle injury leads to these normally quiescent cells becoming molecularly and metabolically activated and embarking on a program of proliferation, migration, differentiation, and fusion culminating in the repair of damaged tissue. These processes are highly coordinated by paracrine signaling events that drive cytoskeletal rearrangement and cell-cell communication. Pannexins are a family of transmembrane channel proteins that mediate paracrine signaling by ATP release. It is known that Pannexin1 (Panx1) is expressed in skeletal muscle, however, the role of Panx1 during skeletal muscle development and regeneration remains poorly understood. Here we show that Panx1 is expressed on the surface of myoblasts and its expression is rapidly increased upon induction of differentiation and that Panx1^-/- mice exhibit impaired muscle regeneration after injury. Panx1^-/- myoblasts activate the myogenic differentiation program normally, but display marked deficits in migration and fusion. Mechanistically, we show that Panx1 activates P2 class purinergic receptors, which in turn mediate a lipid signaling cascade in myoblasts. This signaling induces bleb-driven amoeboid movement that in turn supports myoblast migration and fusion. Finally, we show that Panx1 is involved in the regulation of cell-matrix interaction through the induction of ADAMTS (Disintegrin-like and Metalloprotease domain with Thrombospondin-type 5) proteins that help remodel the extracellular matrix. These studies reveal a novel role for lipid-based signaling pathways activated by Panx1 in the coordination of myoblast activities essential for skeletal muscle regeneration.

The role of the extracellular matrix protein TGFBI in cancer

The secreted extracellular protein, transforming growth factor beta induced (TGFBI or βIGH3), has roles in regulating numerous biological functions, including cell adhesion and bone formation, both during embryonic development and during the pathogenesis of human disease. TGFBI has been most studied in the context of hereditary corneal dystrophies, where mutations in TGFBI result in accumulation of TGFBI in the cornea. In cancer, early studies focused on TGFBI as a tumor suppressor, in part by promoting chemotherapy sensitivity. However, in established tumors, TGFBI largely has a role in promoting tumor progression, with elevated levels correlating to poorer clinical outcomes. As an important regulator of cancer progression, TGFBI expression and function is tightly regulated by numerous mechanisms including epigenetic silencing through promoter methylation and microRNAs. Mechanisms to target TGFBI have potential clinical utility in treating advanced cancers, while assessing TGFBI levels could be a biomarker for chemotherapy resistance and tumor progression.

TGFBI modulates tumour hypoxia and promotes breast cancer metastasis

Breast cancer metastasis is a complex process that depends not only on intrinsic characteristics of metastatic stem cells, but also on the particular microenvironment that supports their growth and modulates the plasticity of the system. In search for microenvironmental factors supporting cancer stem cell (CSC) growth and tumour progression to metastasis, we here investigated the role of the matricellular protein transforming growth factor beta induced (TGFBI) in breast cancer. We crossed the MMTV‐PyMT model of mammary gland tumorigenesis with a Tgfbi^Δ/Δ mouse and studied the CSC content of the tumours. We performed RNAseq on wt and ko tumours, and analysed the tumour vasculature and the immune compartment by IHC and FACS. The source of TGFBI expression was determined by qPCR and by bone marrow transplantation experiments. Finally, we performed in silico analyses using the METABRIC cohort to assess the potential prognostic value of TGFBI. We observed that deletion of Tgfbi led to a dramatic decrease in CSC content and lung metastasis. Our results show that lack of TGFBI resulted in tumour vessel normalisation, with improved vessel perfusion and decreased hypoxia, a major factor controlling CSCs and metastasis. Furthermore, human data mining in a cohort of breast cancer patients showed that higher expression of TGFBI correlates with poor prognosis and is associated with the more aggressive subtypes of breast cancer. Overall, these data reveal a novel biological mechanism controlling metastasis that could potentially be exploited to improve the efficacy and delivery of chemotherapeutic agents in breast cancer.

Transformation of the Transcriptomic Profile of Mouse Periocular Mesenchyme During Formation of the Embryonic Cornea

Purpose

Defects in neural crest development are a major contributing factor in corneal dysgenesis, but little is known about the genetic landscape during corneal development. The purpose of this study was to provide a detailed transcriptome profile and evaluate changes in gene expression during mouse corneal development.

Methods

RNA sequencing was used to uncover the transcriptomic profile of periocular mesenchyme (pNC) isolated at embryonic day (E) 10.5 and corneas isolated at E14.5 and E16.5. The spatiotemporal expression of several differentially expressed genes was validated by in situ hybridization.

Results

Analysis of the whole-transcriptome profile between pNC and embryonic corneas identified 3815 unique differentially expressed genes. Pathway analysis revealed an enrichment of differentially expressed genes involved in signal transduction (retinoic acid, transforming growth factor-β, and Wnt pathways) and transcriptional regulation.

Conclusions

Our analyses, for the first time, identify a large number of differentially expressed genes during progressive stages of mouse corneal development. Our data provide a comprehensive transcriptomic profile of the developing cornea. Combined, these data serve as a valuable resource for the identification of novel regulatory networks crucial for the advancement of studies in congenital defects, stem cell therapy, bioengineering, and adult corneal diseases.

RNA-Seq Analysis Reveals Spatial and Sex Differences in Pectoralis Major Muscle of Broiler Chickens Contributing to Difference in Susceptibility to Wooden Breast Disease

Wooden Breast Disease (WBD) is a novel myopathy affecting the pectoralis major muscle of modern broiler chickens. The etiology of WBD is not currently known, but has been linked to increased feed efficiency, growth rate, and muscle yield in broiler chickens. Differential effect of WBD has been detected between regions of the P. major and between sexes of broilers-male birds and the cranial aspect of the muscle tend to be more severely affected by the disease than females and the caudal aspect. This study aimed to characterize biological differences in the P. major between regions of the muscle and sexes of birds. Samples were taken from the cranial and caudal aspects of P. major muscles of 3-week-old, unaffected male and female birds for RNA sequencing. RNA was extracted and used for preparation of cDNA libraries, which were sequenced by the Delaware Biotechnology Institute (DBI) using HiSeq2500. Sequence reads were aligned to the chicken reference genome with HISAT, and genes were analyzed for differential expression between regions of the breast muscle and sexes of birds using CuffDiff. Functional analysis was performed on differentially expressed genes (DEGs) between sex groups using DAVID and Ingenuity Pathway Analysis (IPA). There were 12 DEGs between cranial and caudal samples, and 260 between male and female birds. Out of the 260 genes differentially expressed between sexes, 189 were upregulated in males. Of this subset, 103 genes (55%) were located on the Z-chromosome. There was increased expression of genes involved in fat metabolism and oxidative stress responses in the cranial region of the P. major muscle, as well as increased expression of fat metabolism, oxidative stress response, antiangiogenesis, and connective tissue proliferation genes in male broilers. These results support the hypothesis that there are biological characteristics in male broilers and the cranial region of the breast muscle that may make them more susceptible to WBD, as well as raising the possibility of a metabolic switch in modern broiler chickens that may be more prominent in males.

Subtle morphometric, behavioral and gene expression effects in larval zebrafish exposed to PFHxA, PFHxS and 6:2 FTOH

Recent studies of perfluoroalkylated substances (PFASs) have focused on the toxicity of long chain PFASs, such as PFOS or PFOA, which have been demonstrated to cause an array of developmental and behavioral effects. However, less is known about low molecular weight PFASs and alternatives. This study examined the morphometric and behavioral effects in zebrafish following developmental exposures of C6 PFASs: perfluorohexanoic acid, PFHxA, perfluorohexane sulfonate, PFHxS, and 6:2 fluorotelomer alcohol, 6:2 FTOH. Embryos were exposed to 0.02-20 μM concentrations of these compounds from the high stage (˜3 h post fertilization, hpf) until 120 hpf. Morphometric and gene expression endpoints were examined at 120 hpf. Genes selected for analysis were previously shown to be altered in zebrafish developmentally exposed to PFOS and PFOA. Additionally, exposed larvae were transferred to clean water and reared until 14 days post fertilization, dpf, when behavioral assays were completed and morphometric endpoints examined. While PFHxA was found to be the most acutely toxic at 120 hpf, few morphometric effects were observed. Gene expression was the most sensitive endpoint with significant increased tgfb1a, bdnf, and ap1s1 expression observed with PFHxA exposure. PFHxS exposure produced morphometric effects in the larvae, specifically increased length and yolk sac area at 2 and 20 μM. This phenotype persisted to the 14 dpf time point, where these larvae additionally displayed decreased distance traveled and crosses through the center of the arena of the behavioral assay. Exposure to 6:2 FTOH caused no morphometric effects at 120 hpf, and this compound was the least acutely toxic. However, expression of both tgfb1a and bdnf were increased by greater than 2 fold change at this time point. Effects also persisted to 14 dpf where a significant increase in distance traveled and velocity were observed in the behavioral assay. This study demonstrates effects on behavioral, morphometric and gene expression endpoints with developmental PFHxA, PFHxS, and 6:2 FTOH exposures in zebrafish.

Paxillin genes and actomyosin contractility regulate myotome morphogenesis in zebrafish

Paxillin (Pxn) is a key adapter protein and signaling regulator at sites of cell-extracellular matrix (ECM) adhesion. Here, we investigated the role of Pxn during vertebrate development using the zebrafish embryo as a model system. We have characterized two Pxn genes, pxna and pxnb, in zebrafish that are maternally supplied and expressed in multiple tissues. Gene editing and antisense gene knockdown approaches were used to uncover Pxn functions during zebrafish development. While mutation of either pxna or pxnb alone did not cause gross embryonic phenotypes, double mutants lacking maternally supplied pxna or pxnb displayed defects in cardiovascular, axial, and skeletal muscle development. Transient knockdown of Pxn proteins resulted in similar defects. Irregular myotome shape and ECM composition were observed, suggesting an "inside-out" signaling role for Paxillin genes in the development of myotendinous junctions. Inhibiting non-muscle Myosin-II during somitogenesis altered the subcellular localization of Pxn protein and phenocopied pxn gene loss-of-function. This indicates that Paxillin genes are effectors of actomyosin contractility-driven morphogenesis of trunk musculature in zebrafish. Together, these results reveal new functions for Pxn during muscle development and provide novel genetic models to elucidate Pxn functions.

A Zebrafish Model for a Human Myopathy Associated with Mutation of the Unconventional Myosin MYO18B

Myosin 18B is an unconventional myosin that has been implicated in tumor progression in humans. In addition, loss-of-function mutations of the MYO18B gene have recently been identified in several patients exhibiting symptoms of nemaline myopathy. In mouse, mutation of Myo18B results in early developmental arrest associated with cardiomyopathy, precluding analysis of its effects on skeletal muscle development. The zebrafish, frozen (fro) mutant was identified as one of a group of immotile mutants in the 1996 Tübingen genetic screen. Mutant embryos display a loss of birefringency in their skeletal muscle, indicative of disrupted sarcomeric organization. Using meiotic mapping, we localized the fro locus to the previously unannotated zebrafish myo18b gene, the product of which shares close to 50% identity with its human ortholog. Transcription of myo18b is restricted to fast-twitch myocytes in the zebrafish embryo; consistent with this, fro mutant embryos exhibit defects specifically in their fast-twitch skeletal muscles. We show that sarcomeric assembly is blocked at an early stage in fro mutants, leading to the disorganized accumulation of actin, myosin, and α-actinin and a complete loss of myofibrillar organization in fast-twitch muscles.

Immobilization of Dystrophin and Laminin α2-Chain Deficient Zebrafish Larvae In Vivo Prevents the Development of Muscular Dystrophy

Muscular dystrophies are often caused by genetic alterations in the dystrophin-dystroglycan complex or its extracellular ligands. These structures are associated with the cell membrane and provide mechanical links between the cytoskeleton and the matrix. Mechanical stress is considered a pathological mechanism and muscle immobilization has been shown to be beneficial in some mouse models of muscular dystrophy. The zebrafish enables novel and less complex models to examine the effects of extended immobilization or muscle relaxation in vivo in different dystrophy models. We have examined effects of immobilization in larvae from two zebrafish strains with muscular dystrophy, the Sapje dystrophin-deficient and the Candyfloss laminin α2-chain-deficient strains. Larvae (4 days post fertilization, dpf) of both mutants have significantly lower active force in vitro, alterations in the muscle structure with gaps between muscle fibers and altered birefringence patterns compared to their normal siblings. Complete immobilization (18 hrs to 4 dpf) was achieved using a small molecular inhibitor of actin-myosin interaction (BTS, 50 μM). This treatment resulted in a significantly weaker active contraction at 4 dpf in both mutated larvae and normal siblings, most likely reflecting a general effect of immobilization on myofibrillogenesis. The immobilization also significantly reduced the structural damage in the mutated strains, showing that muscle activity is an important pathological mechanism. Following one-day washout of BTS, muscle tension partly recovered in the Candyfloss siblings and caused structural damage in these mutants, indicating activity-induced muscle recovery and damage, respectively.

Periostin and TGF-β-induced protein: Two peas in a pod?

Periostin (PN) and TGF-β-induced protein (βig-h3) are paralogs that contain a single emilin and four fasciclin-1 modules and are secreted from cells. PN receives attention because of its up-regulation in cancer and degenerative and allergic diseases. βig-h3 is highly enriched in cornea and best known for harboring mutations in humans associated with corneal dystrophies. Both proteins are expressed widely, and many functions, some over-lapping, have been attributed to PN and βig-h3 based on biochemical, cell culture, and whole animal experiments. We attempt to organize this knowledge so as to facilitate research on these interesting and incompletely understood proteins. We focus particularly on whether PN and βig-h3 are modified by vitamin K-dependent γ-glutamyl carboxylation, a question of considerable importance given the profound effects of γ-carboxylation on structure and function of other proteins. We consider the roles of PN and βig-h3 in formation of extracellular matrix and as ligands for integrin receptors. We attempt to reconcile the contradictory results that have arisen concerning the role of PN, which has emerged as a marker of TH2 immunity, in murine models of allergic asthma. Finally, when possible we compare and contrast the structures and functions of the two proteins.

Membrane-myofibril cross-talk in myofibrillogenesis and in muscular dystrophy pathogenesis: lessons from the zebrafish

Striated muscle has a highly ordered structure in which specialized domains of the cell membrane involved in force transmission (costameres) and excitation-contraction coupling (T tubules) as well as the internal membranes of the sarcoplasmic reticulum are organized over specific regions of the sarcomere. Optimal muscle function is dependent on this high level of organization but how it established and maintained is not well understood. Due to its ex utero development and transparency, the zebrafish embryo is an excellent vertebrate model for the study of dynamic relationships both within and between cells during development. Transgenic models have allowed the delineation of cellular migration and complex morphogenic rearrangements during the differentiation of skeletal myocytes and the assembly and organization of new myofibrils. Molecular targeting of genes and transcripts has allowed the identification of key requirements for myofibril assembly and organization. With the recent advances in gene editing approaches, the zebrafish will become an increasingly important model for the study of human myopathies and muscular dystrophies. Its high fecundity and small size make it well suited to high-throughput screenings to identify novel pharmacologic and molecular therapies for the treatment of a broad range of neuromuscular conditions. In this review, we examine the lessons learned from the zebrafish model regarding the complex interactions between the sarcomere and the sarcolemma that pattern the developing myocyte and discuss the potential for zebrafish as a model system to examine the pathophysiology of, and identify new treatments for, human myopathies and muscular dystrophies.

Force measurement during contraction to assess muscle function in zebrafish larvae

Zebrafish larvae provide models of muscle development, muscle disease and muscle-related chemical toxicity, but related studies often lack functional measures of muscle health. In this video article, we demonstrate a method to measure force generation during contraction of zebrafish larval trunk muscle. Force measurements are accomplished by placing an anesthetized larva into a chamber filled with a salt solution. The anterior end of the larva is tied to a force transducer and the posterior end of the larva is tied to a length controller. An isometric twitch contraction is elicited by electric field stimulation and the force response is recorded for analysis. Force generation during contraction provides a measure of overall muscle health and specifically provides a measure of muscle function. Although we describe this technique for use with wild-type larvae, this method can be used with genetically modified larvae or with larvae treated with drugs or toxicants, to characterize muscle disease models and evaluate treatments, or to study muscle development, injury, or chemical toxicity.

Whole Genome Association Studies of Residual Feed Intake and Related Traits in the Pig

BACKGROUND

Residual feed intake (RFI), a measure of feed efficiency, is the difference between observed feed intake and the expected feed requirement predicted from growth and maintenance. Pigs with low RFI have reduced feed costs without compromising their growth. Identification of genes or genetic markers associated with RFI will be useful for marker-assisted selection at an early age of animals with improved feed efficiency.

METHODOLOGY/PRINCIPAL FINDINGS

Whole genome association studies (WGAS) for RFI, average daily feed intake (ADFI), average daily gain (ADG), back fat (BF) and loin muscle area (LMA) were performed on 1,400 pigs from the divergently selected ISU-RFI lines, using the Illumina PorcineSNP60 BeadChip. Various statistical methods were applied to find SNPs and genomic regions associated with the traits, including a Bayesian approach using GenSel software, and frequentist approaches such as allele frequency differences between lines, single SNP and haplotype analyses using PLINK software. Single SNP and haplotype analyses showed no significant associations (except for LMA) after genomic control and FDR. Bayesian analyses found at least 2 associations for each trait at a false positive probability of 0.5. At generation 8, the RFI selection lines mainly differed in allele frequencies for SNPs near (<0.05 Mb) genes that regulate insulin release and leptin functions. The Bayesian approach identified associations of genomic regions containing insulin release genes (e.g., GLP1R, CDKAL, SGMS1) with RFI and ADFI, of regions with energy homeostasis (e.g., MC4R, PGM1, GPR81) and muscle growth related genes (e.g., TGFB1) with ADG, and of fat metabolism genes (e.g., ACOXL, AEBP1) with BF. Specifically, a very highly significantly associated QTL for LMA on SSC7 with skeletal myogenesis genes (e.g., KLHL31) was identified for subsequent fine mapping.

CONCLUSIONS/SIGNIFICANCE

Important genomic regions associated with RFI related traits were identified for future validation studies prior to their incorporation in marker-assisted selection programs.

βig-h3 potentiates the profibrogenic effect of TGFβ signaling on connective tissue progenitor cells through the negative regulation of master chondrogenic genes

Tendons and cartilage are specialized forms of connective tissues originated from common progenitor cells. Initial stages of differentiation of these tissues are characterized by the formation of cell aggregates, which share many molecular markers. Once differentiated, these cells retain considerable plasticity, and chondral metaplasia of tendon and fibrous connective tissues and eventual ossification often accompany degenerative diseases in the adult musculoskeletal system. While this fact is of great relevance for regenerative medicine and aging biology, its molecular basis remains to be elucidated. Gene expression analysis in several physiological and experimental paradigms suggests that differentiation of tendon and cartilage is regulated by a balance in the expression of chondrogenic versus tenogenic genes in the connective tissue cell precursors. Transforming growth factor β (TGFβ) may function both as a profibrogenic or as a prochondrogenic factor for embryonic limb mesoderm and mesenchymal stem cell cultures, but mice that are null for TGFβ 2 and 3 lack tendons. Here, we identify βig-h3 as a factor downstream TGFβ signaling regulated by Smad 2 and 3, which is highly expressed in the differentiating tendons and joint capsules. Furthermore, gain- and loss-of-function experiments using limb mesoderm micromass cultures show that βig-h3 downregulates the expression of cartilage master genes, including Sox9, type II collagen, and Hif-1α. Positive regulation of Sox9 and type II Collagen observed in micromass cultures grown under hypoxic conditions is prevented by exogenous administration of βIG-H3, and the antichondrogenic influence of βIG-H3 is lost after Hif-1α silencing with shRNA. Collectively, our findings indicate that βig-h3 promotes the fibrogenic influence of TGFβ signaling, neutralizing the prochondrogenic influence of the hypoxic-inducible factor 1 activated by the hypoxic microenvironment characteristic of limb mesenchymal aggregates.

Obscurin depletion impairs organization of skeletal muscle in developing zebrafish embryos

During development, skeletal myoblasts differentiate into myocytes and skeletal myotubes with mature contractile structures that are precisely oriented with respect to surrounding cells and tissues. Establishment of this highly ordered structure requires reciprocal interactions between the differentiating myocytes and the surrounding extracellular matrix to form correctly positioned and well-organized attachments from the skeletal muscle to the bony skeleton. Using the developing zebrafish embryo as a model, we examined the relationship between new myofibril assembly and the organization of the membrane domains involved in cell-extracellular matrix interactions. We determined that depletion of obscurin, a giant muscle protein, resulted in irregular cell morphology and disturbed extracellular matrix organization during skeletal muscle development. The resulting impairment of myocyte organization was associated with disturbance of the internal architecture of the myocyte suggesting that obscurin participates in organizing the internal structure of the myocyte and translating those structural cues to surrounding cells and tissues.

Rbfox-regulated alternative splicing is critical for zebrafish cardiac and skeletal muscle functions

Rbfox RNA binding proteins are implicated as regulators of phylogenetically-conserved alternative splicing events important for muscle function. To investigate the function of rbfox genes, we used morpholino-mediated knockdown of muscle-expressed rbfox1l and rbfox2 in zebrafish embryos. Single and double morphant embryos exhibited changes in splicing of overlapping sets of bioinformatically-predicted rbfox target exons, many of which exhibit a muscle-enriched splicing pattern that is conserved in vertebrates. Thus, conservation of intronic Rbfox binding motifs is a good predictor of Rbfox-regulated alternative splicing. Morphology and development of single morphant embryos were strikingly normal; however, muscle development in double morphants was severely disrupted. Defects in cardiac muscle were marked by reduced heart rate and in skeletal muscle by complete paralysis. The predominance of wavy myofibers and abnormal thick and thin filaments in skeletal muscle revealed that myofibril assembly is defective and disorganized in double morphants. Ultra-structural analysis revealed that although sarcomeres with electron dense M- and Z-bands are present in muscle fibers of rbfox1l/rbox2 morphants, they are substantially reduced in number and alignment. Importantly, splicing changes and morphological defects were rescued by expression of morpholino-resistant rbfox cDNA. Additionally, a target-blocking MO complementary to a single UGCAUG motif adjacent to an rbfox target exon of fxr1 inhibited inclusion in a similar manner to rbfox knockdown, providing evidence that Rbfox regulates the splicing of target exons via direct binding to intronic regulatory motifs. We conclude that Rbfox proteins regulate an alternative splicing program essential for vertebrate heart and skeletal muscle functions.