Mechanism Analysis of Metabolic Fatty Liver on Largemouth Bass (Micropterus salmoides) Based on Integrated Lipidomics and Proteomics

Metabolic fatty liver disease caused by high-starch diet restricted the intensive and sustainable development of carnivorous fish such as largemouth bass. In this study, the combination liver proteomic and lipidomic approach was employed to investigate the key signaling pathways and identify the critical biomarkers of fatty liver in largemouth bass. Joint analysis of the correlated differential proteins and lipids revealed nine common metabolic pathways; it was determined that FABP1 were significantly up-regulated in terms of transporting more triglycerides into the liver, while ABCA1 and VDAC1 proteins were significantly down-regulated in terms of preventing the transport of lipids and cholesterol out of the liver, leading to triglyceride accumulation in hepatocyte, eventually resulting in metabolic fatty liver disease. The results indicate that FABP1, ABCA1 and VDAC1 could be potential biomarkers for treating metabolic fatty liver disease of largemouth bass.


Introduction
Largemouth bass (Micropterus salmoides), originally a popular sport fish in North America, is now reared in many other countries worldwide and also is an important freshwater carnivorous species in China for its high market value as a food [1][2][3][4]. At present, more than 80% of cultured largemouth bass are fed on chilled fish [5]. However, the high costs of chilled fish as well as the potential role as disease vectors hinder the healthy development of largemouth bass aquaculture worldwide [2].
Previous studies discovered that fishes (especially carnivorous ones such as largemouth bass) generally have a low ability to use carbohydrate in compound feed [6][7][8].
Starch is the cheapest energy source in practical diet ingredients, and it contributes to good binders for aquatic feeds [9]. Traditional puffed floating feeds usually require at least 20% starch to obtain adequate swelling and floating properties. However, researchers demonstrated that excess dietary carbohydrate could cause fat deposition, glucose and lipid metabolism disorder, hepatopathy, apoptosis and finally result in liver damage [8,[10][11][12][13]. Previous studies suggested that the starch content in the diet should be lower than 10% to ensure the liver health of largemouth bass [14,15]. Therefore, high carbohydrate is recognized as the primary factor that induced metabolic liver disease (MLD) in largemouth bass [8,9]. MLD of largemouth bass was the first limiting factor for its sustainable development worldwide [7]. Thus, it is crucial to understand the processes and progression of MLD, as well as the key factors in largemouth bass.
Proteomic technologies could be used for identifying and quantifying target proteins to comprehensively explore the etiology of liver disease, due to their high-resolution and highthroughput advantages [16]. Currently, it has become a powerful tool for understanding pathology of liver disease [17]. By identifying groups of changed proteins, we could gain insight into potential pathways and regulatory networks that might contribute to the development of MLD. Many studies have found that proteomics could be used for biomarker identification for early diagnosis in MLD [18][19][20]. To date, no study has ever comprehensively investigated the protein expression changes in fishes of MLD. Lipids are essential components for maintaining various homeostasis, physiological and cellular processes in animals [21]. Lipid metabolism disorders would lead to many major health problems, such as obesity and nonalcoholic fatty liver disease (NAFLD) [22]. Lipidomics, based on UPLC-MS/MS technology, is a novel omics strategy for investigating lipid metabolism and identifying lipid biomarkers [23,24], which enables large-scale and comprehensive studies of lipids [25]. However, there is a lack of studies that systematically assess the protein and lipid changes in the diseased liver of fish for the discovery of the biomarkers of metabolic liver disease.
In this study, we constructed a fatty liver model of largemouth bass by feeding high levels of dietary starch (16.2%) to investigate the complex biological processes and pathogenesis of MLD with proteomic and lipidomic techniques. The mass spectrometry technique was used for proteome-level quantification analyses. Some characteristic changes to different proteins and lipid metabolites are documented for the first time in this study and provide a comprehensive multi-omics framework for MLD biomarker discovery in largemouth bass.

Materials and Methods
During the experiment period, all fishes were maintained in compliance with the Laboratory Animal Welfare Guidelines of China (Decree No. 2 of Ministry of Science and Technology, issued in 2021).

Experimental Diets
Two experimental diets were formulated to be isonitrogenous and isoenergetic. A basal diet was used as the control containing 10.8% starch (named the Normal group), whereas another diet was prepared with 16.2% starch (named the MLD group). Each diet was extruded into 2 mm diameter pellets using a twin-screwed extruder (EXT50A, Yang gong Machine, Beijing, China). The diet formulation and analyzed chemical composition are shown in Table 1.

Experimental Fish, Feeding and Sampling
Largemouth bass were obtained from the Tangshan Aquafarm (Tangshan, Hebei, China). The experiment was conducted in the indoor circulating water system at the National Aquatic Feed Safety Evaluation Base (Nan Kou, Beijing, China). Prior to the formal experiment, the fish were acclimatized for 4 weeks by being fed the Normal diet. The water temperature was maintained at 21-25 • C, pH = 7.2-8.0, dissolved oxygen (DO) >6.0 mg/L, ammonia nitrogen content <0.3 mg/L and NO 2 − < 0.1 mg/L. The fish were fed until apparent satiation twice a day (8:00 am and 4:00 pm) for 10 weeks. Food intake was measured daily.
At the end of the growth trial, fishes from each group were randomly selected and anaesthetized with chlorobutanol (300 mg/mL) after 24 h starvation. The specific growth rate (SGR), final body weight (FBW) and feed conversion ratio (FCR) were detected by weighing the fish at the end of the 10 weeks. The body weight, body length, liver and viscera weight of the fish were recorded individually to calculate condition factor (CF), hepatosomatic index (HSI), viscerosomatic index (VSI) and hepatic liquid (HL), respectively. Blood samples were drawn from the caudal part of the sedated fish using anticoagulant syringes with 2% NaF and 4% potassium oxalate and centrifuged at 4000× g rpm for 10 min at 4 • C to obtain serum. Two parts of liver samples near to the bile duct were collected for histology examination (fixed in 4% paraformaldehyde solution) and biochemical criterion analysis (frozen in liquid nitrogen). All samples (except for histological samples) were stored at −80 • C until analysis.

Histopathological Examination of the Liver Tissue
After 24 h of fixation, all liver samples were dehydrated by the standard procedures, and the samples were embedded in paraffin and cut to 6 µm sections. Liver sections were stained following the protocols of hematoxylin and eosin (H&E) staining and observed with light microscopy (Leica DM2500, Leica, Solms, Germany).

Protein Extraction
The sample (3 mm × 3 mm) was grinded with liquid nitrogen into cell powder and then transferred to a 5 mL centrifuge tube. After that, four volumes of lysis buffer (8 M urea, 1% Triton-100, 10 mM dithiothreitol and 1% Protease Inhibitor Cocktail) were added to the cell powder, followed by sonication three times on ice using a high-intensity ultrasonic processor. The remaining debris was removed by centrifugation at 20,000× g at 4 • C for 10 min. Finally, the protein was precipitated with cold 20% TCA for 2 h at −20 • C. After centrifugation at 12,000× g 4 • C for 10 min, the supernatant was discarded. The remaining precipitate was washed with cold acetone three times. The protein was redissolved in 8 M urea and the protein concentration was determined with BCA kit according to the manufacturer's instructions.

Trypsin Digestion
For digestion, the protein solution was reduced with 5 mM dithiothreitol for 30 min at 56 • C and alkylated with 11 mM iodoacetamide for 15 min at room temperature in darkness. The protein sample was then diluted to urea concentration less than 2 M. Finally, trypsin was added at 1:50 trypsin-to-protein mass ratio for the first digestion overnight and 1:100 trypsin-to-protein mass ratio for a second 4 h digestion.

TMT Labeling
After trypsin digestion, peptide was desalted with Strata X C18 SPE column (Phenomenex) and vacuum dried. Peptide was reconstituted in 0.5 M TEAB and processed according to the manufacturer's protocol for TMT kit. Briefly, one unit of TMT reagent was thawed and reconstituted in acetonitrile. The peptide mixtures were then incubated for 2 h at room temperature and pooled, desalted and dried by vacuum centrifugation.

LC-MS/MS Analysis
The hydrolysated peptides were separated by HPLC connected with a reverse capillary column. The specific steps and related conditions were as follows [19,26,27]: (1) sample desalination: 0.1% TFA cleaning for 5 min (20 µL/min); (2) gradient elution: 2-35% ACN gradient elution for 45 min (350 nL/min), 80% ACN elution for 15 min; (3) the eluent composition of HPLC: A, 0.1% formic acid solution (v/v). B, 0.1% formic acid acetonitrile solution (v/v); (4) the peptides were subjected to NSI source followed by tandem mass spectrometry (MS/MS) in Q Exactive TM Plus (Thermo Fisher Scientific Co., LTD, Shanghai, China) coupled online with the UPLC. The electrospray voltage applied was 2.0 kV. The m/z scan range was 350 to 1000 for full scan, and intact peptides were detected in the Orbitrap at a resolution of 35,000. Peptides were then selected for MS/MS using NCE setting as 27 and the fragments were detected in the Orbitrap at a resolution of 17,500. A data-independent procedure alternated between one MS scan followed by 20 MS/MS scans. Automatic gain control (AGC) was set at 3E6 for full MS and 1E5 for MS/MS. The maximum IT was set at 20 ms for full MS and auto for MS/MS. The isolation window for MS/MS was set at 2.0 m/z.

Protein Identification
Protein identifications were performed by using the ProteinPliot software (AB Sciex). For protein quantitation, proteins were required to contain at least two unique peptides.
To demonstrate the reproducibility of the replicates, protein differential expressions between various biological replicates were compared. Then p-value was calculated by using the two-sample and two-tail t-test. When p-value < 0.05, the change of differential expression exceeding 1.3 was regarded as the threshold for significant up-regulation, and that less than 0.77 was regarded as the threshold for significant down-regulation.

Lipid Extraction
For the lipidomics analysis, the total lipids were extracted from the liver in the Normal and MLD groups. Liver tissues (20 mg) were homogenized in a 2 mL centrifuge tube with 1 mL extracting solution (methyl tertiary butyl ether: methyl alcohol = 3:1, v/v, including internal standard mixture). After homogenization, 200 µL deionized water was added and then centrifuged for 10 min at 12,000× g r/min, 4 • C. The supernatant solution was dried under nitrogen and redissolved in 200 µL mobile phase B for further analysis.  Table 2 lists the lipid types and information of internal and external standards.

Statistical Analysis
Independent t-test was performed by data of growth performance. All data are presented as the mean value ± standard error of the mean (S.E.M). p < 0.05 was considered significantly different. The correlation and linear regression analysis of biological repetitions were carried out using GraphPad Prism 8.0 (Graph-Pad Software Inc., San Diego, CA, USA).
Functional annotations of the proteins and genetic ontology were conducted using the Blast2GO program. The related pathways of each protein and bioinformation mining were analyzed by searching the KEGG pathway database.
R software was used for Principal Component Analysis (PCA), cluster analysis, grouping Principal Component Analysis, differential metabolite screening and KEGG functional annotation of different samples.

High-Starch Diet (HSD) Induced Fatty Liver of Largemouth
Fish liver sections were examined after H&E staining and Sirius red staining for collagen. Typical phenotypes of normal liver and fatty liver are shown in Figure 1. The results showed that HSD led to enlarged lipid droplets and vacuolated cells.

Statistical Analysis
Independent t-test was performed by data of growth performance. All data are presented as the mean value ± standard error of the mean (S.E.M). p < 0.05 was considered significantly different. The correlation and linear regression analysis of biological repetitions were carried out using GraphPad Prism 8.0 (Graph-Pad Software Inc. San Diego, CA, USA).
Functional annotations of the proteins and genetic ontology were conducted using the Blast2GO program. The related pathways of each protein and bioinformation mining were analyzed by searching the KEGG pathway database.
R software was used for Principal Component Analysis (PCA), cluster analysis, grouping Principal Component Analysis, differential metabolite screening and KEGG functional annotation of different samples.

High-Starch Diet (HSD) Induced Fatty Liver of Largemouth
Fish liver sections were examined after H&E staining and Sirius red staining for collagen. Typical phenotypes of normal liver and fatty liver are shown in Figure 1. The results showed that HSD led to enlarged lipid droplets and vacuolated cells. Both Normal and MLD groups showed high survival (99%) and there was no significant difference between the two groups. Compared with the Normal group, fish in the MLD group showed significantly higher HSI and HL (p < 0.05) but significantly lower SGR and FBW (p < 0.05, Table 3).
LDL-C indicators (in plasma and liver, Tables 4 and 5) and the hepatic function indicators (AKP, AST and ALT) increased significantly in the MLD group (p < 0.05, Table 4). Higher ROS level and lower CAT content were recorded in the MLD group (Table 6).

Proteomics Analysis
Liver proteome profiles of largemouth bass with fatty liver (n = 8) and normal liver were compared (n = 8) (Table S1). This revealed an obvious shift in the proteome composition, For the KEGG pathway analysis, differential proteins were primarily associated with enrichment in pathways related to the fat digestion and absorption, gastric acid secretion, ABC transport, cholesterol metabolism, glycolysis/gluconeogenesis, PPAR and MAPK ( Figure 2E,F). The differential pathways that may be related to fatty liver metabolism and the important differentially expressed proteins in these pathways are listed in Table 7. To investigate the functions of the identified and quantified proteins, we annotated their functions and features, gene ontology (GO), subcellular localization, GO-based en-richment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, protein domain and involvement in a protein complex of the characteristics we annotated [26].
For the biological process category, the differential proteins between the Normal and MLD groups were related to cellular processes, biological regulation processes, metabolic processes, response to stimulus processes, multicellular organismal processes, localization processes and development processes. For the cellular component category, these differential proteins were mainly enriched in the organelle, membrane, macromolecular complex and membrane-enclosed lumen. For the molecular function category, fatty liver induced by high-starch diet was primarily enriched in functional clusters, such as catalytic activity, transporter activity and structural molecule activity. (Figure 2B,C).
Subcellular structure localizations of differentially expressed proteins were predicted and classified, the results showed that most of these differential proteins were located in the nucleus, cell membrane, cytoplasmic cells and external matrix, and a few were located in mitochondria and endoplasmic reticulum. Forty were located in the cytoplasm of the 99 differential expressed proteins, with 23 in the extracellular matrix, 14 in the nucleus, 7 in the cell membrane, 6 in the cytoplasm and nucleus, 5 in the mitochondria, 2 in the endoplasmic reticulum and 2 in other cellular structures ( Figure 2D).
For the KEGG pathway analysis, differential proteins were primarily associated with enrichment in pathways related to the fat digestion and absorption, gastric acid secretion, ABC transport, cholesterol metabolism, glycolysis/gluconeogenesis, PPAR and MAPK ( Figure 2E,F). The differential pathways that may be related to fatty liver metabolism and the important differentially expressed proteins in these pathways are listed in Table 7.   (Table 8 and Figure 3). PCA analysis showed discrimination between the MLD group and the Normal group ( Figure 4A), indicating that there was difference in the lipid composition between the two groups. OPLS-DA analysis was applied to screen out the differential metabolites between the two groups. Using the criteria of VIP ≥ 1 and FC ≥ 1.5 or ≤ 0.67, 164 significantly differential lipid species were identified in the Normal versus the MLD group ( Figure 4B). As shown in the heatmap (Figure 4C), the differential lipid species were mainly enriched in TC, PC, DPPE, Cer, SM, PI and DG classes, and all of them were upregulated in the MLD group. Notably, the top 20 differential upregulated lipids were mainly TGs, including TG (51:0), TG (51:1), TG (54:1), TG (50:0), TG (52:0), TG (49:1) and TG (54:0). However, the three downregulated lipids were mainly LPC (20:2), Cer (d34:1) and PE (P-34:2) ( Table 9). Table 8. Major lipid composition in the liver of largemouth bass.
A total of 29 differential metabolic pathways were found with KEGG pathway enrichment analysis. These metabolic pathways mainly involved fat digestion and absorption, cholesterol metabolism, vitamin digestion and absorption, sphingolipid metabolism and insulin resistance, all of which are related to lipid metabolism (Table 10).
The above correlated differential lipids and differential proteins were simultaneously mapped to the KEGG pathway, and were enriched into nine common metabolic pathways, as shown in Table 11. Four KEGG pathways with the highest concentration were selected for subsequent analysis, including thermogenesis (ko04714), fat digestion and absorption (ko04975), cholesterol metabolism (ko04979) and metabolic pathways (ko01100).

Discussion
Carnivorous fish have a low starch utilization rate, and high digestible starch intake leads to accumulation of liver glycogen and persistent high blood sugar, finally causing fish liver disease [28,29]. Generally, the level of digestible carbohydrate above 10% reduced the growth performance of largemouth [9,15,30,31], induced glycogen and lipid accumulation and dysfunction of antioxidant capabilities, leading to MLD in largemouth [9,32]. In this study, compared with the Normal group (starch content: 10.8%), the MLD group (starch content: 16.2%) decreased the growth performance of largemouth bass and induced liver lipid accumulation, liver function injury, oxidative stress and higher hepatosomatic index, which is similar to the non-alcoholic fatty liver symptom. Therefore, we used HSD to construct the fatty liver phenotype in largemouth bass successfully.

Proteomic Analysis
The differentially expressed proteins induced by a high-starch diet of largemouth bass were mainly involved in endoplasmic reticulum protein processing, fat digestion and absorption, cholesterol metabolism, phosphatidylinositol signal system and insulin signal transduction, which affected the occurrence and development of fatty liver through the comprehensive action of multiple metabolic pathways.
Endoplasmic reticulum (ER) is an important site involved in the regulation of substance transport, metabolism and protein synthesis in eukaryotic cells [33]; ER stress caused by various reasons is closely related to the occurrence and development of diseases [34]. The evidence implied that ER played a role in the development of steatosis and nonalcoholic steatohepatitis. ER stress occurred in liver and adipose tissue in patients with non-alcoholic fatty liver disease [35,36]. RRBP1 (reticulum ribosome-binding protein 1) is one of the important proteins involved in the ER unfolded protein reaction, and abnormal up-regulated expression of RRBP1 has been found in the study of human cancer-related diseases [37][38][39]. Xiong found that the expression of RRBP1 in liver cancer tissues was significantly higher than that in normal liver tissues, and the expression of RRBP1 in liver cancer tissues with metastasis was significantly higher than that in liver cancer tissues without metastasis, indicating that RRBP1 is closely related to liver cancer metastasis in the state of liver injury [40]. In this study, RRBP1 was found to be highly expressed in liver in the MLD group, indicating that it may have an important relationship with liver injury and the development of MLD. In addition, TRAPα (translocon-associated protein subunit alpha) is a typical glycosylated membrane protein which plays an important role in the process of signal modification recognition and the transport of new peptide chains. It synergizes with AMPK (adenosine monophosphate activated protein kinase) and ROS pathway to cause resistance to oxidative stress-induced apoptosis [41]. TRAPα also plays an important role in maintaining cell homeostasis and may become a target drug for intervention in related diseases [42]. The low expression of TRAPα may be due to the abnormality in the modification and transport of the peptide chain, resulting in liver damage, which affects the normal progress of subsequent life activities and leads to the occurrence and development of fatty liver. At present, the TRAPα sequence of fish is relatively conserved, but the research on its related functions is still in the blank stage and needs to be further studied.
The pathway of fat digestion and absorption has the important function of maintaining the homeostasis of lipid metabolism. In this study, we found the differential expression of three important proteins FABP1, MTTP (microsomal triglyceride transfer protein) and ABCA1 in its pathway. FABP1 is expressed mainly in the liver and could combine with longchain fatty acid to regulate lipid absorption and fatty acid metabolism in the cytoplasm [43]. FABP1 activity changed significantly in the occurrence and development of fatty liver, liver cirrhosis, liver cancer and other liver diseases [44][45][46]. In this study, HSD induced fatty liver of largemouth bass, leading to significant upregulation of FABP1. Similar results also showed that FABP1 levels were higher in patients with non-alcoholic fatty liver disease [43,47]. Therefore, FABP1 is expected to become a diagnostic marker of liver injury [48][49][50]. The MTTP is a key protein for lipid excretion in liver and affects the metabolism of lipids and lipoproteins, leading to the occurrence and development of many diseases [51][52][53][54]. Specific knockout of the MTTP gene in hepatocytes leads to the accumulation of a large amount lipids in mouse hepatocytes, and eventually results in hepatic steatosis and fatty liver [55]. Some evidence showed that the variability of MTTP is associated with the development of NAFLD, raised lipid and risk of atherosclerotic cardiovascular disease [56,57]. A previous study indicated that in the NAFLD of rat induced by a high-fat diet, the expression level of MTTP is markedly decreased [58]. Similarly, HSDinduced fatty liver of largemouth bass decreased MTTP level in liver in our results. ABCA1 played an essential role in the regulation of high-density lipoproteins (HDLs) and reverse cholesterol transport [59]. ABCA1 overexpression increased the concentration of cholesterol and high-density lipoprotein and decreased hepatic lipid contents [60,61]. The deposition of cholesterol in tissues induced the absence of ABCA1 [62,63]. A previous study reported that decreased hepatic ABCA1 expression can cause steatohepatitis in obese adult patients [61]. It also presented the evidence that the expression level of ABCA1 was significantly lower in the NAFLD patients than in the healthy controls [64]. Consistent with these observations, ABCA1 was down-regulated in the MLD group induced by HSD. These findings suggest that the expression level of ABCA1 was associated with the development of fatty liver disease. In this study, the differential expression of FABP1, ABCA1 and MTTP in the liver of largemouth bass in the MLD group showed a high-starch diet increased hepatic lipid contents and affected fat digestion and absorption pathways, and eventually promoted the development of metabolic fatty liver.

Lipidomic Analysis
Significant changes in hepatic lipids are important pathophysiological markers of fatty liver [65]. Lipidomic studies of the mouse liver have shown that NAFLD is associated with an increase in TG [66]. Similar results were observed in our HSD treatment largemouth bass. In this study, most TG species were increased by HSD. During the occurrence of fatty liver, the up-regulation of TGs is attributed to the rate of TG synthesis in the liver exceeding the catabolic rate of TGs [67,68]. Therefore, preventing TG accumulation in liver may contribute to the attenuating effect of fatty liver disease [66,69].
KEGG pathway enrichment analysis indicated that the insulin-resistance pathway was significantly enriched. Many studies believe that NAFLD is closely related to insulin resistance [70][71][72]. The relationship between NAFLD and insulin resistance is bidirectional. Insulin resistance may promote the development of NAFLD, while NAFLD may promote the development of insulin resistance [71,73]. The related research results of our team also found that the persistent hyperglycemia of largemouth bass after a meal is caused by insufficient insulin secretion (unpublished), which is also the reason for the starch intolerance of largemouth bass. Therefore, insulin resistance is a key factor in the development and progression of NAFLD [74] and may be the direction of targeted intervention of MLD in largemouth bass in the future.
Interestingly, in this study we also found that differential lipids were mainly enriched in the sphingolipid metabolism signaling pathway. Liver is the center of ceramide production and usually contains much higher sphingolipid than all other tissues, especially ceramide and sphingomyelin [75][76][77]. Therefore, liver is prone to sphingolipid toxicity [78]. For example, sphingolipid content was significantly increased in the liver of rats fed with a high-fat diet [79]. Mice exposure to exogenous microbiome sphingolipids no longer had lipid accumulation and rescued diet-induced hepatic steatosis [80]. Evidence suggested sphingolipid levels contribute to the development of NAFLD in a variety of ways, including inflammation, insulin resistance and oxidative stress [78].

Integrative Proteomic and Lipidomic Analysis
In this study, using integrative proteomic and lipidomic analysis, we demonstrated the hepatic TG accumulation was important for the development progress of fatty liver induced by HSD. Theoretically, the primary cause of lipid accumulation is the increased uptake of free fatty acids (FFAs) by the liver [56,81]. Undoubtedly, regulation of TG accumulation is important for improving fatty liver disease [82]. Many studies have shown that overexpression of FABP1 enhances fatty acid uptake [83][84][85]. Researchers also found FABP1-deficient mice were protected from TG accumulation in liver induced by a high-fat diet [86,87]. It is expected that reducing the expression or function of FABP1 could inhibit TG accumulation in liver and ameliorate NAFLD [85]. In this study, we similarly observed HSD increased the expression level of FABP1, resulting in more TG accumulation in liver, ultimately leading to fatty liver in largemouth bass.
Meanwhile, accumulated evidence suggested cholesterol homeostasis and hepatic cholesterol accumulation associated with the risk of liver cirrhosis and other liver diseases [88][89][90]. ABCA1 could prevent cholesterol accumulation by transporting cholesterol and phospholipids from cells to apolipoproteins [91][92][93]. Recently, a report showed that ABCA1 were tightly associated with the lipid variables and lipid metabolism [94]. Human studies also suggest the elevation of TG levels usually develop fatty liver disease, and there is an inverse association between TG concentrations and dysfunctional ABCA1 [95].
In addition, a study also reported when the function of ABCA1 is diminished, the TG secretion from liver increased, suggesting there is a strong association between ABCA1 and fatty liver disease [59]. Consistent with this finding, this study also showed the content of ABCA1 was significantly reduced, which exacerbated the accumulation of cholesterol in the liver and worsened metabolic fatty liver disease.
Furthermore, VDAC1 was involved in cholesterol transport and was thought to be part of a complex process that mediates fatty acid transport [96][97][98]. A lack of VDAC1 caused mitochondria to stop oxidizing fatty acids [99]. It is also indicated that the absence of ABCA1 lead to steatosis in metabolic liver disease [100,101]. Recent study had found hepatocyte fat droplet and glycogen accumulation in high-fat diet-treated mice. However, the VDAC1-based peptide decreased serum TG content and hepatocyte fat droplets and increased hepatocyte glycogen accumulation [98]. Moreover, VDAC1 could interact with many proteins, such as metabolism, apoptosis and anti-oxidation-related proteins [102,103], so it could be regarded as a potential therapeutic for metabolic liver disease, such as NAFLD [98]. Here, we demonstrated that the reduction in VDAC1 expression led to cholesterol and lipid accumulation in liver of largemouth bass induced by HSD. Importantly, in addition to hepatic TG accumulation, abnormal transport of cholesterol also plays a critical role in the occurrence and development of fatty liver disease.

Conclusions
We successfully induced metabolic fatty liver phenotype of largemouth bass with high-starch diet. We found that three differential proteins (FABP1, ABCA1 and VDAC1) in fat digestion and absorption and cholesterol metabolism pathways played important roles in the occurrence and development of MLD in largemouth bass after the combination of targeted lipidomics and proteomics analyses. The three differential proteins may be promising targets for therapeutic intervention in fatty liver disease. Our findings propose the possible key regulatory factors in the occurrence and development of metabolic fatty liver, which provide a theoretical basis for the efficient use of compound feed and targeted nutritional regulation in commercial aquaculture.