RNA sequencing provides insights into the toxicogenomic response of ZF4 cells to methyl methanesulfonate

Abstract

Whole genome transcriptomic studies are powerful for characterizing the molecular mechanisms underlying the physiological effects of chemicals, and are informative for environmental health risk assessment. Alkylating agents are an abundant class of chemicals that can damage DNA in the environment, and are used for anticancer treatments. Currently, little is known regarding the molecular mechanisms of toxic alkylating agents in zebrafish cell lines. In this study, RNA-sequencing was used to investigate the transcriptomic responses of zebrafish ZF4 cells following exposure to the model genotoxicant methyl methanesulfonate (MMS). The half-maximal inhibitory concentration (IC50 ) of MMS was 639.16 ± 61.8 µm, and apoptosis was induced within 24 h of exposure. RNA sequencing identified 3601 differentially expressed genes (DEGs) that were upregulated and 3037 that were downregulated. Gene ontology enrichment analysis revealed that most DEGs belonged to synthesis and metabolism categories. RNA-associated processes were the most upregulated, while cell cycle and adhesion were the most repressed processes, and neuron-related processes were the most downregulated developmental process. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis identified DNA damage repair, cell cycle, apoptosis and spliceosome as overrepresented terms. Six types of alternative splicing were detected. In total, 1156 alternative splicing DEGs were specifically expressed following MMS treatment, many of which belonged to metabolism and catabolic process categories. Cluster analysis of orthologs was able to extrapolate toxicotranscriptomic data between zebrafish and yeast. These results provide insight into the genome-wide response of ZF4 cells following exposure to MMS, and this knowledge will inform future toxicogenomic data analysis and environmental health risk assessment.