The Effects of Rumen-Protected Choline and Rumen-Protected Nicotinamide on Liver Transcriptomics in Periparturient Dairy Cows

To investigate the effects of rumen-protected choline (RPC) and rumen-protected nicotinamide (RPM) on liver metabolic function based on transcriptome in periparturient dairy cows, 10 healthy Holstein dairy cows with similar parity were allocated to RPC and RPM groups (n = 5). The cows were fed experimental diets between 14 days before and 21 days after parturition. The RPC diet contained 60 g RPC per day, and the RPM diet contained 18.7 g RPM per day. Liver biopsies were taken 21 days after calving for the transcriptome analysis. A model of fat deposition hepatocytes was constructed using the LO2 cell line with the addition of NEFA (1.6 mmol/L), and the expression level of genes closely related to liver metabolism was validated and divided into a CHO group (75 μmol/L) and a NAM group (2 mmol/L). The results showed that the expression of a total of 11,023 genes was detected and clustered obviously between the RPC and RPM groups. These genes were assigned to 852 Gene Ontology terms, the majority of which were associated with biological process and molecular function. A total of 1123 differentially expressed genes (DEGs), 640 up-regulated and 483 down-regulated, were identified between the RPC and RPM groups. These DEGs were mainly correlated with fat metabolism, oxidative stress and some inflammatory pathways. In addition, compared with the NAM group, the gene expression level of FGF21, CYP26A1, SLC13A5, SLCO1B3, FBP2, MARS1 and CDH11 in the CHO group increased significantly (p < 0.05). We proposed that that RPC could play a prominent role in the liver metabolism of periparturient dairy cows by regulating metabolic processes such as fatty acid synthesis and metabolism and glucose metabolism; yet, RPM was more involved in biological processes such as the TCA cycle, ATP generation and inflammatory signaling.


Introduction
The peripartum period, lasting from 3 weeks before and 3 weeks after parturition, induces remarkable physiological and metabolic adaptations in mammals which are crucial for a good reproductive performance and to ensure the suitable development of the fetus and to provide adequate substrates that are needed in utero and following birth [1]. Simultaneously, cows may face immune system dysfunction and increased inflammatory status, which are important factors limiting their ability to achieve optimal production performance and immune metabolic status and may have a negative impact on the successful transition of offspring from intrauterine to extrauterine [2]. Perinatal cows are subject to the imbalance between energy demand and dietary energy intake, resulting in stimulation of fat mobilization [3] and an elevated level of non-esterified fatty acids (NEFA) in the body. When the elevated NEFA excesses the clearance limit of the body, there will be excessive Metabolites 2023, 13, 594 2 of 16 accumulation of fat and its metabolites in the liver, such as ketone bodies, leading to fatty liver and ketosis [4]. In addition, the excessive NEFA will also lead to oxidative stress [5], rumen acidosis, protein metabolism imbalance and other metabolic syndromes [6]. These metabolic disorders can increase the incidence of other diseases, such as placental retention and endometritis, which will impair the immune system of cows, reduce the milk yield and reproductive efficiency of cows, and cause serious economic losses [7][8][9]. To counteract these negative effects, rumen-protected products are frequently added to the diet of dairy cows during the periparturient period to maintain the physical condition of dairy cows, enhance productivity performance and reduce economic losses.
Choline (2-hydroxyethyl trimethyl ammonium salt) is an essential nutrient for humans and animals which needs to be obtained from the diet in addition to endogenous synthesis. Choline plays a significant role in cell maintenance and growth at various stages of life, including neurotransmission [10], membrane synthesis, lipid transport and one-carbon metabolism [11]. In the liver, choline can be used to synthesize phosphatidylcholine through the cytidine diphosphate-choline pathway and can also be oxidized to betaine, which, through the phosphatidylethanolamine N-methyl transfer pathway, is used to synthesize very-low-density lipoprotein (VLDL) to participate in lipid transport. It has been found that the imbalance between the proportion of phosphatidylcholine and phosphatidylethanolamine in the liver will affect the integrity of the cell membrane, leading to cell damage and inflammatory reactions such as steatohepatitis [12]. In recent years, accumulating evidence suggests that choline can regulate the metabolism and transport of glucose and lipid, as well as the occurrence of oxidative stress in the body by regulating the expression of related genes and activating signal pathways such as AMPK [13][14][15]. Moreover, a large number of animal experiments have shown that adding choline to the diet of periparturient dairy cows could increase the level of VLDL [16,17], regulate fat metabolism [18], relieve metabolic diseases such as fatty liver and ketosis [19,20], increase milk yield [21,22], improve milk quality [23,24], maintain an animal's physical condition [25,26], boost immune system [27,28], and also promote the growth performance of offspring [29].
Nicotinamide is a member of water-soluble B vitamins, which can be served as a precursor to synthesize nicotinamide adenine dinucleotide (NAD) and nicotinamide adenine dinucleotide phosphate [30,31]. In the liver, nicotinamide can be converted and synthesized from nicotinic acid. In addition, nicotinamide can also be produced through the hydrolysis of nicotinamide adenine dinucleotide (NADH) and NAD+ by coenzyme β-Nicotinamide adenine dinucleotide (NAD+). These metabolites can be directly used for ATP generation, DNA synthesis and repair, as well as other metabolic pathways and cell signal transductions. Moreover, nicotinamide also plays an important role in the regulation of lipid and energy metabolism, oxidative stress, mitochondrial dysfunction and other processes [32,33]. Plenty of studies indicate that the addition of nicotinamide and niacin can reduce the levels of total cholesterol, triglyceride, high-density lipoprotein (HDL) cholesterol, LDL cholesterol and VLDL cholesterol as well as improve dyslipidemia, alleviate adipose tissue inflammation [34][35][36], and increase milk yield and the milk protein content of lactating dairy cows [37,38]. However, the difference between choline and nicotinamide in liver metabolism and its mechanism in periparturient dairy cows is still unclear. Therefore, in this study, we added rumen-protected choline (RPC) and rumen-protected nicotinamide (RPM) to the diet of periparturient dairy cows. We then performed a liver biopsy and explored the effects of choline and nicotinamide on the liver metabolism of periparturient dairy cows and the differences in their mechanisms based on transcriptome.

Animals and Sample Acquisition
Ten healthy Chinese Holstein pregnant cows with similar parity were selected and randomly divided into 2 groups: the rumen-protected choline group (60 g/d RPC, active Metabolites 2023, 13, 594 3 of 16 ingredient 15 g/d, purity 65%, the rumen bypass rate 80.3%) and the rumen-protected nicotinamide group (18.7 g/d RPM, active ingredient 9.6 g/d, purity 60%, the rumen bypass rate 85.6%). The experimental period lasted for 35 days, from 14 days before parturition to 21 days after calving. During the monitoring period, the cows in health status and 1-2 parity delivered healthy, viable full-term calves.
The feeding and management of cattle were carried out in accordance with the pasture management system. During the experimental period, two diets were prepared, namely, the pre-partum diet and the post-partum diet. The diets were prepared in the form of total mixed rations (TMR) and offered twice a day, at 07:00 and 16:00, respectively. All experimental animals had free access to food and water. Liver biopsy was performed on the 21st day after delivery (collected at 13:00 each time) by using a biopsy needle, and the liver samples collected were quickly stored at −80 • C. All the experimental procedures in the present study were conducted according to the Animal Protection Law based on the Guide for the Care and Use of Laboratory Animals approved by the Ethics Committee of Northwest A&F University.

Culture of Human Hepatocyte Line LO2
LO2 cells were cultured in Roswell Park Memorial Institute (RPMI) 1640 medium (HyClone, Logan, UT, USA) containing 10% fetal bovine serum (Gibco BRL, Grand Island, NY, USA), and the culture was performed at 37 • C under 5% CO2. To establish a cell model of fatty liver in periparturient dairy cows, LO2 cells were plated into 6-well culture plates to which was added 1.6 mmol/L NEFA subject to two treatments: CHO and NAM. The CHO group contained 75 µmol/L choline in choline-free RPMI 1640 medium, and the NAM group contained 2 mmol/L nicotinamide in nicotinamide-free medium. Before treatment, cells were starved using choline-free RPMI 1640 medium (Gibco BRL, Grand Island, NY, USA) and nicotinamide-free medium (Coolaber, Beijing, China), respectively, without adding fetal bovine serum for 6 h. After starving, LO2 cells with experimental treatment were incubated for 24 h. NEFA was formulated from palmitic acid (3.19 mmol), palmitoleic acid (0.53 mmol), stearic acid (1.44 mmol), linoleic acid (0.49 mmol) and oleic acid (4.35 mmol). All reagents without instructions were purchased from Sigma-Aldrich.

RNA Isolation and Quantitative Real-Time Polymerase Chain Reaction
Total RNA in cells was isolated with Trizol (Invitrogen). An amount of 1-2 µg of total RNA was treated with DNAse (Invitrogen) and reverse transcribed into cDNA using M-MLV Reverse Transcriptase (Accubate Biology, Changsha, China). Quantitative Real-Time PCR (RT-PCR) was performed using SYBR ® Select Master Mix (Accubate Biology, Changsha, China) and the Roche LightCycler ® 480 Real-Time PCR System. The 2 −∆∆Ct method was used to determine the relative quantitative gene expression levels, normalized by β-actin (Tables 1 and 2) [39].

RNA Isolation and cDNA Library Construction
Total RNA in liver biopsies was extracted, and the RNA concentration and purity were determined using a Nanodrop2000 (Thermo Fisher Scientific, US). Agarose gel electrophoresis was used to detect the integrity of RNA. An Agilent 2100 bioanalyzer was used to measure the RNA value. The enriched mRNA was segmented then reverse synthesized to cDNA, which completed the construction of the library to be sequenced. The following conditions must be met for the creation of a single library to be sequenced: the total amount of RNA ≤ 1 µg, the concentration ≥ 35 ng/µL, OD260/280 ≥ 1.8, OD260/230 ≥ 1.0. After passing the quality inspection, the library was sequenced on the IlluminaNovaseq6000 platform, and the database was established by Shanghai Meiji Biomedical Technology Co.

Analysis of Differential Expression of Genes
The software DESeq2 based on negative binomial distribution was used to analyze the differential expression of genes (DEGs). According to p-value and |log 2 FC| as screening conditions, log 2 FC > 0 indicates gene expression is up-regulated and log 2 FC < 0 indicates gene expression is down-regulated.

Analysis of DEG Enrichment
The software Goatools was used for GO (Gene Ontology) enrichment analysis of genes/transcripts, and R script was used for KEGG (Kyoto encyclopedia of genes and genomes) pathway enrichment analysis of genes/transcripts. Fisher's Exact Test and the Benjamin-Hochberg method was used to correct the p-value. Significance of GO enrichment and KEGG enrichment was considered as a p-value < 0.05.

Statistical Analysis
MetaboAnalyst 5.0 (Wishart Research Group, Edmonton, AB, Canada) software was used to perform orthogonal partial least squares-discriminant analysis (OPLS-DA) on the samples from the RPC group and RPM group. A volcano plot was used to analyze DEGs and their changes in each group, and a cluster heatmap was used to analyze the gene expression patterns in each sample. Statistical significance was determined with IBM SPSS Statistics 26.0. Data were represented as an average ± SD and treated with normal distribution test. Student's t-test was used for two-group comparisons. ** indicates p < 0.01, *** indicates p < 0.001.

Sequencing Data Quality Control
After sequencing data generation, the original sequencing data for each sample was subjected to a quality assessment such as the base group's error rate and content distribution, and quality control was carried out to obtain high-quality control data (clean data). The Q20 and Q30 quality values were 97.51% and 93.02%, respectively. The GC content ratio was within an acceptable range, which indicated that the sequencing results were reliable and could be used for subsequent results analysis (Supplementary Table S1).

Analysis of DEGs
Transcripts with an absolute value of |log 2 FC| ≥ 1 and adjusted p-value < 0.05 were defined as significant DEGs. A total of 11,023 genes were expressed in all 10 samples, while 828 genes were specifically expressed in the RPC group and 481 were specifically expressed in the RPM group (Supplementary Figure S1). From the comparisons between these two groups, 1123 DEGs were identified, 640 up-regulated and 483 down-regulated ( Figure 1; Supplementary Table S2). The cluster heatmap analysis was performed on the identified DEGs, as shown in Figure 2, which shows that the samples from the RPC and RPM groups had obvious clustering. OPLS-DA analysis showed that the samples from the RPC and RPM groups were significantly separated ( Figure 3). The Q20 and Q30 quality values were 97.51% and 93.02%, respectively. The GC content ratio was within an acceptable range, which indicated that the sequencing results were reliable and could be used for subsequent results analysis (Supplementary Table S1).

Analysis of DEGs
Transcripts with an absolute value of |log2 FC| ≥ 1 and adjusted p-value < 0.05 were defined as significant DEGs. A total of 11,023 genes were expressed in all 10 samples, while 828 genes were specifically expressed in the RPC group and 481 were specifically expressed in the RPM group (Supplementary Figure S1). From the comparisons between these two groups, 1123 DEGs were identified, 640 up-regulated and 483 down-regulated ( Figure 1; Supplementary Table S2). The cluster heatmap analysis was performed on the identified DEGs, as shown in Figure 2, which shows that the samples from the RPC and RPM groups had obvious clustering. OPLS-DA analysis showed that the samples from the RPC and RPM groups were significantly separated ( Figure 3).

Analysis of Differential Gene Enrichment
The GO analysis was performed to detect the function of the identified DEGs. The 268 DEGs were enriched in 44 GO terms, including 20 terms in biological process (BP), 13 terms in cellular component (CC), and 10 terms in molecular function (MF) (Figure 4). Table 4 shows the top 20 items under GO enrichment which are mainly involved with terms relating to fat anabolic processes, including lipid metabolic process, small molecule metabolic process, cellular hormone metabolic process, xenobiotic metabolic process, extracellular space, monooxygenase activity, organic acid metabolic process, oxoacid metabolic process, hormone metabolic process, carboxylic acid metabolic process, monocarboxylic acid metabolic process, oxidoreductase activity, tetrapyrrole binding, chemokine activity, steroid metabolic process, detoxification, fatty acid omega-hydroxylase activity, steroid hydroxylase activity, etc.

Effects of Choline and Nicotinamide on Gene Expression Related to Liver Metabolism
Based on the above analysis results, in order to further verify the effect of choline and nicotinamide on liver metabolic function, the mRNA expression of some DEGs (FGF21, CYP26A1, SLC13A5, SLCO1B3, FBP2, MARS1, CDH11) associated with liver metabolic function was tested. As shown in Figure 6, compared with the NAM group, the gene expression of FGF21, CYP26A1, SLC13A5, SLCO1B3, FBP2, MARS1, CDH11 in the CHO group increased significantly (p < 0.05).

Discussion
The severe negative energy balance in periparturient dairy cows stimulates a large amount of body fat mobilized to support the energy demand of animals, resulting in the disorder of fat metabolism, particularly in the liver. Numerous studies have found that dietary choline and dietary nicotinamide supplementation can effectively alleviate the occurrence of metabolic diseases such as fatty liver, oxidative stress and inflammation [40][41][42]. This study compared the transcriptomic differences in liver metabolism in periparturient dairy cows with the addition of RPC and RPM, identified differential genes and

Discussion
The severe negative energy balance in periparturient dairy cows stimulates a large amount of body fat mobilized to support the energy demand of animals, resulting in the disorder of fat metabolism, particularly in the liver. Numerous studies have found that dietary choline and dietary nicotinamide supplementation can effectively alleviate the occurrence of metabolic diseases such as fatty liver, oxidative stress and inflammation [40][41][42]. This study compared the transcriptomic differences in liver metabolism in periparturient dairy cows with the addition of RPC and RPM, identified differential genes and pathways, and discussed the differences between choline and nicotinamide on liver metabolism at the transcriptome level.
In the present study, 640 genes up-regulated by the addition of RPC or RPM and 483 genes down-regulated by the addition of RPC or RPM were identified; the top 20 GO items were mainly involved in fat anabolic processes. The top two DEGs were phenylalanyl-tRNA synthetase 2 and fibroblast growth factor 21. Fibroblast growth factor 21 (FGF21), a member of the fibroblast 19 subfamily, is synthesized by the liver and regulated by peroxisome proliferators-activated receptors-α, carbohydrate response element binding protein and other reaction elements. It can act directly on liver cells as a paracrine or autocrine hormone and can also enter the systemic circulation and act on other tissues and organs to regulate glucose, lipid and energy metabolism [43,44]. Some researchers fed FGF21 knockout mice with a ketogenic diet and found that fat levels in the plasma of the mice were increased, obvious steatosis occurred, and the regulatory function of ketogenesis and sugar metabolism was impaired [45,46]. The chronic accumulation of triglycerides in the liver improves the β-oxidation rate of fatty acids in liver cells. When the energy supply from fatty acid oxidation exceeds the energy demand, it leads to the formation of reactive oxygen species, which, in turn, cause oxidative stress, thus inducing an increase in the synthesis of FGF21 [47,48]. It is worth noting that FGF21 can respond to low protein levels in the body to increase its own circulation and increase energy consumption to cope with the increase in calorie intake [44,49]. The results of the gene expression levels of FGF21 between two treatment groups showed that choline probably produces an effect on liver fat metabolism and reduces body energy consumption by up-regulating the level of FGF21. However, the reason for the reverse change of FGF21 expression under the influence of choline and nicotinamide needs to be further explored.
The top KEGG pathway of the DEGs in response to the RPC and RPM addition was for retinol metabolism in this study. Vitamin A can be stored in hepatic stellate cells in the form of retinol and provide vitamin A for various tissues of the body. All-trans-retinoic acid is an important active metabolite of vitamin A, and excessive retinoic acid levels in the body will cause changes in physiological functions of the body [50,51]. Retinoic acid hydroxylase CYP26A1 is very important for decomposing exogenous retinoic acid [52]. The results of the cell experiment show that the gene level of Cyp26a1 in the CHO group was significantly higher than that in the NAM group, indicating that choline can regulate retinoic acid hydroxylase to regulate liver retinol metabolism and maintain its level balance.
Subsequently, the genes closely related to liver metabolism among the DEGs to clarify the differences in the regulation of liver metabolism by choline and nicotinamide were further discussed. The solute carrier family 13 member 5 (SLC13A5) used to mediate citrate transport is mainly expressed in the liver of mammals. As a key intermediate of the tricarboxylic acid (TCA) cycle, citrate plays a key role in the cell metabolism of carbohydrates and fatty acids, as well as energy generation in mitochondria, and can also be served for the biosynthesis of triglycerides, fatty acids, cholesterol and low-density lipoprotein [53,54]. Brachs et al. (2016) found that in the Slc13a5 knockout of mice fed with a high-fat diet, tolerance to insulin resistance and fat deposition induced by the high-fat diet was increased [55]. The expression of Slc13a5 was also increased in the model of type-2-diabetes rats [56]. In addition, some researchers found that the level of intracellular citrate and phospholipid was significantly reduced after silencing the Slc13a5 in HepG2 and Huh7 cells [57]. Solid carrier organic animion transporter family member 1B3 (SLCO1B3, also known as OATP1B3) is a member of the solute transporter superfamily, which is mainly expressed in the basolateral membrane of hepatocytes. It can transport a variety of endogenous and xenobiotic compounds (such as prostaglandins, steroid hormone conjugates, bilirubin, bile acid, and thyroid hormone) to liver cells for metabolism [58,59]. It has been found that the new splice of OATP1B3, cancer-type OATP1B3, can up-regulate the expression of carnitine palmitoyl transfer and NADH: ubiquinone oxidoreductase subunit A2 through interaction with insulin-like growth factor 2 mRNA-binding protein 2, thereby promoting fatty acid β-oxidative and mitochondrial oxidative phosphorylation activities and increasing ATP production. This results in the formation of plate-like pseudopodia and migration and the invasion of high-grade serous ovarian cancer cells [60]. Sec14-like proteins belong to typical class III phosphatidylinositol transfer proteins, which can transfer phosphatidylinositol and PC between different biological membranes, exchange phosphatidylinositol to PC, and vice versa, so as to maintain membrane-lipid balance [61,62]. This study found that the expression levels of SLC13A5 and SLCO1B3 in the CHO group were significantly higher than that in the NAM group, showing that choline may regulate energy metabolism and material transport of the liver by up-regulation of SLCO1B3, thereby improving fat metabolite accumulation and energy imbalance. Nicotinamide may reduce the accumulation of citrate by promoting the TCA cycle, allowing the metabolites of citrate, acetyl coenzyme A and oxaloacetate to flow towards ATP synthesis rather than fatty acid synthesis, resulting in the lower gene expression of SLC13A5.
Analysis of DEGs revealed a significant increase in the expression levels of FBP and transaminase in the RPC group. Fructose diphosphatase (FBP) can not only regulate the glucose/glycogen synthesis of carbohydrate precursors but can also interact with proteins (such as ATP synthase, HIF1-α and NF-κB) to influence cell-cycle-dependent events, mitochondrial biogenesis and membrane polarization, expression of glycolytic enzymes, induction of synaptic plasticity, and even cancer progression [63]. Transketolase, a key enzyme of the pentose phosphate pathway that is ubiquitous in all organisms, provides a special connection between glycolysis and the non-oxidation period of the pentose phosphate pathway [64]. The study of LO2 cells suggests that the expression level of FBP significantly increased under the regulation of choline, indicating that choline may regulate the process of liver glucose metabolism through FBP and transketolase.
Methyl tRNA synthetase 1 (MARS1) and glycyl tRNA synthetase 1 (GARS1) belong to aminoacyl tRNA synthetases (ARSs). Aminoacyl tRNA synthetases, one of the potential markers of pneumonia, are generally used for protein synthesis and can interact with proteins in the signal pathway of the mammalian target of rapamycin 1, cyclin-dependent kinases 4 and the vascular endothelial growth factor receptor [65][66][67]. Several studies have found that inhibitors of aminoacyl tRNA synthetases have been used in the treatment of many diseases [68][69][70]. Fibroblasts will lose the function of matrix remodeling under pathological conditions, leading to tissue destruction and fibrosis [71,72]. Cadherin (cadherin-11, CDH11) is a mesenchymal cadherin that can be expressed in fibroblasts of various tissues and is also a marker and functional regulator of fibroblasts and can mediate homotypic cell adhesion that is important in histomorphogenesis and structure [73,74]. In CDH11knocked diet-induced obese mice, inflammation was reduced and blood glucose levels were controlled [75]. The result of cell assays found that, compared with the CHO group, the levels of formyl tRNA synthetase 1 and CDH11 genes in the NAM group decreased significantly. Thus, it can be concluded that nicotinamide may affect liver inflammation by regulating the gene expression of MARS1, CDH11 and their related signal pathways, thereby preventing the occurrence of liver diseases.
The GO pathway enrichment analysis of the differential genes between RPC and RPM showed that the enriched pathways were mainly associated with the metabolic processes of lipid, cytokine, carboxylic acid, organic acid, oxoacid, and the activity of fatty acid omega-hydroxylase, steroid hydroxylase, oxidoreductase, chemokine, etc. The different metabolic pathways enriched by the KEGG pathway indicated that the differences between choline and nicotinamide in the regulation of liver metabolism were mainly manifested in fat metabolism, biosynthesis, drug metabolism, etc. This may be due to the fact that choline can also regulate the metabolism process of fat and energy in the liver as well as related signals of oxidative stress and inflammation as a signal, in addition to catabolizing into other substances and then participating in membrane biosynthesis, fat metabolism and transportation, and other substance synthesis [14,75]. Nicotinamide is mostly used to synthesize NAD for energy supply after entering the liver, to maintain the level of reduced glutathione and thioredoxin in the antioxidant system, and to participate in cell energy metabolism and other processes.

Conclusions
In this study, the effects of RPC and RPM on the liver transcriptome of perinatal cows were investigated, and 1123 DEGs and 27 different KEGG pathways were identified, mainly including retinol metabolism, metabolism of xenobiotics by cytochrome P450, steroid hormone biosynthesis, pyrimidine metabolism, bile secretion, linoleic acid metabolism, arachidonic acid metabolism, chemokine signaling pathway, aminoacyl-tRNA biosynthesis, purine metabolism, glutathione metabolism, etc. The expression levels of genes related to liver metabolism were identified, and it was found that, compared with the niacinamide group, the choline group could up-regulate the expression levels of FGF21, CYP26A1, SLC13A5, SLCO1B3, FBP2, MARS1, and CDH11. Overall, this study shows that RPC plays a prominent role in liver metabolism by regulating metabolic processes such as fatty acid synthesis and metabolism and glucose metabolism. Yet, RPM is more involved in biological processes such as the TCA cycle, ATP generation and inflammatory signaling.

Institutional Review Board Statement:
In this study, all animal procedures were approved by the Institutional Animal Care and Use Committee of Northwest A&F University (Yangling, Shaanxi, China). The protocols of this study were specifically approved with the protocol number NWAFAC1008.

Data Availability Statement:
The data presented in this study are available in the article and the supplementary materials.

Conflicts of Interest:
The authors declare no conflict of interest.