Impact and Molecular Mechanism of Microplastics on Zebrafish in the Presence and Absence of Copper Nanoparticles

The adverse effects of microplastics (MPs) in aquatic environments have attracted increasing attention and posed health risks along with nanomaterials. Therefore, the toxic effects of polystyrene microplastics (PS-MPs) with different particle sizes (0.07, 0.7 and 7 μm) on zebrafish in the presence and absence of copper nanoparticles (Cu-NPs) were evaluated. The acute toxicity of MPs on zebrafish was 7 μm > 0.07 μm > 0.7 μm. Both 0.07 and 7 μm MPs acted on chromosomes and significantly affected cell cycle process by affecting palmitoyl hydrolase activity; while 0.7 μm MPs acted on extracellular space and significantly affected the activity of endopeptidase inhibitor to affect the cholesterol transport. And 0.07 and 7 μm MPs dominantly affected “cell cycle” pathway by inhibiting DNA replication, delaying the progression of S phase and G2/M phase, and affecting the accurate arrangement and separation of chromosomes; while the 0.7 μm MPs activated numerous platelets to aggregate and adhere in damaged parts, enhanced the coagulation function of platelets, and promoted the formation of fibrin clots, thus abnormally activating the “hemostasis” pathway. The presence of Cu-NPs significantly changed the toxicity-related pathways induced by 7 μm MPs from “cell cycle” into “hemostasis,” but not for the smaller-sized MPs (0.07 and 0.7 μm). The combined exposure of Cu-NPs and 7 μm MPs acted on the extracellular region and significantly affected cholesterol transport by affecting the activity of cholesterol transporters. This study provides theoretical insights for the health risks of MPs to aquatic species and even humans in the actual ecosystem.


INTRODUCTION
Microplastics (MPs), generally defined as plastic fragments, fibers and particles smaller than 5 mm in size, have attracted widespread attention due to their ubiquity and persistence (Machado et al., 2018;Hartmann et al., 2019). It is of great concern that aquatic organisms with low feeding selectivity are likely to ingest these plastic particles because of their small size (Geyer et al., 2017), accompanying with adverse effects (Wright et al., 2013;Lehner et al., 2019). It has been demonstrated that MPs can be ingested by various aquatic species, including zooplankton (Cole et al., 2019), crustaceans (Yu et al., 2018), bivalves , fish (Foley et al., 2018), mammals (Sanganyado et al., 2018), etc.; previous reports also have confirmed that MPs uptake can cause a series of adverse effects on aquatic organisms, such as reduced survival rate (Lee et al., 2013), intestinal damage (Gu et al., 2020), loss of energy reserves (Watts et al., 2015), oxidative stress (Choi et al., 2018;Ding et al., 2020), reduced food intake (Cole et al., 2019;Lehner et al., 2019), and growth retardation (Xia et al., 2020). More serious is the transfer of MPs at different trophic levels and that may pose higher associated risks (Bouwmeester et al., 2015). Therefore, it is necessary to effectively identify the behavior of MPs in aquatic environments and their adverse effects on organisms.
It should be mentioned that the aquatic environment in natural state is a complex system containing different types of microparticles (or other pollutants), and thus, the toxicity study of only one type of pollutant may not accurately report its environmental impact. Therefore, more studies have taken into account that MPs can absorb many hazardous contaminants due to their hydrophobic surface and large specific surface area, resulting in combined toxicity (Brennecke et al., 2016;Batel et al., 2018). It has been currently demonstrated that MPs have high adsorption capacity to heavy metals (Brennecke et al., 2016), persistent organic pollutants (POPs) (Batel et al., 2018), antibiotics (Prata et al., 2018) and engineered nanomaterials (NMs) (Schirinzi et al., 2017), and may increase their bioavailability after desorption in organisms (Bakir et al., 2014). Much of the recent concern has concentrated upon the combined effects resulting from MPs and organic pollutants or heavy metals. Yet, little is known about the impact of MPs in the presence of NMs. Both MPs and NMs have accumulated in rivers and coastal areas as many similar processes in the migration and transformation exist like sewage treatment and landfill, thereby providing them with the opportunity to co-exist in aquatic systems (Bouwmeester et al., 2015;Lebreton et al., 2017;Pico et al., 2018). Current studies also have provided evidence that MPs can be used as carriers of silver, TiO 2 and other nanoparticles (NPs) in aquatic environments, which are partly attributed to the monolayer adsorption process driven by electrostatic interaction (Fries et al., 2013;Li et al., 2020). Additionally, the particle toxicity effect of MPs on organisms is similar to that of some NMs, so concern is raised that the combined effects of MPs and NMs (Hueffer et al., 2017). Prior work exhibited that when two or more kinds of microparticles were mixed exposed, the effect of interaction may be synergistic, additive or antagonistic. For example, it was revealed that combined exposure of Au-NPs and MPs at a high concentration was more toxic to Tetraselmis chuii than single exposures (Davarpanah and Guilhermino, 2019). In another study, it has already reported that polystyrene microplastics (PS-MPs) could enhance the accumulation of Cu in zebrafish tissues and exacerbate Cutoxicity as they acted as carriers for Cu ions (Qiao et al., 2019;Bernardes Roda et al., 2020). On the contrary, microsized polyvinyl chloride (mPVC) particles reduced the toxicity of Cu-NPs to the microalgae Skeletonema costatum by adsorption of Cu ions and mutual agglomeration with Cu-NPs (Zhu et al., 2020). However, conflicting results have been found for the combined toxicity of MPs and NMs, and the mechanism of toxicity has not been clearly explained at the molecular level. Hence, it warrants further study the potential toxicity mechanism of MPs and speculate related ecological impacts.
In the present research, PS-MPs and Cu-NPs have been selected because of their extensive commercial and industrial use and high accumulation in aquatic environments (Qiao et al., 2019;Wu et al., 2020). The purposes of this study were (1) to verify the impact of PS-MPs with different particle sizes (0.07, 0.7, and 7 µm) on zebrafish (Danio rerio) in the presence and absence of Cu-NPs; (2) to reveal their molecular toxicity mechanisms. In general, this study investigated the combined toxicity mechanism of MPs and NPs to aquatic vertebrates at the transcriptome level, which can better understand the health risks of MPs to aquatic ecosystems.

Zebrafish Husbandry and Acclimation
The zebrafish (0.19 ± 0.03 g in weight, 27 ± 2 mm in body length) were maintained in tap water with ultraviolet (UV) sterilized and fully aerated (culture water) under the following conditions: temperature 24 ± 1 • C, photoperiod 14: 10 h (light: dark), pH 7.3 ± 0.5, dissolved oxygen 6.7 ± 0.4 mg/L, electrical conductivity 0.287 ± 0.004 mS/cm, water hardness 188 ± 5 mg/L CaCO 3 , and acclimated for 2 weeks before the exposure experiment. The fish were fed twice a day and dead fish were removed in time. All experiments were performed in accordance with the laboratory animal guidelines.

Preparation of PS-MPs and Cu-NPs
The non-functionalized PS microbeads (2.5% w/v, 10 mL) with three different sizes (0.07, 0.7, 7 µm) were purchased from Tianjin BaseLine ChromTech Research Centre (Tianjin, China). The 50 nm Cu-NPs was obtained from Beijing Deke Daojin Science and Technology Co., Ltd. (Beijing, China). The dispersion solutions of PS-MPs (0.07 µm MPs, 0.7 µm MPs, 7 µm MPs) and Cu-NPs were diluted with deionized water and sonicated for 30 min as the stock solutions. The specific experimental solutions were obtained by diluting the stock solutions with culture water. The morphology and size of PS-MPs and Cu-NPs were characterized by transmission electron microscope (TEM) (Hitachi, H-7500, Japan) at an accelerating voltage of 80 kV. Dynamic light scattering (DLS; Zetasizer nano series, Malvern Instruments Ltd., United Kingdom) was used to determine the hydrodynamic diameters and zeta potentials of PS-MPs and Cu-NPs in culture water after 0 and 96 h. The composition of PS-MPs was determined by Fourier transform infrared (FTIR) spectroscopy with a Nicolet 330 spectrophotometer (Thermo Electron Corp., United States).

RNA-Sequencing
Total RNA was extracted from each zebrafish sample by TRIzol R Reagent (Invitrogen, California, United States) according to the manufacturer's instructions, including the following operations: homogenizing samples, phase separation, RNA precipitation, RNA wash and RNA resuspension. Qualified Total RNA was quantified with Qubit RNA Assay Kit (Invitrogen, California, United States) and tested with Agilent 2100 Bioanaylzer (Agilent Technologies, California, United States). The mRNA was enriched using beads with oligod (T) and then mixed with the fragment buffer to fragment the mRNA. Two strands of cDNA were then synthesized successively using random primers and a series of substrates and polymerases. After purifying the second strand with AMPure XP Beads, the ds cDNA repair reaction was performed. Then, the sequencing adapter was ligated to the cDNA fragment for PCR amplification on S1000 TM Thermal Cycler (Bio-Rad Laboratories, California, United States). And the Agilent 2100 Bioanaylzer and Qubit DNA HS Assay Kit (Invitrogen, California, United States) were used for quality control and quantification of the constructed library. The flowcell and the prepared library were cluster generated (bridge PCR amplification) on the cBot system and subsequently sequenced on the Illumina Novaseq 6000. After the results pass the quality control, the differentially expression genes (DEGs) were compared and analyzed according to the following standards by DESeq2 method: fold change =2 or =0.5 and a false discovery rate (FDR) adjusted p-value < 0.05.

Gene Ontology Enrichment Analysis
The functional significance enrichment analysis of DEGs was performed. All DEGs and background genes were mapped to each term in Gene Ontology (GO) database, 1 and the number of genes in each term was calculated. The significantly enriched GO terms were selected according to the p-value < 0.05 (hypergeometric test). 1 http://www.geneontology.org/

Pathway Annotation and Enrichment Analysis
The KEGG database, 2 Reactome database 3 and Panther database 4 were used to conduct pathway enrichment analysis on the selected DEGs. The significantly enriched pathways were selected adjusting for p-value < 0.05 via hypergeometric test.

Statistical Analysis
Data were expressed as means ± standard deviations (SD). Oneway analysis of variance (ANOVA) test was used to evaluate the statistical differences of physical and chemical parameters and experimental results. The statistical analysis of all experimental data was based on SPSS 21.0 software, and the LC 50 values of PS-MPs and Cu-NPs on zebrafish were obtained by the probit analysis. The level of statistical significance was set at p-value < 0.05. The experimental data were mapped using Origin 2019 software.

Microbead Characterization
FT-IR analysis demonstrated the same chemical composition of the three MPs powders (i.e., polystyrene) ( Figure 1A). The representative TEM images exhibited that the three sizes of virgin MPs (0.07, 0.7, and 7 µm) were spherical in shape and uniform in size ( Figure 1B). Supplementary Table 1 showed that there was no significant change in hydrodynamic sizes and zeta potentials of the three MPs in culture water at 0 h and 96 h, suggesting that no agglomeration of MPs occurred over time. These results revealed that the three MPs have similar morphology, same chemical composition and little difference in zeta potentials. Taking into consideration all aforementioned, it is reasonable to analyze the impact of particle size on the acute toxicity of PS-MPs.

Acute Toxicity
The physical and chemical parameters of the exposure solution recorded daily during the experiment were shown in Supplementary File 1. As shown in Figure 2, the toxicity of Cu-NPs to zebrafish was significantly concentration-dependent, and the survival rate was negatively correlated with the concentration. The survival rate of PS-MPs decreased sharply with the increase of concentration at higher concentrations. The 96 h LC 50 values for Cu-NPs, 0.07 µm MPs, 0.7 µm MPs, and 7 µm MPs were 1.951, 93.971, 104.352, 91.379 mg/L, respectively. On the basis of available literature, the impact of MPs on organisms is size-dependent. Generally, in contrast to micron-sized particles, nano-sized particles exhibit the increased toxic effect since nanoplastics have larger specific surface area, stronger ability to penetrate cells, and higher bioavailability (Bouwmeester et al., 2015;Trevisan et al., 2019). For instance, it has been reported that the 0.05 µm MPs showed a more deleterious impact on the survival and reproduction of Tigriopus japonicus and Brachionus koreanus than 0.5 µm and 6 µm MPs (Lee et al., 2013;Jeong et al., 2016). The analogous size-dependent toxicity of 0.1 µm and 2 µm MPs was noted by Rist et al. (2017), that is, 0.1 µm MPs posed greater potential threats like decreasing feeding rate and digestibility, and increasing teratogenicity. Nevertheless, our results indicated that the toxicity of MPs to zebrafish was not monotonously negatively correlated with size. Specifically, the acute toxicity of MPs with three sizes was 7 µm MPs > 0.07 µm MPs > 0.7 µm MPs. Previous studies have also suggested that micro-sized particles can travel into organisms through ingestion or other means to cause biochemical perturbations, and are large enough to induce mechanical damage accompanied by inflammatory responses (Ding et al., 2020). Consistent with our findings, it was found that 5 µm MPs provoked higher SOD and CAT activities than 0.07 µm MPs in Lu et al. (2016)'s oxidative stress analysis on zebrafish. Additionally, in another study, the 5 µm MPs caused more pronounced changes in the expression of metabolism-related genes in the mice offspring with respect to 0.5 µm MPs (Luo et al., 2019). Compared with 0.05 and 0.1 µm MPs, 1 µm MPs induced higher immune responses of hemocyte subpopulations in Mytilus galloprovincialis (Sendra et al., 2020). Perhaps it may be attributed to that organism has a species-specific selection mechanism for micro-sized particles, just as in terms of Sussarellu et al. (2016), oysters preferentially ingested 6 µm MPs rather than smaller MPs, and the selection efficiency for uptake of 5 to 6 µm MPs could reach 100%. These studies demonstrated that the size and feeding selection of test organisms need to be considered when evaluating the adverse effects of MPs with different sizes, which was approved to be the reason for the greater toxicity triggered by microparticles (Le et al., 2018). Accordingly, the toxicity of 0.07 µm MPs may be due to their larger specific surface area and ability to penetrate cells, while the toxicity of 7 µm MPs may contribute to the traumatic inflammatory response triggered by feeding selection of zebrafish. In view of this, more researches are needed to enrich and verify the size-dependent toxicity of MPs.

Differentially Expressed Genes Identification
We chose 1/2 LC 50 as the experimental concentration for transcriptome sequencing since this concentration can be used as the concentration causing the expected genotoxic response (Hansen et al., 2010). A total of 408, 915 and 530 DEGs were responsive to 0.07 µm MPs, 0.7 µm MPs and 7 µm MPs exposures in the absence of Cu-NPs, while 419, 863, and 934 DEGs were identified in the presence of Cu-NPs (versus the control group), respectively (Figures 3A-F and Supplementary  File 2). The majority of the DEGs in 0.07 µm MPs and 7 µm FIGURE 2 | Survival rate curve of zebrafish exposed to Cu-NPs, 0.07 µm MPs, 0.7 µm MPs and 7 µm MPs for 96 h.
MPs were down-regulated, whereas most of the DEGs were upregulated in 0.7 µm MPs treatment. Interestingly, the addition of Cu-NPs had little effect on the gene expression profiles (including the number and proportion of up-regulated and down-regulated DEGs) induced by 0.07 and 0.7 µm MPs treatments, but significantly changed the gene expression profile induced by 7 µm MPs treatment, making it dominated by up-regulated DEGs instead of down-regulated DEGs. Furthermore, clustering of each treatment group indicated that the transcription profiles of six treatment groups were roughly two patterns. 0.07 µm MPs, 7 µm MPs and Cu-NPs + 0.07 µm MPs treatments clustered closely, demonstrating that they cause similar forms of toxicity; however, 0.7 µm MPs, Cu-NPs + 0.7 µm MPs and Cu-NPs + 7 µm MPs treatments shared a cluster and were close to control group, which suggested the effects of the three treatments did not differ much (Supplementary Figure 1).

Gene Ontology Enrichment Analysis of Differentially Expression Genes
In order to further understand the potential functional relevance affected by exposures of PS-MPs with and without Cu-NPs, GO enrichment analysis was performed on DEGs of each treatment. The results were assigned to three categories: biological process, cellular component, and molecular function, and the top 10 most significant terms of each category were shown in Supplementary  Figure 2 and Supplementary File 3. 0.07 and 7 µm MPs acted on the chromosome and significantly affected the cell cycle process by affecting palmitoyl hydrolase activity, and the normal progress of cell cycle was considered to be closely related to of palmitoyl hydrolase activity ; while 0.7 µm MPs acted on the extracellular space and significantly affected the activity of endopeptidase inhibitors to affect the cholesterol transport process, which is determined the directional movement of sterols through or between cells (Clay et al., 2015). Interestingly, the presence of Cu-NPs only changed the GO enrichment results of 7 µm MPs exposure, but did not significantly change the GO enrichment results affected by 0.07 µm MPs and 0.7 µm MPs. That is, the combined exposure of Cu-NPs and 7 µm MPs acted on the extracellular region, and significantly affected the cholesterol transport process by affecting the activity of cholesterol transporters. Worthy to note, the activation of cholesterol transport in three treatments (0.7 µm MPs, Cu-NPs + 0.7 µm MPs, and Cu-NPs + 7 µm MPs) were consistent with those in the juvenile redclaw crayfish exposed to PS-MPs, where the level of total cholesterol in hemolymph increased significantly . (F) Cu-NPs + 7 µm MPs. The x-axis represents the fold change of each gene on a log 2 scale (fold change < 0.5 or fold change > 2). The y-axis represents the significance level of each gene on a -log 10 scale (p-value < 0.05). Each dot represents a gene. Red dots represent up-regulated genes. Green dots represent down-regulated genes.

Pathway-Based Analysis of Differentially Expression Genes
By profiling the functional involvement of DEGs of each group in the KEGG, Reactome and Panther pathways, we can further gain insights into the association between DEGs and biological pathways. Based on p-value, gene number and rich factor, we screened the top 15 pathways that significantly enriched in each treatment (Figure 4 and Supplementary File 4). Interestingly, the results of pathway enrichment analysis have similar sizedependent effects to those of acute toxicity. Specifically, most of the DEGs induced by 0.07 µm MPs and 7 µm MPs exposures were both involved in cell cycle (p = 6.06 × 10 −7 for 0.07 µm MPs, p = 1.37 × 10 −9 for 7 µm MPs) and cell cycle, mitotic (p = 6.83 × 10 −6 for 0.07 µm MPs, p = 4.82 × 10 −9 for 7 µm MPs), while the DEGs induced by 0.7 µm MPs exposure were enriched in common pathway of fibrin clot formation (p = 1.01 × 10 −7 ) and platelet degranulation (p = 1.88 × 10 −7 ), which are both crucial sub-processes of hemostasis. The presence of Cu-NPs only significantly changed the enriched biological pathways of 7 µm MPs exposure, i.e., from cell cycle-related pathways to hemostasis-related pathways, while there was no significant change in the enriched pathways of 0.07 µm MPs and 0.7 µm MPs. This is because the addition of Cu-NPs significantly changed the gene expression profile induced by 7 µm MPs, and changed the cellular components, molecular functions and biological processes associated with DEGs. Therefore, we predicted the molecular toxicity mechanisms caused by PS-MPs exposure in the presence or absence of Cu-NPs by explaining the two most significant enrichment pathways "hemostasis" and "cell cycle" screened from six exposure treatments.

Hemostasis
In our study, the activation of cholesterol transport in three treatments (0.7 µm MPs, Cu-NPs + 0.7 µm MPs, and Cu-NPs + 7 µm MPs) based on GO analysis may affect hemostasis as cholesterol plays a crucial role in the activation and adhesion of platelets and the interaction of cells on the vascular walls (Kannel, 2005). Hemostasis is an intricate and effective physiological defense mechanism initiated in response to injured blood vessels (Figure 5) (Jagadeeswaran et al., 2005;Lasne et al., 2006). It is mainly composed of platelet adhesion, coagulation, and fibrinolysis, wherein coagulation and fibrinolysis are two interconnected reaction networks, modulating the formation and dissolution of fibrin clots, respectively (Jagadeeswaran et al., 2005;Medvedev et al., 2010). When acute vascular injury occurs, the vasoconstriction mechanism predominates, and activated platelets adhere to the exposed collagen to promote coagulation (Austin et al., 2011;Berna-Erro et al., 2013). At first, circulating platelets can identify endothelial lesion areas and adhere to the components of the subendothelium (Kannel, 2005). Here they combine with a variety of thrombogenic substrates and interact with the excitatory platelet agonists to ultimately activate themselves (van der Poll and Parker, 2020). Subsequently, a series of signaling reactions will further enhance the adhesion and procoagulant properties of tethered platelets or nearby circulating platelets. Adhered platelets will degranulate and secrete large amounts of stored effector molecules at the injured site in the form of exocytosis, such as dense granules and alpha granules, to regulate the platelet aggregation process through cell signal transduction (Ge et al., 2010;Golebiewska and Poole, 2014; FIGURE 5 | Hemostasis. The platelet degranulation and formation of fibrin clot are the two most representative sub-processes in hemostasis pathway. Zhang et al., 2018). Our results showed that the abundance of components in platelet dense granules and alpha granules will increase significantly under three exposure conditions (0.7 µm MPs, Cu-NPs + 0.7 µm MPs and Cu-NPs + 7 µm MPs) due to up-regulated expression of genes encoding them (f5, serpina1, hgfa, serping1, aldoaa, igf2b, ahsg2, tor4aa, tfa, etc.), which could activate numerous platelets, aggregate and adhere to the damaged part, greatly enhance the coagulation function of platelets, and provide a positive role for the formation of fibrin clots. Previous research has confirmed our results and shown that MPs and pollutants could cause plasma membrane damage in different organs of fish typified by significant changes in cholesterol, platelets and other biochemical parameters (Hamed et al., 2019); and it was reported that physical stress caused by direct contact of PS-MPs with fibroblasts and red blood cells would lead to cell membrane damage and hemolysis (Choi et al., 2020). Similarly, exposure to NPs such as mesoporous silica NPs could damage the adhesion of platelets to endothelial cells, significantly increase the platelet adhesion events, lead to uncontrolled platelet aggregation actions, and ultimately promote the formation of blood clots in vivo (Kim et al., 2014).
It is currently believed that the products of each reaction in the formation of fibrin clots are enzymes or cofactors, capable of catalyzing or guiding the subsequent reactions, and consequently the entire process is a cascade reaction (Figure 5). Hence, it is reasonable to infer that the changes in the expression levels of some intermediate products have a great impact on clot formation. The common pathway, as a crucial portion of fibrin clot formation, includes a series of activation events from the formation of activated factors X, V, and XIII to the generation of activated thrombin, and finally fibrinogen is cleaved to form a stable cross-linked complex (Johari and Loke, 2012). Factor Xa catalyzes the conversion of prothrombin into activated thrombin, which in turn activates factor V into factor Va. Then factors Va and Xa combine on the membrane surface to form a complex that can rapidly activate a large multitude of thrombin to facilitate the conversion of fibrinogen to fibrin (Krupiczojc et al., 2008). Prior work has documented that fga, fgb, and fgg encoding fibrinogen alpha, beta and gamma chains in zebrafish, respectively, are homologous to human fibrinogen gene sequences (Vo et al., 2013). Moreover, protein C inhibitor (SERPINA5) is also a vital factor during hemostasis that can inhibit some coagulation factors, such as prothrombin and factor Xa (Wagenaar et al., 2010). The three treatments (0.7 µm MPs, Cu-NPs + 0.7 µm MPs and Cu-NPs + 7 µm MPs) resulted in the down-regulation of serpina7, an important component of SERPINA5, which weakened inhibitory effect on coagulation factors and reduced the feedback regulation of the coagulation process, indicating that the coagulation function of the organism has been activated. And the overexpression of regulatory factors X, V, and fibrinogen induced by 0.7 µm MPs significantly increased the abundance of substrates and cofactors in the process of fibrin clot formation, and activated the formation of clots in the blood. In the presence of Cu-NPs, both 0.7 and 7 µm MPs caused the up-regulation of fibrinogen and facilitated the formation of fibrin, suggesting that the blood vessels of zebrafish may be damaged and in selfhealing. To be exact, the overexpression of fibrinogen in the above three treatments is the manifestation of zebrafish's acute stress response because fibrinogen is considered as an acute reactant, whose level increases under inflammation and stress to directly affect blood coagulation and platelet aggregation (Khokhlova et al., 2017). And it was reported that the interaction between NPs and fibrinogen could induce blood coagulation under in vitro conditions (Chen et al., 2011). In addition, the treatment of some NPs such as amorphous silica NPs caused acute lethality and abnormal activation of the formation of fibrin clots in mice; the presence of silica NPs could also accelerate the production of fibrin and facilitate the coagulation system (Kudela et al., 2015;Yoshida et al., 2015). It is therefore plausible that zebrafish were likely to suffer injury after exposure to MPs and Cu-NPs, which would promote the activation and aggregation of numerous platelets due to up-regulated expression of platelet dense particles and α particles; and also lead to up-regulation of fibrinogen and important factors in the process of coagulation, promote the formation of fibrin clots, and ultimately enhance the hemostasis pathway to cope with external stimuli.

Cell Cycle
The DEGs under the three conditions (0.07 µm MPs, 7 µm MPs, and Cu-NPs + 0.07 µm MPs) were classified as cell cycle regulatory genes, checkpoint signal genes and DNA replication genes, which are indispensable for the inspection, regulation and driving of cell cycle. Cell cycle process is mainly composed of DNA replication (S phase) and mitosis (M phase), which are separated by G1 and G2 phases (Figure 6). The minichromosome maintenance protein (MCM) complex is a sixmembered replication helicase that controls DNA replication once in each cell cycle of eukaryotic cells. In the process of DNA replication licensing reaction, it is assembled on chromatin before replication to obtain replication capability (Maiorano, 2000;Tanaka and Araki, 2013). After being activated by S-phase protein kinases, the bound MCM complexes unwind doublestranded DNA at the starting point, recruit DNA polymerases and start DNA synthesis. Due to their crucial roles in the genome replication of proliferating cells, the down-regulation of the expressions of mcm3l, mcm6l, and mcm5 induced by the three treatments would prevent DNA replication by inhibiting the initiation and extension of DNA replication, thereby blocking the cell cycle in S phase. This has also been demonstrated in another study that the decrease of MCM proteins could reduce cell viability by inhibiting DNA replication (Lei, 2005). Furthermore, the timing of mitosis entry and the process of mitosis determine the fate of cells during development. Cyclins drive cell cycle progression by modulating the activity of cyclindependent kinases (CDKs). Cyclin A activates CDK2 and CDK1, promoting the transition from S phase to M phase and the initiation of mitosis, respectively (Malumbres and Barbacid, 2009). And cyclin B combines with CDK1 after the degradation of cyclin A in the prophase of mitosis, then continues to drive cells into mitosis in the form of complexes (Gurunathan et al., 2018). In the present study, the lower expression of components of cyclin A (ccna1 and ccna2) may consequently hinder the normal function of CDK2 and CDK1, and delay the progression of S and G2/M phases. The suppressed expression of cyclin B (ccnb2) in three treatments would block the cell cycle progression in the G2 phase and eventually lead to apoptosis or inhibition of cell proliferation. Kim et al. (2013) further confirmed our results and reported that after exposure to PS-NPs, the cyclin levels were strongly disturbed, and DNA synthesis and cell proliferation were inhibited. By combining numerical simulation and experiment, it was found that the inhibition of cell cycle was due to the combined arrest of G1/S phase and G2/M phase transition. Similarly, it was suggested that Cu-NPs could increase cell death and decrease cell proliferation, and lead to cell cycle arrest by acting on G2 phase; another research has also shown that PS-NPs exerted cytotoxicity by affecting the expression of cell cycle regulators (Thit et al., 2013;Maity et al., 2020).
Additionally, mitotic catastrophe is a type of cell death represented by abnormal mitotic, which is usually associated with spindle formation damage. Polo-like kinase 1 (Plk1) is essential for chromosome homeostasis and bipolar spindle assembly during mitosis. Previous studies have shown that decreased Plk1 activity in cells can induce apoptosis after delayed mitosis; it was also found that inhibition of Plk1 in zebrafish embryos resulted in dysregulation of chromosomes, failure of chromosome separation, formation of unipolar spindles, complete cessation of mitosis, and ultimately widespread cell death (Steegmaier et al., 2007;Jeong et al., 2010). Our results indicated that the down-regulated expression of plk1 induced by 0.07 µm MPs and 7 µm MPs could cause abnormal spindle assembly in the cell and affect the accurate arrangement and division of chromosomes. Besides, BUB1, BubR1, and Mad2 are the main components of the spindle assembly checkpoint (SAC), and they form mitotic checkpoint complex (MCC) targeting cell division cycle 20 homolog (Cdc20), which can control the transition from metaphase to anaphase of mitosis to ensure accurate separation of sister chromatids (Tang et al., 2004;Yang et al., 2012). The down-regulated expression of bub1, mad2l1, and cdc20 in three treatments would weaken the function of SAC and inhibit the SAC signaling pathway, thereby increasing the error rate of chromosome separation and causing the failure of the cell mitosis process and mitotic catastrophe. Our results were corroborated by another study that exposure to MPs could cause cytogenetic abnormalities, manifested by a significant increase in mitotic aberrations (Trifuoggi et al., 2019). Bellingeri et al. (2019) also found that exposure to Cu and PS-MPs resulted in the formation of multinucleated cells and abnormalities in the mitotic process by observing cell morphology. Additionally, the exposure of other NPs (TiO 2 ) inhibited the oogenesis of zebrafish and led to cell death, which was manifested by chromatin condensation, mitochondrial swelling and mitotic catastrophe (Akbulut et al., 2017). Therefore, we believe that the exposure of PS-MPs and Cu-NPs would hinder the transition of different phases in the cell cycle by inhibiting the expression of MCM complex and cyclins, thereby affecting the normal operation of the cell cycle and the failure of cell proliferation; it also inhibited the normal formation of spindles and spindle assembly checkpoints, which would lead to abnormal chromosome arrangement and inaccurate separation of genetic materials, and induce mitotic catastrophe.

CONCLUSION
The toxic effects of PS-MPs varied with particle size and the presence of Cu-NPs. 7 µm MPs showed slightly greater acute toxicity on zebrafish than 0.07 and 0.7 µm MPs. We profiled the gene expression of zebrafish exposed to MPs (in the presence and absence of Cu-NPs), and provided evidence that there were two different toxic modes of exposure to MPs and Cu-NPs on zebrafish. Both 0.07 and 7 µm MPs acted on chromosomes and significantly affected the cell cycle process by affecting palmitoyl hydrolase activity; while 0.7 µm MPs acted on extracellular space and significantly affected the activity of endopeptidase inhibitor to affect the cholesterol transport. And 0.07 and 7 µm MPs dominantly affected the "cell cycle" pathway by inhibiting DNA replication, delaying the progression of S phase and G2/M phase, and affecting the accurate arrangement and separation of chromosomes; while the 0.7 µm MPs activated numerous platelets to aggregate and adhere in damaged parts, enhanced the coagulation function of platelets, and promoted the formation of fibrin clots, thus abnormally activating the "hemostasis" pathway. The presence of Cu-NPs only changed the GO enrichment results of 7 µm MPs exposure, but did not significantly change the GO enrichment results affected by 0.07 µm MPs and 0.7 µm MPs, i.e., the combined exposure of Cu-NPs and 7 µm MPs acted on the extracellular region and significantly affected cholesterol transport by affecting the activity of cholesterol transporters, but did not significantly change the GO enrichment results affected by 0.07 µm MPs and 0.7 µm MPs. And the presence of Cu-NPs significantly changed the toxicity-related pathway induced by 7 µm MPs into "hemostasis, " but not for 0.07 and 0.7 µm MPs. These findings help us to better understand the potential toxicity of the MPs and NPs on aquatic organisms at the transcriptome level, and provide inspiration for the combined effects of MPs and other contaminants.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI (accession: PRJNA758596).

ETHICS STATEMENT
The animal study was reviewed and approved by North China Electric Power University.

AUTHOR CONTRIBUTIONS
NG and JH conceptualized and designed the initial experiments and contributed to the writing, reviewing, and editing. NG performed and analyzed all experiments and conceptualized additional experiments. All authors provided final editing and approval for the manuscript.