RNA Splicing and Disease: Animal Models to Therapies.

Alternative splicing of pre-mRNA increases genetic diversity, and recent studies estimate that most human multiexon genes are alternatively spliced. If this process is not highly regulated and accurate, it leads to mis-splicing events, which may result in proteins with altered function. A growing body of work has implicated mis-splicing events in a range of diseases, including cancer, neurodegenerative diseases, and muscular dystrophies. Understanding the mechanisms that cause aberrant splicing events and how this leads to disease is vital for designing effective therapeutic strategies. In this review, we focus on advances in therapies targeting splicing, and highlight the animal models developed to recapitulate disease phenotypes as a model for testing these therapies.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium.

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Splicing factor SRSF1 negatively regulates alternative splicing of MDM2 under damage.

Genotoxic stress induces alternative splicing of the oncogene MDM2 generating MDM2-ALT1, an isoform attributed with tumorigenic properties. However, the mechanisms underlying this event remain unclear. Here we explore MDM2 splicing regulation by utilizing a novel minigene that mimics endogenous MDM2 splicing in response to UV and cisplatinum-induced DNA damage. We report that exon 11 is necessary and sufficient for the damage-specific alternative splicing of the MDM2 minigene and that the splicing factor SRSF1 binds exon 11 at evolutionarily conserved sites. Interestingly, mutations disrupting this interaction proved sufficient to abolish the stress-induced alternative splicing of the MDM2 minigene. Furthermore, SRSF1 overexpression promoted exclusion of exon 11, while its siRNA-mediated knockdown prevented the stress-induced alternative splicing of endogenous MDM2. Additionally, we observed elevated SRSF1 levels under stress and in tumors correlating with the expression of MDM2-ALT1. Notably, we demonstrate that MDM2-ALT1 splicing can be blocked by targeting SRSF1 sites on exon 11 using antisense oligonucleotides. These results present conclusive evidence supporting a negative role for SRSF1 in MDM2 alternative splicing. Importantly, we define for the first time, a clear-cut mechanism for the regulation of damage-induced MDM2 splicing and present potential strategies for manipulating MDM2 expression via splicing modulation.

Regulation of photoreceptor gene transcription via a highly conserved transcriptional regulatory element by vsx gene products.

PURPOSE: The photoreceptor conserved element-1 (PCE-1) sequence is found in the transcriptional regulatory regions of many genes expressed in photoreceptors. The retinal homeobox (Rx or Rax) gene product functions by binding to PCE-1 sites. However, other transcriptional regulators have also been reported to bind to PCE-1. One of these, vsx2, is expressed in retinal progenitor and bipolar cells. The purpose of this study is to identify Xenopus laevis vsx gene products and characterize vsx gene product expression and function with respect to the PCE-1 site. METHODS: X. laevis vsx gene products were amplified with PCR. Expression patterns were determined with in situ hybridization using whole or sectioned X. laevis embryos and digoxigenin- or fluorescein-labeled antisense riboprobes. DNA binding characteristics of the vsx gene products were analyzed with electrophoretic mobility shift assays (EMSAs) using in vitro translated proteins and radiolabeled oligonucleotide probes. Gene transactivation assays were performed using luciferase-based reporters and in vitro transcribed effector gene products, injected into X. laevis embryos. RESULTS: We identified one vsx1 and two vsx2 gene products. The two vsx2 gene products are generated by alternate mRNA splicing. We verified that these gene products are expressed in the developing retina and that expression resolves into distinct cell types in the mature retina. Finally, we found that vsx gene products can bind the PCE-1 site in vitro and that the two vsx2 isoforms have different gene transactivation activities. CONCLUSIONS: vsx gene products are expressed in the developing and mature neural retina. vsx gene products can bind the PCE-1 site in vitro and influence the expression of a rhodopsin promoter-luciferase reporter gene. The two isoforms of vsx have different gene transactivation activities in this reporter gene system.

PDPR Gene Variants Predisposing to Papillary Thyroid Cancer.

Background: Papillary thyroid cancer (PTC) is the predominant subtype of thyroid cancer (THCA), and it can cluster in families with an autosomal dominant (AD) inheritance pattern. The aim of this study was to identify novel genes and mechanisms underlying PTC susceptibility. Methods: Our previous investigation of 17 AD PTC families led us to conduct a deeper analysis on one family (Family Q) with whole-genome sequencing data from 3 PTC-affected individuals. In addition, 323 sporadic THCA cases from Avatar data and 12 familial adenomatous polyposis (FAP) individuals with secondary THCA were screened for pyruvate dehydrogenase phosphatase regulatory (PDPR) variants. CRISPR-Cas9 was used to create PDPR-deficient THCA (TPC1) and transformed normal thyroid cell lines (N-Thyori3-1) to study the metabolic consequences of PDPR loss. Results: We found truncating PDPR splice donor variants (NM_017990.4:c.361 + 1G>C) in all affected PTC Family Q members, and another PDPR splice donor variant (NM_017990.4:c.443 + 1G>C) in a sporadic PTC case. In addition, an ultra-rare missense variant was found in an FAP-PTC patient. The PDPR-deficient cells presented with elevated phosphorylation of pyruvate dehydrogenase and altered glucose metabolism, implying that PDPR plays an essential part in regulating glucose metabolism in thyroid cells. Conclusions: Our finding of novel truncating germline variants in PDPR in Family Q and additional cohorts suggests a role for PDPR loss in PTC predisposition. Also, somatic and RNA sequencing from the thyroid carcinoma (Firehouse Legacy) data showed that PDPR gene expression is much lower in THCA tumor tissue compared with matching normal tissue. Thus, PDPR appears to have a loss of function effect on THCA tumorigenesis.

Splice-switching of the insulin receptor pre-mRNA alleviates tumorigenic hallmarks in rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is an aggressive pediatric tumor with a poor prognosis for metastasis and recurrent disease. Large-scale sequencing endeavors demonstrate that Rhabdomyosarcomas have a dearth of precisely targetable driver mutations. However, IGF-2 signaling is known to be grossly altered in RMS. The insulin receptor (IR) exists in two alternatively spliced isoforms, IR-A and IR-B. The IGF-2 signaling molecule binds both its innate IGF-1 receptor as well as the insulin receptor variant A (IR-A) with high affinity. Mitogenic and proliferative signaling via the canonical IGF-2 pathway is, therefore, augmented by IR-A. This study shows that RMS patients express increased IR-A levels compared to control tissues that predominantly express the IR-B isoform. We also found that Hif-1α is significantly increased in RMS tumors, portraying their hypoxic phenotype. Concordantly, the alternative splicing of IR adapts to produce more IR-A in response to hypoxic stress. Upon examining the pre-mRNA structure of the gene, we identified a potential hypoxia-responsive element, which is also the binding site for the RNA-binding protein CUG-BP1 (CELF1). We designed Splice Switching Oligonucleotides (SSO) against this binding site to decrease IR-A levels in RMS cell lines and, consequently, rescue the IR-B expression levels. SSO treatment resulted in a significant reduction in cell proliferation, migration, and angiogenesis. Our data shows promising insight into how impeding the IGF-2 pathway by reducing IR-A expression mitigates tumor growth. It is evident that Rhabdomyosarcomas use IR alternative splicing as yet another survival strategy that can be exploited as a therapeutic intervention in conjunction with already established anti-IGF-1 receptor therapies.

Transcriptome analysis discloses dysregulated genes in normal appearing tumor-adjacent thyroid tissues from patients with papillary thyroid carcinoma.

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. The molecular characteristics of histologically normal appearing tissue adjacent to the tumor (NAT) from PTC patients are not well characterized. The aim of this study was to characterize the global gene expression profile of NAT and compare it with those of normal and tumor thyroid tissues. We performed total RNA sequencing with fresh frozen thyroid tissues from a cohort of three categories of samples including NAT, normal thyroid (N), and PTC tumor (T). Transcriptome analysis shows that NAT presents a unique gene expression profile, which was not associated with sex or the presence of lymphocytic thyroiditis. Among the differentially expressed genes (DEGs) of NAT vs N, 256 coding genes and 5 noncoding genes have been reported as cancer genes involved in cell proliferation, apoptosis, and/or tumorigenesis. Bioinformatics analysis with Ingenuity Pathway Analysis software revealed that "Cancer, Organismal Injury and Abnormalities, Cellular Response to Therapeutics, and Cellular Movement" were major dysregulated pathways in the NAT tissues. This study provides improved insight into the complexity of gene expression changes in the thyroid glands of patients with PTC.

A novel essential splice site variant in SPTB in a large hereditary spherocytosis family.

BACKGROUND: We studied a large family with 22 individuals affected with autosomal dominant hereditary spherocytosis (HS). METHODS: Genome-wide linkage, whole-genome sequencing (WGS), Sanger sequencing, RT-PCR, and ToPO TA cloning analyses were performed. RESULTS: We revealed a heterozygous G>A transition in the 14q23 locus, at position +1 of the intron 8 donor splice site of the spectrin beta, erythrocytic (SPTB) gene. This splice variant (SPTB c.1064+1G>A) was confirmed by Sanger sequencing and showed complete co-segregation with HS in the family. Further RT-PCR reactions and sequencing analysis indicated that the variant leads to the exclusion of exon 8 and subsequent frameshift in exon 9 and a premature stop codon in SPTB. Translation of the altered allele would lead to a truncation with a loss of all spectrin repeat domains in SPTB protein. CONCLUSION: This variant is novel and has not been found in any databases. We propose that this splice variant explains the spherocytosis phenotype observed in this large family.

SRSF2 Regulation of MDM2 Reveals Splicing as a Therapeutic Vulnerability of the p53 Pathway.

MDM2 is an oncogene and critical negative regulator of tumor suppressor p53. Genotoxic stress causes alternative splicing of MDM2 transcripts, which leads to alterations in p53 activity and contributes to tumorigenesis. MDM2-ALT1 is one of the alternatively spliced transcripts predominantly produced in response to genotoxic stress, and is comprised of terminal coding exons 3 and 12. Previously, we found that SRSF1 induces MDM2-ALT1 by promoting MDM2 exon 11 skipping. Here we report that splicing regulator SRSF2 antagonizes the regulation of SRSF1 by facilitating the inclusion of exon 11 through binding at two conserved exonic splicing enhancers. Overexpression of SRSF2 reduced the generation of MDM2-ALT1 under genotoxic stress, whereas SRSF2 knockdown induced the expression of MDM2-ALT1 in the absence of genotoxic stress. Blocking the exon 11 SRSF2-binding sites using oligonucleotides promoted MDM2-ALT1 splicing and induced p53 protein expression, and apoptosis in p53 wild-type cells. The regulation of MDM2 splicing by SRSF2 is also conserved in mice, as mutation of one SRSF2-binding site in Mdm2 exon 11, using CRISPR-Cas9, increased the expression of the MDM2-ALT1 homolog Mdm2-MS2. IMPLICATIONS: Taken together, the data indicate that modulating MDM2 splicing may be a useful tool for fine-tuning p53 activity in response to genotoxic stress.

Characterizing the function of EPB41L4A in the predisposition to papillary thyroid carcinoma.

Papillary thyroid carcinoma (PTC) is the most common histotype of thyroid carcinoma. The heritability of PTC is high compared to other cancers, but its underlying causes are unknown. A recent genome-wide association study revealed the association of a variant at the 5q22 locus, rs73227498, with PTC predisposition. We report that rs17134155, a variant in high linkage disequilibrium with rs73227498, is located in an enhancer region downstream of coding transcripts of EPB41L4A. Rs17134155 showed significant enhancer activity in luciferase assays, and haplotypes containing the protective allele of this variant conferred a significantly lower risk of PTC. While the index SNP, rs73227498, acted as a significant cis-eQTL for expression of EPB41L4A, rs17134155 was a significant cis-sQTL for the alternative splicing of a non-coding transcript of EPB41L4A, called EPB41L4A-203. We also performed knockdown of EPB41L4A followed by microarray analysis. Some of the top differentially-expressed genes were represented among regulators of the WNT/ß-catenin signaling pathway. Our results indicate that an enhancer region at 5q22 regulates the expression and splicing of EPB41L4A transcripts. We also provide evidence that EPB41L4A expression is involved in regulating growth and differentiation pathways, suggesting that decreased expression of EPB41L4A is a mechanism in the predisposition to PTC.

A Truncating Germline Mutation of TINF2 in Individuals with Thyroid Cancer or Melanoma Results in Longer Telomeres.

Background: Our genome sequencing analysis revealed a frameshift mutation in the shelterin gene TINF2 in a large family with individuals affected with papillary thyroid carcinoma (PTC) and melanoma. Here, we further characterized the mutation and screened for coding variants in the 6 shelterin genes in 24 families. Methods: Sanger sequencing was performed to screen for the TINF2 mutation in the key family. Quantitative reverse transcription-polymerase chain reaction (PCR) was used for TINF2 gene expression analysis. Exogenous expression and co-immunoprecipitation techniques were used for assessing TINF2 binding to TERF1. Relative telomere length (RTL) was quantified in DNAs from lymphocytes by using quantitative real-time PCR. Whole exome sequencing (WES) was performed in seven families with individuals affected with PTC and other cancer types. Screening for DNA variants in shelterin genes was performed by using whole genome sequencing data from 17 families and WES data from 7 further families. Results: The TINF2 mutation (TINF2 p.Trp198fs) showed complete co-segregation with PTC and melanoma in the key family. The mutation is not reported in databases and not identified in 23 other families we screened. The expression of TINF2 was borderline reduced in individuals with the mutation. The truncated TINF2 protein showed abolished binding to TERF1. The RTL in the individuals with the mutation was significantly longer when compared with those without the mutation from the same family as well as compared with 62 healthy controls. Among the 24 families, we identified 3 missense and 1 synonymous variant(s) in 2 shelterin genes (TINF2 and ACD). Conclusions: The rare frameshift mutation in the TINF2 gene and the associated longer telomere length suggest that dysregulated telomeres could be a mechanism predisposing to PTC and melanoma. DNA coding variants in shelterin genes are rare. Further studies are required to evaluate the roles of variants in shelterin genes in thyroid cancer and melanoma.

Clinical implications of GWAS variants associated with differentiated thyroid cancer.

The genetic risk of differentiated thyroid cancer (DTC) probably consists of multiple low-penetrance, single-nucleotide polymorphisms (SNP). Such markers are difficult to uncover by linkage analysis but can be revealed by association studies. Genome-wide association studies (GWASs) have uncovered 31 SNPs associated with DTC. These markers carry a low to moderate risk for DTC, but their cumulative effect increases with each successive risk allele. These data support the important contribution of low penetrance variants in the pathogenesis of DTC. Contrary to somatic mutations such as BRAFV600E, germline variants can be ascertained prior to surgical treatment. Therefore, we hypothesise that GWAS SNPs might impact the clinical course of DTC and we can benefit from this knowledge in choosing a treatment strategy. Several associations between clinical factors and GWAS markers have been reported so far. The most important are associations between rs966423 and mortality (HR = 1.60, p = 0.038), extrathyroidal extension (ETE) (OR = 1.57, p = 0.019); rs965513 and tumour diameter (slope of regression 0.14, p = 0.025), lymph node metastasis (OR = 1.59, p = 0.030) and ETE (OR = 1.29, p = 0.045); rs944289 and distant metastasis (OR = 0.58, p = 0.042); and rs116909374 and lymph node metastasis (OR = 0.61, p = 0.016). These findings show that GWAS SNPs are not only the ignition factors (together with environmental factors) for malignant transformation of thyrocytes but might also impact the clinical course of DTC. Surprisingly, it is not always the risk allele for DTC that is associated with worse clinical outcome. The second interesting observation is that GWAS SNPs show different associations with DTC clinical features depending on their histological subtypes. These point to the complexity of DTC with putatively different roles of genes at different stages of DTC development. (Endokrynol Pol 2019; 70 (5): 423-429).

Identification of Rare Variants Predisposing to Thyroid Cancer.

Background: Familial non-medullary thyroid cancer (NMTC) accounts for a relatively small proportion of thyroid cancer cases, but it displays strong genetic predisposition. So far, only a few NMTC susceptible genes and low-penetrance variants contributing to NMTC have been described. This study aimed to identify rare germline variants that may predispose individuals to NMTC by sequencing a cohort of 17 NMTC families. Methods: Whole-genome sequencing and genome-wide linkage analysis were performed in 17 NMTC families. MendelScan and BasePlayer were applied to screen germline variants followed by customized filtering. The remaining candidate variants were subsequently validated by Sanger sequencing. A panel of 277 known cancer predisposition genes was also screened in these families. Results: A total of 41 rare coding candidate variants in 40 genes identified by whole-genome sequencing are reported, including 24 missense, five frameshift, five splice change, and seven nonsense variants. Sanger sequencing confirmed all 41 rare variants and proved their co-segregation with NMTC in the extended pedigrees. In silico functional analysis of the candidate genes using Ingenuity Pathway Analysis showed that cancer was the top category of "Diseases and Disorders." Additionally, a targeted search displayed six variants in known cancer predisposition genes, including one frameshift variant and five missense variants. Conclusions: The data identify rare germline variants that may play important roles in NMTC predisposition. It is proposed that in future research including functional characterization, these variants and genes be considered primary candidates for thyroid cancer predisposition.

Stress-induced alternative splice forms of MDM2 and MDMX modulate the p53-pathway in distinct ways.

MDM2 and MDMX are the chief negative regulators of the tumor-suppressor protein p53 and are essential for maintaining homeostasis within the cell. In response to genotoxic stress and also in several cancer types, MDM2 and MDMX are alternatively spliced. The splice variants MDM2-ALT1 and MDMX-ALT2 lack the p53-binding domain and are incapable of negatively regulating p53. However, they retain the RING domain that facilitates dimerization of the full-length MDM proteins. Concordantly, MDM2-ALT1 has been shown to lead to the stabilization of p53 through its interaction with and inactivation of full-length MDM2. The impact of MDM2-ALT1 expression on the p53 pathway and the nature of its interaction with MDMX remain unclear. Also, the role of the architecturally similar MDMX-ALT2 and its influence of the MDM2-MDMX-p53 axis are yet to be elucidated. We show here that MDM2-ALT1 is capable of binding full-length MDMX as well as full-length MDM2. Additionally, we demonstrate that MDMX-ALT2 is able to dimerize with both full-length MDMX and MDM2 and that the expression of MDM2-ALT1 and MDMX-ALT2 leads to the upregulation of p53 protein, and also of its downstream target p21. Moreover, MDM2-ALT1 expression causes cell cycle arrest in the G1 phase in a p53 and p21 dependent manner, which is consistent with the increased levels of p21. Finally we present evidence that MDM2-ALT1 and MDMX-ALT2 expression can activate subtly distinct subsets of p53-transcriptional targets implying that these splice variants can modulate the p53 tumor suppressor pathway in unique ways. In summary, our study shows that the stress-inducible alternative splice forms MDM2-ALT1 and MDMX-ALT2 are important modifiers of the p53 pathway and present a potential mechanism to tailor the p53-mediated cellular stress response.

Stress-induced isoforms of MDM2 and MDM4 correlate with high-grade disease and an altered splicing network in pediatric rhabdomyosarcoma.

Pediatric rhabdomyosarcoma (RMS) is a morphologically and genetically heterogeneous malignancy commonly classified into three histologic subtypes, namely, alveolar, embryonal, and anaplastic. An issue that continues to challenge effective RMS patient prognosis is the dearth of molecular markers predictive of disease stage irrespective of tumor subtype. Our study involving a panel of 70 RMS tumors has identified specific alternative splice variants of the oncogenes Murine Double Minute 2 (MDM2) and MDM4 as potential biomarkers for RMS. Our results have demonstrated the strong association of genotoxic-stress inducible splice forms MDM2-ALT1 (91.6% Intergroup Rhabdomyosarcoma Study Group stage 4 tumors) and MDM4-ALT2 (90.9% MDM4-ALT2-positive T2 stage tumors) with high-risk metastatic RMS. Moreover, MDM2-ALT1-positive metastatic tumors belonged to both the alveolar (50%) and embryonal (41.6%) subtypes, making this the first known molecular marker for high-grade metastatic disease across the most common RMS subtypes. Furthermore, our results show that MDM2-ALT1 expression can function by directly contribute to metastatic behavior and promote the invasion of RMS cells through a matrigel-coated membrane. Additionally, expression of both MDM2-ALT1 and MDM4-ALT2 increased anchorage-independent cell-growth in soft agar assays. Intriguingly, we observed a unique coordination in the splicing of MDM2-ALT1 and MDM4-ALT2 in approximately 24% of tumor samples in a manner similar to genotoxic stress response in cell lines. To further explore splicing network alterations with possible relevance to RMS disease, we used an exon microarray approach to examine stress-inducible splicing in an RMS cell line (Rh30) and observed striking parallels between stress-responsive alternative splicing and constitutive splicing in RMS tumors.