Temporal Modulation of Differential Alternative Splicing in HaCaT Human Keratinocyte Cell Line Chronically Exposed to Arsenic for up to 28 Wk

Abstract

BACKGROUND: Chronic arsenic exposure via drinking water is associated with an increased risk of developing cancer and noncancer chronic diseases. Pre-mRNAs are often subject to alternative splicing, generating mRNA isoforms encoding functionally distinct protein isoforms. The resulting imbalance in isoform species can result in pathogenic changes in critical signaling pathways. Alternative splicing as a mechanism of arsenic-induced toxicity and carcinogenicity is understudied. OBJECTIVE: This study aimed to accurately profile differential alternative splicing events in human keratinocytes induced by chronic arsenic exposure that might play a role in carcinogenesis. METHODS: Independent quadruplicate cultures of immortalized human keratinocytes (HaCaT) were maintained continuously for 28 wk with 0 or 100 nM sodium arsenite. RNA-sequencing (RNA-Seq) was performed with poly(A) RNA isolated from cells harvested at 7, 19, and 28 wk with subsequent replicate multivariate analysis of transcript splicing (rMATS) analysis to detect and quantify differential alternative splicing events. Reverse transcriptase-polymerase chain reaction (RT-PCR) for selected alternative splicing events was performed to validate RNA-Seq predictions. Functional enrichment was performed by gene ontology (GO) analysis of the differential alternative splicing event data set at each time point. RESULTS: At least 600 differential alternative splicing events were detected at each time point tested, comprising all the five main types of alternative splicing and occurring in both open reading frames (ORFs) and untranslated regions (UTRs). Based on functional relevance ELK4, SHC1, and XRRA1 were selected for validation of predicted alternative splicing events at 7 wk by RT-PCR. Densitometric analysis of RT-PCR data corroborated the rMATS predicted alternative splicing for all three events. Protein expression validation of the selected alternative splicing events was challenging given that very few isoform-specific antibodies are available. GO analysis demonstrated that the enriched terms in differential alternatively spliced mRNAs changed dynamically with the time of exposure. Notably, RNA metabolism and splicing regulation pathways were enriched at the 7-wk time point, when the greatest number of differentially alternatively spliced mRNAs are detected. Our preliminary proteomic analysis demonstrated that the expression of the canonical isoforms of the splice regulators DDX42, RMB25, and SRRM2 were induced upon chronic arsenic exposure, corroborating the splicing predictions. DISCUSSION: These results using cultures of HaCaT cells suggest that arsenic exposure disrupted an alternative splice factor network and induced time-dependent genome-wide differential alternative splicing that likely contributed to the changing proteomic landscape in arsenic-induced carcinogenesis. However, significant challenges remain in corroborating alternative splicing data at the proteomic level.