Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants.

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Role of Phytochromes in Red Light-Regulated Alternative Splicing in Arabidopsis thaliana: Impactful but Not Indispensable.

Light is both the main source of energy and a key environmental signal for plants. It regulates not only gene expression but also the tightly related processes of splicing and alternative splicing (AS). Two main pathways have been proposed to link light sensing with the splicing machinery. One occurs through a photosynthesis-related signal, and the other is mediated by photosensory proteins, such as red light-sensing phytochromes. Here, we evaluated the relative contribution of each of these pathways by performing a transcriptome-wide analysis of light regulation of AS in plants that do not express any functional phytochrome (phyQ). We found that an acute 2-h red-light pulse in the middle of the night induces changes in the splicing patterns of 483 genes in wild-type plants. Approximately 30% of these genes also showed strong light regulation of splicing patterns in phyQ mutant plants, revealing that phytochromes are important but not essential for the regulation of AS by R light. We then performed a meta-analysis of related transcriptomic datasets and found that different light regulatory pathways can have overlapping targets in terms of AS regulation. All the evidence suggests that AS is regulated simultaneously by various light signaling pathways, and the relative contribution of each pathway is highly dependent on the plant developmental stage.

The paralogues MAGOH and MAGOHB are oncogenic factors in high-grade gliomas and safeguard the splicing of cell division and cell cycle genes.

The exon junction complex (EJC) plays key roles throughout the lifespan of RNA and is particularly relevant in the nervous system. We investigated the roles of two EJC members, the paralogs MAGOH and MAGOHB, with respect to brain tumour development. High MAGOH/MAGOHB expression was observed in 14 tumour types; glioblastoma (GBM) showed the greatest difference compared to normal tissue. Increased MAGOH/MAGOHB expression was associated with poor prognosis in glioma patients, while knockdown of MAGOH/MAGOHB affected different cancer phenotypes. Reduced MAGOH/MAGOHB expression in GBM cells caused alterations in the splicing profile, including re-splicing and skipping of multiple exons. The binding profiles of EJC proteins indicated that exons affected by MAGOH/MAGOHB knockdown accumulated fewer complexes on average, providing a possible explanation for their sensitivity to MAGOH/MAGOHB knockdown. Transcripts (genes) showing alterations in the splicing profile are mainly implicated in cell division, cell cycle, splicing, and translation. We propose that high MAGOH/MAGOHB levels are required to safeguard the splicing of genes in high demand in scenarios requiring increased cell proliferation (brain development and GBM growth), ensuring efficient cell division, cell cycle regulation, and gene expression (splicing and translation). Since differentiated neuronal cells do not require increased MAGOH/MAGOHB expression, targeting these paralogs is a potential option for treating GBM.

Computational approaches for circRNAs prediction and in silico characterization.

Circular RNAs (circRNAs) are single-stranded and covalently closed non-coding RNA molecules originated from RNA splicing. Their functions include regulatory potential over other RNA species, such as microRNAs, messenger RNAs and RNA binding proteins. For circRNA identification, several algorithms are available and can be classified in two major types: pseudo-reference-based and split-alignment-based approaches. In general, the data generated from circRNA transcriptome initiatives is deposited on public specific databases, which provide a large amount of information on different species and functional annotations. In this review, we describe the main computational resources for the identification and characterization of circRNAs, covering the algorithms and predictive tools to evaluate its potential role in a particular transcriptomics project, including the public repositories containing relevant data and information for circRNAs, recapitulating their characteristics, reliability and amount of data reported.

Evaluating the Contribution of Cell Type-Specific Alternative Splicing to Variation in Lipid Levels.

BACKGROUND: Genome-wide association studies have identified hundreds of loci associated with lipid levels. However, the genetic mechanisms underlying most of these loci are not well-understood. Recent work indicates that changes in the abundance of alternatively spliced transcripts contribute to complex trait variation. Consequently, identifying genetic loci that associate with alternative splicing in disease-relevant cell types and determining the degree to which these loci are informative for lipid biology is of broad interest. METHODS: We analyze gene splicing in 83 sample-matched induced pluripotent stem cell (iPSC) and hepatocyte-like cell lines (n=166), as well as in an independent collection of primary liver tissues (n=96) to perform discovery of splicing quantitative trait loci (sQTLs). RESULTS: We observe that transcript splicing is highly cell type specific, and the genes that are differentially spliced between iPSCs and hepatocyte-like cells are enriched for metabolism pathway annotations. We identify 1384 hepatocyte-like cell sQTLs and 1455 iPSC sQTLs at a false discovery rate of <5% and find that sQTLs are often shared across cell types. To evaluate the contribution of sQTLs to variation in lipid levels, we conduct colocalization analysis using lipid genome-wide association data. We identify 19 lipid-associated loci that colocalize either with an hepatocyte-like cell expression quantitative trait locus or sQTL. Only 2 loci colocalize with both a sQTL and expression quantitative trait locus, indicating that sQTLs contribute information about genome-wide association studies loci that cannot be obtained by analysis of steady-state gene expression alone. CONCLUSIONS: These results provide an important foundation for future efforts that use iPSC and iPSC-derived cells to evaluate genetic mechanisms influencing both cardiovascular disease risk and complex traits in general.

Potential role of exitron-containing homeobox genes in cancer.

Homeobox genes are protagonists in developmental and cancer biology, making comprehending their regulation pivotal in multiple molecular pathways. Exitrons, also known as intronic exons, are new players in the transcriptional organization, providing additional splicing variants whose functions are still vastly unknown. Exitron splicing sites were identified in eight homeobox genes, which has not been yet debated in the scientific literature. Due to the intimate connection between homeobox genes and tumorigenesis, it is worth investing more time in understanding how these less explored exitron-containing transcriptional isoforms could play a role in modulating the homeobox gene's biological functions. The perspectives devised in this article are meant to instigate fresh debates on how the transcriptional variants retaining exitrons identified in the human homeobox genes HOXA1, HOXA9, HOXD8, NKX3.1, and DLX6 can be examined in the context of tumorigenesis. This article is categorized under: Cancer > Genetics/Genomics/Epigenetics.

UHMK1 is a novel splicing regulatory kinase.

The U2AF Homology Motif Kinase 1 (UHMK1) is the only kinase that contains the U2AF homology motif, a common protein interaction domain among splicing factors. Through this motif, UHMK1 interacts with the splicing factors SF1 and SF3B1, known to participate in the 3' splice site recognition during the early steps of spliceosome assembly. Although UHMK1 phosphorylates these splicing factors in vitro, the involvement of UHMK1 in RNA processing has not previously been demonstrated. Here, we identify novel putative substrates of this kinase and evaluate UHMK1 contribution to overall gene expression and splicing, by integrating global phosphoproteomics, RNA-seq, and bioinformatics approaches. Upon UHMK1 modulation, 163 unique phosphosites were differentially phosphorylated in 117 proteins, of which 106 are novel potential substrates of this kinase. Gene Ontology analysis showed enrichment of terms previously associated with UHMK1 function, such as mRNA splicing, cell cycle, cell division, and microtubule organization. The majority of the annotated RNA-related proteins are components of the spliceosome but are also involved in several steps of gene expression. Comprehensive analysis of splicing showed that UHMK1 affected over 270 alternative splicing events. Moreover, splicing reporter assay further supported UHMK1 function on splicing. Overall, RNA-seq data demonstrated that UHMK1 knockdown had a minor impact on transcript expression and pointed to UHMK1 function in epithelial-mesenchymal transition. Functional assays demonstrated that UHMK1 modulation affects proliferation, colony formation, and migration. Taken together, our data implicate UHMK1 as a splicing regulatory kinase, connecting protein regulation through phosphorylation and gene expression in key cellular processes.

The Intracellular and Secreted Sides of Osteopontin and Their Putative Physiopathological Roles.

Classically, osteopontin (OPN) has been described as a secreted glycophosprotein. Indeed, most data concerning its physiological and pathological roles are mainly related to the secreted OPN (sOPN). However, there are several instances in which intracellular OPN (iOPN) has been described, presenting some specific roles in distinct experimental models, such as in the immune system, cancer cells, and neurological disorders. We herein aimed to highlight and discuss some of these secreted and intracellular roles of OPN and their putative clinical and biological impacts. Moreover, by consolidating data from the OPN protein database, we also analyzed the occurrence of signal peptide (SP) sequences and putative subcellular localization, especially concerning currently known OPN splicing variants (OPN-SV). Comprehending the roles of OPN in its distinct cellular and tissue environments may provide data regarding the additional applications of this protein as biomarkers and targets for therapeutic purposes, besides further describing its pleiotropic roles.

Further delineation of the CWC27-associated spliceosomeopathy: Case report and review of the literature.

Pre-mRNA splicing factors are crucial in regulating transcript diversity, by removing introns from eukaryotic transcripts, an essential step in gene expression. Splicing of pre-mRNA is catalyzed by spliceosomes. CWC27 is a cyclophilin associated with spliceosome, in which genetic defects of its components have been linked to spliceosomopathies with clinical phenotypes including skeletal developmental defects, retinitis pigmentosa (RP), short stature, skeletal anomalies, and neurological disorders. We report two siblings (male and female) of Mexican descent with a novel homozygous frameshift variant in CWC27 and aim to highlight the cardinal features among the previously described 12 cases as well as expand the currently recognized phenotypic spectrum. Both siblings presented with a range of ocular and extraocular manifestations including novel features such as solitary kidney and tarsal coalition in the male sibling, together with gait abnormalities, and Hashimoto's thyroiditis in the female sibling. Finally, we highlight ectodermal involvement including sparse scalp hair, eyebrows and lashes, pigmentary differences, nail dysplasia, and dental anomalies as a core phenotype associated with the CWC27 spliceosomopathy.

Trans-splicing in the cestode Hymenolepis microstoma is constitutive across the life cycle and depends on gene structure and composition.

Spliced leader (SL) trans-splicing is a key process during mRNA maturation of many eukaryotes, in which a short sequence (SL) is transferred from a precursor SL-RNA into the 5' region of an immature mRNA. This mechanism is present in flatworms, in which it is known to participate in the resolution of polycistronic transcripts. However, most trans-spliced transcripts are not part of operons, and it is not clear if this process may participate in additional regulatory mechanisms in this group. In this work, we present a comprehensive analysis of SL trans-splicing in the model cestode Hymenolepis microstoma. We identified four different SL-RNAs which are indiscriminately trans-spliced to 622 gene models. SL trans-splicing is enriched in constitutively expressed genes and does not appear to be regulated throughout the life cycle. Operons represented at least 20% of all detected trans-spliced gene models, showed conservation to those of the cestode Echinococcus multilocularis, and included complex loci such as an alternative operon (processed as either a single gene through cis-splicing or as two genes of a polycistron). Most insertion sites were identified in the 5' untranslated region (UTR) of monocistronic genes. These genes frequently contained introns in the 5' UTR, in which trans-splicing used the same acceptor sites as cis-splicing. These results suggest that, unlike other eukaryotes, trans-splicing is associated with internal intronic promoters in the 5' UTR, resulting in transcripts with strong splicing acceptor sites without competing cis-donor sites, pointing towards a simple mechanism driving the evolution of novel SL insertion sites.

IARA: a complete and curated atlas of the biogenesis of spliceosome machinery during RNA splicing.

Splicing is one of the most important post-transcriptional processing systems and is responsible for the generation of transcriptome diversity in all living eukaryotes. Splicing is regulated by the spliceosome machinery, which is responsible for each step of primary RNA processing. However, current molecules and stages involved in RNA splicing are still spread over different studies. Thus, a curated atlas of spliceosome-related molecules and all involved stages during RNA processing can provide all researchers with a reliable resource to better investigate this important mechanism. Here, we present IARA (website access: https://pucpr-bioinformatics.github.io/atlas/), an extensively curated and constantly updated catalog of molecules involved in spliceosome machinery. IARA has a map of the steps involved in the human splicing mechanism, and it allows a detailed overview of the molecules involved throughout the distinct steps of splicing.

Do not panic: An intron-centric guide to alternative splicing.

This review is an attempt to establish concepts of splicing and alternative splicing giving proper relevance to introns, the key actors in this mechanism. It might also work as a guide for those who found their favorite gene undergoes alternative splicing and could benefit from gaining a theoretical framework to understand the possible impacts of this process. This is not a thorough review of all the work in the field, but rather a critical review of some of the most relevant work done to understand the underlying mechanisms of splicing and the key questions that remain unanswered such as: What is the physiological relevance of alternative splicing? What are the functions of the different outcomes? To what extent do different alternative splicing types contribute to the proteome? Intron retention is the most frequent alternative splicing event in plants and, although scientifically neglected, it is also common in animals. This is a heterogeneous type of alternative splicing that includes different sub-types with features that have distinctive consequences in the resulting transcripts. Remarkably, intron retention can be a dead end for a transcript, but it could also be a stable intermediate whose processing is resumed upon a particular signal or change in the cell status. New sequencing technologies combined with the study of intron lariats in different conditions might help to answer key questions and could help us to understand the actual relevance of introns in gene expression regulation.

PICLN modulates alternative splicing and light/temperature responses in plants.

Plants undergo transcriptome reprograming to adapt to daily and seasonal fluctuations in light and temperature conditions. While most efforts have focused on the role of master transcription factors, the importance of splicing factors modulating these processes is now emerging. Efficient pre-mRNA splicing depends on proper spliceosome assembly, which in plants and animals requires the methylosome complex. Ion Chloride nucleotide-sensitive protein (PICLN) is part of the methylosome complex in both humans and Arabidopsis (Arabidopsis thaliana), and we show here that the human PICLN ortholog rescues phenotypes of Arabidopsis picln mutants. Altered photomorphogenic and photoperiodic responses in Arabidopsis picln mutants are associated with changes in pre-mRNA splicing that partially overlap with those in PROTEIN ARGININE METHYL TRANSFERASE5 (prmt5) mutants. Mammalian PICLN also acts in concert with the Survival Motor Neuron (SMN) complex component GEMIN2 to modulate the late steps of UsnRNP assembly, and many alternative splicing events regulated by PICLN but not PRMT5, the main protein of the methylosome, are controlled by Arabidopsis GEMIN2. As with GEMIN2 and SM PROTEIN E1/PORCUPINE (SME1/PCP), low temperature, which increases PICLN expression, aggravates morphological and molecular defects of picln mutants. Taken together, these results establish a key role for PICLN in the regulation of pre-mRNA splicing and in mediating plant adaptation to daily and seasonal fluctuations in environmental conditions.

An Unusual U2AF2 Inhibits Splicing and Attenuates the Virulence of the Human Protozoan Parasite Entamoeba histolytica.

E. histolytica is the etiological agent of intestinal amebiasis and liver abscesses, which still poses public health threat globally. Metronidazole is the drug of choice against amebiasis. However, metronidazole-resistant amoebic clinical isolates and strains have been reported recently, challenging the efforts for amebiasis eradication. In search of alternative treatments, E. histolytica transcriptomes have shown the association of genes involved in RNA metabolism with the virulence of the parasite. Among the upregulated genes in amoebic liver abscesses are the splicing factors EhU2AF2 and a paralog of EhSF3B1. For this reason and because EhU2AF2 contains unusual KH-QUA2 (84KQ) motifs in its lengthened C-terminus domain, here we investigated how the role of EhU2AF2 in pre-mRNA processing impacts the virulence of the parasite. We found that 84KQ is involved in splicing inhibition/intron retention of several virulence and non-virulence-related genes. The 84KQ domain interacts with the same domain of the constitutive splicing factor SF1 (SF1KQ), both in solution and when SF1KQ is bound to branchpoint signal RNA probes. The 84KQ-SF1KQ interaction prevents splicing complex E to A transition, thus inhibiting splicing. Surprisingly, the deletion of the 84KQ domain in EhU2AF2 amoeba transformants increased splicing and enhanced the in vitro and in vivo virulence phenotypes. We conclude that the interaction of the 84KQ and SF1KQ domains, probably involving additional factors, tunes down Entamoeba virulence by favoring intron retention.

Counteracting chromatin effects of a splicing-correcting antisense oligonucleotide improves its therapeutic efficacy in spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a motor-neuron disease caused by mutations of the SMN1 gene. The human paralog SMN2, whose exon 7 (E7) is predominantly skipped, cannot compensate for the lack of SMN1. Nusinersen is an antisense oligonucleotide (ASO) that upregulates E7 inclusion and SMN protein levels by displacing the splicing repressors hnRNPA1/A2 from their target site in intron 7. We show that by promoting transcriptional elongation, the histone deacetylase inhibitor VPA cooperates with a nusinersen-like ASO to promote E7 inclusion. Surprisingly, the ASO promotes the deployment of the silencing histone mark H3K9me2 on the SMN2 gene, creating a roadblock to RNA polymerase II elongation that inhibits E7 inclusion. By removing the roadblock, VPA counteracts the chromatin effects of the ASO, resulting in higher E7 inclusion without large pleiotropic effects. Combined administration of the nusinersen-like ASO and VPA in SMA mice strongly synergizes SMN expression, growth, survival, and neuromuscular function.

CircRNAs: promising factors for regulating angiogenesis in colorectal cancer.

Colorectal cancer (CRC) is one of the most common cancers in the world. The incidence rate of cancer is high. The overall response to traditional treatment methods such as surgery, radiotherapy, and chemotherapy is not very satisfactory. Therefore, finding new therapeutic targets is very important for improving CRC treatment. In recent reports, the role of circRNAs in regulating colorectal angiogenesis has been gradually revealed. CircRNAs can indirectly act on angiogenesis pathways and regulate the expression of growth factors such as vascular endothelial growth factor (VEGF). CircRNAs are endogenous noncoding RNAs formed by pre-mRNAs through exon circular splicing. The covalent closed-loop structure makes these RNAs highly conserved and stable. CircRNAs have been found in human plasma, serum, urine, and other body fluids. Their highly conserved characteristics play important roles in many biological activities. CircRNAs can participate in the progression of many diseases by sponging miRNAs, interacting with proteins, and regulating transcription. Angiogenesis can provide nutrients and oxygen for tumour proliferation and metastasis. Angiogenesis is an important sign of the formation of the tumour microenvironment. Here, we will summarize the role of the latest circRNAs in the mechanism of angiogenesis in CRC and provide potential therapeutic targets for clinical treatment.

Proteomics Uncovers Novel Components of an Interactive Protein Network Supporting RNA Export in Trypanosomes.

In trypanosomatids, transcription is polycistronic and all mRNAs are processed by trans-splicing, with export mediated by noncanonical mechanisms. Although mRNA export is central to gene regulation and expression, few orthologs of proteins involved in mRNA export in higher eukaryotes are detectable in trypanosome genomes, necessitating direct identification of protein components. We previously described conserved mRNA export pathway components in Trypanosoma cruzi, including orthologs of Sub2, a component of the TREX complex, and eIF4AIII (previously Hel45), a core component of the exon junction complex (EJC). Here, we searched for protein interactors of both proteins using cryomilling and mass spectrometry. Significant overlap between TcSub2 and TceIF4AIII-interacting protein cohorts suggests that both proteins associate with similar machinery. We identified several interactions with conserved core components of the EJC and multiple additional complexes, together with proteins specific to trypanosomatids. Additional immunoisolations of kinetoplastid-specific proteins both validated and extended the superinteractome, which is capable of supporting RNA processing from splicing through to nuclear export and cytoplasmic events. We also suggest that only proteomics is powerful enough to uncover the high connectivity between multiple aspects of mRNA metabolism and to uncover kinetoplastid-specific components that create a unique amalgam to support trypanosome mRNA maturation.

A novel strategy to uncover specific GO terms/phosphorylation pathways in phosphoproteomic data in Arabidopsis thaliana.

BACKGROUND: Proteins are the workforce of the cell and their phosphorylation status tailors specific responses efficiently. One of the main challenges of phosphoproteomic approaches is to deconvolute biological processes that specifically respond to an experimental query from a list of phosphoproteins. Comparison of the frequency distribution of GO (Gene Ontology) terms in a given phosphoproteome set with that observed in the genome reference set (GenRS) is the most widely used tool to infer biological significance. Yet, this comparison assumes that GO term distribution between the phosphoproteome and the genome are identical. However, this hypothesis has not been tested due to the lack of a comprehensive phosphoproteome database. RESULTS: In this study, we test this hypothesis by constructing three phosphoproteome databases in Arabidopsis thaliana: one based in experimental data (ExpRS), another based in in silico phosphorylation protein prediction (PredRS) and a third that is the union of both (UnRS). Our results show that the three phosphoproteome reference sets show default enrichment of several GO terms compared to GenRS, indicating that GO term distribution in the phosphoproteomes does not match that of the genome. Moreover, these differences overshadow the identification of GO terms that are specifically enriched in a particular condition. To overcome this limitation, we present an additional comparison of the sample of interest with UnRS to uncover GO terms specifically enriched in a particular phosphoproteome experiment. Using this strategy, we found that mRNA splicing and cytoplasmic microtubule compounds are important processes specifically enriched in the phosphoproteome of dark-grown Arabidopsis seedlings. CONCLUSIONS: This study provides a novel strategy to uncover GO specific terms in phosphoproteome data of Arabidopsis that could be applied to any other organism. We also highlight the importance of specific phosphorylation pathways that take place during dark-grown Arabidopsis development.

Chronic Granulomatous Disease and Myelodysplastic Syndrome in a Patient with a Novel Mutation in CYBB.

Chronic Granulomatous Disease (CGD) is an inborn error of immunity characterized by impaired phagocyte function, recurrent fungal and bacterial infections and granuloma formation in multiple organs. Pediatric myelodysplastic Syndrome (MDS) is a rare hematological stem cell disease that leads to an ineffective hematopoiesis with variable risk of evolution to acute leukemias. Both disorders are rare and have distinct pathophysiologic mechanisms, with no known association. A 7-month-old boy presenting with recurrent infections and anemia at age 2 months underwent immunological, hematological and genetic investigation that culminated in the diagnosis of both CGD and MDS. Next generation sequencing was performed and identified a silent variant predicted as of Uncertain Significance, located in the splicing site at the end of exon 5 in CYBB. CYBB variants account for at least two thirds of CGD cases, but no previous descriptions of this variant were found in ClinVar or The Human Gene Mutation Database (HGMD) databases. We were able to demonstrate an exon 5 skipping on the proband's cDNA, which strongly suggests the disruption of the NADPH oxidase complex, abrogating the formation of reactive oxygen species from neutrophils. Moreover, erythroid cell lineage could be also affected by NADPH oxidase complex damages. Further investigation is needed to evaluate the potential effect of CYBB gene alterations in hematopoiesis, as well as in MDS and CGD association.

Intronic variants of MITF (rs7623610) and CREB1 (rs10932201) genes may enhance splicing efficiency in human melanoma cell line.

We previously reported that intronic single nucleotide variations (SNVs) in MITF (c.938-325G>A, rs7623610) and CREB1 (c.303+373G>A, rs10932201) genes were associated with risk, aggressiveness, and prognosis of cutaneous melanoma (CM). In this study, we investigated the influence of the above SNVs in splicing patterns and efficiency. We constructed minigenes with wild type and variant alleles from MITF and CREB1 to assess the effect of the SNVs on splicing. The minigenes were transfected in the human melanoma cell line (SK-MEL-28). RT-PCR and DNA sequencing investigated the constructs' splicing patterns. Minigenes constructs' splicing efficiency and HNRNPA1 and SF1 splicing genes' expression were investigated by qPCR. We found that MITF and CREB1 SNVs did not alter the splicing pattern, but they influenced the splicing efficiency. MITF-A (p= 0.03) and CREB1-A (p= 0.005) variant minigenes yielded an increase of mRNA generated from the constructions. Additionally, lower mRNA levels of HNRNPA1 and SF1 were seen in the variant minigenes MITF-A (p= 0.04) and CREB1-A (p= 0.005). We described for the first time the potential importance of MITF rs7623610 and CREB1 rs10932201 SNVs in splicing efficiency and its relationship with CM.