Integrated Proteomic and Mirna Transcriptional Analysis Reveals the Hepatotoxicity Mechanism of PFNA Exposure in Mice
Total Page:16
File Type:pdf, Size:1020Kb
Integrated proteomic and miRNA transcriptional analysis reveals the hepatotoxicity mechanism of PFNA exposure in mice Jianshe Wang, Shengmin Yan, Wei Zhang, Jiayin Dai * Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, 100101, P.R. China *Corresponding author: Dai JY; Tel. +86-10-64807185; E-mail: [email protected] Supplementary materials and methods PFNA concentration detection in the serum and liver PFNA in liver and serum samples were quantified using HPLC-MS/MS. Briefly, samples with internal standard were extracted with 5 mL of acetonitrile in 15 mL polypropylene (PP) tubes, and all tubes were placed on a mechanical shaker for 20 min followed by centrifugation at 3,000 × g for 10 min. The top layers, which contained PFNA (analytes and internal standard), were transferred into new PP tubes. The extraction procedure was repeated and a final solution of 10 mL of acetonitrile was combined and concentrated to 0.5 mL under nitrogen gas at 40°C. After the addition of 0.5 mL methanol, the final solution was diluted into 10 mL Milli-Q water for SPE cleanup. All samples were then extracted using an Oasis WAX cartridge (Oasis1 HLB; 150 mg, 6 cc; Waters). The cartridge was pre-equilibrated by the addition of a sequence of 4 mL of 0.1% NH 4OH in methanol, 4 mL methanol, and 4 mL water at a rate of 1 drop per second. Samples (11 mL) were then passed through these cartridges at a rate of 1 drop per second. After loading all samples, cartridges were rinsed with 5 mL Milli-Q water and then washed with 4 mL of 25 mM acetate buffer solution (pH 4). Any water remaining in the cartridges was removed by centrifugation at 3,000 rpm for 2 min, and PFNA were eluted by 4 mL of 0.1% NH 4OH in methanol and then concentrated to 1 mL under a stream of nitrogen. The instrumental chromatographic setup consisted of a P680 binary gradient pump, an UltiMate 3000 autosampler, and a Chromeleon 6.70 chromatography workstation (Dionex, USA). Mass spectra were collected using an API 3200 triple quadrupole tandem mass spectrometer, fitted with an electrospray ionization source and operated in negative ionization mode. Quantification using these transitions was performed using Analyst 1.4.1 software. Chromatographic separations were carried out on an Acclaim 120 C18 column (4.6 × 150 mm, 3 µm) (Dionex, USA), with a binary gradient. Methanol (A) and 50 mM ammonium acetate (NH 4Ac) (B) were employed as mobile phases. The flow rate was 1 mL/min and the injection volume was 10 µL. The elution gradient was: 0 to 4 min, from 28 to 5% B linearly; 4 to 7 min, 5% B; 7 to 10 min, 28% B. A calibration curve was prepared from a series of concentrations (0, 10, 50, 100, 500, 1000, 5000, 20000, and 50000 pg/mL), and standard deviations were less than 20%. Blanks and recoveries were assessed following the same procedure as described above with each group of extractions. The blanks were all below the limit of quantifications (LOQs). Quality control is given in the Supplementary Table S1. The concentrations of PFNA in the experimental samples were not corrected for their corresponding recoveries. Supplementary Fig. S1. The PFNA content in mouse liver (A) and serum (B) after treatment. Supplementary Fig. S2. Volcano plots of miRNAs expression in PFNA treatment mice liver tissues. (A) miRNA expression in 1 mg/kg/d dose of PFNA treatment groups compared with control. (B) miRNA expression in 5 mg/kg/d dose of PFNA treatment groups compared with control. On all plots, the X axis represents the log2 ratio of the normalized miRNA signal (mean of the three replicates) of the treatment condition compared with the appropriate control, and the Y axis represents the -log10 P value obtained from the student’s t-test comparing the mean normalized miRNA signal in treatment to that in control. Differentially expressed miRNAs were selected based on fold change at least 2 and student’s t-test ( P < 0.05). Supplementary Fig. S3. Potential binding sites of target genes. (A) Lactate dehydrogenase A (Ldha); (B) fucosyltransferase 8 (Fut8). Supplementary Table S1. Quality control for PFNA concentration detection. Analyte LOQ Instrument Cartridge Procedural Procedural blank blank blank ( n = 3) recovery % ( n = 3) Native Mean Mean Mean Mean (ng/mL) Mean S. D chemical (ng/mL) (ng/mL) (ng/mL) PFNA 0.01 < 0.01 < 0.01 < 0.01 94 8 Surrogate 13 C5-PFNA 0.01 < 0.01 < 0.01 < 0.01 71 1 Supplementary Table S2. Gene-specific primers used in qPCR. Gene symbol Genbank Sequence (5'-3') Amplicon Accession no. size (bp) Me1-F NM_008615 ACAAGGTGACCAAGGGAC Me1-R CAGGGAACACGTAGGAATT 117 Elovl6-F NM_130450 TCAGTTACCTGGGAGTGT Elovl6-R TCTGTATCTGGCAGTTCTT 76 Apoa4-F NM_007468 GGAGCACCTGAAGCCCTAT Apoa4-R ATGCGCTGGATGTATGG 96 Apoc3-F NM_023114 AAGTAGCGTGCAGGAGTC Apoc3-R ATAGCTGGAGTTGGTTGG 150 Gapdh-F NM_008084 CCTTCCGTGTTCCTACCC Gapdh-R GCCTGCTTCACCACCTTC 100 Fut8-F NM_178155 TTCAAGTGGTCGAGCTTCCC 96 Fut8-R GTACAAGTCGATCTGCGAGGT Ldha-F NM_005566 ATGGCAACTCTAAAGGATCAGC 86 Ldha-R CCAACCCCAACAACTGTAATCT Ldlr-F NM_010700 GAGGAACTGGCGGCTGAA Ldlr-R GTGCTGGATGGGGAGGTCT 239 Hmgcr-F NM_008255 GATTCTGGCAGTCAGTGGGAA Hmgcr-R GTTGTAGCCGCCTATGCTCC 214 Fasn-F NM_007988 TGCTCCCAGCTGCAGGC Fasn-R GCCCGGTAGCTCTGGGTGTA 107 Acaca-F NM_133360 TGAAGGGCTACCTCTAATG Acaca-R TCACAACCCAAGAACCAC 182 Hmgcs2-F NM_008256 CCGTGCTCCCTTCTTAG Hmgcs2-R GCTCTTCGTGGGTTCTGT 100 Fdft1-F NM_010191 GCCTTTCTCGTCTATTCTC Fdft1-R GCTGTAGGGTGTTGGTTATG 261 Lpl-F NM_008509 TCTCCTGATGACGCTGAT Lpl-R CTTGGAGTTGCACCTGTAT 288 Gpam-F NM_008149 CCGTGGGCATCTCGTAT Gpam -R CGGACGTAGCCGTAGTTT 139 Acox1-F NM_015729 AGCACTGGTCTCCGTCAT Acox1-R TTCCCTCATCTTCTTCACC 244 Acsl1-F NM_007981 ATTTCAGGTCGCAGATAGA Acsl1-R TGGTACGAGGAGGATTGT 102 Mlycd-F NM_019966 GGCTGTGATGTGGCGTAT Mlycd-R TTGGAGCCCAGGTAGGA 135 Scd1-F NM_009127 GACCTGAAAGCCGAGAA Scd1-R CGTTGAGCACCAGAGTGT 169 Cyp4a10-F NM_010011 CCCAAGTGCCTTTCCTA Cyp4a10-R GCAAACCATACCCAATCC 149 Abbreviations: Me1, malic enzyme 1, NADP(+)-dependent, cytosolic; Elovl6, ELOVL family member 6, elongation of long chain fatty acids; Ldlr, low density lipoprotein receptor; Hmgcr, 3-hydroxy-3-methylglutaryl-Coenzyme A reductase; Fasn, fatty acid synthase; Acaca, acetyl-Coenzyme A carboxylase alpha; Fut8, Fucosyltransferase 8; Ldha, Lactate dehydrogenase; Fdft1, farnesyl diphosphate farnesyl transferase 1; Lpl, lipoprotein lipase; Gpam, glycerol-3-phosphate acyltransferase, mitochondrial; Acsl1, acyl-CoA synthetase long-chain family member 1; Hmgcs2, 3-hydroxy-3-methylglutaryl Coenzyme A synthase 2; Scd1, stearoyl-Coenzyme A desaturase 1; Apoa4, apolipoprotein A-IV; Apoc3, apolipoprotein C-III; Acox1, acyl-Coenzyme A oxidase 1, palmitoyl; Mlycd, malonyl-CoA decarboxylase; Cyp4a10, cytochrome P450, family 4, subfamily a, polypeptide 10 (Cyp4a10); Gapdh, glyceraldehyde-3-phosphate dehydrogenase. Supplementary Table S3. The primers used for 3’ UTR psiCHECK-2 vector preparation. Name Genbank Primer sequence accession no. Ldha-F NM_010699 ACCGCTCGAGACACCCTTCTCGTCTGA L dha-R ATAAGAATGCGGCCGCGACTTAACTGGGAACTCG Fut8-F NM_016893 ACCGCTCGAGAGAAGACGCAGACTAACA Fut8-R ATAAGAATGCGGCCGCGAGGTGGATAAACTGAACT Abbreviations: Ldha, Lactate dehydrogenase A; Fut8, fucosyltransferase 8. The XhoI/NotI sites were underlined. Supplementary Table S4. Significantly altered proteins after 1 mg/kg/d PFNA treatment. Fold change (1 mg/kg/d Unique Accession Description Score Coverage Proteins PFNA treatment P value Peptides group/Ctrl) Q8VCH0 3-ketoacyl-CoA thiolase B, peroxisomal OS=Mus musculus GN=Acaa1b 1264.18 66.51 1 9 2.72 0.0000 PE=2 SV=1 [THIKB_MOUSE] P97412 Lysosomal-trafficking regulator OS=Mus musculus GN=Lyst PE=1 SV=1 0.00 0.32 1 1 2.48 0.0000 [LYST_MOUSE] Q9DBM2 Peroxisomal bifunctional enzyme OS=Mus musculus GN=Ehhadh PE=1 2744.07 70.33 1 61 2.21 0.0000 SV=4 [ECHP_MOUSE] O35728 Cytochrome P450 4A14 OS=Mus musculus GN=Cyp4a14 PE=2 SV=1 54.29 20.32 1 9 1.93 0.0000 [CP4AE_MOUSE] Q9Z0K8 Pantetheinase OS=Mus musculus GN=Vnn1 PE=1 SV=3 33.64 5.27 1 2 1.89 0.0000 [VNN1_MOUSE] P47740 Fatty aldehyde dehydrogenase OS=Mus musculus GN=Aldh3a2 PE=2 412.74 32.02 4 18 1.81 0.0000 SV=2 [AL3A2_MOUSE] Q9EQ06 Estradiol 17-beta-dehydrogenase 11 OS=Mus musculus GN=Hsd17b11 72.92 20.13 1 4 1.78 0.0075 PE=2 SV=1 [DHB11_MOUSE] Q9R0H0 Peroxisomal acyl-coenzyme A oxidase 1 OS=Mus musculus GN=Acox1 588.04 53.4 1 31 1.74 0.0000 PE=1 SV=5 [ACOX1_MOUSE] Q9QYR9 Acyl-coenzyme A thioesterase 2, mitochondrial OS=Mus musculus 489.06 60.71 1 10 1.70 0.0001 GN=Acot2 PE=1 SV=2 [ACOT2_MOUSE] P41216 Long-chain-fatty-acid--CoA ligase 1 OS=Mus musculus GN=Acsl1 PE=1 569.85 46.35 3 35 1.69 0.0000 SV=2 [ACSL1_MOUSE] O55137 Acyl-coenzyme A thioesterase 1 OS=Mus musculus GN=Acot1 PE=1 350.31 62.53 1 6 1.65 0.0002 SV=1 [ACOT1_MOUSE] O88833 Cytochrome P450 4A10 OS=Mus musculus GN=Cyp4a10 PE=2 SV=2 125.23 26.52 1 11 1.64 0.0000 [CP4AA_MOUSE] P47934 Carnitine O-acetyltransferase OS=Mus musculus GN=Crat PE=1 SV=3 49.18 19.33 1 11 1.62 0.0000 [CACP_MOUSE] P12710 Fatty acid-binding protein, liver OS=Mus musculus GN=Fabp1 PE=1 922.06 70.87 1 13 1.61 0.0000 SV=2 [FABPL_MOUSE] Q08857 Platelet glycoprotein 4 OS=Mus musculus GN=Cd36 PE=1 SV=2 0.00 2.54 1 1 1.60 0.0191 [CD36_MOUSE] P10648 Glutathione S-transferase A2 OS=Mus musculus GN=Gsta2 PE=1 SV=2 387.86 37.39 2 6 1.59 0.0004 [GSTA2_MOUSE] P12791 Cytochrome P450 2B10 OS=Mus musculus GN=Cyp2b10 PE=1 SV=1 32.22 6.2 1 3 1.58 0.0001 [CP2BA_MOUSE] Q8JZR0 Long-chain-fatty-acid--CoA ligase 5 OS=Mus musculus GN=Acsl5 PE=2 185.27 33.24 1 21 1.56 0.0000 SV=1 [ACSL5_MOUSE] P55096 ATP-binding cassette sub-family D member