PGC-1α-mediated changes in phospholipid profiles of exercise-trained skeletal muscle[S]

Exercise training influences phospholipid fatty acid composition in skeletal muscle and these changes are associated with physiological phenotypes; however, the molecular mechanism of this influence on compositional changes is poorly understood. Peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α), a nuclear receptor coactivator, promotes mitochondrial biogenesis, the fiber-type switch to oxidative fibers, and angiogenesis in skeletal muscle. Because exercise training induces these adaptations, together with increased PGC-1α, PGC-1α may contribute to the exercise-mediated change in phospholipid fatty acid composition. To determine the role of PGC-1α, we performed lipidomic analyses of skeletal muscle from genetically modified mice that overexpress PGC-1α in skeletal muscle or that carry KO alleles of PGC-1α. We found that PGC-1α affected lipid profiles in skeletal muscle and increased several phospholipid species in glycolytic muscle, namely phosphatidylcholine (PC) (18:0/22:6) and phosphatidylethanolamine (PE) (18:0/22:6). We also found that exercise training increased PC (18:0/22:6) and PE (18:0/22:6) in glycolytic muscle and that PGC-1α was required for these alterations. Because phospholipid fatty acid composition influences cell permeability and receptor stability at the cell membrane, these phospholipids may contribute to exercise training-mediated functional changes in the skeletal muscle.

angiogenesis, and a fi ber-type switch to oxidative fi bers ( 23,24 ). Previously, we generated mice that overexpressed PGC-1 ␣ -b in skeletal muscle but not in heart ( 25,26 ). PGC-1 ␣ -b, whose N terminus is different from that of PGC-1 ␣ -a protein, is the predominant PGC-1 ␣ isoform in skeletal muscles that is expressed in response to exercise ( 26 ). Overexpression of PGC-1 ␣ -b promoted fi ber-type switch, mitochondrial biogenesis, and exercise capacity, increased the expression of fatty acid transporters, and enhanced angiogenesis and oxygen utilization kinetics in skeletal muscle ( 25,27 ). Recently, we also found that PGC-1 ␣ -b induced branched chain amino acid metabolism, which might also be involved in endurance capacity ( 28 ). Because exercise training induces changes in skeletal muscle phospholipid species, it is possible that PGC-1 ␣ -b may be the underlying mechanism of induction. In this study, to understand how exercise training induces changes in phospholipid species, we performed lipidomics of skeletal muscle from genetically modifi ed mice that overexpressed PGC-1 ␣ -b and mice that carry KO alleles of PGC-1 ␣ in skeletal muscle. We found that PGC-1 ␣ is involved in exercise training-induced changes in skeletal muscle phospholipid species.

Genetically modifi ed mice
The methods for generating transgenic mice overexpressing PGC-1 ␣ -b in skeletal muscle (PGC-1 ␣ -Tg mice) were described previously ( 26 ). The promoter for human ␣ -skeletal actin, provided by Drs. E. C. Hardeman and K. L. Guven (Children's Medical Research Institute, Australia) was used to express PGC-1 ␣ -b in skeletal muscle. The transgenic mice (heterozygotes, BDF 1 background) and WT C57BL6 mice were crossed and female 10-13-week-old offspring (heterozygote and WT, from the same litter) were used for the experiments. To generate skeletal muscle-specifi c PGC-1 ␣ KO mice (muscle PGC-1 ␣ -KO mice), we inactivated PGC-1 ␣ expression in skeletal muscles by crossing mice carrying a fl oxed PGC-1 ␣ allele with mice transgenic for the human ␣ -skeletal actin promoter driven-Cre transgenic. PGC-1 ␣ fl ox/fl ox mice were obtained from the Jackson Laboratory (Bar Harbor, ME) ( 29,30 ). Mice were maintained in a 12 h light/dark cycle at 22°C and were fed a normal chow diet ad libitum (CE-2; CLEA Japan, Tokyo, Japan). Mice were cared for in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and our institutional guidelines. All animal experiments were conducted with the approval of the Institutional Animal Care and Use Committee of the University of Shizuoka (number 135024).

Voluntary exercise training
Male (9-week-old) muscle PGC-1 ␣ -KO mice and control PGC-1 ␣ fl ox/fl ox mice were randomly assigned to one of two experimental groups: the sedentary control group or the training group. Mice assigned to training were housed individually in cages (22 × 9 × 8 cm) equipped with a running wheel (20 cm diameter; Shinano Co., Tokyo, Japan) for 5 weeks. The running wheel was equipped with a tachometer to determine the total running distance. Sedentary mice were housed in cages without a running wheel. exercise training (4)(5)(6)(7) can both infl uence the composition of skeletal muscle membrane fatty acids, changes in phospholipid fatty acids may be involved in diet-induced or exercise training-induced physiological adaptation of the skeletal muscle. This effect on skeletal muscle adaptation may ultimately infl uence its function. Previous studies examined the effect of endurance training on the molecular species of skeletal muscle phospholipids. Exercise training increases phosphatidylcholine (PC) in muscle ( 8 ), and the effects of exercise on numerous phospholipid species in skeletal muscle have been examined in rats fed a standard chow diet ( 9 ) or a high-fat diet ( 10 ). With chow feeding, exercise training increased the abundance of two PC species ( 9 ) ( 11 ) performed lipidomic analysis of skeletal muscle from mice subjected to chronic exercise training and high-fat diet using imaging MS (IMS) and TLC-Blot-MALDI-IMS. They found that PC (16:0/18:2), PC (18:0/22:6), and SM (d18:1/16:0) were chronic exercise training-induced lipids and, in contrast, PC (18:0/20:4) and SM (d18:1/24:1) were high-fat diet-induced lipids. The largest differences were also observed when comparing oxidative and glycolytic muscles, such as a decrease in plasmenyl-PE (16:0/20:4) and an increase in PE (18:0/22:6) in the oxidative muscle ( 9 ). Although evidence has shown that exercise training induces changes in skeletal muscle phospholipid species, it is not fully understood how exercise training induces these changes, or what roles these phospholipids play in the functional changes of exercise-trained skeletal muscle.

IMS analysis
The localization of each lipid was analyzed by IMS. The developed TLC plates were transferred to a polyvinylidene difl uoride (PVDF) membrane by the TLC-blot method, as described previously ( 31 ), and transferred PVDF membranes were attached to the MALDI target plate for IMS analyses. For TLC-blot-imaging analyses, we used a QSTAR Elite high-performance hybrid quadrupole TOF mass spectrometer (Applied Biosystems, Foster City, CA). The laser irradiated 500 times per position on the PVDF membrane. We set the spatial resolution to 400 m. All analyses were performed in the positive ion mode within the mass ranges of m/z 400-1,200, with 2,5-dihydroxybenzoic acid at 50 mg/ml as a matrix. The ion images were constructed using BioMap software (Novartis, Basel, Switzerland).
The tissue blocks (gastrocnemius) were rapidly frozen in isopentane cooled by liquid nitrogen. Transverse cross-sections of 10 m were made with a cryostat (Leica, CM1510; Germany) at Ϫ 20°C. For positive-ion mode, 2,5-dihydroxybenzoic acid of 50 mg/ml in methanol/water (7:3, v/v) was uniformly sprayed over the muscle tissue sections with a 0.2 mm nozzle caliber airbrush (Procon Boy FWA Platinum; Mr. Hobby, Tokyo, Japan). We used a MALDI TOF/TOF-type instrument, the Ultrafl ex II (Bruker Daltonics, Billerica, MA). The laser irradiated 200 times per position. All pixel sizes of imaging were 100 m. The MS parameters were set to obtain the highest sensitivity with m/z values in the range of 400-1,000. The ion images were constructed using Flex-Imaging software (Bruker Daltonics). Normalization by total ion current was performed using the same software.

Quantitative RT-PCR
RNA preparation methods and quantitative (q)RT-PCR were performed as described previously ( 26 ). The mouse-specifi c primer pairs used are shown in supplementary Table 2.

Statistical analysis
For lipidomic analyses, the detected peaks were aligned according to the m/z value and normalized retention time using Signpost MS (Reifycs, Tokyo, Japan). After applying autoscaling, mean-centering, and scaling by standard deviation on a per-peak basis as pretreatment, a hierarchical clustering analysis and a principal component analysis (PCA) were conducted using JMP version 11 (SAS Institute, Cary, NC). In hierarchical clustering analysis, the resulting data sets of each genotype were clustered using Euclidean distance with Ward's method ( 32 ). The relative area value of each peak was calculated and used for the comparison between the PGC-1 ␣ -Tg and WT groups. In PCA, a score plot of the fi rst and second principal components was generated. Statistical hypothesis testing of factor loading in PCA was performed to select species that had a statistically signifi cant correlation to the principal component score. The P value was calculated as reported previously ( 33 ). A change trend was defi ned at P < 0.05. Furthermore, a Bonferroni adjustment was applied to determine the level of signifi cance for multiple testing (the adjusted ␣ = 0.05/97 = 0.000515 for retention time 14-29 min peaks and ␣ = 0.05/80 = 0.000625 for retention time 30-41 min peaks ).

LC/MS
Lipids were extracted from the extensor digitorum longus (EDL) and the soleus. Frozen muscle was homogenized and powdered in liquid nitrogen. Total lipids were extracted from homogenates with 1 ml chloroform/methanol (2:1, v/v with 0.2 mg/ml butyl hydroxyl toluene) overnight. For targeted analysis of phospholipids, 0.05 mg/ml PC (17:0/17:0) was added to the chloroform/methanol for use as an internal standard. The lipid fractions were evaporated to dryness under vacuum. Samples were reconstituted in an equal volume of acetonitrile/isopropanol/ water (65:30:5, v/v/v). Ten microliters of samples were injected onto the LC/MS system. Lipidomic analysis was performed using a Q-Exactive TM benchtop orbitrap mass spectrometer (Thermo Fisher Scientifi c, Waltham, MA) equipped with an electrospray source ionization probe and an autosampler, Accela quaternary HPLC pump (Thermo Fisher Scientifi c). For LC analysis, an Acquity UPLC CSH C18 column (1.7 m, 2.1 × 150 mm; Waters, Milford, MA) was used. Mobile phase A consisted of water/acetonitrile (60:40, v/v) and mobile phase B consisted of isopropanol/acetonitrile (90:10, v/v). Both mobile phases, A and B, were supplemented with 10 mM ammonium formate and 0.1% formic acid. The fl ow rate was 0.4 ml/min. The gradient was as follows: 10% B at 0 min, 10% B at 2 min, 50% B at 8 min, 75% B at 20 min, 90% B at 55 min, and 10% B at 60 min. For MS analysis, the spray voltage was set to 3.5 kV, the capillary temperature was set to 350°C, the S-lens radio frequency (RF) level was set to 50 , and heater temperature was set to 300°C. The sheath gas fl ow rate was set to 40, and the auxiliary gas fl ow rate was set to 10. These conditions were applied to both positive and negative ionization modes. All samples were analyzed by both positive and negative ionization mode acquiring full scan MS, and the scan range was between m/z 120 and 1,200. To identify lipid species, m/z and retention time were used to obtain MS/MS data of the target peak by tandem mass spectrometric analysis. Using the LIPID MAPS online MS tool, the MS/MS data obtained were searched against a database of glycerophospholipid (http://www.lipidmaps.org/tools/ms/GP_ prod_search.html) and glycerolipid (http://www.lipidmaps. org/tools/ms/GL_prod_search.html) precursor/product ions.
Targeted LC/MS/MS analysis was performed using an LCMS-8040 triple quadrupole mass spectrometer (Shimadzu, Kyoto, Japan) in positive ionization mode equipped with an electrospray source ionization probe, LC-30AD binary pump (Shimadzu), SIL-30AC auto sampler (Shimadzu), and CTO-20AC column oven (Shimadzu). For HPLC analysis, an Accucore RP-MS C18 column (2.6 m, 2.1 × 50 mm, Thermo Fisher Scientifi c) was used. The composition of the mobile phase and the fl ow rate were the same as above. The gradient was as follows: 40% B at 0 min, 40% B at 2 min, 52% B at 8 min, 60% B at 20 min, 100% B at 25 min, and 40% B at 30 min. For MS analysis, the nebulizer gas fl ow was set to 3.0 l/min, the drying gas fl ow was set to 15.0 l/ min, the desolvation line temperature was set to 250°C, the heat block temperature was set to 400°C, and the collision-induced dissociation gas was set to 17 kPa . Target analysis of phospholipid species was operated in selected reaction monitoring (SRM) mode. The details of the conditions for SRM for each phospholipid are shown in supplementary Table 1. The relative peak area for each species was normalized by the peak area of internal standard and muscle weight.

TLC analysis
For quantitative analysis of each lipid, we used conventional TLC. Briefl y, total lipid extracts from each gastrocnemius dissolved in chloroform/methanol (2:1, v/v) were manually applied as 5 mm wide spots to silica gel 60 HPTLC plates (Merck, Darmstadt, To identify the lipid species that contribute to PGC-1 ␣driven changes in lipid profi le, PCA was performed using lipid species in the phospholipid and TG fractions. In PCA of the phospholipid fraction, the fi rst principal component separated the WT and PGC-1 ␣ -Tg mice (x axis), and the differences were obvious in EDL ( Fig. 1C ). Statistical hypothesis testing for factor loading in the fi rst and second principal components was performed, and the lipid species that showed statistically signifi cant differences ( P < 0.05) are shown in Fig. 1D and Table 1 . The phospholipids, 12 PC and 6 PE, were identifi ed as signifi cant in the fi rst principal component, and 3 PC and 1 SM were identifi ed in the second principal component. In PCA of the TG fraction, the fi rst principal component separated EDL and soleus, and the second principal component separated soleus from WT and PGC-1 ␣ -Tg mice (data not shown). Fifty-fi ve TGs were detected as signifi cant in statistical hypothesis testing for factor loading in the fi rst principal component; however, none were found to be signifi cant in the second principal component ( Table 2 ). The results indicated that TG molecular species might not be major factors contributing to the alteration of lipid profi le by PGC-1 ␣ .

Targeted analysis of phospholipids altered by PGC-1 ␣ overexpression
To measure the relative amount of phospholipids altered by PGC-1 ␣ overexpression in EDL and soleus, 11 PC and 5 PE, identifi ed as signifi cant in the fi rst principal component, were quantitatively measured by SRM using quadrupole MS ( Fig. 2 ). PE (p-16:0/22:6) was not determined in these samples . In EDL, 3.5-and 1.9-fold increases were observed in PC Other data were analyzed by one-way ANOVA. In case of signifi cant differences, each group was compared with the other groups by a Student's t -test (JMP, version 11). Values are shown as the mean ± SE.

Differences between lipid profi les of glycolytic and oxidative fi ber and the impact of PGC-1 ␣
To examine the differences between the lipid profi les of glycolytic and oxidative muscle fi bers such as EDL and soleus, and to determine the impact of PGC-1 ␣ on these profi les, lipidomic analyses were performed using highresolution LC/MS that allows for accurate identifi cation of lipid species. Figure 1A shows PCA scatter plots of the samples. The fi rst principal component effectively and distinctly separated the mice based on muscle fi ber type (x axis), and the second principal component separated the mice based on the genotype (y axis). The results suggested that overexpression of PGC-1 ␣ in the skeletal muscle caused a signifi cant change in the overall lipid profi les of the muscle. However, the lipid profi le of EDL from PGC-1 ␣ -Tg mice, which showed oxidative characteristics, was different from the profi le of the originally oxidative muscle, such as the soleus. In the loading plot of this PCA ( Fig.  1B ), lipid species with chromatographic retention times between 14 and 29 min contributed to PGC-1 ␣ -driven alterations in the lipid profi le. On the other hand, lipid species having chromatographic retention times between 30 and 41 min contributed to the differences in lipid profi les between the EDL and soleus. Because our preliminary study showed that the fractions having chromatographic retention times between 14 and 29 min contained many phospholipid species, this fraction was termed the phospholipid fraction. The fraction having chromatographic retention times between 30 and 41 min contained many TG species and was termed the TG fraction. Exercise training-induced changes in lipid profi les and the role of PGC-1 ␣ PGC-1 ␣ altered lipid profi les in skeletal muscles. Notably, several PC and PE species signifi cantly increased in EDL overexpressing PGC-1 ␣ , as shown in Fig. 2 . To determine whether exercise training-induced changes in the PC and PE profi les were similar to those observed in PGC-1 ␣ -Tg mice, and whether PGC-1 ␣ was involved in these changes, the amounts of PC and PE species in sedentary and trained muscle were measured in muscle PGC-1 ␣ -KO mice and control PGC-1 ␣ fl ox/fl ox mice. Expression of PGC-1 ␣ mRNA was not detected in muscle of PGC-1 ␣ -KO mice (supplementary Fig. 1). The averages of the daily running distances were 6.8 ± 0.6 km/day and 6.4 ± 0.5 km/day in PGC-1 ␣ fl ox/fl ox mice and muscle PGC-1 ␣ -KO mice, respectively ( P = 0.7704). As shown in

Lipid analyses by TLC and determination of PC (18:0/22:6) and PE (18:0/22:6) by TLC-blot-MALDI-IMS and IMS-based histological examination
Changes in the lipid composition of muscle might be due to an increase in membrane fraction, such as plasma membrane and mitochondrial membrane, because PGC-1 ␣ stimulated mitochondrial biogenesis. To determine whether the total amounts of PC and PE were altered by PGC-1 ␣ overexpression, lipids extracted from gastrocnemius of WT and PGC-1 ␣ -Tg mice were separated by TLC ( Fig. 3A ). TG increased in the muscle from PGC-1 ␣ -Tg mice ( P = 0.055), but diacylglycerol (DG) content did not. In addition to PC and PE, SM, PI, and phosphatidylserine (PS) were detected with phospholipid separation. The amount of PI was slightly higher in the muscle from PGC-1 ␣ -Tg mice ( P = 0.065); whereas, there were no signifi cant differences between the genotypes in any other phospholipid classes.
TLC-blot-MALDI-IMS was used to determine the quantities of PC (   from the PGC-1 ␣ fl ox/fl ox mice, but also from the muscle PGC-1 ␣ -KO mice, suggesting that exercise training-induced changes of these lipid species did not involve PGC-1 ␣ . In soleus, exercise-induced increases were not observed in either group of mice.

PGC-1 ␣ -mediated changes in acyltransferase expression in the skeletal muscle
Because glycerophospholipids and TG are formed by de novo synthesis in the Kennedy pathway, and because the former are modifi ed by the Lands cycle remodeling pathway using acyl-CoAs as donors (supplementary Fig. 2) ( 36 ), PGC-1 ␣ -mediated changes in skeletal muscle lipid profi les may be caused by changes in the expression of acyltransferases involved in the Kennedy pathway or the Lands cycle. To test this possibility, we measured the expression of glycerophosphate acyltransferase (GPAT) and lysophosphatidic acid acyltransferase (LPAAT), which are components of the Kennedy pathway, and lysophosphatidylcholine acyltransferase (LPCAT), lysophosphatidylglycerol LPEAT did not change in EDL, suggesting that PGC-1 ␣mediated increase in PC (18:0/22:6) and PE (18:0/22:6) in EDL cannot be explained by the expression of these genes.

DISCUSSION
It has been reported that exercise training modifi es phospholipid species in the skeletal muscles ( 4-11 ); however, the molecular mechanisms involved in these modifi cations acyltransferase, lysocardiolipin acyltransferase, lysophosphatidylethanolamine acyltransferase (LPEAT), lysophosphatidylinositol acyltransferase, and acyl glycerophosphate acyltransferase (AGPAT), which are components of the Lands cycle, using qRT-PCR ( Fig. 5 ). Although the expressions of GPAT1, GPAT3, LPAAT2, and lysophosphatidylglycerol acyltransferase 1 in EDL and GPAT1, GPAT3, and LPCAT4 in soleus signifi cantly increased in the muscle overexpressing PGC-1 ␣ , the expressions of the genes involved in PC and PE remodeling, AGPAT5, LPCAT, and    of TG carrying 22:6 fatty acid moiety to the difference in the muscle from WT and PGC-1 ␣ -Tg mice ( Table 2 ). These results suggest that the Kennedy pathway may not be involved in PGC-1 ␣ -mediated changes of PC (18:0/22:6) and PE (18:0/22:6) content in the muscle.
In addition to the acyl-CoA selectivity of GPAT and LPAAT, the Lands cycle is also involved in fatty acid diversity of glycerophospholipids. The Lands cycle provides a route for acyl remodeling to alter fatty acid composition at sn -2 of glycerophospholipids derived from the Kennedy pathway through the concerted actions of phospholipase A2 and lysophospholipid acyltransferase (LPLAT). Recently, several LPLATs were discovered and the substrate specifi city was uncovered ( 36 ). For example, LPCAT3 showed higher activities with polyunsaturated fatty acyl-CoAs, 20:4-CoA and 18:2-CoA, than with saturated fatty acyl-CoA; whereas LPCAT4 and LPEAT1 had a clear preference for 18:1-CoA ( 47 ). These substrate preferences of LPLATs may explain the diversity in membrane glycerophospholipids, which vary among tissues and can change in response to external stimuli. However, the acyltransferases, which prefer 22:6-CoA as a substrate, have not been discovered in the Lands cycle. In the present study, measurements of LPLAT expression using qRT-PCR did not explain the changes of PC (18:0/22:6) and PE (18:0/22:6) content in the muscle. Further studies are required to determine the mechanisms of PGC-1 ␣ -mediated modifi cation of the phospholipid profi le in skeletal muscle.
In skeletal muscle, several studies have shown the relation between phospholipid species and physiological or pathophysiological phenotypes. For instance, a higher proportion of total n-3 polyunsaturated fatty acid, particularly docosahexaenoic acid [22:6 (n-3)], in skeletal muscle phospholipids has been associated with a lower fasting plasma glucose level in young children ( 48 ) and higher insulin action in rats ( 49 ). Docosahexaenoic acid [22:6 (n-3)] may have another function after being released from membrane glycerophospholipids by phospholipase A2, because D series resolvins and protectins derived from docosahexaenoic acid [22:6 (n-3)] have anti-infl ammatory and immunoregulatory properties (50)(51)(52)(53)(54). Recently, protectin DX, also derived from docosahexaenoic acid [22:6 (n-3)], was reported to stimulate the release of the prototypic myokine, interleukin-6, from skeletal muscle and thereby initiated a myokine-liver signaling axis, which blunted hepatic glucose production ( 55 ). Protectin DX also activates AMP-activated protein kinase ( 55 ), which stimulates glucose and lipid metabolism in skeletal muscle. From these facts, some of the health benefi ts of exercise training may be explained by PGC-1 ␣ -mediated alteration of skeletal muscle phospholipid profi les. On the other hand, a decrease in docosahexaenoic acid [22:6 (n-3)] from the skeletal muscle phospholipid fraction was observed in dystrophic mdx mice ( 56 ), and PGC-1 ␣ ameliorates muscular dystrophy in mdx mice ( 57 ). Although docosahexaenoic acid [22:6 (n-3)] has not been analyzed in the phospholipid fraction of dystrophic mdx mice overexpressing PGC-1 ␣ , it is likely that PGC-1 ␣ infl uences muscular dystrophy by altering phospholipid composition. PGC-1 ␣ expression, like the previously recognized changes in mitochondrial enzymes with training in the PCG-1 ␣ mice ( 37 ).
Previously, exercise training has been shown to increase in proportion to docosahexaenoic acid [22:6 (n-3)] in human skeletal muscle phospholipids ( 4,7 ). For example, Andersson et al. ( 4 ) demonstrated that both the distribution of type I fi bers in the skeletal muscle and the estimated physical activity level of humans were positively correlated with the percentage of docosahexaenoic acid [22:6 (n-3)], suggesting a relation between the proportion of docosahexaenoic acid [22:6 (n-3)] in human skeletal muscle phospholipids and the increase in endurance capacity with fi ber type switch from glycolytic to oxidative. In the rat skeletal muscle, exercise increased the amount of PC (18:0/22:6) in EDL ( 11 ), and the content of PE (18:0/22:6) was much higher in red vastus lateralis than in white vastus lateralis ( 9 ). In the present study, we found that among phospholipids containing 22:6 fatty acids, PC (18:0/22:6) and PE (18:0/22:6) were increased by exercise training and overexpression of PGC-1 ␣ in the skeletal muscle. Furthermore, muscle PGC-1 ␣ was required for exerciseinduced increase in PC (18:0/22:6) and PE (18:0/22:6) in EDL. These results suggested that exercise training-induced expression of PGC-1 ␣ changed the muscle fi ber type from glycolytic to oxidative concomitantly with increase in PC (18:0/22:6) and PE (18:0/22:6). In soleus, the involvement of PGC-1 ␣ in these changes was not observed, probably because soleus originally has characteristics of red muscle.
Both the Kennedy pathway and the Lands cycle contribute to the fatty acid composition of glycerophospholipids (supplementary Fig. 2) ( 38,39 ). The Kennedy pathway begins with the stepwise acylation of glycerol-3-phosphate, which is catalyzed by two acyltransferases, GPAT and LPAAT, and leads to the formation of PA. PA is subsequently dephosphorylated to generate DG. The transfer of phosphocholine from cytidine diphosphate-choline to DG is a fi nal step in the synthesis of PC. TG is also synthesized through acylation of DG. These acyltransferases have substrate specifi city because saturated and monounsaturated fatty acids are normally found at sn -1, whereas polyunsaturated fatty acids are enriched at sn -2 of glycerolipids. GPAT1 prefers 16:0-CoA as a substrate ( 40 ), and GPAT3 and GPAT4 recognize a broad range of substrates, from 12:0-CoA to 18:1-or 18:2-CoA ( 41,42 ). LPAAT1 has high activities with 14:0-, 16:0-, and 18:2-CoAs, and intermediate activities with 18:1-and 20:4-CoAs ( 43 ). LPAAT2 prefers 20:4-CoA over 16:0-or 18:0-CoA ( 43,44 ). Two acyltransferases in the Kennedy pathway, LPAAT3 and LPAAT4, prefer 22:6-CoA as a substrate ( 45,46 ). However, it is diffi cult to explain the PGC-1 ␣ -driven alteration of PC (18:0/22:6) and PE (18:0/22:6) content in muscle by these acyltransferases because expression of LPAAT3 and LPAAT4 was not induced by overexpression of PGC-1 ␣ . Furthermore, although PGC-1 ␣ -mediated increase in TG might occur if the Kennedy pathway is involved in the alteration of PC (18:0/22:6) and PE (18:0/22:6) content, we found no evidence from PCA to support the contribution PGC-1 ␣ regulates several exercise-associated aspects of muscle function. To the best of our knowledge, the present study is the fi rst to demonstrate the involvement of PGC-1 ␣ in the exercise-induced rearrangement of phospholipid molecular species in skeletal muscle. Our fi nding suggests that exercise-mediated increase in PGC-1 ␣ expression alters phospholipid profi les, and that these alterations may be one of the adaptive changes of skeletal muscle in response to endurance training. Although further studies are required, the alteration of the phospholipid profi le may contribute to the endurance capacity of skeletal muscle and may partly explain the benefi cial effects of exercise training.