Integrating Transcriptomics and Free Fatty Acid Profiling Analysis Reveal Cu Induces Shortened Lifespan and Increased Fat Accumulation and Oxidative Damage in C. elegans

Nowadays, human beings are exposed to Cu in varieties of environmental mediums, resulting in health risks needing urgent attention. Our research found that Cu shortened lifespan and induced aging-related phenotypes of Caenorhabditis elegans (C. elegans). Transcriptomics data showed differential expression genes induced by Cu were mainly involved in regulation of metabolism and longevity, especially in fatty acid metabolism. Quantitative detection of free fatty acid by GC/MS further found that Cu upregulated free fatty acids of C. elegans. A mechanism study confirmed that Cu promoted the fat accumulation in nematodes, which was owing to disorder of fatty acid desaturase and CoA synthetase, endoplasmic reticulum unfolded protein response (UPRER), mitochondrial membrane potential, and unfolded protein response (UPRmt). In addition, Cu activated oxidative stress and prevented DAF-16 translocating into nuclear with a concomitant reduction in the expression of environmental stress-related genes. Taken together, the research suggested that Cu promoted aging and induced fat deposition and oxidative damage.


Introduction
In recent decades, environmental pollution has become a globalization problem [1,2]. Environment pollution would induce water quality problems, especially excessive heavy metals [3,4]. Copper (Cu) has attracted widespread attention due to its utilization for water pipes since 19th century [5]. As an essential trace element, Cu is needed in trace amounts in human and animal [6]. The daily intake of Cu ranges from 0.7 to 3.0 mg/day [7]. Several foods are rich in Cu, including meats, seafood, kidney and liver, mushrooms, and nuts, which is a sufficient supply for daily intake. Thus, the Cu content in the environment and the Cu level in drinking water are particularly important.
Due to favorable mechanical properties, Cu was used extensively in drinking water distribution systems [8,9]. However, metals could not escape the effects of corrosion. The release of Cu particles or ions into the drinking water would cause human health risk [10]. Except for Cu pipes, Cu mining and use of Cu-containing fungicides might affect soil Cu concentration, then directly or indirectly alter the Cu content in animals and plants, finally contribute to excessive Cu intake of human [11,12]. As the Paracelsus principle pointed, "the dose makes the poison" [13].
Cu participates in the synthesis of various proteins and enzymes and is very important in tissue structure and various enzymatic reactions. For example, Cu is involved in the formation of connective tissue, myelin, and formation of melanin pigment in the skin, hair, and eyes [14][15][16]. Cu also plays as an enzymatic cofactor, such as bounding to ceruloplasmin for oxidase activity and system circulation [17], and constitutes cytochrome c oxidase to participate in the respiratory chain in mitochondria [18]. In addition, Cu is a component of cytochrome P450 monooxygenase, dopamine-beta-hydroxylase, and copper, zinc-superoxide dismutase (Cu, Zn-SOD) [18][19][20]. Cu is found in cells and tissues throughout the body, especially abundant in the liver and brain [7]. Therefore, Cu overload induced-malfunctions contributes to the development of hepatic, neurological, and other disorders.
The US Environmental Protection Agency (EPA) had proposed the Optimized corrosion control treatment (OCCT) to carry out the lead and copper rule (LCR) [21]. However, some water systems still exceeded the action level of Cu [22]. Not to mention the hazards of exceeding the upper limit, even the upper limit of the drinking water regulations might not be safety. Thus, this study was focus on two standards for drinking water quality L set up by WHO and 1 mg/L set up by China (GB 5749-2006). The longterm chronic effect of Cu on human health has not been well studied. Actually, chronic effect could also occur when exposed to Cu during a period of time.
Caenorhabditis elegans (C. elegans) is an outstanding experimental system, due to its small size, rapid life cycle, transparency, and well-annotated genome [23]. C. elegans could provide data with aging, reproductive, endocrine, sensory, and neuromuscular systems [24]. In addition, C. elegans has the high conservation of genes and signaling pathways with mammals [25]. Thus, this study is aimed at exploring the effects of Cu exposure on the health, lifespan, and the underlying mechanism of model organisms C. elegans.

Cu Preparation and
Treatment. CuSO 4 ·5H 2 O (purity > 98%, CAS No. 7758-99-8, Sigma-Aldrich, USA) was suitable for C. elegans culture. The CuSO 4 ·5H 2 O was dissolved in K-medium and formulated into a 1 g/L stock solution. As reported, guidelines for drinking water quality issued by the World Health Organization (WHO, fourth edition) and China (GB 5749-2006, China) stipulate the guideline value of Cu were 2.0 mg/L and 1.0 mg/L [27,28]. Thus, Cu concentrations were designed as 0 (K-medium control), 1, and 2 mg/L. Synchronized L1 stage larvae were cultivated for 24-28 h to obtain L4 stage. L4 stage larvae were treated with Cu (control, 1 and 2 mg/L) for 24 h in 35 [29]. The worms were confirmed dead if no movements were observed by touching with a platinum wire at the tail and the head. In brief, synchronized N 2 worms and CF1038 worm were raised to the L4 stage and then were treated with Cu and control on 35 mm NGM plates (with 30 nematodes/plate, total 360 worms/condition, three independent experiments were carried out). After Cu treatment until adult, the C. elegans were transferred to new clean plate and for recognizing the viability of worms under a stereomicroscope every day. Floxuridine (FUdR) was not used to avoid progeny production, and the animals were transferred every day during the fertile period.
2.3.2. Body Growth Index Assay. Body length and width were obtained to evaluate the development. Body growth indexes were detected using Zeiss microscope (Zeiss AX10, Germany) to capture the worms at least 20 for each group.

Pharyngeal
Pumping and Defecation Assay. Pharyngeal pumping and defecation interval are age-related changes [29]. Pumping frequency is determined by counting olfactory bulb pumping over a 30 s timespan [30]. Defecation interval is determined as the period between the two posterior intestinal muscles [30]. At least 20 worms were counted for each group.
All experiments were independently repeated three times.
2.4. Sample Preparation, Library Construction, and RNA Sequencing. L4 stage larvae were treated with Cu (control, 1 and 2 mg/L) until adult. A total of 3 treatment groups were obtained (L4-0, L4-1, and L4-2) each with three biological replicates. C. elegans were collected to extract total RNA using TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA) according to the protocol. Then, sample preparation, library construction, and RNA sequencing were conducted according to previous description [31]. RNA sequencing was performed using Illumina HiSeq TM 6000 by Gene Denovo Biotechnology Co., Ltd. (Guangzhou, China). The experimental details were descripted in the supplementary information (SI).
The raw RNA-seq data of the study have been deposited in NCBI SRA (http://www.ncbi.nlm. http://nil.gov/sra) with the accession number SRP323540.

Analysis of the Subcellular
Localization of the DAF-16. As described, TJ356 (the strain expressing DAF-16 fused to GFP) was used to observe the localization the transcription factor DAF-16 [34]. L4 stage TJ356 worms were firstly transferred to NGM plates without OP50 for starving 16 h to induce nuclear DAF-16 localization. Then, this worms were exposed with 0, 1, and 2 mg/L Cu in NGM plates for 24 h at 20°C. The subcellular DAF-16 localization was analyzed by fluorescence microscopy (Zeiss). The localization of the DAF-16::GFP fusion protein of TJ356 worms were classified into nuclear, cytosolic, and both nuclear and cytosolic three categories. Experiments 2.7.1 to 2.7.3 were independently repeated three times, at least 30 nematodes/group were examined.
2.8. Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR). Gene expression was conducted by RT-qPCR using Gradient PCR instrument (Eppendorf, Germany), StepOne Plus (ABI, USA), and SYBR Green Master Mix (Takara, Tokyo, Japan). Briefly, RNA of C. elegans was reversed transcription to cDNA; then, quantitative PCR was performed according to the protocols. Primer sequences were listed in Table S1. Expression levels of target genes were evaluated by the 2 -ΔΔCt method. Experiments were independently repeated three times.
2.9. Metabolite Extraction and GC-MS Free Fatty Acid Analysis. L4 stage larvae were treated with Cu (control, 1 and 2 mg/L) until adult. A total of 3 treatment groups were obtained (L4-0, L4-1, and L4-2) each with five biological replicates. The N 2 worms were washed for 5 times and taken into the 1.5 mL EP tubes; then, worms were conducted according to previously described methods [35]. The GC-MS analysis was performed using an Agilent 7890B gas chromatograph system (Agilent, California, USA) coupled with an Agilent 5977B mass spectrometer (Agilent, California, USA). In this study, 31 peaks were detected and 31 metabolites were left after relative standard deviation denoising. Then, the missing values were filled up by the median value. The final dataset containing the information of peak number, sample name, and normalized peak area was imported to SIMCA16.0.2 software package (Sartorius Stedim Data Analytics AB, Umea, Sweden) for multivariate analysis. Data was scaled and logarithmic transformed to minimize the impact of both noise and high variance of the variables. After these transformations, PCA (principle component analysis), an unsupervised analysis that reduces the dimension of the data, was carried out to visualize the distribution and the grouping of the samples. 95% confidence interval in the PCA score plot was used as the threshold to identify potential outliers in the dataset. The experimental details could be found in the SI.
2.10. Data Analysis. SPSS 23 (SCR: 002865) and the Graph-Pad Prism 9 (SCR: 002798) software were used to conduct data analysis. The statistics were described as mean ± SD from three independent experiments. The univariate analysis of variance (ANOVA) was performed to analyze the differences among multigroup, followed by Dunnett's test to identify significant difference between treatment and control. Kaplan-Meier survival analysis with log-rank test were used to analyze survival of the worms. ImageJ (SCR: 003070) was used for quantitative of fluorescence intensity.

Results
3.1. The Roles of Cu Exposure in Aging, Growth, and Behavior in C. elegans. The L4 stage worms were treated with Cu to adult then were observed the survival until all worms died. The results showed lifespan and average lifespan of N 2 worms were gradually decreased with Cu exposure (Figures 1(a) and 1(b)). Cu could induce negative effect on growth of C. elegans, manifested by shortening of body length and width (Figures 1(c) and 1(d)). In addition, body movement was detected. The results showed that Cu inhibited the pumping rate and extended defecation interval in a dose-dependent manner (Figures 1(e) and 1(f)), implying the promoting aging effect of Cu in C. elegans.

Identification of Differentially Expressed Genes (DEGs).
RNA-sequencing was performed on three groups with 3 biological replicates each group of C. elegans and finally to construct a total of 9 cDNA libraries. Reads were obtained from the sequencing machines and were further filtered by fastp (version 0.18.0) following protocol [36]. Sample expression violin plot showed maximum, minimum, and median of the 9 sample genes were basically the same ( Figure S1). Correlation of two parallel experiments provided the evaluation of the reliability of experimental results as well  Oxidative Medicine and Cellular Longevity as operational stability. Correlation analysis showed that the correlation of each sample was above 0.97, suggesting a good repeatability among samples ( Figure S2, Figure 2(a)). Cluster analysis was performed on the DEG sets in 9 samples, and genes with similar expression patterns were clustered together. These genes might have similarity functions or participate in common metabolic pathways and signal pathways. Red in the hierarchical clustering diagram indicated high expression, and blue indicated low expression (Figure 2(a)). Cluster analysis was performed on the DEG sets in 9 samples, and genes with similar expression patterns were clustered together. These genes might have similarity functions or participate in common metabolic pathways and signal pathways. Red in the hierarchical clustering diagram indicated high expression, and blue indicated low expression (Figure 2(b)). To further analyze the differential expression induced by Cu, the top 20 DEGs of the 0-VS-1 group and 0-VS-2 group were analyzed. As shown in Figure 2(d) and Table 1, there were 3 upregulated genes and 17 downregulated genes in the 0-VS-1 group, among which fmo-2, gst-33, cht-1, gst-21, and gst-35 were participating in metabolism, six (hsp-70, hsp-16.11, hsp-16.2, hsp-16.41, hsp-16.48, and hsp-16.49) were heat shock protein, col-41, dpy-17, and dsl-2 were related with a component of membrane, cav-1 was involved in reproduction, his-66 was related to genetics, and nhr-17 had metal ion binding activity. As shown in Figure 2(d) and Table 2, there were 2 upregulated genes and 18 downregulated genes in the 0-VS-2 group. Among them, fat-5, fat-7, elo-7, and icl-1 were involved in metabolism, abu-1 and abu-7 participated in endoplasmic reticulum stress, adt-1, bli-1, cdh-7, col-17, col-154, col-157, dpy-4, dpy-8, and dpy-14 were structural constituent of cuticle, adt-1 also could bind with metal ion, fis-1 regulated mitochondrial damage, and grd-2 was related to reproduction. The descriptions and functions of genes were shown in Tables 1 and 2. 3.3. Pathway Analysis of DEGs in Cu-Exposed C. elegans. DEGs were annotated through GO terms for gene function analysis. As showed in SI files Figure S3 A and B, there were 25 GO terms in biological process (BP), 12 GO terms in molecular function (MF), and 13 GO terms in biological process (CC) of 0-VS-1 group. As for the 0-VS-2 group, GO terms in BP, MF, and CC of were 24, 16, and 17. Compared to the control, it was obvious that DEGs and GO terms induced by Cu were much more as the exposure concentration increased.
The KEGG pathway database was used to identify differently enriched pathways of DEGs, which was contributed to explicate genes biological functions of C. elegans response to Cu. As shown in Figures 2(e) and 2(f), pathways mainly related to metabolism and organismal systems, including metabolism of xenobiotics by cytochrome P450, arachidonic acid metabolism, oxidative phosphorylation, glutathione metabolism, longevity regulating pathway, circadian rhythm, thermogenesis, and adipocytokine signaling pathway. Compared with the control group, 1 and 2 mg/L Cu induced analogous types of pathways but with different number of DEGs. Most pathways and DEGs were enriched in metabolism-related pathways, including fatty acid metabolism, fatty acid degradation, polyketide sugar unit biosynthesis, valine, leucine and isoleucine degradation, fatty acid biosynthesis, biosynthesis of amino acids, and cysteine and methionine metabolism. Further, top 10 pathway enrichment analyses showed most identified DEGs involved in longevity regulating pathway in the 0-VS-1 group, while most identified DEGs involved in fatty acid metabolism and degradation in the 0-VS-2 group ( Figure S3 C and D). Practically, the twelfth pathway of the 0-VS-2 group was also longevity regulating pathway (Table S2).
3.5. Cu Exposure Promoted Fat Deposition in C. elegans. To explore the effect of Cu exposure on fat deposition in C. elegans, Oil red O staining and Nile Red O staining were conducted. As shown, triglyceride (TG) or fats were obviously enhanced by Cu (Figure 4(a) and 4(b)). Lipid droplets in worms were increased in a dose-dependent manner (Figures 4(c) and 4(d)). To further explain the lipidregulating effect of Cu, lipid metabolism-related genes of C. elegans were analyzed from RNA-seq results and were confirmed by RT-qPCR (Figures 4(e) and 4(f)). The Δ9-desaturase fat-5 and fat-7 gene expressions were decreased at 1 mg/ L concentration and increased at 2 mg/L concentration; genes acs-1 and acs-2 encoding fatty acid β-oxidation enzymes were gradually inhibited by Cu. The polyunsaturated fatty acid (PUFA) synthesis pathway and schematic of lipid metabolism of C. elegans were shown in Figure 4(g), FAT-5 worked on C16:0 to synthesize C16:1n7, and FAT-6 and FAT-7 acted on C18:0 to synthesize C18:1n9 [37]. ACS enzymes, as fatty acid CoA synthetase family, were required for activating fatty acids to acyl CoA [37]. Our GC/MS results found saturated fatty acids (SFA), monounsaturated fatty acids (MUFA), and PUFA were all altered (Table S4 and 5, Figure 4(h)). To further evaluate the activity of fat-5 and fat-7, product C16:1n7/ precursor C16:0 ratio and product C18:1n9/precursor 5 Oxidative Medicine and Cellular Longevity       Oxidative Medicine and Cellular Longevity C18:0 ratio was compared. The C16:1n7/C16:0 ratio was increased significantly by Cu, which also verified the enhanced activity of FAT-5 (Figure 4(i)).
Mitochondrial dysfunction impaired lipid homeostasis. Thus, mitochondrial membrane potential (MMP, Δψm) and mitochondrial unfolded protein response (UPR mt ) were detected. The result of JC-1 probes staining showed that J-   Table S2 showed mRNA expression and biological functions of the above 4 genes.

DAF-16 Mediated Cu-Induced Aging and Oxidative
Stress. IIS pathway was related to regulation of longevity, aging, and fat metabolism. Results of RNA-seq and RT-qPCR revealed that daf-2 and age-1 were promoted (Figure 6(a)) and downstream genes of daf-16 were inhibited by Cu; however, the expression of daf-16 was not changed significantly. Therefore, we speculated whether Cu prevented the daf-16 trans-location to the nucleus through insulin signaling to regulate fat storage and lifespan. DAF-16 was mainly in the cytosolic of worm cells, for inducing nuclear localization of DAF-16, TJ356 worms were starved for 16 h. Then, Cu exposure altered the DAF-16 localization. The result showed (Figures 6(b) and 6(c)) that DAF-16 translocated out of the nucleus, indicating DAF-16 translocation to the nucleus was inhibited by Cu. DAF-16 directly regulated the expression of a variety of genes, such as sod-3, gst-4, ctl-2, and mtl-1, which are involved in oxidative stress [39]. The H2DCFDA probe staining results showed that ROS level was increased with Cu treatment (

12
Oxidative Medicine and Cellular Longevity and 6(g)). To further elucidate the regulation of DAF-16 on lifespan, daf-16 (mu86) strain was used to evaluate survival after Cu treatment. However, Cu did not alter the lifespan of the mutant strain significantly compared with untreated worm (Figure 6(h)). In addition, RNA-seq results revealed that mRNA level of sod-1, sod-2, sod-3, mlt-1, and mlt-2 were inhibited by Cu (Figure 6(i)), the specific FC value and description of the genes were shown in Table S2 of SI.

Discussion
The early-life period is an important window of susceptibility to environmental exposures. Different developmental stages would be critical windows of exposure. Developmental exposure to toxicants may result in an acceleration of age-related decline in function. A study showed different exposure paradigms (i.e., L1 stage treatment or L4 stage treatment) induced differentially expressed genes (DEGs) [40]. Therefore, we conducted two different exposure paradigms of Cu (L1 stage and L4 stage of C. elegans) during RNA sequencing process. Interestingly, Cu exposure promotes aging in C. elegans at different stages, but the involving different mechanisms. Our previous results show that Cu induced 2332 DEGs (567 up-and 1765 downregulated genes) of the 1 mg/L group and 2449 DEGs (724 up-and 1725 downregulated genes) of 2 mg/L group in the L1 stage exposure paradigm [31]. In this study, Cu induced 759 DEGs (280 up-and 479 downregulated genes) of the 1 mg/L group and 1342 DEGs (155 up-and 1187 downregulated genes) of the 2 mg/L group in L4 stage exposure paradigm. Further, the top 20 DEGs were analyzed. In the L1 stage exposure paradigm, the top 20 DEGs were involved in structural constituent of cuticle or integral component of membrane, immune response, lipid transport, epigenetic regulation, and gene expression regulation [31]. However, in the L4 stage exposure paradigm, the top 20 DEGs were participating in component of membrane, metabolism, heat shock stress, endoplasmic reticulum stress, and mitochondrial damage. Further, identification of the pathway showed two concentrations in the L1 stage exposure paradigm both influenced on the lysosome, longevity regulating, and peroxisome pathway [31]. In the L4 stage exposure paradigm, two concentrations both induced longevity regulating pathway and adherens junction. Interest-ingly, Cu mainly affected the lipid metabolism with increasing concentrations. Therefore, the direction of the mechanism verification of the two Cu exposure paradigms is also different. Cu exposure widely affects longevity regulation pathways in the L1 stage exposure paradigm, thereby promoting aging [31]. In the L4 stage exposure paradigm, we focus on fat metabolism, mitochondrial dysfunction, endoplasmic reticulum emergency, and oxidative stress according to clues of RNA sequencing.
This study is aimed at exploring the health damage after L4 stage C. elegans exposed to Cu for 24 h. Transcriptome sequencing and quantitative detection of free fatty acid of Cu-treated nematodes were conducted to analyze the molecular mechanism of Cu-induced abnormal behavior and function of nematodes.

Cu Promoted Aging and Behavioral Disorders of C.
elegans. Aging composites the important risk factor for chronic degenerative diseases, such as cancer and cardiovascular [41]. It was very important to explore the occurrence and mechanism of aging. C. elegans was an excellent model to study with aging due to its physical decline strikingly similar to humans [42]. Typical markers of aging-relative phenotype in C. elegans included the locomotion decline, feeding ability (pharyngeal pumping frequency), defecation ability (defecation interval), morphological changes, and the aging-associated genes [29,43].
First of all, lifespan was shortened by two environmentally relevant doses of Cu. Body length and width were used to evaluate the development of C. elegans, the result of declined body length and body width suggested that C. elegans was sensitivity to growth inhibition by Cu. Pharyngeal pumping was controlled by corpus and isthmus muscles, which was regulated by neurons and muscle-intrinsic  19 Oxidative Medicine and Cellular Longevity molecules [44]. A healthy well-fed adults C. elegans pumped 200-300 times per minute [29,44]. Defecation was a kind of rhythmic exercise that consists a three-part motor program coordinated by neural activity [29]. A healthy well-fed young adult wild-type C. elegans defecated once every 45 seconds [29]. Both pharyngeal pumping frequency and defecation interval were correlated with the lifespan.
In the study, the pharyngeal pumping rate was decreased by Cu, the defecation interval was prolonged. The above research results preliminary indicated that Cu decreased To illustrate the regulatory mechanism behind these recessions induced by Cu, RNA-seq was performed. Transcriptomics data showed cuticle structure related genes included cuticle collagen bli-1; collagen col-17, col-41, col-154 and col-157; and collagen structural dpy-4, dpy-8, dpy-14, and dpy-17 were downregulated; besides, other 54 collagen-encoded genes were also suppressed including col-8, col-14, and col-39. col-8, col-14, and col-39 belonged to C. elegans cuticle collagen cysteine families [45]. These genes encoded external structure, that is, cuticle, and were crucial for development and survival of nematode [46,47], which well explained that Cu inhibited their development. Behavior and lifespan were affected by several pathways, including the inositol phosphate pathway, lipid metabolism pathways, and energy metabolism pathways. Therefore, we would conduct an in-depth discussion in conjunction with subsequent experiments.

Cu Might Regulate
Fat Metabolism via UPR mt and UPR ER . The effect on metabolism of C. elegans seemed to be more pronounced with Cu exposure increases. Analysis of DGEs, GO function annotation, and KEGG pathway were all pointing that abnormal lipid metabolism was produced by Cu. DGEs were involved in fatty acid degradation, metabolism, and biosynthesis. Lipids executed many essential functions, for example, they were efficient energy storage molecules, formed a double-layer membrane structure, and played as signaling molecules [48]. However, incorrectly lipid homeostasis would develop dysregulated neurological and developmental consequences [37]; moreover, lipid overload could reduce lifespan of C. elegans [37,49].
Lipid metabolism was set of process that comprised lipid synthesis at endoplasmic reticulum, storage and transport of lipid droplets, β-oxidation at mitochondria and peroxisomes, and lipid hydrolysis and recycling at lysosomes [50][51][52]. Firstly, we found 2 mg/L Cu upregulated delta (9)-fatty-acid desaturase (Δ9 desaturase) fat-5 and fat-7. Δ9 desaturase was considered to be lipogenic enzymes and active in endoplasmic reticulum membrane both in mammals and C. elegans [53,54]. FAT-5 worked on C16:0 to produce C16:1n-7, and FAT-7 mainly acted on C18:0 to produce C18:1n-9 [55]. Thus, we conducted targeted FFA metabolomic analysis in C. elegans and found FFA was increased after Cu treated. Correspondingly, C16:0, C16:1n-7, C18:0, and C18:1n-9 were all increased. Oil red staining and Nile red staining further confirmed that Cu promoted fat accumulation in N2 nematodes. Further, we found Cu induced ER stress of C. elegans, and multiple proteins associated with endoplasmic reticulum were affected by Cu. The ER was to facilitate the folding and maturation of protein molecules [56,57]. AS reported that misfolded proteins would be activated to cope with ER stress [56]. Thus, Cu stimulated ER stress on the one hand to activate expression of fat-5 and fat-7, and on the other hand, to increase the expression of HSP-4::GFP sequestrated to the misfolded protein.
The end products of lipid mobilization, FFAs, went to further degradation via mitochondrial β-oxidation, which was related to longevity mechanisms [58]. In the C. elegans intestine, asc-1 and acs-2 encoded acyl-CoA synthetases (ACSs), which could activate FFAs β-oxidation at mitochondrial [59,60]. However, asc-1 and acs-2 were downregulated 23 Oxidative Medicine and Cellular Longevity in our study. To explore the mechanism behind, it was found that MMP (Δψm) and UPR mt were both reduced in Cutreated N 2 worms. The loss of MMP (Δψm) was a hallmark of mitochondrial dysfunction [61]. Mitochondria was the core organelles for energy metabolism, programmed cell death, and metabolic pathways involving lipids and amino acids [62]. Environmental stress or dysfunction led to accumulation of unfolded protein would trigger UPR mt , which was a compensatory mechanism to maintain mitochondrial homeostasis [63,64]. Surprisingly, we found that a series of ribosomal proteins located in the mitochondrial were dysfunctional, and some of them were firstly downexpressed then overexpressed as Cu exposure increasing (SI Figure S4). The influence of the mitochondrial ribosomal proteins (MRPs) on lifespan was conserved from mammals to C. elegans that triggers the UPR mt by downexpression of MRPs was a life-extending mechanism [65]. On the contrary, suppressed UPR mt contributed to aging.
Phenotype is the result of a combination of inheritance and gene regulation and involves a myriad of changes at the organ, tissue, and cellular levels. Inactivation of C. elegans desaturase and elongase family members, encoded by fat and elo genes, respectively, causes imbalances in fatty acid composition and is associated with metabolic, physiological, and behavioral phenotypes [66]. These include altered total fat levels, growth retardation, slowed movement, reduction in body size, germ cell maintenance and reproductive defects, aberrations in rhythmic behavior, defects in sensory signaling, defects in neurotransmission, and reduced adult lifespan [37]. Similarly, lipid storage phenotypes in this study was a comprehensive result of Cu-induced oxidative damage, mitochondrial dysfunction, and endoplasmic reticulum stress, even though the gene expression changes of lipid metabolism genes vary. Thus, the exact mechanism of lipid storage induced by Cu was still need discussed in depth.  [67][68][69]. HSPs were regulated by DAF-16 and participating in the longevity regulating pathway [70,71]. DAF-16/FOXO transcription factor, controlled by DAF-2/IGFR activating phosphoinositide 3-kinase AGE-1/PI3K kinase cascade, governed most of the functions of the insulin/IGF-1 signaling (IIS) pathway [72]. Activated DAF-2 receptor could recruit and activate AGE-1, resulting in phosphorylation of and cytoplasmic sequestration of the DAF-16. The C. elegans IIS pathway extensively regulated metabolism, growth, development, behavior, and longevity [73,74].

Cu Prevented
To further confirm the effect of Cu on IIS pathway, genes (daf-2, age-1, and daf-16) based on the RNA-seq data were selected for RT-qPCR validation. It showed that daf-2 and age-1 were activated, and image analysis of TJ356 suggested DAF-16 nuclear localization was decreased. Next, glutathione S-transferase, gst-21 and gst-33 were found downregulated, and gst-35 was upregulated, which were a response to oxidative stress [75]. Accordingly, we further confirmed that Cu caused accumulation of ROS accompanied by decreased expression of antioxidant SOD-3. As shown in Table S2, C. elegans had five SOD genes that sod-1 and sod-5 encoded cytoplasmic CuZnSOD, sod-4 encoded extracellular CuZnSOD, and sod-2 and sod-3 encoded mitochondrial MnSOD. sod-5 was stimulant increased while another sod-1 to sod-4 were inhibited, especially sod-1 and sod-2. Two other genes that could not be ignored, mtl-1 and mtl-2, that encoded metallothionein were involved in stress response to heavy metal in C. elegans [76]. As reported, ROS served as pleiotropic physiological signaling agents which could also stimulate DAF-2/IGFR signaling [77,78]. Additionally, sod-3, mtl-1, and mtl-2 were direct transcript targets of DAF-16 and were related to longevity [79,80].
Taken together, Cu exposure activated ROS; then, oxidative damage signals directly or indirectly sequestrated DAF-16 in the cytoplasm, contributing to the DAf-16 fail of regulation target gene transcription (such as sod-3, MTLs and HSPs), finally decreased stress resistance of C. elegans.

Conclusion
In summary, our study implicated Cu accelerated C. elegans aging accompanied by developmental and behavioral suppression. Further mechanism study showed Cu induced fat deposition and oxidation damage. Cu, on the one hand, triggered UPRER which led to the upregulation of Δ9 desaturase FAT-5 and FAT-7 to increase lipid synthesis, and on the other hand, inhibited MMP and induced an imbalance of MRP, thereby inhibiting UPR mt and FFAs β-oxidation. In addition, Cu activated ROS and at least in part prevented antioxidant capacity through transcription factor DAF-16/ FOXO.
The study illustrated the potential threats of Cu to humans at part degree and provided some theoretical basis for preventing or mitigating aging induced by Cu.

Data Availability
All data included in this study were available upon request by contact with the corresponding author.  Figure S1: violin plot for all samples. Figure S2: correlation analysis for all samples. The correlation of each sample was above 0.97. Figure S3: GO enrichment analysis and top 10 pathways enrichment. Figure S4: mitochondrial ribosomal proteins were transcriptionally repressed by Cu. Table  S1: the gene primer used in the study. Table S2: expression of differentially expressed genes (DEGs) in copper-exposed C. elegans. Table S3: top 12 statistics of pathways enrichment. Table S4: differential free fatty acid induced by 1 mg/L copper. Table S5: differential free fatty acid induced by 2 mg/L copper. (Supplementary Materials)