Alternative Splice Site Research Articles

Cacna1b alternative splicing impacts excitatory neurotransmission and is linked to behavioral responses to aversive stimuli

Presynaptic CaV2.2 channels control calcium entry that triggers neurotransmitter release at both central and peripheral synapses. The Cacna1b gene encodes the α1-pore forming subunit of CaV2.2 channels. Distinct subsets of splice variants of CaV2.2 derived from cell-specific alternative splicing of the Cacna1b pre-mRNA are expressed in specific subpopulations of neurons. Four cell-specific sites of alternative splicing in Cacna1b that alter CaV2.2 channel function have been described in detail: three cassette exons (e18a, e24a, and e31a) and a pair of mutually exclusive exons (e37a/e37b). Cacna1b mRNAs containing e37a are highly enriched in a subpopulation of nociceptors where they influence nociception and morphine analgesia. E37a-Cacna1b mRNAs are also expressed in brain, but their cell-specific expression in this part of the nervous system, their functional consequences in central synapses and their role on complex behavior have not been studied. In this report, we show that e37a-Cacna1b mRNAs are expressed in excitatory projection neurons where CaV2.2 channels are known to influence transmitter release at excitatory inputs from entorhinal cortex (EC) to dentate gyrus (DG). By comparing behaviors of WT mice to those that only express e37b-CaV2.2 channels, we found evidence that e37a-CaV2.2 enhances behavioral responses to aversive stimuli. Our results suggest that alternative splicing of Cacna1b e37a influences excitatory transmitter release and couples to complex behaviors.

Clinical Challenge and Call for Research on Keloid Disorder: Meeting Report from The 3rd International Keloid Research Foundation Symposium, Beijing 2019

Clinical Challenge and Call for Research on Keloid Disorder: Meeting Report from The 3rd International Keloid Research Foundation Symposium, Beijing 2019

Lnc9141-a and -b Play a Different Role in Bovine Myoblast Proliferation, Apoptosis, and Differentiation

Previously, our transcriptome sequencing revealed that lnc9141 was differentially expressed in muscles of fetal bovine, calf, and adult bovine, which is considered to provide the basis for raising the muscle mass. In this study, we identified lnc9141 characters. lnc9141 has different transcription start sites and 3′ alternative splicing sites of exon 1, producing lnc9141-a and lnc9141-b transcripts that were highly expressed in the heart and lung. Moreover, neither lnc9141-a nor lnc9141-b had the ability to encode proteins. The functions of lnc9141-a and lnc9141-b were explored by cell cycle, 5-ethynyl-2'-deoxyuridine (EdU), and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The results showed that lnc9141-a or lnc9141-b overexpression decreased the number of myoblasts in the S phase and increased the proportion of cells in the G0/G1 phase. Furthermore, overexpressing lnc9141-a and lnc9141-b respectively downregulated the expression of Cyclin D1. However, lnc9141-a or lnc9141-b interference was found to increase the number of S-phase myoblasts, and upregulate Cyclin D1 and Cyclin E expression. Through Annexin V-FITC/propidium iodide (PI) double staining and the expression of apoptosis marker genes (Bax, Bcl2, and Caspase-3), it was found that lnc9141-b could regulate the expression of Bax gene. Meantime, high expression of lnc9141-b could decrease MyHC expression. In addition, the intergenic region between lnc9141 and IRX5 was 2.3 kb, with a head-to-head orientation. The study also revealed the core regions of the lnc9141 and IRX5 promoter. Our study demonstrated that both lnc9141-a and -b expression inhibited bovine myoblast proliferation. However, lnc9141-b regulated Bax and MyHC expression. The regulatory mechanism of lnc9141-a and lnc9141-b needs to be further explored.

The Oncogene Metadherin Interacts with the Known Splicing Proteins YTHDC1, Sam68 and T-STAR and Plays a Novel Role in Alternative mRNA Splicing.

Oncogenic metadherin is a key contributor to tumourigenesis with metadherin expression and cytoplasmic localisation previously linked to poor survival. A number of reports have shown metadherin localises specifically to nuclear speckles known to be rich in RNA-binding proteins including the splicing proteins YTHDC1, Sam68 and T-STAR, that have been shown to select alternative splice sites in mRNA of tumour-associated proteins including BRCA, MDM2 and VEGF. Here we investigate the interaction and relationship between metadherin and the splice factors YTHDC1, T-STAR and Sam68. Using a yeast two-hybrid assay and immunoprecipitation we show that metadherin interacts with YTHDC1, Sam68 and T-STAR and demonstrate that T-STAR is significantly overexpressed in prostate cancer tissue compared to benign prostate tissue. We also demonstrate that metadherin influences splice site selection in a dose-dependent manner in CD44v5-luc minigene reporter assays. Finally, we demonstrate that prostate cancer patients with higher metadherin expression have greater expression of the CD44v5 exon. CD44v5 expression could be used to discriminate patients with poor outcomes following radical prostatectomy. In this work we show for the first time that metadherin interacts with, and modulates, the function of key components of splicing associated with cancer development and progression.

Abstract 772: The Highly Prevalent 25bp Intronic Deletion in MYBPC3 is Benign Under Baseline Conditions

Background: Hypertrophic cardiomyopathy (HCM) affects at least 1 in 500 people worldwide, and results in the thickening of the ventricular walls and reduced cardiac function. Mutations in MYBPC3 , encoding cardiac myosin binding protein-C, are the most common cause of HCM. Previously, a highly prevalent 25bp deletion within intron 32 of MYBPC3 was described in the South Asian population. The MYBPC3 d25bp variant is present in approximately 100 million people, and encompasses a splicing branch point predicted to result in abnormal splicing of exon 33. Thus, there is a critical need to understand the mechanism by which MYBPC3 d25bp may cause cardiomyopathy. Methods: To determine the role of the 25bp deletion in vivo , knock-in humanized mice were created in which intron 32 was replaced with the human intron 32, with or without the MYBPC3 d25bp mutation. Mice were characterized at 3- and 6-months of age by echocardiography, histological, and protein analysis. The presence of aberrant exon splicing was also determined in mice carrying the MYBPC3 d25bp variant through RT-PCR and mini-gene assays. Finally, exon trapping experiments were performed to understand the mechanism behind exon skipping. Results: Under baseline conditions, MYBPC3 d25bp displayed no changes in cardiac function or morphology as measured by echocardiography (FS (%): NTG 35.3%, WT 32.8%, Het 33.7%), heart weight to body weight ratio, or histology. While exon 33 skipping was not detected by RT-PCR, the presence of an alternative splice site within exon 33 was identified in MYBPC3 d25bp mice. However, this did not affect the protein levels of cMyBP-C. Furthermore, mini-gene experiments demonstrated that the MYBPC3 d25bp mutation significantly reduced the percentage of correctly spliced transcripts (86.2% vs. 77.5%). Conclusions: These data demonstrate that the presence of the highly prevalent 25bp deletion is not sufficient to cause disease under baseline conditions. However, it is possible that the increased levels of aberrant splicing may increase the risk for developing HCM.

Deep Splicing Code: Classifying Alternative Splicing Events Using Deep Learning

Alternative splicing (AS) is the process of combining different parts of the pre-mRNA to produce diverse transcripts and eventually different protein products from a single gene. In computational biology field, researchers try to understand AS behavior and regulation using computational models known as “Splicing Codes”. The final goal of these algorithms is to make an in-silico prediction of AS outcome from genomic sequence. Here, we develop a deep learning approach, called Deep Splicing Code (DSC), for categorizing the well-studied classes of AS namely alternatively skipped exons, alternative 5’ss, alternative 3’ss, and constitutively spliced exons based only on the sequence of the exon junctions. The proposed approach significantly improves the prediction and the obtained results reveal that constitutive exons have distinguishable local characteristics from alternatively spliced exons. Using the motif visualization technique, we show that the trained models learned to search for competitive alternative splice sites as well as motifs of important splicing factors with high precision. Thus, the proposed approach greatly expands the opportunities to improve alternative splicing modeling. In addition, a web-server for AS events prediction has been developed based on the proposed method.

Epigenetic marks on the chromatin are associated with RNA splicing in human leukemia cells

In this work we estimated associations between distribution patterns of several epigenetic marks and splicing events on the level of full genome and transcriptome in the cells of two leukemic cell lines containing two different reciprocal chromosome translocations. Significant difference in distribution of epigenetic marks was found, contributing to more opened or more closed chromatin in loci of donor vs acceptor and canonical vs alternative splice sites in expressing genes. Marks of the opened chromatin are significantly more often present in the genomic regions with alternative splicing events than in regions with canonical splicing, while for the mark of the histone 3 trimethylation at lysine 36, the opposite trend is observed. The obtained results reveal the presence of an additional, still very poorly studied layer in the regulation of alternative splicing in human cells.

Improving Silkworm Genome Annotation Using a Proteogenomics Approach.

The silkworm genome has been deeply sequenced and assembled, but accurate genome annotation, which is important for modern biological research, remains far from complete. To improve silkworm genome annotation, we carried out a proteogenomics analysis using 9.8 million mass spectra collected from different tissues and developmental stages of the silkworm. The results confirmed the translational products of 4307 existing gene models and identified 1701 novel genome search-specific peptides (GSSPs). Using these GSSPs, 74 novel gene-coding sequences were identified, and 121 existing gene models were corrected. We also identified 1182 novel junction peptides based on an exon-skipping database that resulted in the identification of 973 alternative splicing sites. Furthermore, we performed RNA-seq analysis to improve silkworm genome annotation at the transcriptional level. A total of 1704 new transcripts and 1136 new exons were identified, 2581 untranslated regions (UTRs) were revised, and 1301 alternative splicing (AS) genes were identified. The transcriptomics results were integrated with the proteomics data to further complement and verify the new annotations. In addition, 14 incorrect genes and 10 skipped exons were verified using the two analysis methods. Altogether, we identified 1838 new transcripts and 1593 AS genes and revised 5074 existing genes using proteogenomics and transcriptome analyses. Data are available via ProteomeXchange with identifier PXD009672. The large-scale proteogenomics and transcriptome analyses in this study will greatly improve silkworm genome annotation and contribute to future studies.

Evolution of Alternative Splicing in Eudicots.

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon–exon junctions were the major determinants for both alternative 5′ donor site and 3′ acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Factors affecting HLA expression: A review.

The detection and semiquantitative measurement of circulating human leucocyte antigen (HLA)-specific antibodies is essential for the management of patients before and after transplantation. In addition, the pretransplant cross-match to assess the reactivity of recipient HLA antibody against donor lymphocytes has long been the gold standard to prevent hyperacute rejection. Whilst both of these tests assume that recipient HLA-specific antibody is the only variable in the assessment of transplant risk, this is not the case. Transplant immunologists recognize that some HLA antigens are expressed at levels a magnitude lower than others (e.g., HLA-C, HLA-DQ), but within loci, and between different cell types there are many factors that influence HLA expression in both resting and activated cells. HLA is not usually expressed without the specific promoter proteins NLRC5, for HLA class I, and CIITA, for class II. The quantity of HLA protein production is then affected by factors including promoter region polymorphisms, alternative exon splice sites, methylation and microRNA-directed degradation. Different loci are influenced by multiple combinations of these control mechanisms making prediction of HLA regulation difficult, but an ability to measure the cellular expression of each HLA antigen, in conjunction with knowledge of circulating HLA-specific antibody, would lead to a more informed algorithm to assess transplant risk.

The arginine and serine-rich domains of Acinus modulate splicing.

Apoptotic chromatin condensation inducer in the nucleus (Acinus) is an RNA-binding protein that has a functional role in inducing apoptotic chromatin condensation and regulating messenger RNA (mRNA) processing. Acinus interacts with the spliceosomal machinery and is a member of the ASAP (apoptosis and splicing-associated protein complex) as well as the EJC (exon junction complex), which gets deposited onto mRNA during splicing. In this study, we have used in vivo splicing assays to characterize the function of Acinus in pre-mRNA splicing more closely. We show that full-length Acinus-S', an isoform of Acinus, does not have a role in modulating splice site selection in human immunodeficiency virus 1 minigene reporter system. In contrast, we observed that the tethering of arginine/serine (RS) and RNPS1-SAP18-binding (RSB) domains of Acinus could regulate the selection of alternative splice sites, thereby revealing the potential of Acinus in stimulating alternative splicing. Altogether, our data suggest that the RS and RSB domains play a critical role in regulating splicing activity via selection of distinct splice sites during pre-mRNA splicing.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents.

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B.carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.

The Role of Transposable Elements in the Differentiation of Stem Cells

The activation of mobile genetic elements is a prerequisite for controlling the expression of genes in successive cell divisions, with their change specifically depending on the tissue and developmental stage. These activation patterns are characteristic of totipotent and pluripotent embryo cells, as well as for stem cells in the postnatal period. We have suggested that in the evolution of multicellular eukaryotes, optimal ratios are selected in the arrangement of transposons relative to exons and introns of host genes. These ratios, specific to each species, can be the basis for controlling the sequential differentiation of stem cells necessary for the development of the whole organism. Despite the fact that cell differentiation in ontogenesis is controlled by a very conservative set of genes, mobile genetic elements are involved in the fine-tuning of regulatory networks that control the expression of these genes, which is reflected in the phenotypic characteristics of each species. Transposons are important sources of genome structures that are actively used to regulate the multicellular embryonic development. These structures include binding sites with transcription factors, enhancers and silencers, promoters, insulators, alternative splicing sites, non-coding RNA. Moreover, transposons are involved in the emergence and evolution of new protein-coding genes through exonization, domestication, and the formation of retrogenes. The activation of transposons is needed to regulate the differentiation and reproduction of cells in the body; however, in terminally differentiated cells, upon reaching predetermined sizes of organs, molecular systems are activated that block a further cascade of transposon activation. We suggest that the imperfection of systems aimed at specific suppression of transposon activity in mature cells may be the cause of aging and age-related diseases due to the pathological activation of mobile genetic elements. Identifying the tissue-specific mechanisms of inherited transposon activation in stem cells, as well as their pathological activation in terminally differentiated cells, may be the basis for finding ways to fight aging.

Differential Alternative Splicing Genes in Response to Boron Deficiency in Brassica napus.

Alternative splicing (AS) can increase transcriptome diversity, protein diversity and protein yield, and is an important mechanism to regulate plant responses to stress. Oilseed rape (Brassica napus L.), one of the main oil crops in China, shows higher sensitivity to boron (B) deficiency than other species. Here, we demonstrated AS changes that largely increased the diversity of the mRNA expressed in response to B deficiency in B. napus. Each gene had two or more transcripts on average. A total of 33.3% genes in both Qingyou10 (QY10, B-efficient cultivar) and Westar10 (W10, B-inefficient cultivar) showed AS in both B conditions. The types of AS events were mainly intron retention, 3′ alternative splice site, 5′ alternative splice site and exon skipping. The tolerance ability of QY10 was higher than that of W10, possibly because there were far more differential alternative splicing (DAS) genes identified in QY10 at low B conditions than in W10. The number of genes with both DAS and differentially expressed (DE) was far lower than that of the genes that were either with DAS or DE in QY10 and W10, suggesting that the DAS and DE genes were independent. Four Serine/Arginine-rich (SR) splicing factors, BnaC06g14780D, BnaA01g14750D, BnaA06g15930D and BnaC01g41640D, underwent differentially alternative splicing in both cultivars. There existed gene–gene interactions between BnaC06g14780D and the genes associated with the function of B in oilseed rape at low B supply. This suggests that oilseed rape could regulate the alterative pre-mRNA splicing of SR protein related genes to increase the plant tolerance to B deficiency.

Gymnotic Delivery of LNA Mixmers Targeting Viral SREs Induces HIV-1 mRNA Degradation.

Transcription of the HIV-1 provirus generates a viral pre-mRNA, which is alternatively spliced into more than 50 HIV-1 mRNAs encoding all viral proteins. Regulation of viral alternative splice site usage includes the presence of splicing regulatory elements (SREs) which can dramatically impact RNA expression and HIV-1 replication when mutated. Recently, we were able to show that two viral SREs, GI3-2 and ESEtat, are important players in the generation of viral vif, vpr and tat mRNAs. Furthermore, we demonstrated that masking these SREs by transfected locked nucleic acid (LNA) mixmers affect the viral splicing pattern and viral particle production. With regard to the development of future therapeutic LNA mixmer-based antiretroviral approaches, we delivered the GI3-2 and the ESEtat LNA mixmers “nakedly”, without the use of transfection reagents (gymnosis) into HIV-1 infected cells. Surprisingly, we observed that gymnotically-delivered LNA mixmers accumulated in the cytoplasm, and seemed to co-localize with GW bodies and induced degradation of mRNAs containing their LNA target sequence. The GI3-2 and the ESEtat LNA-mediated RNA degradation resulted in abrogation of viral replication in HIV-1 infected Jurkat and PM1 cells as well as in PBMCs.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Human Antithrombin III Minigene with an Optimized Splicing Pattern

Antithrombin III (AT3) belongs to the superfamily of serine protease inhibitors (serpins) and is a major anticoagulant in physiological conditions. Based on SERPINC1 gene, a minigene coding for human AT3, which is valuable for medicine and biotechnology, was constructed by minimizing the size of lengthy introns and preserving the splicing site-flanking sequences. An analysis of the minigene splicing pattern identified one correct AT3 transcript and two alternatively spliced transcripts, which formed either due to minigene exons 2 and 3 skipping or an aberrant exon insertion via splicing at cryptic splicing sites in intron 1 of the minigene. Site-directed mutagenesis of the cryptic splicing sites successfully optimized the splicing pattern of the AT3 minigene to completely prevent the generation of the alternative transcripts. The presence of the cryptic splicing sites in intron 1 of the minigene was confirmed with Human Splicing Finder v. 3.1 software, thus demonstrating that putative alternative splicing sites are possible to identify in minimized or hybrid introns of minigenes and to eliminate via mutagenesis before experimentally testing the minigene splicing patterns. The approach to the design of minigenes together with the bioinformatical analysis of the nucleotide sequences of minigene introns can be used to construct minigenes in order to generate transgenic animals producing economically valuable proteins in the milk.

Tissue-specific alternative splicing of pentatricopeptide repeat (PPR) family genes in Arabidopsis thaliana.

Alternative splicing is a post- and co-transcriptional regulatory mechanism of gene expression. Pentatricopeptide repeat (PPR) family proteins were recently found to be involved in RNA editing in plants. The aim of this study was to investigate the tissue-specific expression and alternative splicing of PPR family genes and their effects on protein structure and functionality. Of the 27 PPR genes in Arabidopsis thaliana, we selected six PPR genes of the P subfamily that are likely alternatively spliced, which were confirmed by sequencing. Four of these genes show intron retention, and the two remaining genes have 3' alternative-splicing sites. Alternative-splicing events occurred in the coding regions of three genes and in the 3' UTRs of the three remaining genes. We also identified five previously unannotated alternatively spliced isoforms of these PPR genes, which were confirmed by PCR and sequencing. Among these, three contain 3' alternative-splicing sites, one contains a 5' alternative-splicing site, and the remaining gene contains a 3'-5' alternative-splicing site. The new isoforms of two genes affect protein structure, and three other alternative-splicing sites are located in 3' UTRs. These findings suggest that tissue-specific expression of different alternatively spliced transcripts occurs in Arabidopsis, even at different developmental stages.

Identification of deleterious and regulatory genomic variations in known asthma loci

BackgroundCandidate gene and genome-wide association studies have identified hundreds of asthma risk loci. The majority of associated variants, however, are not known to have any biological function and are believed to represent markers rather than true causative mutations. We hypothesized that many of these associated markers are in linkage disequilibrium (LD) with the elusive causative variants.MethodsWe compiled a comprehensive list of 449 asthma-associated variants previously reported in candidate gene and genome-wide association studies. Next, we identified all sequence variants located within the 305 unique genes using whole-genome sequencing data from the 1000 Genomes Project. Then, we calculated the LD between known asthma variants and the sequence variants within each gene. LD variants identified were then annotated to determine those that are potentially deleterious and/or functional (i.e. coding or regulatory effects on the encoded transcript or protein).ResultsWe identified 10,130 variants in LD (r2 > 0.6) with known asthma variants. Annotations of these LD variants revealed that several have potentially deleterious effects including frameshift, alternate splice site, stop-lost, and missense. Moreover, 24 of the LD variants have been reported to regulate gene expression as expression quantitative trait loci (eQTLs).ConclusionsThis study is proof of concept that many of the genetic loci previously associated with complex diseases such as asthma are not causative but represent markers of disease, which are in LD with the elusive causative variants. We hereby report a number of potentially deleterious and regulatory variants that are in LD with the reported asthma loci. These reported LD variants could account for the original association signals with asthma and represent the true causative mutations at these loci.

EDA Fibronectin-TLR4 Axis Sustains Megakaryocyte Expansion and Inflammation during Bone Marrow Fibrosis Progression

Fibronectin (FN) is a glycoprotein of approximately 220 kDa, whose mRNA has three alternative splicing sites (termed EDA, EDB and IIICS) that allow 20 different isoforms of FN mRNA. Circulating plasma FN lacks both EDA and EDB segments and is a soluble form secreted by the hepatocytes, while cellular FN contains variable proportions of EDA and EDB segments and is organized as fibrils in the tissue matrix. FN containing EDA segment (EDA FN) presents unique biochemical properties as compared to the isoform lacking this domain, specifically: 1) pro-fibrotic, as the EDA-containing isoform becomes predominantly expressed during wound healing to sustain myofibroblast differentiation; 2) pro-inflammatory, as inclusion of EDA domain activates Toll Like Receptor 4 (TLR4), involved in the defense of the innate immune response following recognition of pathogens associated molecular patterns resulting in NF-κβ-dependent cytokine release; 3) pro-thrombotic, as EDA FN was shown to promote agonist-induced platelet aggregation and thrombus formation in vivo through TLR4 dependent mechanisms. EDA FN is instrumental in fibrogenesis but, to date, its expression and function in bone marrow (BM) fibrosis have not been explored. We identified the up-regulation, at both molecular and protein levels, of EDA FN in an experimental mouse model of BM fibrosis treated with supra-pharmacological doses of the thrombopoietin (TPO) mimetic Romiplostim (TPOhigh). Higher expression levels of EDA FN were detected in endothelial and stromal cells after inducing experimental fibrosis in wild type (WT) mice as compared to untreated controls. To unravel the role of EDA segment on the progression of BM fibrosis in vivo, we exploited two unique transgenic mouse models that present an aberrant splicing of the EDA exon. Mice with constitutive inclusion of EDA exon (EIIIA+/+), in tissues, were more prone to BM fibrosis development as compared to knockout mice (EIIIA-/-). Upon fibrosis induction with TPO mimetic treatment, EIIIA+/+ mice presented extensive reticulin deposition, thrombocytosis, splenomegaly and increased extra-medullary hematopoiesis.Applying numerous in vitro methods, we demonstrated that EDA FN as opposed to plasma FN lacking the EDA domain, induces megakaryocyte (Mk) proliferation in a TPO-independent fashion through engagement of TLR4 and activation of STAT5 and ERK 1/2 signalling pathways. Additionally, in vitro activation of EDA FN/TLR4 axis on Mks resulted in lipopolysaccharide-like responses, such as NF-kb activation, the release of pro-fibrotic Interleukin-6 (IL-6), and no anti-fibrotic TNF-a. Pharmacological inhibition of TLR4 in TPOhigh and EIIIA+/+/TPOhigh mice or TLR4 deletion in TPOhigh mice abrogated Mk hyperplasia, BM fibrosis, IL-6 release, extramedullary hematopoiesis and splenomegaly. Finally, developing a novel ELISA assay, we analysed samples from patients affected by primary myelofibrosis (PMF). PMF is one of the Philadelphia-negative Myeloproliferative Neoplasms (MPNs), a group of heterogeneous clonal disorders affecting the hematopoietic stem cell. MPNs are characterised by neoplastic proliferation of the myeloid lineage leading to an abnormal increase of platelet count in essential thrombocythemia (ET), red blood cells in polycythemia vera (PV), or Mks with BM fibrosis in primary myelofibrosis (PMF). To date, three major driver mutations have been identified affecting, respectively, the TPO receptor gene (MPL), the intracellular Janus Kinase 2 gene (JAK2) and the latest discovered mutations affecting the endoplasmic reticulum (ER) molecular chaperone calreticulin gene (CALR). All of these mutations induce a permanent activation of the JAK/STAT signalling pathways. We found that the EDA FN is increased in plasma and BM biopsies of PMF patients, as compared to healthy controls and ET patients, correlating with the fibrotic phase. In conclusion, we identified EDA FN/TLR4 as a new pathological axis that sustains Mk expansion and inflammation during BM fibrosis progression. EDA FN and TLR4 targeting has very promising potential for new therapeutic approaches in PMF patients with high phase BM fibrosis. DisclosuresNo relevant conflicts of interest to declare.