Relationship between nucleosome positioning and progesterone-induced alternative splicing in breast cancer cells

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

Alternative splicing is coupled to gene expression in a subset of variably expressed genes

Numerous factors regulate alternative splicing of human genes at a co-transcriptional level. However, how alternative splicing depends on the regulation of gene expression is poorly understood. We leveraged data from the Genotype-Tissue Expression (GTEx) project to show a significant association of gene expression and splicing for 6874 (4.9%) of 141,043 exons in 1106 (13.3%) of 8314 genes with substantially variable expression in ten GTEx tissues. About half of these exons demonstrate higher inclusion with higher gene expression, and half demonstrate higher exclusion, with the observed direction of coupling being highly consistent across different tissues and in external datasets. The exons differ with respect to sequence characteristics, enriched sequence motifs, RNA polymerase II binding, and inferred transcription rate of downstream introns. The exons were enriched for hundreds of isoform-specific Gene Ontology annotations, suggesting that the coupling of expression and alternative splicing described here may provide an important gene regulatory mechanism that might be used in a variety of biological contexts. In particular, higher inclusion exons could play an important role during cell division.

Histone Marks-Dependent Effect on Alternative Splicing: New Perspectives for Targeted Splicing Modulation in Cancer?

Alternative splicing (AS) is a tightly regulated mechanism that generates the complex human proteome from a small number of genes. Cis-regulatory RNA motifs in exons and introns control AS, recruiting positive and negative trans-acting splicing regulators. At a higher level, chromatin affects splicing events. Growing evidence indicates that the popular histone code hypothesis can be extended to RNA-level processes, such as AS. In addition to nucleosome positioning, which can generate transcriptional barriers to shape the final splicing outcome, histone post-translational modifications can contribute to the detailed regulation of single exon inclusion/exclusion. A histone-based system can identify alternatively spliced chromatin stretches, affecting RNAPII elongation locally or recruiting splicing components via adaptor complexes. In tumor cells, several mechanisms trigger misregulated AS events and produce cancer-associated transcripts. On a genome-wide level, aberrant AS can be the consequence of dysfunctional epigenetic splicing code, including altered enrichment in histone post-translational modifications. This review describes the main findings related to the effect of histone modifications and variants on splicing outcome and how a dysfunctional epigenetic splicing code triggers aberrant AS in cancer. In addition, it highlights recent advances in programmable DNA-targeting technologies and their possible application for AS targeted epigenetic modulation.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Histone marks regulate the epithelial-to-mesenchymal transition via alternative splicing

Histone modifications impact final splicing decisions. However, there is little evidence of the driving role of these marks in inducing cell-specific splicing changes. Using CRISPR epigenome editing tools, we show in an epithelial-to-mesenchymal cell reprogramming system (epithelial-to-mesenchymal transition [EMT]) that a single change in H3K27ac or H3K27me3 levels right at the alternatively spliced exon is necessary and sufficient to induce a splicing change capable of recapitulating important aspects of EMT, such as cell motility and invasiveness. This histone-mark-dependent splicing effect is highly dynamic and mediated by direct recruitment of the splicing regulator PTB to its RNA binding sites. These results support a role for H3K27 marks in inducing a change in the cell's phenotype via regulation of alternative splicing. We propose the dynamic nature of chromatin as a rapid and reversible mechanism to coordinate the splicing response to cell-extrinsic cues, such as induction of EMT.

Integrative epigenomic profiling reveal AP-1 is a key regulator in intrahepatich cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is a malignant tumor with poor prognosis while its mechanisms of pathogenesis remain elusive. In this study, we performed systemic epigenomic and transcriptomic profiling via MNase-seq, ChIP-seq and RNA-seq in normal cholangiocyte and ICC cell lines. We showed that active histone modifications (H3K4me3, H3K4me1 and H3K27ac) were less enriched on cancer-related genes in ICC cell lines compared to control. The region of different histone modification patterns is enrichment in sites of AP-1 motif. Subsequent analysis showed that ICC had different nucleosome occupancy in differentially expressed genes compared to a normal cell line. Furthermore, we found that AP-1 plays a key role in ICC and regulates ICC-related genes through its AP-1 binding site. This study is the first report showing the global features of histone modification, transcript, and nucleosome profiles in ICC; we also show that the transcription factor AP-1 might be a key target gene in ICC.

Progesterone receptor antagonists reverse stem cell expansion and the paracrine effectors of progesterone action in the mouse mammary gland

Background

The ovarian hormones estrogen and progesterone (EP) are implicated in breast cancer causation. A specific consequence of progesterone exposure is the expansion of the mammary stem cell (MSC) and luminal progenitor (LP) compartments. We hypothesized that this effect, and its molecular facilitators, could be abrogated by progesterone receptor (PR) antagonists administered in a mouse model.

Methods

Ovariectomized FVB mice were randomized to 14 days of treatment: sham, EP, EP + telapristone (EP + TPA), EP + mifepristone (EP + MFP). Mice were then sacrificed, mammary glands harvested, and mammary epithelial cell lineages separated by flow cytometry using cell surface markers. RNA from each lineage was sequenced and differential gene expression was analyzed using DESeq. Quantitative PCR was performed to confirm the candidate genes discovered in RNA seq. ANOVA with Tukey post hoc analysis was performed to compare relative expression. Alternative splicing events were examined using the rMATs multivariate analysis tool.

Results

Significant increases in the MSC and luminal mature (LM) cell fractions were observed following EP treatment compared to control (p < 0.01 and p < 0.05, respectively), whereas the LP fraction was significantly reduced (p < 0.05). These hormone-induced effects were reversed upon exposure to TPA and MFP (p < 0.01 for both). Gene Ontology analysis of RNA-sequencing data showed EP-induced enrichment of several pathways, with the largest effect on Wnt signaling in MSC, significantly repressed by PR inhibitors. In LP cells, significant induction of Wnt4 and Rankl, and Wnt pathway intermediates Lrp2 and Axin2 (confirmed by qRTPCR) were reversed by TPA and MFP (p < 0.0001). Downstream signaling intermediates of these pathways (Lrp5, Mmp7) showed similar effects. Expression of markers of epithelial-mesenchymal transition (Cdh1, Cdh3) and the induction of EMT regulators (Zeb1, Zeb2, Gli3, Snai1, and Ptch2) were significantly responsive to progesterone. EP treatment was associated with large-scale alternative splicing events, with an enrichment of motifs associated with Srsf, Esrp, and Rbfox families. Exon skipping was observed in Cdh1, Enah, and Brd4.

Conclusions

PR inhibition reverses known tumorigenic pathways in the mammary gland and suppresses a previously unknown effect of progesterone on RNA splicing events. In total, our results strengthen the case for reconsideration of PR inhibitors for breast cancer prevention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01455-2.

Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana

Alternative splicing (AS) is a major gene regulatory mechanism in plants. Recent evidence supports co-transcriptional splicing in plants, hence the chromatin state can impact AS. However, how dynamic changes in the chromatin state such as nucleosome occupancy influence the cold-induced AS remains poorly understood. Here, we generated transcriptome (RNA-Seq) and nucleosome positioning (MNase-Seq) data for Arabidopsis thaliana to understand how nucleosome positioning modulates cold-induced AS. Our results show that characteristic nucleosome occupancy levels are strongly associated with the type and abundance of various AS events under normal and cold temperature conditions in Arabidopsis. Intriguingly, exitrons, alternatively spliced internal regions of protein-coding exons, exhibit distinctive nucleosome positioning pattern compared to other alternatively spliced regions. Likewise, nucleosome patterns differ between exitrons and retained introns, pointing to their distinct regulation. Collectively, our data show that characteristic changes in nucleosome positioning modulate AS in plants in response to cold.

Splicing-associated chromatin signatures: a combinatorial and position-dependent role for histone marks in splicing definition

Alternative splicing relies on the combinatorial recruitment of splicing regulators to specific RNA binding sites. Chromatin has been shown to impact this recruitment. However, a limited number of histone marks have been studied at a global level. In this work, a machine learning approach, applied to extensive epigenomics datasets in human H1 embryonic stem cells and IMR90 foetal fibroblasts, has identified eleven chromatin modifications that differentially mark alternatively spliced exons depending on the level of exon inclusion. These marks act in a combinatorial and position-dependent way, creating characteristic splicing-associated chromatin signatures (SACS). In support of a functional role for SACS in coordinating splicing regulation, changes in the alternative splicing of SACS-marked exons between ten different cell lines correlate with changes in SACS enrichment levels and recruitment of the splicing regulators predicted by RNA motif search analysis. We propose the dynamic nature of chromatin modifications as a mechanism to rapidly fine-tune alternative splicing when necessary.

Chromatin is known to regulate splicing by modulating recruitment of splicing factors. Using machine learning approaches, the authors have underlined a chromatin code for alternative splicing regulation that is conserved amongst cell lines.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

dSreg: a Bayesian model to integrate changes in splicing and RNA-binding protein activity

MOTIVATION

Alternative splicing (AS) is an important mechanism in the generation of transcript diversity across mammals. AS patterns are dynamically regulated during development and in response to environmental changes. Defects or perturbations in its regulation may lead to cancer or neurological disorders, among other pathological conditions. The regulatory mechanisms controlling AS in a given biological context are typically inferred using a two-step framework: differential AS analysis followed by enrichment methods. These strategies require setting rather arbitrary thresholds and are prone to error propagation along the analysis.

RESULTS

To overcome these limitations, we propose dSreg, a Bayesian model that integrates RNA-seq with data from regulatory features, e.g. binding sites of RNA-binding proteins. dSreg identifies the key underlying regulators controlling AS changes and quantifies their activity while simultaneously estimating the changes in exon inclusion rates. dSreg increased both the sensitivity and the specificity of the identified AS changes in simulated data, even at low read coverage. dSreg also showed improved performance when analyzing a collection of knock-down RNA-binding proteins' experiments from ENCODE, as opposed to traditional enrichment methods, such as over-representation analysis and gene set enrichment analysis. dSreg opens the possibility to integrate a large amount of readily available RNA-seq datasets at low coverage for AS analysis and allows more cost-effective RNA-seq experiments.

AVAILABILITY AND IMPLEMENTATION

dSreg was implemented in python using stan and is freely available to the community at https://bitbucket.org/cmartiga/dsreg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ARGLU1 is a transcriptional coactivator and splicing regulator important for stress hormone signaling and development

Stress hormones bind and activate the glucocorticoid receptor (GR) in many tissues including the brain. We identified arginine and glutamate rich 1 (ARGLU1) in a screen for new modulators of glucocorticoid signaling in the CNS. Biochemical studies show that the glutamate rich C-terminus of ARGLU1 coactivates multiple nuclear receptors including the glucocorticoid receptor (GR) and the arginine rich N-terminus interacts with splicing factors and binds to RNA. RNA-seq of neural cells depleted of ARGLU1 revealed significant changes in the expression and alternative splicing of distinct genes involved in neurogenesis. Loss of ARGLU1 is embryonic lethal in mice, and knockdown in zebrafish causes neurodevelopmental and heart defects. Treatment with dexamethasone, a GR activator, also induces changes in the pattern of alternatively spliced genes, many of which were lost when ARGLU1 was absent. Importantly, the genes found to be alternatively spliced in response to glucocorticoid treatment were distinct from those under transcriptional control by GR, suggesting an additional mechanism of glucocorticoid action is present in neural cells. Our results thus show that ARGLU1 is a novel factor for embryonic development that modulates basal transcription and alternative splicing in neural cells with consequences for glucocorticoid signaling.

The Intricate Interplay between Epigenetic Events, Alternative Splicing and Noncoding RNA Deregulation in Colorectal Cancer

Colorectal cancer (CRC) results from a transformation of colonic epithelial cells into adenocarcinoma cells due to genetic and epigenetic instabilities, alongside remodelling of the surrounding stromal tumour microenvironment. Epithelial-specific epigenetic variations escorting this process include chromatin remodelling, histone modifications and aberrant DNA methylation, which influence gene expression, alternative splicing and function of non-coding RNA. In this review, we first highlight epigenetic modulators, modifiers and mediators in CRC, then we elaborate on causes and consequences of epigenetic alterations in CRC pathogenesis alongside an appraisal of the complex feedback mechanisms realized through alternative splicing and non-coding RNA regulation. An emphasis in our review is put on how this intricate network of epigenetic and post-transcriptional gene regulation evolves during the initiation, progression and metastasis formation in CRC.

Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

How Are Short Exons Flanked by Long Introns Defined and Committed to Splicing?

The splice sites (SSs) delimiting an intron are brought together in the earliest step of spliceosome assembly yet it remains obscure how SS pairing occurs, especially when introns are thousands of nucleotides long. Splicing occurs in vivo in mammals within minutes regardless of intron length, implying that SS pairing can instantly follow transcription. Also, factors required for SS pairing, such as the U1 small nuclear ribonucleoprotein (snRNP) and U2AF65, associate with RNA polymerase II (RNAPII), while nucleosomes preferentially bind exonic sequences and associate with U2 snRNP. Based on recent publications, we assume that the 5' SS-bound U1 snRNP can remain tethered to RNAPII until complete synthesis of the downstream intron and exon. An additional U1 snRNP then binds the downstream 5' SS, whereas the RNAPII-associated U2AF65 binds the upstream 3' SS to facilitate SS pairing along with exon definition. Next, the nucleosome-associated U2 snRNP binds the branch site to advance splicing complex assembly. This may explain how RNAPII and chromatin are involved in spliceosome assembly and how introns lengthened during evolution with a relatively minimal compromise in splicing.

Interconnections Between RNA-Processing Pathways Revealed by a Sequencing-Based Genetic Screen for Pre-mRNA Splicing Mutants in Fission Yeast

Pre-mRNA splicing is an essential component of eukaryotic gene expression and is highly conserved from unicellular yeasts to humans. Here, we present the development and implementation of a sequencing-based reverse genetic screen designed to identify nonessential genes that impact pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, an organism that shares many of the complex features of splicing in higher eukaryotes. Using a custom-designed barcoding scheme, we simultaneously queried ∼3000 mutant strains for their impact on the splicing efficiency of two endogenous pre-mRNAs. A total of 61 nonessential genes were identified whose deletions resulted in defects in pre-mRNA splicing; enriched among these were factors encoding known or predicted components of the spliceosome. Included among the candidates identified here are genes with well-characterized roles in other RNA-processing pathways, including heterochromatic silencing and 3ʹ end processing. Splicing-sensitive microarrays confirm broad splicing defects for many of these factors, revealing novel functional connections between these pathways.

Chromatin structure and pre-mRNA processing work together

Chromatin is the natural context for transcription elongation. However, the elongating RNA polymerase II (RNAPII) is forced to pause by the positioned nucleosomes present in gene bodies. Here, we briefly discuss the current results suggesting that those pauses could serve as a mechanism to coordinate transcription elongation with pre-mRNA processing. Further, histone post-translational modifications have been found to regulate the recruitment of factors involved in pre-mRNA processing. This view highlights the important regulatory role of the chromatin context in the whole process of the mature mRNA synthesis.

Coupling of RNA Polymerase II Transcription Elongation with Pre-mRNA Splicing

Pre-mRNA maturation frequently occurs at the same time and place as transcription by RNA polymerase II. The co-transcriptionality of mRNA processing has permitted the evolution of mechanisms that functionally couple transcription elongation with diverse events that occur on the nascent RNA. This review summarizes the current understanding of the relationship between transcriptional elongation through a chromatin template and co-transcriptional splicing including alternative splicing decisions that affect the expression of most human genes.

A deformation energy-based model for predicting nucleosome dyads and occupancy

Nucleosome plays an essential role in various cellular processes, such as DNA replication, recombination, and transcription. Hence, it is important to decode the mechanism of nucleosome positioning and identify nucleosome positions in the genome. In this paper, we present a model for predicting nucleosome positioning based on DNA deformation, in which both bending and shearing of the nucleosomal DNA are considered. The model successfully predicted the dyad positions of nucleosomes assembled in vitro and the in vitro map of nucleosomes in Saccharomyces cerevisiae. Applying the model to Caenorhabditis elegans and Drosophila melanogaster, we achieved satisfactory results. Our data also show that shearing energy of nucleosomal DNA outperforms bending energy in nucleosome occupancy prediction and the ability to predict nucleosome dyad positions is attributed to bending energy that is associated with rotational positioning of nucleosomes.

Defective histone supply causes changes in RNA polymerase II elongation rate and cotranscriptional pre-mRNA splicing

RNA polymerase II (RNAPII) transcription elongation is a highly regulated process that greatly influences mRNA levels as well as pre-mRNA splicing. Despite many studies in vitro, how chromatin modulates RNAPII elongation in vivo is still unclear. Here, we show that a decrease in the level of available canonical histones leads to more accessible chromatin with decreased levels of canonical histones and variants H2A.X and H2A.Z and increased levels of H3.3. With this altered chromatin structure, the RNAPII elongation rate increases, and the kinetics of pre-mRNA splicing is delayed with respect to RNAPII elongation. Consistent with the kinetic model of cotranscriptional splicing, the rapid RNAPII elongation induced by histone depletion promotes the skipping of variable exons in the CD44 gene. Indeed, a slowly elongating mutant of RNAPII was able to rescue this defect, indicating that the defective splicing induced by histone depletion is a direct consequence of the increased elongation rate. In addition, genome-wide analysis evidenced that histone reduction promotes widespread alterations in pre-mRNA processing, including intron retention and changes in alternative splicing. Our data demonstrate that pre-mRNA splicing may be regulated by chromatin structure through the modulation of the RNAPII elongation rate.

Promoter-like epigenetic signatures in exons displaying cell type-specific splicing

BACKGROUND

Pre-mRNA splicing occurs mainly co-transcriptionally, and both nucleosome density and histone modifications have been proposed to play a role in splice site recognition and regulation. However, the extent and mechanisms behind this interplay remain poorly understood.

RESULTS

We use transcriptomic and epigenomic data generated by the ENCODE project to investigate the association between chromatin structure and alternative splicing. We find a strong and significant positive association between H3K9ac, H3K27ac, H3K4me3, epigenetic marks characteristic of active promoters, and exon inclusion in a small but well-defined class of exons, representing approximately 4 % of all regulated exons. These exons are systematically maintained at comparatively low levels of inclusion across cell types, but their inclusion is significantly enhanced in particular cell types when in physical proximity to active promoters.

CONCLUSION

Histone modifications and other chromatin features that activate transcription can be co-opted to participate in the regulation of the splicing of exons that are in physical proximity to promoter regions.

Connecting the dots: chromatin and alternative splicing in EMT

Nature has devised sophisticated cellular machinery to process mRNA transcripts produced by RNA Polymerase II, removing intronic regions and connecting exons together, to produce mature RNAs. This process, known as splicing, is very closely linked to transcription. Alternative splicing, or the ability to produce different combinations of exons that are spliced together from the same genomic template, is a fundamental means of regulating protein complexity. Similar to transcription, both constitutive and alternative splicing can be regulated by chromatin and its associated factors in response to various signal transduction pathways activated by external stimuli. This regulation can vary between different cell types, and interference with these pathways can lead to changes in splicing, often resulting in aberrant cellular states and disease. The epithelial to mesenchymal transition (EMT), which leads to cancer metastasis, is influenced by alternative splicing events of chromatin remodelers and epigenetic factors such as DNA methylation and non-coding RNAs. In this review, we will discuss the role of epigenetic factors including chromatin, chromatin remodelers, DNA methyltransferases, and microRNAs in the context of alternative splicing, and discuss their potential involvement in alternative splicing during the EMT process.

Role of six single nucleotide polymorphisms, risk factors in coronary disease, in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.