Abstract
Parasitic infection improves metabolic homeostasis in “western diet”-induced obesity through the regulation of adipogenesis. However, the underlying mechanism is not yet fully understood. Using microarray analysis, this study investigated the long non-coding RNA (lncRNA) and messenger RNA (mRNA) profiles of subcutaneous adipose tissues from mice infected with Echinococcus granulosus protoscoleces. A total of 1052 mRNA (541 upregulated, 511 downregulated) and 220 lncRNA (126 upregulated, 94 downregulated) transcripts were differentially expressed (fold change ≥2, P < 0.05) in the infected subcutaneous adipose tissues. When compared with the control group, the infected mice showed a decrease in adipose tissue mass and a reduction in adipocyte size. Indirect calorimetry revealed the change in the energy metabolism after infection, characterized by a lower CO2 production and O2 consumption, a sharp decline in carbohydrate oxidation, and a slight increase in fat oxidation. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the parasitic infection reprogrammed a complex metabolic network. Notably, “lipoxygenase” and “arginine and proline metabolism” pathways were significantly upregulated while “glycolysis,” “tricarboxylic acid cycle,” “de novo lipogenesis,” and “lipid droplet” pathways were dramatically downregulated. In addition, several key lncRNAs were associated with insulin resistance and adipocyte differentiation. Overall, the present study suggested that E. granulosus infection could enhance lipolysis. Thus, our findings provide novel insights into parasite-mediated metabolic homeostasis, and into the mechanism of hypertrophic adipocytes in obesity.
Introduction
Parasites pose a threat to human health since they are widely considered as agents of disease and death (). Following infection, these organisms can compete with the host for nutrients, while also regulating its metabolism and immune system. Parasites have developed refined machinery to evade the host’s immune system by manipulating its response (). It is well known that T-helper type 2 (Th2) responses are central to parasitic infections, and can participate in regulation of host’s immune system (). Th2 responses can create an anti-inflammatory immune microenvironment and improve host metabolic disorders (). Such immunological characteristics make parasites and their derived molecules attractive models for the treatment of inflammatory and metabolic disorders. For example, Trichuris muris infection was shown to protect mice deficient in NOD2 (susceptibility gene for Crohn’s disease) from Th2 response-dependent intestinal abnormalities (). Additionally, the excreted or secreted molecules-62 (ES-62) derived from Acanthocheilonema viteae (a filarial nematode parasite of gerbils) are anti-inflammatory by virtue of covalently attached phosphorylcholine moieties (). The subcutaneous administration of this immunomodulator was recently reported to protect against collagen-induced arthritis ().
Currently, obesity is a main public health concern since obese individuals are at a greater risk of developing chronic disease and often present clinical parameters of metabolic syndrome (MetS), insulin resistance, and systemic markers of chronic low-grade inflammation (; ; ). Adipose tissue accumulation is a common clinical phenomenon in obese individuals (). Adipose tissue is the main lipid storage depot in the human body and, as such, plays a crucial role in the regulation of dietary fat daily influx into circulation (; ). However, in obesity, the adipose tissue may not expand properly when there is a need to store the energy surplus. In turn, this may lead to ectopic fat deposition in tissues involved in metabolic homeostasis and, consequently, insulin resistance and MetS (). Thus, adipogenesis is an essential event in obesity and MetS development.
An increasing number of studies have shown that parasitic infection can regulate adipogenesis. The intestinal nematode parasite Heligmosomoides polygyrus was reported to attenuate obesity in high-fat diet-induced obese mice (). The infected individuals exhibited a markedly upregulated expression of uncoupling protein 1 (a key protein involved in energy expenditure) in the adipose tissue, suppression of glucose and triglyceride levels, and alteration in the expression of key genes involved in lipid metabolism (). Trypanosome brucei is a parasite that uses the adipose tissue as a niche, enhancing fatty acid β-oxidation and causing weight loss in model mice (). These studies suggested that helminths may play a protective role against MetS () and that exploring the underlying mechanism of adipogenesis in a parasitic infection context may provide a novel strategy for the treatment of obesity and MetS.
Echinococcus granulosus (E. granulosus) is a known cestode with high host adaptability and worldwide distribution. After oral uptake of the eggs by intermediate hosts, the activated oncospheres hatch in the gastrointestinal tract, where they penetrate the intestinal wall and enter the bloodstream. Eventually, the parasites reach the visceral organs (such as the liver, lung, heart, and kidney), where they develop into hydatid cysts. E. granulosus can infect many intermediate hosts and go unnoticed for several decades, since it has evolved strategies to subvert the immune responses (). Previous studies have focused on the immunoregulatory role of E. granulosus and its derived molecules (, ; ). However, their metabolic regulatory activity of E. granulosus and its derived molecules has been largely neglected.
Long non-coding RNAs (lncRNAs) are emerging as important regulators of cellular signaling and gene expression in numerous cell types. These molecules are long RNA transcripts (>200 bp) that lack protein-coding ability (). lncRNAs can regulate cell function through a variety of mechanisms (). For instance, they can function as scaffolds to bring two or more proteins into a functional ribonucleoprotein complex. lncRNAs can also act as decoys to titrate a protein away from its original target, as guides to recruit chromatin modification enzymes to specific loci on chromosomes, and as microRNA sponges to buffer the inhibitory functions of microRNAs on gene expression. Recent studies have revealed an extensive list of regulators that govern adipocyte differentiation, and the scope of these regulators has been greatly expanded to lncRNAs. For example, the very recently identified adipogenic lncRNA Plnc1 was found to promote differentiation into mature adipocytes (). However, the understanding of the underlying mechanisms of lncRNA in adipocyte regulation is still in its infancy.
The present study investigated whether the larval E. granulosus infection in mice could regulate adipogenesis and determined the differential expression profiles of both mRNA and lncRNA. The results showed that the fat mass and adipocytes’ size were significantly decreased, following infection. Moreover, several differentially expressed lncRNAs and mRNAs were found to regulate the metabolic reprogramming and differentiation of adipocytes. These results provided a novel insight into host-parasite interplay and several clues for the treatment of obesity and its associated Mets.
Materials and Methods
Ethics Approval
This study was carried out in strict accordance with the recommendations of the Guide for the Care and Use of Laboratory Animals of the Ministry of Health, China. The protocol was approved by the Laboratory Animal Welfare and Ethics Committee (LAWEC) of Xuzhou Medical University (Xuzhou, China, SCXK (Su) 2015-0009). All surgery was performed under sodium pentobarbital (100 mg/kg) anesthesia and all efforts were made to minimize suffering.
Mice and Parasites
Female C57BL/6J mice (6–8 weeks old) were obtained from Shanghai Laboratory Animal Center (SLAC, Shanghai, China) and were bred in the university facilities. All mice were fed by the experimental animal growth and reproduction diet (Amufi company, Suzhou, China) containing 20% protein, 4% lipid and 5% carbohydrate, and were housed in an air-conditioned room at 24°C with a 12 h dark/light cycle and permitted free access to food and water.
The protoscoleces of E. granulosus (EgPSC) were obtained from the hydatid cysts of naturally infected sheep livers under aseptic conditions, and washed three times using 0.9% NaCl containing 100 mg/ml penicillin and 100 U/ml streptomycin. Female C57BL/6J mice were intraperitoneally injected with a 200 μl suspension containing 2000 live EgPSC in 0.9% NaCl to establish the Eg model, and the control mice were received 200 μl 0.9% NaCl (, ; ). The infected and control mice were sacrificed for experiment after 7 months of infection.
Indirect Calorimetry
The profile of systemic energy metabolism was determined by indirect calorimetry using the Oxymax CLAMS system (Columbus Instruments) at 7th month post infection. Mice were housed individually in CLAMS cages and allowed to acclimate for 24 h (12 h day and 12 h night) in 25°C with unrestricted excess to food and water. O2 consumption (VO2), CO2 release (VCO2), Respiratory Exchange Ratio (RER) and food intake were recorded for 24 h spanning a single light-dark cycle. Carbohydrate oxidation was calculated using the formula [(4.585∗VCO2)–(3.226∗VO2)]∗4, in which the 4 represents the conversion from mass per time unit to kcal per time unit (). Similarly, fat oxidation was calculated using the formula [(1.695∗VO2)-(1.701∗VCO2)]∗9 ().
Glucose Tolerance Test (GTT)
Glucose tolerance test was performed at 7th month post-infection. After fasting for 16 h, mice were given a glucose solution in the abdominal cavity at a dose of 1 g/kg. Blood samples were collected at 0, 30, 60, 90, and 120 min after injection, and blood glucose test paper was used to detect blood glucose. The curve of blood glucose over time was drawn by Graphpad 6 software and the total area under the curve (AUC) was calculated.
Serum Biochemical Indices
After sacrifice, the blood samples from mice eyeballs were centrifuged in 12,000 rpm, 4°C for 20 min. The serum concentration of aspartate aminotransferase (AST), alanine aminotransferase (ALT), triglyceride, cholesterol, high density lipoprotein (HDL) and low density lipoprotein (LDL) was measured by automatic biochemical analyzer (Roche Cobas 701, Roche Diagnostics Ltd.).
Histopathological Examination
Adipose tissue was immersed in 4% paraformaldehyde for 4 h, and transferred to 70% ethanol. Individual lobes of adipose tissue were placed in processing cassettes, dehydrated through a serial alcohol gradient, and embedded in paraffin wax blocks. Before immunostaining, 5 μm-thick adipose tissue sections were dewaxed in xylene, rehydrated through decreasing concentrations of ethanol, and washed in PBS. Stain the slide with 1% eosin Y solution for 10–30 s with agitation. Dehydrate the sections with two changes of 95% alcohol and two changes of 100% alcohol for 30 s each. At last, extract the alcohol with two changes of xylene (). After staining, sections were dehydrated through increasing concentrations of ethanol and xylene and cover with a coverslip.
Microarray Profiling
The lncRNA and mRNA expression patterns of three EgPSC infected mice and three control mice were detected. The Agilent-085631 was used in this experiment and data analysis of the 6 samples were conducted by OE Biotechnology Co., Ltd. (Shanghai, China). Total RNA was quantified by the NanoDrop ND-2000 (Thermo Fisher Scientific, Waltham, MA, United States) and the RNA integrity was assessed using Agilent Bioanalyzer 2100 (Agilent Technologies). The sample labeling, microarray hybridization and washing were performed based on the manufacturer’s standard protocols. Briefly, total RNAs were transcribed to double strand cDNA, then synthesized into cRNA and labeled with Cyanine-3-CTP. The labeled cRNAs were hybridized onto the microarray. After washing, the arrays were scanned by the Agilent Scanner G2505C (Agilent Technologies).
Differential Expression Analysis
The raw data were analyzed using Feature Extraction software (version 10.7.1.1; Agilent Technologies) and then normalized using percentile normalization. Probes with least one of two conditions flagged in “P” were chosen for further data analysis. Calculated with the t-test, differentially expressed lncRNAs were identified through fold change and the P-value. Aberrantly expressed lncRNAs and mRNAs were defined as fold change ≥2.0 and P < 0.05. Hierarchical clustering was performed on six mouse adipose tissue samples using Cluster 3.0 (Stanford University School of Medicine, California, United States) and TreeView 2.0 (Baryshnikova Lab, Princeton University, New Jersey, United States) to distinguish the distinguishable gene expression pattern among the samples.
GO and KEGG Enrichment Analysis
The Gene Ontology (GO) analysis was executed to investigate the biological functions of mRNAs and lncRNAs. Biological process (BP), cellular component (CC), and molecular function (MF) were involved in the GO terms. The enrichment analysis of GO terms with P < 0.05 was considered significantly. Moreover, to predict the possible pathways, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was adopted to map the possible pathways of these differentially expressed genes.
Co-expression Network Analysis
Cytoscape version 3.5.1 (US National Institute of General Medical Sciences) was used to describe the co-expression of lncRNAs and the protein-coding genes. Correlations with P < 0.05 were considered statistically significant.
Quantitative Reverse Transcription PCR
Quantitative reverse transcription-PCR (qRT-PCR) was performed to validate the microarray results. Quantification was performed with a two-step reaction process: reverse transcription (RT) and PCR. Each RT reaction consisted of 0.5 μg RNA, 2 μl of 5 × TransScript All-in-one SuperMix for qPCR and 0.5 μl of gDNA Remover, in a total volume of 10μl. Reactions were performed in a GeneAmp® PCR System 9700 (Applied Biosystems, United States) for 15 min at 42°C, 5 s at 85°C. The 10 μl RT reaction mix was then diluted × 10 in nuclease-free water and held at −20°C. Real-time PCR was performed using LightCycler® 480 II Real-time PCR Instrument (Roche, Swiss) with 10 μl PCR reaction mixture that included 1 μl of cDNA, 5 μl of 2 × PerfectStartTM Green qPCR SuperMix, 0.2 μl of forward primer, 0.2 μl of reverse primer and 3.6 μl of nuclease-free water. Reactions were incubated in a 384-well optical plate (Roche, Swiss). The exact thermal cycler conditions were described in a previous study (). Each sample was run in triplicate for analysis. At the end of the PCR cycles, melting curve analysis was performed to validate the specific generation of the expected PCR product. The primer sequences are shown in Supplementary Table S1, which were designed in the laboratory and synthesized by TsingKe Biotech based on the mRNA sequences obtained from the NCBI database. The expression levels were calculated using the 2–ΔΔCt method ().
Statistical Analysis
All statistical analysis in this study was performed using the software Graphpad Prism 6.0. All measurement data were presented as the mean with standard error of the mean (SEM). Differentially expressed mRNA and lncRNAs were identified using t-tests. P < 0.05 considered as statistically significant.
Results
Enhanced Lipolysis After EgPSC Infection
There is accumulating evidence showing that parasitic infection causes pathophysiological changes in host adipose tissues (). However, whether this occurs following EgPSC infection remains unknown. EgPSC-infected mice gained body weight due to the cysts in their inner organs, which demonstrated the successful establishment of the infection model. After removing the cysts, the EgPSC-infected mice showed a decrease in body weight (P = 0.041, Figure 1A). Moreover, the mass of subcutaneous fat were significantly lower in the infected mice than in the control group (P < 0.001, Figure 1B). Consistent with this observation, the percentages of subcutaneous fat mass/body weight (without cysts) decreased remarkably (P < 0.001, Figure 1C). Hematoxylin-eosin (H&E) staining revealed shrunken adipocytes in post-infection subcutaneous fat tissue (Figure 1D). The diameter and superficial area of subcutaneous adipocytes in infected mice also decreased (P < 0.001, Figures 1E,F). These results suggested that EgPSC could promote lipolysis in the host. In addition, the relative expression of M1 type macrophage markers (tumor necrosis factor alpha, interleukin-6 and inducible nitric oxide synthase) was significantly downregulated (P < 0.05, Supplementary Figure S1A), while the relative expression of M2 macrophage markers (arginase-1 and interleukin-10) was remarkably upregulated (P < 0.01, Supplementary Figure S1B), which was consistent with the fact that M2 macrophages promote lipolysis ().
FIGURE 1
In contrast, there was no obvious change in liver mass, triglyceride, cholesterol, low- and high-density lipoproteins after infection (P > 0.05, Supplementary Figures S2A–E), while serum levels of aspartate aminotransferase (AST) and alanine aminotransferase (ALT) were decreased (Supplementary Figures S2F–H). However, the parasitic infection significantly improved the mice glucose tolerance (P < 0.001, Supplementary Figure S3), indicating an alternative host metabolism.
Characterization of Systemic Energy Metabolism After EgPSC Infection
The adipose tissue is the main lipid storage depot in the body and has a crucial role in burning lipids for thermogenesis (). As a result, this tissue contributes to the total energy expenditure and can modify the energy metabolism of the whole body. We hypothesized that EgPSC infection could alter the systemic energy metabolism in mice because of the evident decrease in adipose tissue mass, following infection. To corroborate this, an indirect calorimetry assay was conducted with mice housed in metabolic cages for 24 h (1 day and 1 night).
Compared with the control group, the infected mice showed a comparable level of Respiratory Exchange Ratio (RER) throughout the experimental period (P > 0.05, Figure 2A). However, each EgPSC-infected mouse released less CO2 (P < 0.05, Figure 2B) and consumed less O2 (P < 0.05, Figure 2C). The results indicated a lower energy metabolism following infection. Moreover, a sharp decline in carbohydrate oxidation (P < 0.01, Figure 3A) and increased fat oxidation per mouse (P > 0.05, Figure 3B) were observed, suggesting that EgPSC infection could force the mice to use fatty acids as a carbon source and enhance lipid combustion (), which is consistent with fat loss following infection. However, there was no difference in food intake between the EgPSC-infected and control mice (Supplementary Figure S4), indicating that low appetite was not the cause of the observed metabolic changes.
FIGURE 2
FIGURE 3
Differentially Expressed mRNAs and lncRNAs in Adipose Tissues After EgPSC Infection
Microarray analysis was performed to identify the pathophysiological events of EgPSC infection and the mechanism of adipogenesis inhibition. An Agilent Mouse lncRNA Microarray (design ID: 085631) was used in this experiment. Differentially expressed mRNAs or lncRNAs were identified through fold change and the P-values were calculated using the t-test. The threshold between upregulated and downregulated RNAs was set as fold change ≥2.0 and P < 0.05.
Among all the mRNAs and lncRNAs detected in the subcutaneous adipose tissues of infected mice, 1052 mRNAs (Supplementary Table S3) and 220 lncRNAs (Supplementary Table S4) were differentially expressed (upregulated: 541 mRNAs, 126 lncRNAs; downregulated: 511 mRNAs, 94 lncRNAs). The volcano plots of differentially expressed mRNAs and lncRNAs are shown in Figures 4A,B, respectively. Red and blue dots indicated significantly upregulated and downregulated transcripts, respectively. The different expression profiles of lncRNAs and mRNAs between control and infected mice indicated that these transcripts may be involved in the differential process of adipose tissues. Moreover, the scatter plots of differentially expressed mRNAs and lncRNAs are shown in Supplementary Figure S5.
FIGURE 4
To validate the results obtained with microarray analysis, 7 differentially expressed mRNAs (Acod1, Alox12e, Acss2, Fabp5, Fgf18, Lrp2bp, GAPDH) and 5 lncRNAs (Gm36475, Olfr1372-ps1, Gm34973, Lmcbate1, and H19) were randomly selected for qRT-PCR. A consistent expression trend of these genes was observed (Figure 4C).
GO and KEGG Pathway Analyses of Differentially Expressed mRNAs
To evaluate the enrichment of mRNAs in biological processes, cellular components, and molecular functions, Gene Ontology (GO) analysis was performed on 1052 significantly expressed mRNAs noted from the microarray outcome. Based on the biological processes in the gene ontology classification, the differentially expressed mRNAs (fold change ≥2) were classified into different functional categories. Figure 5 shows the GO annotation of upregulated (Figure 5A) and downregulated (Figure 5B) mRNAs with the top 10 enrichment scores covering the three domains. According to GO analysis, upregulated genes were mainly involved in a variety of biological processes, including cell signaling transduction and lipid metabolism, such as “G-protein coupled receptor activity (GO: 0004930)” and “lipoxygenase pathway (GO: 0019372)” Conversely, “lipid metabolic process (GO: 0006629),” “cholesterol biosynthetic process (GO: 0006695),” and “lipid droplet (GO: 0005811)” were downregulated (also listed in Supplementary Table S5). These terms may play a key role in the disorder of lipid metabolism.
FIGURE 5
To further understand the predicted biological functions of the aforementioned differentially expressed mRNAs, their pathways and molecular interactions were then predicted by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The top 30 enriched pathways are shown in Figure 6. Among others (Supplementary Table S2), “Arginine and proline metabolism (mmu00330),” “cyclic adenosine monophosphate (cAMP) signaling pathway (mmu04024),” “carbon metabolism (mmu01200),” and “insulin resistance (mmu04931)” were thought to be closely associated with adipocyte metabolism and deserved further study. Notably, in comparison to the control mice, the expression of enzymes in “glycolysis,” “TCA,” and “de novo lipogenesis” pathways were significantly downregulated, while “fatty acid oxidation” was activated after infection (Figure 7 and Supplementary Table S6). This suggested that the EgPSC infection reprogrammed the adipose tissue metabolism.
FIGURE 6
FIGURE 7
lncRNA Transcription Factor Network Analysis
Hypergeometric distribution was used to characterize the relationship between transcription factors and differentially expressed lncRNAs. The frequency of each transcription factor (TF) was counted, and the one with more functional annotations was recorded to reflect the overall functional distribution of different lncRNAs. A core network of the top 100 lncRNA-TF pairs was constructed by ranking the P-value (Figure 8). These lncRNAs were mostly regulated by contactin 2, metal-regulatory transcription factor 1, aristaless like homeobox 4 (ALX4), aryl hydrocarbon receptor nuclear translocator 2 (ARNT2), and nuclear receptor 1 type H2 (NR1H2). Among these TFs, ARNT2 has been demonstrated to be associated with the prevention of obesity and obesity-related diseases (). Studies have shown that ARNT2 is necessary for the production of secretory hormones, and an induced missense mutation in ARNT2 developed hepatic steatosis and abnormalities in glucose homeostasis ().
FIGURE 8
lncRNA Target-TF Network Analysis
To further identify the functions of each dysregulated lncRNA, the top 50 differentially expressed lncRNAs and their co-expressed mRNAs pairs were analyzed to conduct the lncRNA-target-TF network (Figure 9). LncRNA NR130973 was predicted to be mainly regulated by these TFs and participates in the Wnt signaling pathway, which has been reported to be related to inflammation and adipose cell function ().
FIGURE 9
Discussion
The adipose tissue is an important endocrine organ that can respond differentially to physiological cues or metabolic stress by releasing endocrine factors that regulate diverse processes (). Abnormal accumulation of body fat is associated with insulin resistance and results in obesity, which is the most frequently observed component of metabolic syndrome (). Thus, adipocyte differentiation regulators have been considered as possible therapeutic interventions for metabolic diseases (; ). Although the role of adipose tissue in parasite infection has been neglected for a long time, some progress has been made in recent years. Trypanosoma cruzi () and Nippostrongylus brasiliensis () were successively shown to inhibit adipocyte differentiation in obese mice that were fed a high-fat diet, providing new insight into the treatment of metabolic diseases. Moreover, a recent study showed that adipose tissue is a major niche for T. brucei, which can adapt its metabolism to the fatty acids that are present in the adipose tissue and improve the host’s fatty acid metabolism (). Since the adipose tissue lies in a lipid-rich environment and most drugs are hydrophilic, it may pose some limitations to the efficacy of drug treatments for parasitic infections (). Therefore, it is paramount to understand the role and mechanisms of parasitic infection in adipogenesis regulation.
This study has examined the effect of infection with the EgPSC on metabolic measures. EgPSC infection resulted in improved glucose tolerance and showed no obvious change of triglyceride, cholesterol, LDL, HDL, and AST/ALT. When compared with the control group, the infected mice showed a decrease in subcutaneous adipose tissue mass and a reduction in adipocyte size, however, the mass of brown adipose tissue had no obvious change post-infection. As a result, this study focused on the subcutaneous adipose tissues. Meanwhile, there was no statistically significant difference of food intake between control and EgPSC infected mice, which indicated that the mass loss of fat tissue should not result from sickness and poor appetite. The results also showed the existence of a M2-like macrophage in infected adipose tissues, which may promote fat lipolysis ().
The present study was the first to address the mRNA and lncRNA expression profiles of adipogenesis in EgPSC-infected mice. Parasitic infection can inhibit the differentiation of subcutaneous adipocytes, which is accompanied by a sharp decline in carbohydrate oxidation and increased fat oxidation. Moreover, microarray analysis revealed that 1052 mRNAs and 220 lncRNAs were differentially expressed in the subcutaneous adipose tissue of infected mice. GO and KEGG analyses showed a metabolic switch in the adipose tissue, which was in agreement with the indirect calorimetry results. Additionally, several differentially expressed mRNAs and lncRNAs were involved in insulin resistance and adipocyte differentiation. Thus, the results provided new insights into the metabolic homeostasis mediated by parasitic infection.
The expression of metabolic enzymes in adipose tissue is strongly related to adipocyte morphology and function () because these proteins control the cellular energy demand by regulating the metabolic pathways. This study found that the expression of key enzymes in “glycolysis,” “TCA cycle,” and “de novo lipogenesis” was markedly downregulated, while the expression of enzymes in the “lipoxygenase” pathway was significantly upregulated after EgPSC infection. Acetyl coenzyme A carboxylase 1 (ACC1) and fatty acid synthase (FAS) are key enzymes in adipose synthesis, and play a crucial role in regulating lipid homeostasis (; ; ; ). Sterol regulatory element-binding protein is considered to be a master regulator of lipogenesis () and is able to increase the activity of ACC1 and FAS (, ). The present study showed that the expression of these proteins was significantly downregulated, which indicated that lipogenesis was inhibited. In contrast, the fatty acid oxidation pathway could be enhanced, since the key enzyme carnitine palmitoyltransferase-1 (CPT1) (, ; ) was significantly downregulated. The downregulation of citrate synthase (CS) and malate dehydrogenase (MDH) indicated that the TCA cycle was weakened, which could contribute to an anti-inflammatory microenvironment in the adipose tissue (). Through indirect calorimetry, we observed a sharp decline in carbohydrate oxidation (P < 0.01) and increased fat oxidation (P > 0.05), which agreed with metabolic reprogramming. Additionally, the change of carbohydrate and fat oxidations may be explained by the newly discovered myokine β-aminoisobutyric acid (BAIBA), which induces the adipose tissue transition to a “beige” phenotype and enhances fatty acids oxidation (). Overall, EgPSC infection might regulate gene expression, thereby reprogramming the host’s metabolism to control adipocyte differentiation.
The lncRNAs are important members of the non-coding RNA (ncRNA) family and have been considered as a new class of non-protein regulators in adipocyte biology (). A total of 220 lncRNAs were found to be differentially expressed, some of which were believed to be closely associated with adipocyte differentiation. For example, WNT5b is targeted in the function of signaling pathways that regulate pluripotency in stem cells (), and it was downregulated in our study. The overexpression of the WNT5b gene was previously shown to promote adipogenesis and enhance adipocytokine gene expression (). Previous studies also substantiated the possible roles of WNT genes in conferring susceptibility to type 2 diabetes (; ). Moreover, the differentially expressed transcript AKT-2, which is required for adipocyte lipid filling (), is involved in apoptotic pathways (). After apoptosis, the adipocytes’ debris are reabsorbed by macrophages, inducing chronic low-grade inflammation and potentially contributing to insulin resistance and type 2 diabetes (). In addition, the downregulated lncRNA H19 (gene ID: 14955) was of particular interest. It has been reported that the silencing of H19 inhibits adipogenesis and inflammation response by upregulating miRNA-130b, and the expression of H19 was shown to be upregulated in atherosclerotic patients (; ). However, the function of these lncRNAs in EgPSC-infected adipose tissue requires further investigation.
Overall, this study proposed that EgPSC infection inhibits adipogenesis along with metabolic reprogramming in the subcutaneous adipose tissue of mice. It has been increasingly recognized that parasite derived excretory-secretory products (ESPs), the key substances employed by helminths to downregulate the host’s anti-infection immunity (), are able to regulate adipogenesis and improve insulin sensitivity in a high fat diet-induced obesity model (). Thus, in the future, we will use the ESPs of EgPSC to mock the infection to study whether and how EgPSC alters host metabolism.
Conclusion
The present study showed that EgPSC infection resulted in a reduction of subcutaneous fat along with a decrease in the size of adipocytes. Consistently, indirect calorimetry revealed the change in energy metabolism that was characterized by a lower CO2 production and O2 consumption, a sharp decline in carbohydrate oxidation, and a slight increase in fat oxidation. Moreover, several differentially expressed lncRNAs and mRNAs were identified to be associated with adipocyte differentiation, insulin resistance, and metabolic reprogramming. Overall, these findings deepened the understanding of parasite-host interplay, metabolic homeostasis mediated by parasitic infection, and provided novel insight into the mechanism of hypertrophic adipocytes in obesity.
Statements
Data availability statement
The microarray data has been deposited into GEO (accession: GSE154979, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? acc=GSE154979).
Ethics statement
The animal study was reviewed and approved by the Laboratory Animal Welfare and Ethics Committee of Xuzhou Medical University, China [SCXK (Su) 2015-0009].
Author contributions
WP, HL, and XY conceived and designed the experiments. YL, JW, TL, YG, and FS performed the experiments. YL, JL, and MD analyzed the data. WP and HL contributed reagents, materials, and analysis tools. YL, WP, and FS wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
Project support was provided in part by the National Natural Science Foundation of China (Nos. 81871670 and 81800718), the Jiangsu Planned Projects for Postdoctoral Research Funds (No. 2019K063), the General Project of Natural Science Foundation of Jiangsu Province, and the Priority Academic Program Development of Jiangsu Higher Education Institutions, the Training Programs of Innovation and Entrepreneurship for College Students in Jiangsu Province (Nos. 202010313035Z and 202010313035). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.01078/full#supplementary-material
References
1
AffinassN.ZhangH.LöhningM.HartmannS.RauschS. (2018). Manipulation of the balance between Th2 and Th2/1 hybrid cells affects parasite nematode fitness in mice.Eur. J. Immunol.481958–1964. 10.1002/eji.201847639
2
AlbertiK.ZimmetP.ShawJ. (2005). The metabolic syndrome—a new worldwide definition.Lancet3661059–1062. 10.1016/S0140-6736(05)67402-8
3
AndersenC. J.MurphyK. E.FernandezM. L. (2016). Impact of obesity and metabolic syndrome on immunity.Adv. Nutr.766–75. 10.3945/an.115.010207
4
CaoJ.SinghS. P.McClungJ. A.JosephG.VanellaL.BarbagalloI.et al (2017). EET intervention on Wnt1, NOV, and HO-1 signaling prevents obesity-induced cardiomyopathy in obese mice.Am. J. Physiol. Heart Circ. Physiol.313H368–H380. 10.1152/ajpheart.00093
5
CintiS. (2018). Adipose organ development and remodeling.Compr. Physiol.81357–1431. 10.1002/cphy.c170042
6
CroweJ.LumbF. E.HarnettM. M.HarnettW. (2017). Parasite excretory-secretory products and their effects on metabolic syndrome.Parasite Immunol.39:e12410. 10.1111/pim.12410
7
DahlmanI.RydénM.BrodinD.GrallertH.StrawbridgeR. J.ArnerP. (2016). Numerous genes in loci associated with body fat distribution are linked to adipose function.Diabetes65433–437. 10.2337/db15-0828
8
DanforthE.Jr. (2000). Failure of adipocyte differentiation causes type II diabetes mellitus?Nat. Genet.26:13. 10.1038/79111
9
De NizM.MeehanG. R.BrancucciN. M. B.MartiM.RotureauB.FigueiredoL. M.et al (2019). Intravital imaging of host-parasite interactions in skin and adipose tissues.Cell Microbiol.21:e13023. 10.1111/cmi.13023
10
DíazÁ (2017). Immunology of cystic echinococcosis (hydatid disease).Br. Med. Bull.124121–133. 10.1093/bmb/ldx033
11
DoonanJ.TarafdarA.PinedaM. A.LumbF. E.CroweJ.KhanA. M.et al (2019). The parasitic worm product ES-62 normalises the gut microbiota bone marrow axis in inflammatory arthritis.Nat. Commun.10:1554. 10.1038/s41467-019-09361-0
12
EberléD.HegartyB.BossardP.FerréP.FoufelleF. (2004). SREBP transcription factors:master regulators of lipid homeostasis.Biochimie86839–848. 10.1016/j.biochi.2004.09.018
13
EnginA. (2019). The definition and prevalence of obesity and metabolic syndrome.Adv. Exp. Med. Biol.9601–17. 10.1007/978-3-319-48382-5_1
14
FeakinsR. M. (2016). Obesity and metabolic syndrome: pathological effects on the gastrointestinal tract.Histopathology68630–640. 10.1111/his.12907
15
FerrèF.ColantoniA.Helmer-CitterichM. (2016). Revealing protein-lncRNA interaction.Brief Bioinform.17106–116. 10.1093/bib/bbv031
16
FischerA. H.JacobsonK. A.RoseJ.ZellerR. (2008). Hematoxylin and eosin staining of tissue and cell sections.CSH Protoc.2008:pdb.prot4986. 10.1101/pdb.prot4986
17
FullertonM. D.GalicS.MarcinkoK.SikkemaS.PulinilkunnilT.ChenZ. P.et al (2013). Single phosphorylation sites in Acc1 and Acc2 regulate lipid homeostasis and the insulin-sensitizing effects of metformin.Nat. Med.191649–1654. 10.1038/nm.3372
18
GijsH. (2017). The metabolic phenotype in obesity: fat mass, body fat distribution, and adipose tissue function.Obes Facts10207–215. 10.1159/000471488
19
GomesP. S.BhardwajJ.Rivera-CorreaJ.Freire-De-LimaC. G.MorrotA. (2016). Immune escape strategies of malaria parasites.Front. Microbiol.7:1617. 10.3389/fmicb.2016.01617
20
GonzálezF. B.VillarS. R.ToneattoJ.PaciniM. F.MárquezJ.D’AttilioL.et al (2019). Immune response triggered by Trypanosoma cruzi infection strikes adipose tissue homeostasis altering lipid storage, enzyme profile and adipokine expression.Med. Microbiol. Immunol.208651–666. 10.1007/s00430-018-0572-z
21
HanY.MaJ.WangJ.WangL. (2018). Silencing of H19 inhibits the adipogenesis and inflammation response in ox-LDL-treated Raw264.7 cells by up-regulating miR-130b.Mol. Immunol.93107–114. 10.1016/j.molimm.2017.11.017
22
HarnettW. (2014). Secretory products of helminth parasites as immunomodulators.Mol. Biochem. Parasitol.195130–136. 10.1016/j.molbiopara.2014.03.007
23
HellardE.FouchetD.VavreF.PontierD. (2015). Parasite-parasite interactions in the wild: how to detect them.Trends Parasitol.31640–652. 10.1016/j.pt.2015.07.005
24
HortonJ. D.ShahN. A.WarringtonJ. A.AndersonN. N.ParkS. W.BrownM. S.et al (2003a). Combined analysis of oligonucleotide microarray data from transgenic and knockout mice identifies direct SREBP target genes.Proc. Natl. Acad. Sci. U.S.A.10012027–12032. 10.1073/pnas.1534923100
25
HortonJ. D.ShimomuraI.IkemotoS.BashmakovY.HammerR. E. (2003b). Overexpression of sterol regulatory element-binding protein-1a in mouse adipose tissue produces adipocyte hypertrophy, increased fatty acid secretion, and fatty liver.J. Biol. Chem.27836652–36660. 10.1074/jbc.M306540200
26
HussaartsL.García-TardónN.van BeekL.HeemskerkM. M.HaeberleinS.van der ZonG. C.et al (2015). Chronic helminth infection and helminth-derived egg antigens promote adipose tissue M2 macrophages and improve insulin sensitivity in obese mice.FASEB J.293027–3039. 10.1096/fj.14-266239
27
JaisA.BrüningJ. C. (2017). Hypothalamic inflammation in obesity and metabolic disease.J. Clin. Invest.12724–32. 10.1172/JCI88878
28
KanazawaA.TsukadaS.SekineA.TsunodaT.TakahashiA.KashiwagiA.et al (2004). Association of the gene encoding wingless-type mammary tumor virus integration-site family member 5B (WNT5B) with type 2 diabetes.Am. J. Hum. Genet.75832–843. 10.1086/425340
29
LallyJ. S. V.GhoshalS.DePeraltaD. K.MoavenO.WeiL.MasiaR.et al (2019). Inhibition of acetyl-CoA carboxylase by phosphorylation or the inhibitor ND-654 suppresses lipogenesis and hepatocellular carcinoma.Cell Metab.29174–182.e5. 10.1016/j.cmet.2018.08.020
30
LiC.YangD.CaoX.WangF.JiangH.GuoH.et al (2016). LFG-500, a newly synthesized flavonoid, attenuates lipopolysaccharide-induced acute lung injury and inflammation in mice.Biochem. Pharmacol.11357–69. 10.1016/j.bcp.2016.05.007
31
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method.Methods25402–408. 10.1006/meth.2001.1262
32
MaegawaH.MaedaS. (2006). Role of TFAP2B and Wnt5B in adipose tissue.Nihon Rinsho.64(Suppl. 9)244–248.
33
MaierT.LeibundgutM.BoehringerD.BanN. (2010). Structure and function of eukaryotic fatty acid synthases.Q. Rev. Biophys.43373–422. 10.1017/S0033583510000156
34
MárquezS.FernándezJ. J.ManceboC.Herrero-SánchezC.AlonsoS.SandovalT. A.et al (2019). Tricarboxylic acid cycle activity and remodeling of glycerophosphocholine lipids support cytokine induction in response to fungal patterns.Cell Rep.27525–536.e4. 10.1016/j.celrep.2019.03.033
35
McGarryJ. D.MillsS. E.LongC. S.FosterD. W. (1983). Observations on the affinity for carnitine, and malonyl-CoA sensitivity, of carnitine palmitoyltransferase I in animal and human tissues.Biochem. J.21421–28. 10.1042/bj2140021
36
McGarryJ. D.WoeltjeK. F.KuwajimaM.FosterD. W. (1989). Regulation of ketogenesis and the renaissance of carnitine palmitoyltransferase.Diabetes Metab. Rev.5271–284. 10.1002/dmr.5610050305
37
Mota de SáP.RichardA. J.HangH.StephensJ. M. (2017). Transcriptional regulation of adipogenesis.Compr. Physiol.7635–674. 10.1002/cphy.c160022
38
NguyenK. D.QiuY.CuiX.GohY. P. S.MwangiJ.DavidT.et al (2011). Alternatively activated macrophages produce catecholamines to sustain adaptive thermogenesis.Nature480104–108. 10.1038/nature10653
39
PanW.XuH. W.HaoW. T.SunF. F.QinY. F.HaoS. S.et al (2018). The excretory-secretory products of Echinococcus granulosus protoscoleces stimulated IL-10 production in B cells via TLR-2 signaling.BMC Immunol.19:29. 10.1186/s12865-018-0267-7
40
PanW.ZhouH. J.ShenY. J.WangY.XuY. X.HuY.et al (2013). Surveillance on the status of immune cells after Echinnococcus granulosus protoscoleces infection in Balb/c mice.PLoS One8:e59746. 10.1371/journal.pone.0059746
41
PeónA. N.Ledesma-SotoY.TerrazasL. I. (2016). Regulation of immunity by Taeniids: lessons from animal models and in vitro studies.Parasite Immunol.38124–135. 10.1111/pim.12289
42
PéronnetF.MassicotteD. (1991). Table of nonprotein respiratory quotient: an update.Can. J. Sport Sci.1623–29.
43
PurdomT.KravitzL.DokladnyK.MermierC. (2018). Understanding the factors that effect maximal fat oxidation.J. Int. Soc. Sports Nutr.15:3. 10.1186/s12970-018-0207-1
44
RamananD.BowcuttR.LeeS. C.TangM. S.KurtzZ. D.DingY.et al (2016). Helminth infection promotes colonization resistance via type 2 immunity.Science352608–612. 10.1126/science.aaf3229
45
Sanchez-GurmachesJ.Martinez CalejmanC.JungS. M.LiH.GuertinD. A. (2019). Brown fat organogenesis and maintenance requires AKT1 and AKT2.Mol. Metab.2360–74. 10.1016/j.molmet.2019.02.004
46
SchejaL.HeerenJ. (2019). The endocrine function of adipose tissues in health and cardiometabolic disease.Nat. Rev. Endocrinol.15507–524. 10.1038/s41574-019-0230-6
47
SchilperoortM.van DamA. D.HoekeG.ShabalinaI. G.OkoloA.HanyalogluA. C.et al (2018). The GPR120 agonist TUG-891 promotes metabolic health by stimulating mitochondrial respiration in brown fat.EMBO Mol. Med.10:e8047. 10.15252/emmm.201708047
48
SuC. W.ChenC. Y.LiY.LongS. R.MasseyW.KumarD. V.et al (2018). Helminth infection protects against high fat diet-induced obesity via induction of alternatively activated macrophages.Sci. Rep.8:4607. 10.1038/s41598-018-22920-7
49
SunL.LinJ. D. (2019). Function and mechanism of long noncoding RNAs in adipocyte biology.Diabetes68887–896. 10.2337/dbi18-0009
50
SunY.ZhongL.HeX.WangS.LaiY.WuW.et al (2019). LncRNA H19 promotes vascular inflammation and abdominal aortic aneurysm formation by functioning as a competing endogenous RNA.J. Mol. Cell. Cardiol.13166–81. 10.1016/j.yjmcc.2019.04.004
51
SwarbrickM. M.EvansD. S.ValleM. I.FavreH.WuS. H.NjajouO. T.et al (2011). Replication and extension of association between common genetic variants in SIM1 and human adiposity.Obesity (Silver Spring).192394–2403. 10.1038/oby.2011.79
52
TanianskiiD. A.JarzebskaN.BirkenfeldA. L.O’SullivanJ. F.RodionovR. N. (2019). Beta-aminoisobutyric acid as a novel regulator of carbohydrate and lipid metabolism.Nutrients11:E524. 10.3390/nu11030524
53
TanowitzH. B.SchererP. E.MotaM. M.FigueiredoL. M. (2017). Adipose tissue: a safe haven for parasites?Trends Parasitol.33276–284. 10.1016/j.pt.2016.11.008
54
TrindadeS.Rijo-FerreiraF.CarvalhoT.Pinto-NevesD.GueganF.Aresta-BrancoF.et al (2016). Trypanosoma brucei parasites occupy and functionally adapt to the adipose tissue in mice.Cell Host Microbe19837–848. 10.1016/j.chom.2016.05.002
55
TurerE. E.San MiguelM.WangK. W.McAlpineW.OuF.LiX.et al (2018). A viable hypomorphic Arnt2 mutation causes hyperphagic obesity, diabetes and hepatic steatosis.Dis. Model Mech.11:dmm035451. 10.1242/dmm.035451
56
van der ZandeH. J. P.Zawistowska-DeniziakA.GuigasB. (2019). Immune regulation of metabolic homeostasis by helminths and their molecules.Trends Parasitol.35795–808. 10.1016/j.pt.2019.07.014
57
WallaceM.GreenC. R.RobertsL. S.LeeY. M.McCarvilleJ. L.Sanchez-GurmachesJ.et al (2018). Enzyme promiscuity drives branched-chain fatty acid synthesis in adipose tissues.Nat. Chem. Biol.141021–1031. 10.1038/s41589-018-0132-2
58
YangZ.GrinchukV.SmithA.QinB.BohlJ. A.SunR.et al (2013). Parasitic nematode-induced modulation of body weight and associated metabolic dysfunction in mouse models of obesity.Infect. Immun.811905–1914. 10.1128/IAI.00053-13
59
YuA.WangY.YinJ.ZhangJ.CaoS.CaoJ.et al (2018). Microarray analysis of long non-coding RNA expression profiles in monocytic myeloid-derived suppressor cells in Echinococcus granulosus-infected mice.Parasit Vectors11:327. 10.1186/s13071-018-2905-6
60
ZhaoX. Y.LiS.DelPropostoJ. L.LiuT.MiL.PorscheC.et al (2018). The long noncoding RNA Blnc1 orchestrates homeostatic adipose tissue remodeling to preserve metabolic health.Mol. Metab.1460–70. 10.1016/j.molmet
Summary
Keywords
Echinococcus granulosus, protoscoleces, adipocyte, adipogenesis, long non-coding RNA, microarray analysis
Citation
Lu Y, Liu H, Yang X, Liu J, Dai M, Wu J, Guo Y, Luo T, Sun F and Pan W (2020) Microarray Analysis of lncRNA and mRNA Reveals Enhanced Lipolysis Along With Metabolic Remodeling in Mice Infected With Larval Echinococcus granulosus. Front. Physiol. 11:1078. doi: 10.3389/fphys.2020.01078
Received
20 May 2020
Accepted
05 August 2020
Published
21 August 2020
Volume
11 - 2020
Edited by
Naim Akhtar Khan, Université de Bourgogne, France
Reviewed by
Scott Widenmaier, University of Saskatchewan, Canada; P. Trayhurn, University of Liverpool, United Kingdom
Updates
Copyright
© 2020 Lu, Liu, Yang, Liu, Dai, Wu, Guo, Luo, Sun and Pan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Pan, panwei525@126.com
†These authors have contributed equally to this work
This article was submitted to Lipid and Fatty Acid Research, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.