Franc ¸ois Roudier

Lionel Gissot

Fre ´de ´ric Beaudoin

Richard Haslam

Louise Michaelson

Jessica Marion

Diana Molino

Amparo Lima

ˆn Lie

Bach

Halima Morin

Fre ´de ´rique Tellier

Jean-Christophe Palauqui

Yannick Bellec

Charlotte Renne

Martine Miquel

Marco Dacosta

Julien Vignard

Christine Rochat

Jonathan E Markham

Patrick Moreau

Johnathan Napier

Jean-Denis Faure
email: faure@versailles.inra.fr

Very-Long-Chain Fatty Acids Are Involved in Polar Auxin Transport and Developmental Patterning in Arabidopsis W

Keywords:

Very-long-chain fatty acids (VLCFAs) are essential for many aspects of plant development and necessary for the synthesis of seed storage triacylglycerols, epicuticular waxes, and sphingolipids. Identification of the acetyl-CoA carboxylase PASTICCINO3 and the 3-hydroxy acyl-CoA dehydratase PASTICCINO2 revealed that VLCFAs are important for cell proliferation and tissue patterning. Here, we show that the immunophilin PASTICCINO1 (PAS1) is also required for VLCFA synthesis. Impairment of PAS1 function results in reduction of VLCFA levels that particularly affects the composition of sphingolipids, known to be important for cell polarity in animals. Moreover, PAS1 associates with several enzymes of the VLCFA elongase complex in the endoplasmic reticulum. The pas1 mutants are deficient in lateral root formation and are characterized by an abnormal patterning of the embryo apex, which leads to defective cotyledon organogenesis. Our data indicate that in both tissues, defective organogenesis is associated with the mistargeting of the auxin efflux carrier PIN FORMED1 in specific cells, resulting in local alteration of polar auxin distribution. Furthermore, we show that exogenous VLCFAs rescue lateral root organogenesis and polar auxin distribution, indicating their direct involvement in these processes. Based on these data, we propose that PAS1 acts as a molecular scaffold for the fatty acid elongase complex in the endoplasmic reticulum and that the resulting VLCFAs are required for polar auxin transport and tissue patterning during plant development.

INTRODUCTION

Very-long-chain fatty acids (VLCFAs) are defined in plants as fatty acids with an acyl chain of at least 20 carbons in length. VLCFAs are components of seed storage triacylglycerols, cutic-ular and epicuticular lipids, and sphingolipids. VLCFA synthesis requires four endoplasmic reticulum (ER)-bound enzymes constituting the elongase complex that carry out four sequential reactions: first, the condensation of a substrate acyl-CoA with malonyl-CoA catalyzed by ketoacyl-CoA synthase; second, the reduction of 3-ketoacyl-CoA by a ketoacyl-CoA reductase (KCR) followed by dehydration of the resulting 3-hydroxy acyl-CoA by the 3-hydroxy acyl-CoA dehydratase (PASTICCINO2 [PAS2]). The final enzymatic step is the reduction of the enoyl acyl-CoA by an enoyl-CoA reductase (ECR), resulting in an acyl-CoA that is two carbons longer [START_REF] Zheng | Disruptions of the Arabidopsis Enoyl-CoA reductase gene reveal an essential role for very-long-chain fatty acid synthesis in cell expansion during plant morphogenesis[END_REF][START_REF] Joubes | The VLCFA elongase gene family in Arabidopsis thaliana: Phylogenetic analysis, 3D modelling and expression profiling[END_REF]. Malonyl-CoA is synthetized from acetyl-CoA by the cytosolic isoform of the acetyl-CoA carboxylase PAS3/GU ¨RKE [START_REF] Baud | Multifonctional acetyl-CoA carboxylase 1 is essential for very long chain fatty acids elongation and embryo development in Arabidopsis[END_REF][START_REF] Baud | gurke and pasticcino3 mutants affected in embryo development are impaired in acetyl-CoA carboxylase[END_REF]. Loss of function of the 3-hydroxy acyl-CoA dehydratase PASTICCINO2 leads to embryo lethality, demonstrating that VLCFAs are essential for plant development [START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF]. Reduction of VLCFA levels resulting from weak pas3 or pas2 alleles is also associated with abnormal development and ectopic cell proliferation [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF]. Loss of Arabidopsis thaliana ECR CER10 function results in reduced cell expansion and eventually in reduced size of aerial organs [START_REF] Zheng | Disruptions of the Arabidopsis Enoyl-CoA reductase gene reveal an essential role for very-long-chain fatty acid synthesis in cell expansion during plant morphogenesis[END_REF]. Similar observations were recently made upon RNA interference silencing of Arabidopsis KCR [START_REF] Beaudoin | Functional characterization of the Arabidopsis beta-ketoacyl-coenzyme A reductase candidates of the fatty acid elongase[END_REF]. Overexpression of the condensing enzyme FATTY ELONGASE1 (FAE1) was found to alter chloroplast structure and cell shape, whereas an increase of the dehydratase activity leads to cell division and expansion defects [START_REF] Millar | Very-long-chain fatty acid biosynthesis is controlled through the expression and specificity of the condensing enzyme[END_REF][START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF]. Interestingly, ectopic expression of FAE1 with the epidermal-specific promoter of FIDDLEHEAD results in trichome cell death [START_REF] Voisin | Misexpression of FATTY ACID ELONGATION1 in the Arabidopsis epidermis induces cell death and suggests a critical role for phospholipase A2 in this process[END_REF]. Altogether these studies show that VLCFA homeostasis is essential and limiting for plant growth and development.

While the analysis of the role of VLCFAs in cell expansion could be explained by a deficit of essential lipids at the plasma membrane, the involvement of VLCFAs in cell proliferation and differentiation remains unclear. The pas mutants were initially identified as three complementation groups, Pas1, Pas2, and Pas3, showing ectopic cell proliferation in hypocotyls, cotyledons, and leaves. The pas mutations result in profound developmental alterations, which are already detectable during embryogenesis and lead to dwarf seedlings with short and thick hypocotyls and fused and misshapen leaves [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF]. The pas seedlings also show an increased competency for cell division in meristematic and differentiated cells that is enhanced by exogenous cytokinins and leads to callus-like development [START_REF] Haberer | The Arabidopsis gene PEPINO/PASTICCINO2 is required for proliferation control of meristematic and non-meristematic cells and encodes a putative anti-phosphatase[END_REF][START_REF] Harrar | Hormonal control of cell proliferation requires PASTICCINO genes[END_REF]. Moreover, in these mutants, the shoot apical meristem is enlarged and shows an increased expression of meristematic markers [START_REF] Haberer | The Arabidopsis gene PEPINO/PASTICCINO2 is required for proliferation control of meristematic and non-meristematic cells and encodes a putative anti-phosphatase[END_REF][START_REF] Harrar | Hormonal control of cell proliferation requires PASTICCINO genes[END_REF]. This higher proliferative activity is illustrated by the ability of the pas2 mutation to rescue the SHOOT MERISTEMLESS mutant that is devoid of meristematic activity [START_REF] Harrar | Hormonal control of cell proliferation requires PASTICCINO genes[END_REF]. Finally, PAS2 was found to interact in vitro with Cyclin-Dependent Kinase 1, suggesting that VLCFA synthesis might have a direct link with cell cycle progression [START_REF] Da Costa | Arabidopsis PASTICCINO2 is an antiphosphatase involved in regulation of cyclin-dependent kinase A[END_REF]. Cell fate acquisition is also altered in the pas mutants. The absence of cotyledon development in pas3 (gu¨rke) embryos is associated with overlapping cell fate territories within the apex [START_REF] Kajiwara | The GURKE gene encoding an acetyl-CoA carboxylase is required for partitioning the embryo apex into three subregions in Arabidopsis[END_REF].

The identification of two PAS genes as essential enzymes of VLCFA biosynthesis suggest that the pleiotropic pas phenotype is directly associated with lipids that contain such acyl chains. However, the role of VLCFAs in cell fate acquisition and tissue patterning remains unclear, as does the involvement of PAS1 in VLCFA synthesis. PAS1 encodes a large molecular weight member of the immunophilin type FK506 binding protein family that have been reported in other eukaryotes to regulate the activity of many proteins involved in signaling pathways [START_REF] Vittorioso | Mutation in the Arabidopsis PASTICCINO1 gene, which encodes a new FK506-binding protein-like protein, has a dramatic effect on plant development[END_REF][START_REF] Harrar | FKBPs: At the crossroads of folding and transduction[END_REF]. In previous work, we have shown that PAS1 interacts in vitro and in vivo through its C terminus with a member of the NAC family of transcription factors and regulates its subcellular localization [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. Here, we show that the pas1 mutant, similar to the two other pas mutants, has reduced levels of VLCFAs and that PAS1 interacts in the ER with the core elongase components, suggesting a new role of the immunophilin in scaffolding the fatty acid elongase complex. We also show that the impairment of PAS1 function results, like in pas3, in defective cotyledon formation associated with an altered patterning of the embryo apex. We provide evidence that this morphological defect and the absence of lateral root initiation later in development are associated with alterations in auxin response, resulting from a lack of polar targeting of the auxin efflux carrier PIN1. In pas1, defective PIN1 polar targeting is correlated with hypersensitivity to the trafficking inhibitor brefeldin A. Fatty acids are directly involved in polar auxin transport since exogenous application of lipids restore auxin distribution in lateral roots. Based on these data, we propose that PAS1 is involved in fatty acid elongation and that the resulting VLCFAs are essential regulators of cell differentiation during development by regulating polar auxin distribution.

RESULTS

PAS1 Is Required for VLCFA Accumulation in Triacylglycerols, Free Fatty Acids, and Sphingolipids

The fact that the three pas mutants share similar phenotypes and that both PAS2 and PAS3 are essential for VLCFA synthesis [START_REF] Baud | Multifonctional acetyl-CoA carboxylase 1 is essential for very long chain fatty acids elongation and embryo development in Arabidopsis[END_REF][START_REF] Baud | gurke and pasticcino3 mutants affected in embryo development are impaired in acetyl-CoA carboxylase[END_REF][START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF] prompted us to determine whether the pas1-3 mutant was also altered in lipid metabolism. Quantification of VLCFA content in pas1-3 mature seeds revealed that the levels of 20:2 and 22:1 fatty acids were reduced by 50% compared with wild-type levels and that, conversely, the level of short-chain 16:0 and 18:1 fatty acids showed a 40% increase (Figure 1A). Similar changes in the levels of 20:1 and 22:1 but also 16:0 were observed in weak pas2 and pas3 alleles [START_REF] Baud | gurke and pasticcino3 mutants affected in embryo development are impaired in acetyl-CoA carboxylase[END_REF][START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF]. VLCFA levels in pas1 roots were also reduced by 38, 42, and 60% for 22:0, 22:1, and 24:0, respectively (Figure 1B). VLCFAs are also found in sphingolipids [START_REF] Dunn | A post-genomic approach to understanding sphingolipid metabolism in Arabidopsis thaliana[END_REF]. Contrary to seed triacylglycerols, sphingolipids are essential for embryo development [START_REF] Chen | The essential nature of sphingolipids in plants as revealed by the functional identification and characterization of the Arabidopsis LCB1 subunit of serine palmitoyltransferase[END_REF]. Sphingolipids are composed of a long-chain base (amino alcohol) and an amide linked with a fatty acyl chain [START_REF] Dunn | A post-genomic approach to understanding sphingolipid metabolism in Arabidopsis thaliana[END_REF]. These lipids are characterized by the acyl chain length and its degree of unsaturation, as well as the nature of the polar head group, which can be modified by glycosylation and phosphorylation. Like pas2-1, both pas1-3 and pas3-1 mutants were affected in sphingolipid content showing a lower amount of very-long-chain (C 24 and C 26 ) sphingolipids, which are the predominant forms in Arabidopsis (Figure 1C) [START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF].

Specific Sphingolipid Modifications Are Associated with Altered PAS1 Function

Reduction of C 24 and C 26 sphingolipids in pas1-3 and pas3-1 was associated with lower levels of D 8 -unsaturated longchain-base moieties 8-sphingenine (d18:1 8 ) and 4-hydroxy-8sphingenine (t18:1 8 ), likely due to modifications of the different classes of sphingolipids [START_REF] Borner | Analysis of detergent-resistant membranes in Arabidopsis. Evidence for plasma membrane lipid rafts[END_REF][START_REF] Markham | Plant sphingolipids: Separation and identification of major sphingolipid classes from leaves[END_REF][START_REF] Chen | Sphingolipid long-chain base hydroxylation is important for growth and regulation of sphingolipid content and composition in Arabidopsis[END_REF] (see Supplemental Figure 1A online). To assess whether PAS1 targets specific classes of sphingolipids, we performed a complete sphingolipid profiling of both mutants. The levels of simple sphingolipids like Ceramide (Cer) were increased by ;30 and 20% for pas1-3 and pas3-1, respectively, whereas hydroxyceramide levels increased (by ;30%) only in pas1-3 (see Supplemental Figure 1B online). Global increase of Cer in both mutants was mainly caused by the accumulation of long-chain 16:0-Cer (see Supplemental Figure 2 online), whereas the very-long-chain 24:0-and 26:0-Cer levels were substantially reduced in pas3-1 and decreased only modestly in pas1-3 (see Supplemental Figure 2 online). We then examined the complex glycosylated sphingolipids glucosyl inositol phosphoryl ceramide (GIPC) and glucosyl ceramide (GluCer). GIPCs levels were not significantly modified in pas3-1 and were slightly reduced in pas1-3 (;25% decrease). Contrary to GIPCs, GluCer levels were reduced by ;50 and 63%, respectively, for pas1-3 and pas3-1 mutants (see Supplemental Figure 1C online). Very-long (C 20+ ) and long-chain (C 18-20 ) GluCers were specifically reduced (see Supplemental Figure 2 online). Reduction of GluCer levels could result from either a decrease in synthesis or a higher turnover rate (by specific hydrolases). To discriminate between these two possibilities, GluCer synthesis was measured by 14 C-acetate incorporation followed by extraction of this sphingolipid fraction. GluCer labeling in pas1-3 and pas3-1 was strongly reduced compared with the wild type, demonstrating that synthesis of GluCer was impaired in both mutants (see Supplemental Figure 1D online). Mutations affecting biosynthesis of sterols, which are associated with sphingolipids in membrane domains, also lead to strong developmental modifications especially during embryogenesis [START_REF] Schaller | New aspects of sterol biosynthesis in growth and development of higher plants[END_REF]. Contrary to VLCFAs, the levels of sterol glycosides were not modified and the levels of sterol were even increased in both pas1-3 and pas3-1 mutants compared with the wild type (see Supplemental Table 1 online). In conclusion, the pas1-3 mutant is characterized by a specific reduction of VLCFAs like for the other pas mutants, the 3-hydroxyacyl-CoA dehydratase pas2, and the cytosolic acetyl CoA carboxylase pas3.

PAS1 Associates with the VLCFA Elongase Complex in the ER PAS1 is a member of the immunophilin family known in animals and plants to target protein complexes and to regulate their assembly or activities [START_REF] Barik | Immunophilins: for the love of proteins[END_REF][START_REF] Bouchard | Immunophilin-like TWISTED DWARF1 modulates auxin efflux activities of Arabidopsis P-glycoproteins[END_REF][START_REF] Fu | A chloroplast cyclophilin functions in the assembly and maintenance of photosystem II in Arabidopsis thaliana[END_REF]. An attractive hypothesis to explain the decrease in VLCFAs would be that PAS1 interacts directly with the VLCFA elongase complex in the ER. PAS1 was originally described to accumulate in the cytosolic space but was also found in the nucleus upon auxin-induced dedifferentiation [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. Detailed analysis of the distribution of a green fluorescent protein (GFP)-PAS1 fusion showed that it was mainly associated with the ER as confirmed by its colocalization with the CD3-959 marker (Figure 2A). The four enzymes forming the elongase complex are all localized in the ER [START_REF] Zheng | Disruptions of the Arabidopsis Enoyl-CoA reductase gene reveal an essential role for very-long-chain fatty acid synthesis in cell expansion during plant morphogenesis[END_REF][START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF][START_REF] Joubes | The VLCFA elongase gene family in Arabidopsis thaliana: Phylogenetic analysis, 3D modelling and expression profiling[END_REF]. To probe direct interactions between PAS1 and members of the elongase complex, we used bimolecular fluorescence complementation (BiFC). The different genes were cloned in frame with the N-or C-terminal half of yellow fluorescent protein (YFP), and the different combinations were then tested in transient assays [START_REF] Desprez | Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana[END_REF][START_REF] Marion | Systematic analysis of protein subcellular localization and interaction using high-throughput transient transformation of Arabidopsis seedlings[END_REF]. YFP fluorescence was detected for all the combinations involving PAS1 and the elongase complex core members KCR, PAS2, and ECR (Figures 2B and2D). In this assay, the interaction between the elongase complex core members KCR, PAS2, and CER10 was used as a positive control [START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF] (Figure 2B,right). Conversely, no interaction was detected between PAS1, KCR, or CER10 and the ER-localized enzymes DPL1 and SBH2, involved in sphingolipid long chain base metabolism (Figure 2C; see Supplemental Figure 3C online; [START_REF] Tsegaye | Arabidopsis mutants lacking long chain base phosphate lyase are fumonisin-sensitive and accumulate trihydroxy-18:1 long chain base phosphate[END_REF][START_REF] Chen | Sphingolipid long-chain base hydroxylation is important for growth and regulation of sphingolipid content and composition in Arabidopsis[END_REF], demonstrating the specificity of interaction between PAS1 and the elongase core enzymes. PAS1 interaction with the core elongase component PAS2 was confirmed by pull down of either PAS1-GFP expressed in BY2 cell culture using a His 6 -PAS2 affinity column or radiolabeled in vitrotranslated PAS2 using a His 6 -PAS1 affinity column (see Supplemental Figures 3A and 3B online). Interaction of PAS1 with the elongase enzymes in the ER was shown by colocalizing the YFP N -CER10/YFP C -PAS1 fusion protein with the ER marker CD3-959 (Figure 2D). Altogether, these data demonstrate that PAS1 is associated in the ER with the VLCFA elongase complex and thus is directly linked to VLCFA synthesis, as PAS2 and PAS3 are.

PAS1 Is Involved in Cotyledon Development during Embryogenesis

The fact that the three pas mutants are severely altered in their development raises the question of the role of VLCFAs in cell proliferation and differentiation. We first characterized the cellular basis of the pas1 phenotype during early stages of development. Mutations in the PAS1 gene lead to seedlings with short hypocotyls and reduced cotyledon growth. These phenotypic alterations were already observed in 25% of the mature embryos present in siliques from heterozygous pas1/+ plants [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF] and ranged from aborted growth (fox-head shape) to complete loss of cotyledons (almond shape) (Figure 3A; see Supplemental Figure 4 online). The remaining embryos (+/+ and pas1/+) were phenotypically indistinguishable, confirming that pas1-3 was fully recessive [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF] and will be referred to phenotypically as wild type. Morphological analysis of embryogenesis revealed that at the heart stage, mutant embryos are characterized by a flattened apex with larger hypocotyls (Figure 3A) and the accumulation of starch granules (Figure 3A; see Supplemental Figure 4 online). In comparison with the wild type, the protodermis of the mutant apex showed uncoordinated division planes resulting in severe loss of planar polarity. Interestingly, L1 cells displayed heterogeneous shapes, suggesting that epidermal polarity is altered at the single cell level (Figure 3A). We previously showed that PAS1 is expressed mostly in cotyledons of mature embryos [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. In situ hybridizations revealed that PAS1 showed enhanced expression in cotyledons during late developmental stages (Figure 3B). The PAS1 expression pattern is therefore in agreement with the phenotypic defect observed during pas1 embryo development.

PAS1 Is Required for Apex Patterning of the Embryo

To understand better the phenotypic alterations observed in mutant embryos, different cellular markers were introduced into the pas1-3 mutant. Expression analysis of the mitotic marker pCYCB1;1:db-b-glucuronidase (GUS) [START_REF] Colon-Carmona | Technical advance: Spatio-temporal analysis of mitotic activity with a labile cyclin-GUS fusion protein[END_REF] revealed no obvious differences in mitotic patterns between wild-type and pas1-3 embryos (see Supplemental Figure 5 online). However, several cell identity markers representative of specific domains in the embryo apex showed abnormal expression patterns. The expression domain of WUSCHEL (WUS), which characterizes the organizing center of the apical meristem and is normally restricted to two cells in wild-type embryos, was clearly enlarged, encompassing up to 10 cells in the pas1 embryo. This is in agreement with the enhanced expression of the whole KNAT gene family that was observed in the different pas mutants [START_REF] Harrar | Hormonal control of cell proliferation requires PASTICCINO genes[END_REF]. The extension of the WUS expression in pas1 was generated by contiguous cells still defining a stem cell niche (Figure 3C) or by cells scattered into the apical domain of the embryo (see Supplemental Figure 6 online). The WUS domain is normally limited by the expression of CUP SHAPE COTYLEDON2 (CUC2), which defines the boundaries between the cotyledons and the central domain of the meristem (Figure 3C). In pas1-3 embryos, CUC2 expression was not restricted to boundary domains but expanded throughout the entire apical domain. Similarly, AINTEGUMENTA, which marks cotyledon primordia, showed an enlarged expression zone spreading throughout the entire apical part of pas1-3 embryos (Figure 3C). Thus, in the apex of pas1-3 embryos, the different cell territories required for proper organogenesis strongly overlap, leading to the coexistence of otherwise mutually exclusive cell fates. In the most severe cases (5% of homozygous pas1 embryos), this aberrant apex patterning resulted in the total absence of cotyledon development and the lack of bilateral symmetry.

Defective Cotyledon Development in pas1 Is Associated with Modified Polar Auxin Patterning

Cotyledon development is associated with dynamic auxin mobilization with an auxin response restricted to the globular embryo hypophysis at first, followed by a redistribution of the response toward the apical domain and in the cotyledon domains at the triangular and heart stages [START_REF] Aida | Roles of PIN-FORMED1 and MONOPTEROS in pattern formation of the apical region of the Arabidopsis embryo[END_REF][START_REF] Benkova | Local, efflux-dependent auxin gradients as a common module for plant organ formation[END_REF][START_REF] Furutani | PIN-FORMED1 and PINOID regulate boundary formation and cotyledon development in Arabidopsis embryogenesis[END_REF][START_REF] Treml | The gene ENHANCER OF PINOID controls cotyledon development in the Arabidopsis embryo[END_REF]; Izhaki and Bowman, 2007). To monitor auxin response during cotyledon development, we introduced the synthetic auxinresponsive promoter DR5 driving the GFP reporter gene in pas1. As reported previously, wild-type torpedo embryos showed DR5 expression at the root pole, at the tip of the cotyledons, and in the vascular tissues [START_REF] Friml | Efflux-dependent auxin gradients establish the apical-basal axis of Arabidopsis[END_REF]. By contrast, DR5 expression was detected only at the root pole in pas1 (Figure 3D). No DR5 expression could be detected in the apical domain or the vascular tissues of pas1, suggesting that auxin was not properly redistributed. Auxin distribution in the embryo involves specific polar influx and efflux carriers [START_REF] Jenik | Surge and destroy: The role of auxin in plant embryogenesis[END_REF]. In particular, the PIN1 efflux carrier is involved in downward auxin transport through the vascular tissues to the root pole and in upward auxin transport through epidermal cells from the hypophysis to the cotyledon tips. In order to determine whether the altered distribution of auxin observed in pas1-3 could be associated with a defect in polar auxin transport, PIN1 expression and subcellular localization were investigated with a PIN1-GFP fusion driven by the PIN1 promoter. In pas1-3, expression of the pPIN1:PIN1-GFP construct was observed in both epidermal and vascular cells of the embryo, like in the wild type (Figure 3E). However, whereas PIN1-GFP was localized in the lateral membranes of epidermal cells in wild-type embryos (Figure 3F,left), its subcellular distribution was often diffused or aggregated inside pas1 epidermal cells (Figure 3F, right). Immunolocalization analysis confirmed the defective polar distribution of PIN1 in pas1-3 apical cells, where it was found aggregated in the cell and showed no clear polar organization when present at the plasma membrane (Figure 3G, right). For instance, PIN1 was observed in contiguous lateral membranes of two adjacent epidermal cells or in the outward membrane facing the inner seed cavity (Figures 3H and 3I,arrows). The alteration of PIN1 polar localization in pas1-3 could already be observed at the late globular stage before any detectable changes in morphology, suggesting that impaired auxin distribution is likely the cause of the subsequent patterning defects observed in the pas1-3 embryo (see Supplemental Figure 7 online). In agreement with the specific expression of PAS1 at the embryo apex, loss of polar distribution of PIN1 was restricted to the apical domain of pas1-3 embryo but was not modified in vascular cells of the hypocotyl and did not affect DR5 expression at the root pole (Figures 3D and3J). Moreover, similar alterations were observed in embryos of the VLCFA biosynthetic mutant pas3-1 (see Supplemental Figure 8 online), indicating that pas1-3 auxin transport defects were also linked with reduced VLCFA levels. In conclusion, the patterning defects observed in the apical domain of pas1 embryos were associated with abnormal polar distribution of PIN1 and compromised auxin response due to defective VLCFA synthesis.

PAS1 Is Involved in Lateral Root Development

Since PAS1 is also expressed during postembryonic development, we investigated whether other developmental defects observed in pas1 seedlings could also be related to altered polar auxin distribution. The pas1-3 mutant was characterized by the absence of lateral root development (see Supplemental Figure 9A online; [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF][START_REF] Vittorioso | Mutation in the Arabidopsis PASTICCINO1 gene, which encodes a new FK506-binding protein-like protein, has a dramatic effect on plant development[END_REF], the initiation of which is strongly dependent on local auxin accumulation in pericycle, cortex, and epidermal cells [START_REF] Benkova | Local, efflux-dependent auxin gradients as a common module for plant organ formation[END_REF][START_REF] Tanaka | Spatiotemporal asymmetric auxin distribution: A means to coordinate plant development[END_REF]. In primary and lateral roots, PAS1 was found to be specifically expressed in the meristem and the columella cells but not in the elongation or differentiation zones (see Supplemental Figure 9B online). PAS1 was expressed at the initial stage of lateral root formation where periclinal division of pericycle cells is observed, and its expression was maintained mostly at the tip of the root primordia until meristem formation (see Supplemental Figure 9C online). Analysis of DR5:GFP expression indicated that auxin response in pas1-3 primary roots and during the early steps of lateral root initiation was similar to the wild type (Figure 4A; see Supplemental Figure 10A online). However, DR5 expression was rarely detected in the outer cortical cells and never observed in epidermal cells during lateral root formation (Figure 4B). In some instances, pas1-3 mutants produced a few degenerative lateral root primordia with weak or no DR5 expression (Figures 4C and4K) that never developed into lateral roots with proper columella and vascular DR5 expression (see Supplemental Figure 11 online). Thus, as in embryo apex patterning, de novo organogenesis is compromised in pas1 roots.

Defective Lateral Root Development in pas1 Is Associated with Altered Polar Auxin Transport

Auxin distribution during primary and lateral root development is associated with specific expression patterns of PIN1 [START_REF] Benkova | Local, efflux-dependent auxin gradients as a common module for plant organ formation[END_REF][START_REF] Tanaka | Spatiotemporal asymmetric auxin distribution: A means to coordinate plant development[END_REF]. In the primary root, PIN1 distribution was similar between the wild type and the mutant, except that pas1-3 showed a higher occurrence of cytosolic PIN1:GFP aggregates (Figure 5; see Supplemental Figure 10B online). During the initial steps of lateral root initiation, PIN1 expression was clearly visible in the newly divided pericycle cells for both wild-type and pas1-3 roots (Figure 4D, right). At the stage of the first periclinal divisions of the primordia, PIN1:GFP was not sharply localized to the plasma membrane as in the wild type (Figures 4E and4F, right) but showed diffused or aggregated accumulation in the cell (Figure 4G,right). Interestingly, the influx carrier AUX1, which is also expressed during the early steps of lateral root development, was properly targeted to the plasma membrane in the pas1-3 mutant (see Supplemental Figures 10C to 10G online). In conclusion, defective organ formation in pas1 embryos and roots was associated with specific PIN1 mistargeting leading to altered auxin distribution.

BFA Sensitivity of PIN1 Distribution Is Enhanced in pas1

Polar localization of PIN1 results from a continuous turnover of the protein from the plasma membrane to endosomal compartments [START_REF] Geldner | Auxin transport inhibitors block PIN1 cycling and vesicle trafficking[END_REF][START_REF] Geldner | The Arabidopsis GNOM ARF-GEF mediates endosomal recycling, auxin transport, and auxin-dependent plant growth[END_REF][START_REF] Grebe | Arabidopsis sterol endocytosis involves actin-mediated trafficking via ARA6-positive early endosomes[END_REF][START_REF] Muday | Vesicular cycling mechanisms that control auxin transport polarity[END_REF][START_REF] Boutte | The plasma membrane recycling pathway and cell polarity in plants: studies on PIN proteins[END_REF][START_REF] Jaillais | AtSNX1 defines an endosome for auxin-carrier trafficking in Arabidopsis[END_REF]. The lower VLCFA levels in pas1-3 mutant cells could result in an altered trafficking of PIN1 and its subsequent accumulation in the cytosol. In order to test this hypothesis, we used brefeldin A (BFA), which inhibits anterograde (toward the plasma membrane) trafficking of PIN1 by affecting the activity of GNOM, a small G protein member of the GTP/GDP exchange factor family [START_REF] Steinmann | Coordinated polar localization of auxin efflux carrier PIN1 by GNOM ARF GEF[END_REF]. In the wild type, BFA treatment results in the accumulation of PIN1 in large endosomal and Golgi compartments called BFA compartments [START_REF] Boutte | The plasma membrane recycling pathway and cell polarity in plants: studies on PIN proteins[END_REF] and thus provides an interesting tool to monitor PIN1 trafficking. The sensitivity of PIN1 distribution to increasing concentrations of BFA was investigated in wild-type and pas1-3 primary roots. BFA was applied to wild-type and pas1-3 seedlings for 30 min, and the occurrence of PIN1:GFP in large BFA compartments was monitored in stele cells (Figure 5A). BFA (100 mM) led to the formation of BFA compartments of PIN1:GFP in most cells in wild-type seedlings, whereas 25 mM affected ;50% of the cells (Figure 5B). On the contrary, pas1-3 mutant cells showed saturating BFA responses already at 25 mM, and >80% cells were BFA responsive at 17.8 mM (compared with 30% in the wild type), indicating that PIN1 localization was more sensitive to BFA in the pas1-3 background. This increased sensitivity of PIN1 localization to BFA suggests that PAS1, probably via the synthesis of VLCFAs, modifies protein recycling.

VLCFAs Are Required for Auxin Distribution and for Lateral Root Development

Finally, to demonstrate the role of VLCFAs in both auxin distribution and lateral root development, pas1 mutants expressing pDR5:GFP or pPIN1:PIN1-GFP were directly grown in the presence or absence of a mixture of VLCFAs (18:0, 20:0, 22:0, and 24:0). After 4 d of treatment, no significant differences could be observed for wild-type roots (data not shown). However, VLCFA treatment enhanced lateral root development in pas1-3 compared with the wild type (Figures 4I and4J). While only 2% of pas1 seedlings (n = 74) showed lateral root initiations in the control condition, the application of VLCFAs led to a fourfold increase in lateral root outgrowth (9% of pas1 seedlings, n = 88). This result was also confirmed by the direct application of acyl-CoA to pas1 seedlings. Whereas 16:0-CoA had no effect on the number of seedlings showing lateral root development, and 18:0-CoA had only a slight effect (13 and 20%, respectively, of seedlings compared with 10% for the control, n = 30), a mixture of 20:0-, 22:0-, and 24:0-CoA increased the number of seedlings initiating lateral root development (42% of seedlings, n = 30). The effect of VLCFA treatment on pDR5:GFP and pPIN1:PIN1-GFP distribution was also analyzed in wild-type and pas1 roots. As for the root growth assay, the different treatment did not affect DR5 expression and PIN1-GFP localization in the wild type (data not shown). However in fatty acid-treated pas1, DR5 was expressed at the tips of lateral roots (Figures 4K and4L) and PIN1-GFP was targeted to the plasma membrane (Figures 4M and4N), demonstrating that VLCFAs are directly involved in the distribution of auxin during lateral root development.

DISCUSSION

PAS1 Is Involved in VLCFA Synthesis

Several lines of evidence suggest that altered lipid metabolism is directly involved in the pas1 phenotype. First, pas1, like the pas2 and pas3 mutants, is characterized by lower levels of VLCFAs. Second, most of the developmental defects observed in pas1 were also found in the pas3 mutant, which is deficient in the synthesis of malonyl-CoA, a precursor of VLCFA synthesis [START_REF] Baud | gurke and pasticcino3 mutants affected in embryo development are impaired in acetyl-CoA carboxylase[END_REF]. Third, we showed that PAS1 is localized in the ER, at the site of VLCFA synthesis and directly interacts with at least three members of the VLCFA elongase complex, KCR, ECR (CER10), and PAS2. Finally, application of exogenous VLCFAs was sufficient to rescue lateral root development in pas1.

VLCFAs are specific components of seed storage triacylglycerols, cuticular waxes, suberins, some phospholipids, and sphingolipids. Lack of aliphatic suberin in the root led to retarded root growth but did not change the root architecture and, in particular, lateral root development [START_REF] Franke | The DAISY gene from Arabidopsis encodes a fatty acid elongase condensing enzyme involved in the biosynthesis of aliphatic suberin in roots and the chalaza-micropyle region of seeds[END_REF]. Defective accumulation of VLCFAs in seed triacylglycerols or cuticle did not result in embryonic and postembryonic phenotypes equivalent to those observed in pas1 or pas3 mutants. However, sphingolipids are known to regulate cell polarity and differentiation [START_REF] Hoekstra | Membrane dynamics and cell polarity: The role of sphingolipids[END_REF]. Polarized epithelial cells showed specific sphingolipid enrichments in apical and basolateral membranes, moreover perturbing the association of polar protein with sphingolipid domains and preventing their sorting to the plasma membrane [START_REF] Nyasae | Transcytotic efflux from early endosomes is dependent on cholesterol and glycosphingolipids in polarized hepatic cells[END_REF]. Interestingly, sterols and sphingolipids are both involved in detergent-resistant membrane structure, and it is noteworthy that the observed GluCer reduction in pas1 is likely associated with a compensatory increase in sterol and sterol glycoside levels [START_REF] Mongrand | Lipid rafts in higher plant cells: Purification and characterization of Triton X-100-insoluble microdomains from tobacco plasma membrane[END_REF][START_REF] Borner | Analysis of detergent-resistant membranes in Arabidopsis. Evidence for plasma membrane lipid rafts[END_REF][START_REF] Laloi | Insights into the role of specific lipids in the formation and delivery of lipid microdomains to the plasma membrane of plant cells[END_REF]. A similar observation was described in yeast where synthetic lethality resulting from multiple mutations in the sterol biosynthetic pathway can be suppressed by mutations in the sphingolipid biosynthesis pathway [START_REF] Valachovic | Cumulative mutations affecting sterol biosynthesis in the yeast Saccharomyces cerevisiae result in synthetic lethality that is suppressed by alterations in sphingolipid profiles[END_REF]. Interestingly, yeast does not seem to have a PAS1 ortholog, suggesting that the plant elongase complex, which harbors a large variety of components, such as the FAE1-like ketoacyl-CoA synthase family and ELO homologs, requires an extra enzymatic chaperone to carry out its activity.

PAS1-Dependent Synthesis of VLCFAs Is Required for Polar Auxin Distribution

In animals, part of protein sorting to apical and basolateral membranes relies on protein partitioning into sphingolipidenriched microdomains in the trans-Golgi network and in endosomes [START_REF] Hoekstra | Membrane dynamics and cell polarity: The role of sphingolipids[END_REF]. Our work shows that a defect in VLCFA synthesis results in altered subcellular distribution of the polar auxin influx carrier PIN1. PIN1 polarity can be modified by several factors, including the Ser/Thr protein kinase PINOID and the lipid content of membranes. Indeed, a mutation in the STEROL METHYL TRANSFERASE1 gene, which caused a marked decrease in sitosterol and an increase in cholesterol, leads to the redistribution of PIN1 to lateral instead of basal membranes of elongated vascular cells [START_REF] Willemsen | Cell polarity and PIN protein positioning in Arabidopsis require STEROL METHYLTRANSFERASE1 function[END_REF].

Involvement of sterols in the polar distribution of PIN proteins is supported by the fact that the PIN2 efflux carrier colocalizes with the sterol marker filipin [START_REF] Grebe | Arabidopsis sterol endocytosis involves actin-mediated trafficking via ARA6-positive early endosomes[END_REF]. PAS1 is thus not a general regulator of polar trafficking but instead displays some specificity toward protein cargo in the anterograde pathway, which is most likely not limited to PIN1.

Defective Auxin Distribution in pas1 Is Associated with Abnormal Patterning during Organogenesis Defective patterning in the apical domain of pas1 embryos was associated with round shaped cells and random division planes, confirming that cell dedifferentiation and proliferation were associated with modified cell polarity. Alteration of cell polarity is often observed in highly proliferative cells like in mammalian cancers [START_REF] Bissell | Putting tumours in context[END_REF]. In Drosophila melanogaster, loss of cell polarity is associated with Ras-induced tumor progression and invasiveness [START_REF] Igaki | Loss of cell polarity drives tumor growth and invasion through JNK activation in Drosophila[END_REF]. In plants, the involvement of cell polarity in cell differentiation was clearly demonstrated for polar auxin transport during organogenesis. The sequential reorientation of epidermal PIN proteins in the late globular embryo was found to be involved in the generation of new auxin maxima at the site of cotyledon initiation [START_REF] Benkova | Local, efflux-dependent auxin gradients as a common module for plant organ formation[END_REF][START_REF] Friml | Efflux-dependent auxin gradients establish the apical-basal axis of Arabidopsis[END_REF]. Thus, the defective polarity of PIN1 in the apical region of pas1-3 embryos is involved in the alteration of polar auxin transport and the lack of local redistribution of auxin accumulation. The altered auxin gradient at the embryo apex in pas1-3 would then prevent the establishment of bilateral symmetry during embryo development as well as the formation of cotyledons. The absence of cotyledons in pas1 would thus result from the absence of proper determination of the different apical territories that are characterized by mutually exclusive cell fates [START_REF] Kajiwara | The GURKE gene encoding an acetyl-CoA carboxylase is required for partitioning the embryo apex into three subregions in Arabidopsis[END_REF]. Consistent with this, PAS1 expression is enhanced in the cotyledons of late developing embryos. This result confirms the comparative transcriptome analysis of the different domains of the embryo, in which PAS1 appears to be one of the most differentially expressed gene in the apical domains of the embryos at globular and heart stages [START_REF] Casson | Laser capture microdissection for the analysis of gene expression during embryogenesis of Arabidopsis[END_REF]. A similar association between embryo patterning and PIN1 polarity was observed for the KANADI gene family during the establishment of embryo bilateral symmetry [START_REF] Izhaki | KANADI and class III HD-Zip gene families regulate embryo patterning and modulate auxin flow during embryogenesis in Arabidopsis[END_REF]. Indeed, embryos of kan1 kan2 kan4 triple mutants showed ectopic outgrowths on the hypocotyl that were caused by the reversal of PIN polarity and localized accumulation of auxin. Similarly, changes of PIN1 polarity are also known to be essential to relocate auxin into lateral root primordia. As expected, we found that the absence of lateral root development in pas1-3 was also associated with the abnormal distribution of PIN1 but not of AUX1, demonstrating that PAS1 regulates PIN1 polarity in specific cell types throughout plant development, from early embryonic stages to postgerminative de novo organogenesis. Interestingly, the related immunophilin TWISTED DWARF has been characterized as a regulator of auxin carrier activities of the PGP1/PGP19 ABCB transporters [START_REF] Bouchard | Immunophilin-like TWISTED DWARF1 modulates auxin efflux activities of Arabidopsis P-glycoproteins[END_REF]. Moreover, PIN1 stability at the plasma membrane is stabilized by the presence of PGP1/PGP19 [START_REF] Titapiwatanakun | ABCB19/PGP19 stabilises PIN1 in membrane microdomains in Arabidopsis[END_REF]. Altogether, these results support the fact that membranelocalized immunophilins are involved in the regulation of polar auxin transport.

Like many immunophilins, PAS1 likely interacts with different substrates [START_REF] Barik | Immunophilins: for the love of proteins[END_REF] and most probably acts as chaperone or scaffold factor for several ER proteins. We previously showed that PAS1 interacts with a transcription factor of the NAC family and regulates its subcellular localization within the nucleus [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. Overexpression of this transcription factor in pas1 led to a partial rescue of the shoot apical growth [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF], suggesting that PAS1 may interact with other protein substrates. Here, we showed that PAS1 interacts with several enzymes of the elongase complex in the ER and is required for VLCFA elongation (Figure 6). Proteins of the elongase complex likely represent major PAS1 substrates since mutations in PAS2 and PAS3 genes, which are directly involved in VLCFA synthesis, are epistatic to pas1 and lead to similar developmental phenotypes [START_REF] Faure | The PASTICCINO genes of Arabidopsis thaliana are involved in the control of cell division and differentiation[END_REF]. We showed that VLCFAs could rescue lateral root formation in pas1, demonstrating that PAS1-dependent fatty acid elongation is necessary for lateral root development. Identification of the VLCFA-derived molecules and their precise cellular functions will be essential to understand the role of lipids in the polar targeting of proteins.

METHODS

Plant Material and Growth Conditions

The pas1-3 and pas3-1 mutants are ethyl methanesulfonate alleles in the Columbia-0 (Col-0) background that were maintained as heterozygous stocks [START_REF] Baud | Multifonctional acetyl-CoA carboxylase 1 is essential for very long chain fatty acids elongation and embryo development in Arabidopsis[END_REF][START_REF] Baud | gurke and pasticcino3 mutants affected in embryo development are impaired in acetyl-CoA carboxylase[END_REF][START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. Lines expressing pDR5:GFP [START_REF] Friml | Efflux-dependent auxin gradients establish the apical-basal axis of Arabidopsis[END_REF], pPIN1:PIN1-GFP [START_REF] Benkova | Local, efflux-dependent auxin gradients as a common module for plant organ formation[END_REF], pCycB1;1:db-GUS, pWUS:GFP, and pAUX1:AUX1-YFP [START_REF] Swarup | Localization of the auxin permease AUX1 suggests two functionally distinct hormone transport pathways operate in the Arabidopsis root apex[END_REF] were crossed with the pas1-3 or pas3-1 mutants (female donor). PAS1 promoter sequence corresponding to the 2407-bp sequence upstream of the ATG was amplified (59-AAAAAAGCAGGCTCCCACAT-CAGGTGGATATAT-39 and 59-AAGAAAGCTGGGTCCGCCGATCAAAT-CCAGAGT-39), cloned into pDONOR207 and recombined into pMDC162 to generate the pPAS1:GUS construct [START_REF] Curtis | A gateway cloning vector set for high-throughput functional analysis of genes in planta[END_REF]. Several independent Col-0 lines expressing pPAS1:GUS were selected, and a representative line was selected for further studies. All the embryo analyses were performed on siliques from heterozygous pas1-3/+ or pas3-1/+ lines homozygous for the GFP/GUS marker construct to ensure similar developmental stages and maternal environment for wild-type and mutant embryos. Plants were grown in the greenhouse in soil (Tref Substrates) and watered with Plant-Prod nutritive solution (Fertil). Seedlings were germinated on Arabidopsis thaliana agar medium [START_REF] Estelle | Auxin-resistant mutants of Arabidopsis thaliana with an altered morphology[END_REF] and grown for 10 d under constant temperature 188C, 16-h/8-h light/dark cycle, and with 60% humidity. VLCFA complementation was performed by germinating seedlings on horizontal Arabidopsis medium plates then transferring them to a vertical position for 10 d. Seedlings were then transferred to horizontal plates supplemented with 200 mM free fatty acids for 4 to 6 d. Free fatty acids and acyl-CoA (18:0, 20:0, 22:0, and 24:0) were prepared as 20 mM stock solutions in toluene and DMSO, respectively. Control experiments were performed with toluene or DMSO supplemented plates.

Cytological Analyses

Embryos were fixed, stained with propidium iodide, and finally cleared with chloral hydrate [START_REF] Truernit | High-resolution wholemount imaging of three-dimensional tissue organization and gene expression enables the study of phloem development and structure in Arabidopsis[END_REF]. In situ hybridizations were performed on serial sections of siliques at different stages as described previously [START_REF] Vernoux | PIN-FORMED 1 regulates cell fate at the periphery of the shoot apical meristem[END_REF][START_REF] Nikovics | The balance between the MIR164A and CUC2 genes controls leaf margin serration in Arabidopsis[END_REF]. For immunolocalization, samples were fixed for 1 h in 4% (w/v) paraformaldehyde, embedded in paraplast, and sectioned as previously described [START_REF] Smyczynski | The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor[END_REF]. Epitope demasking was performed by incubating the slides in 10 mM citrate buffer, pH 6, in a microwave oven until boiling and then washing in PBS. Samples were then incubated for 2 h at 378C in primary antibody diluted 1/250 in PBS containing 2% (w/v) BSA, washed with PBS (three changes, 3 min per wash), and then incubated for 2 h at 378C in secondary antibody, goat anti-rabbit-IgG antibody conjugated to Alexa Fluor 488 (Molecular Probes, Invitrogen), diluted 1:500 in PBS. After washing in PBS (three changes, 3 min per wash), the samples were mounted in antifading agent (Citifluor; Oxford Instruments). PIN1 was probed with an anti-PIN1 polyclonal serum described previously [START_REF] Boutte | The plasma membrane recycling pathway and cell polarity in plants: studies on PIN proteins[END_REF].

Observations were performed using an inverted Leica TCS SP2-AOBS spectral confocal laser microscope (Leica Microsystems) using either a PL APO 20X0.70 NA or 63X1.20 NA water immersion objective. GFP and Alexa 488 fluorescence were recorded after an excitation at 488 nm (argon laser) and a selective emission band of 505 to 525 nm. Propidium iodide fluorescence was recorded with an excitation at 488 nm and emission band of 600 to 700 nm. GUS staining was performed as described previously [START_REF] Harrar | Hormonal control of cell proliferation requires PASTICCINO genes[END_REF]. BiFC experiments were performed as described earlier [START_REF] Marion | Systematic analysis of protein subcellular localization and interaction using high-throughput transient transformation of Arabidopsis seedlings[END_REF]. PAS2, KCR, CER10, and DPL1 BiFC constructs were described by [START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF], and SBH2 cloning was reported previously [START_REF] Chen | Sphingolipid long-chain base hydroxylation is important for growth and regulation of sphingolipid content and composition in Arabidopsis[END_REF][START_REF] Marion | Systematic analysis of protein subcellular localization and interaction using high-throughput transient transformation of Arabidopsis seedlings[END_REF]. PAS1

Figure 1 .

 1 Figure 1. Lipids with Very Long Acyl Chains Are Altered in pas1-3 and pas3-1 Mutants. (A) Triacylglycerol composition of wild-type (Col-0; black bars) and pas1-3 seeds (white bars). (B) Total fatty acid composition of wild-type (Col-0; black bars) and pas1-3 mutant 10-d-old roots (white bars). (C) Sphingolipid composition of wild-type (Col-0; black bars), pas1-3 (white bars), and pas3-1 (gray bars) seedlings. Data represents the mean of three independent analyses. Error bar indicates SE.

Figure 2 .

 2 Figure 2. PAS1 Interacts with Elongase Enzymes in the ER of Arabidopsis Epidermal Cells. (A) Coexpression of 35S:GFP-PAS1 (left) and ER marker CD3-959: mCherry (middle) showed colocalization (merged, right). (B) Coexpression of the split YFP pairs YFP N -PAS1/KCR-YFP C , YFP N -PAS2/ YFP C -PAS1, and YFP N -CER10/KCR-YFP C led to BiFC of YFP (green). Chloroplast autofluorescence is red. (C) Coexpression of split YFP pairsYFP N -PAS1/ YFP C -DPL1, YFP N -DPL1/ YFP C -PAS1, and YFP N -CER10/YFP C -DPL1 did not produce any YFP fluorescence. (D) Coexpression of YFP N -CER10/YFP C -PAS1 GFP (left) and ER marker CD3-959:mCherry (middle) showed BiFC colocalization in the ER (merged, right). Bars = 10 mm in (A) to (C) and 40 mm in (D).

Figure 3 .

 3 Figure 3. PAS1 Is Required for Cell Patterning and Polarity in the Embryo Apex.(A) Development of wild-type (top row) embryo at the dermatogen, globular, heart, torpedo, and late torpedo stage, respectively (from left to right). The first phenotypic alteration in pas1-3 embryos is visible at the heart stage with the absence of cotyledon formation (bottom row). Mutant embryos (bottom) were taken at the same respective stage (i.e., from the same silique) as the wild type (top row). Apical cells of pas1 embryos have lost their polar growth (inset, bottom) compared with the wild type (inset, top). Embryos were fixed and stained with propidium iodide. (B) In situ hybridization of PAS1 mRNA during embryo development in the wild type. Embryos were taken at globular, young heart, late heart, young torpedo, and late torpedo stages (left to right). (C) In situ hybridization of WUS, CUC2, and ANTEGUMENTA (ANT) mRNA in wild-type (left panel of each pair) and pas1-3 (right panel of each pair) embryos. Bars = 40 mm except for the inset in (A), which is10 mm. (D) pDR5-GFP distribution in the wild type (left) and pas1-3 mutant (right). (E) pPIN1:PIN1-GFP distribution in the wild type (left) and pas1-3 mutant (right). (F) Detail of pPIN1:PIN1-GFP distribution in the tip of a wild-type cotyledon (left) and the apex of pas1-3 embryo (right) at heart stage. (G) to (J) Immunolocalization of PIN1 in wild-type and pas1-3 embryos. (G) Detail of PIN1 distribution in the tip of a wild-type cotyledon (left) and in the apex of the pas1-3 mutant where aggregates are visible (right). (H) and (I) Altered polar distribution of PIN1 in the apex of the pas1-3 mutant with PIN1 localizations facing each other in adjacent cells ([H], arrow) or facing outward ([I], arrow). (J) PIN1 polarity is normal in pas1-3 provascular and root pole cells. Bars = 40 mm in (D) and (E), 20 mm in (F), 10 mm in (G) to (I), and 30 mm in (J).

Figure 4 .

 4 Figure 4. VLCFAs Are Involved in Lateral Root Development and Auxin Polar Distribution.

  (A) to (G) pDR5:GFP ([A] to [C]) and pPIN1:PIN1-GFP ([D] to [G]) expression during sequential steps of lateral root development in the wild type (left) and the pas1-3 mutant (right). In the pas1 mutant, PIN1-GFP was found accumulated inside primordia cells ([E] to [G], right) often in aggregates ([G], arrows). (H) to (J) Exogenous application of VLCFAs restored lateral root development in pas1-3 mutants. Seedlings were grown in presence ([H], right) or absence ([H], left) of 200 mM fatty acids (18:0, 20:0, 22:0, and 24:0). Details of control ([I], green bracket in [H]) or treated roots ([J], red bracket in [H]) are shown. Arrows point to lateral root outgrowth in the pas1-3 mutant ([H] and [J]). (K) to (N) VLCFA application restores polar auxin transport in pas1-3 lateral roots. Normal pDR5:GFP ([K] and [L]) and pPIN1:PIN1-GFP ([M] and [N]) expression patterns were observed in treated pas1-3 lateral root tips ([L] and [N]) but not in untreated mutant roots ([K] and [M]). Bars = 45 mm in (A) to (F), 30 mm in (G) (left), 20 mm in (G) (right), 1 mm in (H), 300 mm in (I) and (J), and 20 mm in (K) to (N).

Figure 5 .

 5 Figure 5. BFA-Dependent PIN1-GFP Aggregation Is Enhanced in pas1 Mutant Cells. (A) BFA induction of PIN1-GFP aggregation in pas1. In contrast with the wild type (left), PIN1-GFP aggregation can be observed in the pas1-3 mutant (right) even in the absence of BFA (0 mM) and is strongly enhanced with 25 mM BFA treatment. (B) The pas1-3 mutation enhances sensitivity to BFA. The relative number of vascular cells showing BFA compartments as illustrated in (A) at 100 mM was monitored in the wild type and pas1 according to the BFA concentration. Data are the mean of three replicates of 25 roots 6 SE.

Figure 6 .

 6 Figure 6. A Model for the Role of PAS Proteins and VLCFAs in Plant Cell Differentiation. Fatty acid elongation requires long-chain fatty acyl-CoA (C n LCFA-CoA) of n carbons and malonyl-CoA produced by the acetyl-CoA carboxylase PAS3. Elongation occurs in the ER membrane with four sequential reactions (boxed enzymes) to eventually produce very-long-chain fatty acyl-CoA of n+2 carbons (C n+2 VLCFA-CoA). PAS1, by its association with elongase enzymes (dashed arrows), is required for fatty acid elongation. The role of VLCFA-CoA on plant development is most probably associated with the synthesis of membrane sphingolipids with long-chain bases (LCBs). Sphingolipids have been described to be involved in membrane trafficking and cell polarity, which are key determinants of cell differentiation.

Very-Long-Chain Fatty Acids and Cell Patterning

ACKNOWLEDGMENTS F.R. was funded by a "Haignere ´" Institut National de la Recherche Agronomique postdoctoral fellowship. A.L. was funded by the European Union Marie Curie short-stay training program ( Versailles-Evry Research Training program). J.M. and C.R. were funded by the 6th European Integrated Project AGRON-OMICS (Grant LSHG-CT-2006-037704). L.B. was funded by the Bourses Cance ´ropo ˆle Ile-de-France. D.M. was funded by the European Union Versailles-Evry Research Training program and Agence Nationale pour la Recherche (07-BLAN-202). We thank Jan Traas (Ecole Normale Supe ´rieure, Lyon), Malcolm Bennett (University of Nottingham), Patrick Laufs (Institut National de la Recherche Agronomique, Versailles) for the gifts of pPIN1-PIN1:GFP, pAUX1-AUX1:YFP, and pWUS:GFP constructs, respectively. We thank Olivier Grandjean and the Plateforme Cytologie et d'Imagerie du Ve ´ge ´tal de Versailles for their great help in confocal microscopy as well as Bruno Letarnec for taking care of the plants. The work was partially funded by the ANR blanc SphingopolaR (ANR-07-BLAN-0202). Rothamsted Research receives grant-aided support from the Biotechnology and Biological Sciences Research Council. Received September 8, 2009; revised January 21, 2010; accepted January 27, 2010; published February 9, 2010.

The Plant Cell coding sequence was amplified (59-AAAAAAGCAGGCTTCATGGCGG-TAGGCGATCAGACG-39 and 59-CAAGAAAGCTGGGTCTGTAAATTTG-GCGCTCACAAA-39) and cloned in pDONOR207. Briefly, all the different open reading frames were cloned from pDONOR207 by GATEWAY recombination (Invitrogen) in pBIFC vectors in both N and C configuration with either N-ter YFP or C-ter YFP and under the control of the 35S promoter [START_REF] Desprez | Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana[END_REF][START_REF] Marion | Systematic analysis of protein subcellular localization and interaction using high-throughput transient transformation of Arabidopsis seedlings[END_REF]. The four different constructs expressed under the control of the 35S promoter were systematically tested leading to eight interaction assays per pair of proteins. The ER marker CD3-959:mCherry was used as previously described [START_REF] Marion | Systematic analysis of protein subcellular localization and interaction using high-throughput transient transformation of Arabidopsis seedlings[END_REF].

BFA sensitivity was measured on 7-d-old seedlings (pPIN1:PIN1-GFP in Col-0 or in pas1-3 background) grown on Arabidopsis agar medium and treated with 10, 17.8, 25, 50, or 100 mM BFA for 30 min. Seedlings were mounted in the BFA solution, and GFP fluorescence was recorded for 30 min. Percentages of cells carrying BFA compartments were estimated in an average of 25 roots (700 cells, three replicates) per treatment.

Lipid Analyses

Mass Quantification of Lipid Species

Triacylglycerol analysis was performed as previously described [START_REF] Bach | The plant very long chain hydroxy fatty Acyl-CoA dehydratase PASTICCINO2 is essential and limiting for plant development[END_REF]. GIPC, GluCer, and Cer from Arabidopsis wild-type and mutant seedlings were extracted, isolated, and quantified as detailed by [START_REF] Markham | Rapid measurement of sphingolipids from Arabidopsis thaliana by reversed-phase highperformance liquid chromatography coupled to electrospray ionization tandem mass spectrometry[END_REF]. Long-chain bases of sphingolipids (LCB) were determined as previously described by [START_REF] Borner | Analysis of detergent-resistant membranes in Arabidopsis. Evidence for plasma membrane lipid rafts[END_REF]. To isolate neutral lipids, total lipids were loaded onto HPTLC plates developed in hexane/ethylether/acetic acid (90:15:2, v/v) and separated into diacylglycerols (R F 0.08), sterols (R F 0.17), fatty alcohols (R F 0.22), and free fatty acids (R F 0.29). Lipids were identified by comigration with known standards and quantified by densitometry analysis [START_REF] Macala | Analysis of brain lipids by high performance thin-layer chromatography and densitometry[END_REF]) using a TLC scanner 3 (CAMAG).

De Novo Synthesis of Ceramides and Glucosylceramides

Arabidopsis wild-type and mutant seedlings were placed in small tubes and incubated with 10 mCi of [1-14 C]acetate (54 Ci/mol) in 1 mL of distilled water at 248C for 0 to 240 min. After incubations, seedlings were washed with distilled water three times to recover the remaining acetate. Lipids were extracted by chloroform:methanol (2:1, v/v) for 60 min at 60 to 708C and then washed three times with 9% NaCl. The solvent was evaporated and lipids were dissolved in an appropriate volume of chloroform/ methanol (1:1, v/v). Polar lipids were analyzed by loading total lipids onto HPTLC plates (60F254; Merck), which were developed in methyl acetate/n-propanol/chloroform/methanol/0.25% aqueous KCl (25:25: 25:10:9, v/v) [START_REF] Heape | Improved one-dimensional thin-layer chromatographic technique for polar lipids[END_REF]. Under such conditions, glucosylceramides and sterol-glycosides were not separated. To allow their separation, the corresponding spot was further applied onto an HPTLC plate developed with chloroform/methanol (85:15, v/v) [START_REF] Hillig | Formation of glucosylceramide and sterol glucoside by a UDP-glucose-dependent glucosylceramide synthase from cotton expressed in Pichia pastoris[END_REF]. Ceramides were isolated on HPTLC plates as described previously [START_REF] Heape | Sphingolipid metabolic disorders in Trembler mouse peripheral nerves in vivo result from an abnormal substrate supply[END_REF]. Labeling of the different lipids was determined and quantified with a phosphor imager (Molecular Dynamics).

Accession Numbers

Sequence data from this article can be found in the Arabidopsis Genome Initiative or GenBank/EMBL databases under the following accession numbers: PAS1, At3g54010; PAS2, At5g10480; PAS3, At1g36160; DPL1, At1g27980; KCR, At1g67730; CER10, At3g55360; SBH2, At1g14290; PIN1, At1g73590; AUX1, At2g38120; WUS, At2g17950; and ANT, At4g37750.

Supplemental Data

The following materials are available in the online version of this article. Sterols and Sterol-glycosides are expressed as % of total lipid mass. Each value is the mean of several experiments (n) ±SEM.


A Rosendo

T Druet

† C Péry

J P Bidanel

Correlative responses for carcass and meat quality traits to selection for ovulation rate or prenatal survival in French Large White pigs 1

Keywords: carcass composition, meat quality, ovulation rate, pig, prenatal survival, selection experiment

Correlated effects of selection for components of litter size on carcass and meat quality traits were estimated using data from 3 lines of pigs derived from the same Large White base population. Two lines were selected for 6 generations on high ovulation rate at puberty (OR) or high prenatal survival corrected for ovulation rate in the first 2 parities (PS). The third line was an unselected control (CON). The 3 lines were kept for a 7th generation, but without any selection. Carcass and meat quality traits were recorded on the 5th to 7th generation of the experiment. Carcass traits included dressing percentage, carcass length (LGTH), average backfat thickness (ABT), estimated lean meat content, and 8 carcass joint weight traits. Meat quality traits included pH recorded 24 h after slaughter (pH24) of LM, gluteus superficialis (GS), biceps femoris (BF), and adductor femoris (AD) muscles, as well as reflectance and water-holding capacity (WHC) of GS and BF muscles. Heritabilities of carcass and meat quality traits and their genetic correlations with OR and PS were estimated using REML methodology applied to a multiple trait animal model. Correlated responses to selection were then estimated by computing differences between OR or PS and CON lines at generations 5 to 7 using least squares and mixed model methodology. Heritability (h 2 ) estimates were 0.08 ± 0.04, 0.58 ± 0.10, 0.70 ± 0.10, and 0.74 ± 0.10 for dressing percentage, LGTH, ABT, and lean meat content, respectively, ranged from 0.28 to 0.72 for carcass joint traits, from 0.28 to 0.45 for pH24 and reflectance measurements, and from 0.03 to 0.11 for WHC measurements. Both OR and PS had weak genetic correlations with carcass (r G = -0.09 to 0.17) and most meat quality traits. Selection for OR did not affect any carcass composition or meat quality trait. Correlated responses to selection for PS were also limited, with the exception of a decrease in pH24 of GS and BF muscles (-0.12 to -0.14 after 6 generations; P < 0.05), in WHC of GS muscle (-18.9 s after 6 generations; P < 0.05) and a tendency toward an increase in loin weight (0.44 kg after 6 generations; P < 0.10) .

INTRODUCTION

Litter size is a major component of the breeding goal in pig dam lines [START_REF] Estany | Correlated response to selection for litter size in pigs: II. Carcass, meat, and fat quality traits[END_REF][START_REF] Holl | Result from nine generations of selection for increased litter size in swine[END_REF][START_REF] Hamann | Estimation of genetic parameters for litter size as a sow and boar trait in German herdbook Landrace and Pietrain swine[END_REF]. There is experimental evidence that prolificacy has been substantially improved by applying high selection intensities in pig maternal lines (Bidanel et al., 1994a;[START_REF] Tribout | Estimation, par utilisation de semence congelée, du progrès génétique réalisé en France entre 1977 et 1998 dans la race porcine Large White: Résultats pour quelques caractères de reproduction femelle[END_REF]. However, direct selection for litter size remains difficult in closed maternal lines where extreme selection intensities cannot be applied. Selecting on ovulation rate and prenatal survival has been proposed as an alternative to direct selection to improve litter size at birth [START_REF] Johnson | Selection for components of reproduction in swine[END_REF]. A selection experiment on ovulation rate or prenatal survival was set up at the Institut National de la Recherche Agronomique (INRA) to estimate their genetic parameters and check their potential as selection criteria for litter size [START_REF] Blasco | A Bayesian analysis of genetic parameters and selection response for litter size components in pigs[END_REF]Rosendo et al., 2007b). Before commercial application, direct and correlated responses to selection for the main traits of interest in pig production should be estimated. Correlated responses for reproductive and growth traits were presented in Rosendo et al. (2007a,c). The objective of the present study is to analyze correlated responses to selection for ovulation rate or prenatal survival on various carcass and meat quality traits.

MATERIALS AND METHODS

Animal care followed the general guidelines outlined in the European welfare regulation (directive 91/630/ EC).

Animals and Experimental Design

The experiment was carried out over 7 generations (6 generations of selection and a seventh generation of random mating) at the INRA experimental herd of Galle (Avord, France).

Two lines of pigs were selected for high ovulation rate at puberty (OR) or high prenatal survival (PS) over the first 2 parities corrected for ovulation rate at fertilization (ORF). Prenatal survival was computed as (total number born/ORF) + (0.018 × ORF) (Rosendo et al., 2007b). The correction term was introduced to avoid trends in PS associated with variation in OR. The term 0.018 represented an average literature value for the phenotypic regression coefficient of PS on ORF. A third line was kept as an unselected control (CON line). At each generation, approximately 50 gilts and 6 to 8 boars from first litters were kept for breeding. Boars were chosen on a within-sire family basis in the 3 lines, whereas sows were selected on a population basis in the 2 selected lines and within-dam families in the CON line (Rosendo et al., 2007b). A mating plan was established that minimized inbreeding at each generation. Additional details on the experimental design are given by Rosendo et al. (2007b). Average reproductive performance of the 3 lines on the 5th, 6th, and 7th generations are given in Table 1.

The sow herd was managed under a batch farrowing system. Females were distributed into 7 farrowing batches, which then became postweaning and performance test batches of their progeny. Seven gilts from each line were introduced in each farrowing batch. Females produced 2 experimental litters. Replacement animals were chosen in first parity litters, whereas pigs from second parity litters were all slaughtered. All piglets were weaned at 4 wk of age and moved to a postweaning unit until 10 wk of age. They were then allotted to a performance test building in which they were housed in pens of 10 to 12 animals of the same line, where they stayed until the end of the test period when the average BW within a pen reached 90 kg. Animals were fed ad libitum with a single pelleted diet based on cereals and soybean meal containing 3,100 kcal of DE/kg and 170 g of CP/kg during the whole test period and until the day before slaughter. Additional details on the performance test can be found in Rosendo et al. (2007a).

Carcass and meat quality traits were measured only on second litters from the 5th, 6th, and 7th generation of the experiment. On average, 2 pigs per litter (i.e., 1 gilt and 1 castrate) were chosen at random within each litter and were slaughtered in a commercial abattoir at an average BW of 100.6 ± 2.6 kg. All animals were fasted for 16 h before they were transported for 2 h to the abattoir where they were allowed to rest for an additional 18 h before they were killed by electrical stunning and immediate exsanguination.

Measurements

Carcass measurements were recorded on the day after slaughter. Dressing percentage (DP) was calculated as the ratio of cold carcass (with head and feet) to unfasted BW. Carcass length (LGTH), measured from the cervical vertebra to the anterior edge of the pubic symphysis, and average backfat thickness (ABT), mea-Table 1. Average reproductive performance of the 3 lines 1 (and their SE) on the 5th, 6th, and 7th generations (from Rosendo et al., 2007a,b) sured at the levels of first rib, last rib, and last lumbar vertebrae, were recorded on the right half-carcass. The kidney and leaf fat were removed, and the right halfcarcass was divided into 7 cuts. The front and back feet were separated from the limbs at the levels of the carpal bones and tibio-tarsal joint, respectively. The ham was isolated along a first line parallel to the general direction of the sacrum and a second line perpendicular to the long axis of the carcass between the last lumbar and the first sacral vertebrae. The loin was separated from the belly and the shoulder with a cut starting under the psoas muscle at the level of the last lumbar vertebra and ending under the blade bone. The belly and the shoulder were separated with a cut perpendicular to the long axis of the carcass between the 5th and the 6th rib. The backfat was dissected from the loin, but the other cuts remained untrimmed. The weight of 8 cuts [i.e., backfat (BFWT), belly (BEWT), feet (FWT), ham (HAWT), head (HEWT), leaf fat (LFWT), loin (LOWT), and shoulder (SWT)] were recorded, and carcass lean meat content (LMC) was estimated from the 3 joint weights using the following equation (Larzul et al., 1999b); LMC = 16.56 + (71.6 HAWT + 83.0 LOWT -76.2 BFWT)/half-carcass weight.

Meat quality traits were also measured 24 h postmortem. Ultimate pH (pH24) was taken directly on different types of muscles [i.e., on LM, gluteus superficialis (GS), biceps femoris (BF), and adductor femoris (AD)] muscles using a combined glass electrode (Ingold, Metter Toledo, Switzerland) and a portable pH meter (CG818, Schott Geräte, Mainz, Germany). A Minolta Chromameter CR-300 (Minolta Camera Co. Ltd., Osaka, Japan) with a 8-mm aperture, a pulsed xenon arc lamp, CIE illuminant, and calibration to the white calibration plate, was used to measure lightness (L*), redness/greenness (a*), and yellowness/blueness (b*) on BF and GS muscles. Water-holding capacity was measured on GS and BF muscles by the filter paper imbibition time method [START_REF] Charpentier | Correlations between carcass characteristics and meat quality in Large White pigs[END_REF]. This method consists of measuring the time required for the complete wetting of a piece of filter paper (around 1 cm 2 ) put on the freshly cut surface of the muscle. The time of the observation is limited to 3 min. Measurements on ham muscles (i.e., GS, BF, and AD muscles) were performed on the cut surface of the ham at the intersection of the ham and loin joints. Measurements on LM muscle were carried out at the level of the last rib. The number of records for each line × generation subclass ranged from 40 to 59.

Statistical Analyses

A total of 12 carcass traits were analyzed (i.e., DP, LGTH, ABT, LMC, and the BFWT, HAWT, LOWT, SWT, FWT, LFWT, BEWT, and HEWT). Twelve meat quality traits were considered (i.e., the 4 pH24 measurements, denoted by pH24_LM, pH24_GS, pH24_BF, and pH24_AD), 6 color measurements [i.e., a*GS and a*BF (redness/greenness), b*GS and b*BF (yellowness/blueness), and L*GS and L*BF (lightness)], and 2 water-holding capacity measures denoted by WHC_GS and WHC_BF. Elementary statistics for the traits studied are given in Table 2.

Least Squares Analyses of Line Differences. The data were first analyzed using least squares methodology with the GLM procedure (SAS Inst. Inc., Cary, NC). Least squares means for each line-generation subclass were computed using a linear model. Fixed effects for all traits included line (OR, PS, or CON), generation number (5, 6, or 7), their interaction, sex (female or castrate), and contemporary group. Contemporary group was defined by the fattening batch for carcass traits and slaughter date for meat quality traits. Parity of dam was included as a fixed effect in the model for carcass traits, and BW at slaughter and animal inbreeding coefficient were included in all models as linear covariates.

Mixed Model Analyses. Variance components were first estimated using REML methodology [START_REF] Patterson | Recovery of inter-block information when block sizes are unequal[END_REF] applied to 3 trait mixed linear animal models including the selection criterion (OR and PS) and an additional trait. The models used depended on the trait (Table 3), but all derived from the following base model:

y = Xβ + Z a a + W c c + e,
where y represents the vector of observations; X, Z a , and W c are known incidence matrices relating observations to fixed and random effects; β = vector of fixed effects (i.e., generation, line, sex, parity number of the dam, and contemporary group); a is the vector of direct genetics effects of the pig; c is a vector of random common litter effects of the animals; and e is the vector of random residual effects. The models used to estimate the variance components for OR and SP are detailed in Rosendo et al. (2007b). The following means and (co) variance structures were assumed across random effects in the model: where σ a 2 is the additive genetic variance; A is the additive relationship matrix; σ c 2 and σ e 2 are the random common litter and residual variances, respectively; and I are the identity matrices of appropriate dimension. The analyses were performed using VCE [START_REF] Neumaier | Restricted maximum likelihood estimation of covariances in sparse linear models[END_REF] and ASREML [START_REF] Gilmour | ASReml User Guide Release 1.0[END_REF] computer packages. Additive genetic effects were then estimated as back-solutions from REML analyses at convergence and used to compute average values for each line × generation combination.

E(y) = Xβ

RESULTS

Significance of fixed effects, covariates, random effects, and total variation explained by fixed effects (R 2 ) are given in Table 3. A moderate proportion (24 to 29%) of the total variation was explained by the fixed effect model for LGTH, ABT, and LMC. In contrast, the fixed effect model explained 30 to 50% of the total variation for DP and the weight of joints traits, except for LFWT. Inbreeding significantly decreased performance level for BFWT, BEWT, FWT, and LOWT. Body weight affected (P < 0.05) all traits except DP and LMC. Date of slaughter was the main fixed effect influencing meat quality characteristics: it explained from 36 to 50% of the total variation of ultimate pH and color measurements. Its effect was of less importance for muscle lightness and water-holding capacity, explaining 13 to 21% of the phenotypic variation for these traits.

Genetic parameter estimates for OR and PS were remarkably stable across the different 3-trait analyses and were the same as those reported by Rosendo et al. (2007b). The REML estimates of heritability and common litter effects for carcass and meat quality traits are given in Table 4. Estimates had rather large SE due to the rather limited amount of data. With the exception of the small value obtained for DP (0.08) and the moderate estimates for SWT and HEWT (0.28 and 0.38, respectively), carcass traits had high heritability values (from 0.48 to 0.72). Common litter effects were close to zero for some traits, but were far from negligible for other traits (up to 0.09). Estimates of heritability were moderate for pH measurements, for a* and b* color parameters and the lightness of GS traits (from 0.26 to 0.45). Conversely, very low h 2 values were obtained for L*BF and the 2 water-holding capacity measurements. The magnitude of common litter effects was in the same range as that obtained for carcass traits (0.00 to 0.09).

Estimates of phenotypic and genetic correlations between the selection criteria and carcass traits and meat quality traits are shown in Table 5. Phenotypic correlations between both OR or SP and carcass composition traits were all very weak (-0.10 to 0.06). Genetic correlations between OR and carcass traits were also very low; only one value had an absolute value above 0.10 (-0.17 ± 0.14 with LGTH). The situation was not much different for PS; only 4 estimates were outside the interval -0.10 to 0.10. Though not significant, these values (i.e., -0.19 with LMC; -0.13 with LOWT; 0.16 with ABT) might indicate a small antagonism between PS and carcass leanness.

Similarly, phenotypic correlations with meat quality traits were all small (less than 0.10 in absolute values). With a single exception, genetic correlations of OR and PS with pH measurements were nonsignificant (P > 0.05), but consistently positive (0.25 to 0.37). Consistently positive genetic correlations were also obtained between PS and a* and b* measurements, whereas positive and negative values were obtained with OR. Conversely, genetic correlations of muscle lightness and water-holding capacity with OR and PS were all negative.

First analyses of line differences on generations 5 to 7 did not show any significant (P > 0.10) line × generation interaction or generation effect. Average line effects and differences over the 3 generations were hence computed. Estimates for LC and mixed model of contrasts between these overall line effects are shown in Table 6. They were consistent in sign and magnitude. Line differences were not significant (P > 0.10) for the vast majority of the traits analyzed. The only exceptions for carcass traits were loin weight, which tended to increase in the PS line (P < 0.10), and SWT, which tended to increase in the OR line (P < 0.10). No difference (P > 0.10) was observed for meat quality traits in the OR line. Conversely, pH values tended to increase in the PS, with significant differences (P < 0.05) for 2 of the 4 pH measurements. There were also tendencies toward a greater b* value and a reduced water-holding capacity on GS muscle.

DISCUSSION

For most composition carcass traits, the heritability values found are in agreement with the average literature values reported by [START_REF] Sellier | Genetics of meat and carcass traits[END_REF] and with the estimates obtained by Larzul et al. (1999b) using data that were partly common (CON line) to the current data. Heritability of DP (0.08 ± 0.04) is close to previous estimates (0.11 ± 0.04) of Larzul et al. (1999b), but is noticeably less than the estimates reported by [START_REF] Ducos | Multivariate restricted maximum likelihood estimation of genetic parameters for growth, carcass and meat quality traits in French Large White pigs[END_REF] and [START_REF] Tribout | Paramètres génétiques de quelques caractères de qualité de la viande dans les races porcines Large White et Landrace Français[END_REF] and the average literature value of 0.30 reported by [START_REF] Ducos | Paramètres génétiques des caractères de production chez le porc. Mise au point bibliographique[END_REF]. As suggested by Larzul et al. (1999b), the low heritability estimate might be due to the fact that slaughter BW was recorded before feed withdrawal in the present study, whereas it was recorded on fasted animals in most other studies. Conversely, the heritability value obtained for ABT is greater than most estimates reported in the literature [START_REF] Lo | Genetic analyses of growth, real-time ultrasound, carcass, and pork quality traits in Duroc and Lan-drace Pigs: II. Heritabilities and correlations[END_REF]Bidanel et al., 1994b;[START_REF] Serenius | Genetic associations of prolificacy with performance, carcass, meat quality, and leg conformation traits in the Finnish Landrace and Large White pig populations[END_REF] and close to the estimates reported by [START_REF] Larzul | Phenotypic and genetic parameters for longissimus muscle fiber characteristics in relation to growth, carcass, and meat quality traits in Large White pigs[END_REF].

In this study, the ultimate pH traits had high heritability values. Estimates were 0.35 on average, superior to the average literature value of 0.21 reported in [START_REF] Sellier | Genetics of meat and carcass traits[END_REF]. They were also greater than the estimates reported by Larzul et al. (1999a) and van Wijk et al. ( 2005), but were close to the value reported by Nguyen GS, BF, and AF = gluteus superficialis, biceps femoris, and adductor femoris muscles, respectively. pH24 = pH recorded 24 h after slaughter. a*, b*, and L* = reflectance measurements: redness/greenness, yellowness/ blueness, and lightness, respectively. WHC = water-holding capacity.

2

The dash (-) indicates that the common litter effect was nonsignificant and was removed from final analyses. et al. (2006). The heritability estimates for reflectance are within the range found in the literature [START_REF] Lo | Genetic analyses of growth, real-time ultrasound, carcass, and pork quality traits in Duroc and Lan-drace Pigs: II. Heritabilities and correlations[END_REF][START_REF] De Vries | Genetic parameters of pork quality and production traits in Yorkshire populations[END_REF][START_REF] Knapp | Genetic parameters for lean meat content and meat quality traits in different pigs breeds[END_REF][START_REF] Sellier | Genetics of meat and carcass traits[END_REF]. Among meat quality traits, heritability values for WHC_GS and WHC_BF were in agreement with those reported by [START_REF] Tribout | Paramètres génétiques de quelques caractères de qualité de la viande dans les races porcines Large White et Landrace Français[END_REF] and [START_REF] Sonesson | Genetic parameters and trends of meat quality, carcass composition and performance traits in two selected lines of large white pigs[END_REF], but less than the estimates presented by the most recent literature reviews (Larzul et al., 1999a;Hermesch et al., 2000a). With estimates ranging from 2 to 9% of the phenotypic variance, common litter effects fell in the range of values reported in the literature.

Overall, genetic correlations between OR or SP and carcass composition traits were not significantly different from zero. Yet, correlation estimates with fatness/ leanness traits tended to show a small antagonism between carcass leanness and OR or SP. Literature estimates of the genetic correlation between ABT and OR are close to zero [START_REF] Rothschild | Biology and Genetics of Reproduction[END_REF][START_REF] Ruiz-Flores | Direct and correlated responses to two-stage selection for ovulation rate and number of fully formed pigs at birth in swine[END_REF]). An absence of genetic association between production and reproduction traits has been reported in numerous studies (e.g., Noguera et al., 2002;[START_REF] Serenius | Genetic associations of prolificacy with performance, carcass, meat quality, and leg conformation traits in the Finnish Landrace and Large White pig populations[END_REF][START_REF] Arango | Threshold-linear estimation of genetic parameters for farrowing mortality, litter size, and test performance of Large White sows[END_REF]. Yet other recent studies have found a small unfavorable genetic correlation between carcass composition and litter size (e.g., Hermesch et al., 2000b;[START_REF] Zhang | Genetic parameters and genetic trends in the Chinese x European Tiameslan composite line. I. Genetic parameters[END_REF][START_REF] Chen | Genetic correlations between lean growth and litter traits in U.S. Yorkshire, Duroc, Hampshire, and Landrace pigs[END_REF][START_REF] Holm | Genetic correlations between reproduction and production traits in swine[END_REF][START_REF] Tribout | Paramètres génétiques du nombre de porcelets nés vivants par portée et relations avec les caractères du contrôle en ferme dans les populations porcines Large White type femelle et Landrace Français[END_REF]. This unfavorable relationship may result from the detrimental effects on placental development of increased uterine crowding during early gestation, which is a critical period for the muscle development of embryos [START_REF] Foxcroft | The biological basis for prenatal programming of postnatal performance in pigs[END_REF]. Embryos are then more likely to suffer from intrauterine growth retardation and have reduced numbers of muscle fibers and result in pigs with reduced growth performance and lean content at slaughter [START_REF] Foxcroft | The biological basis for prenatal programming of postnatal performance in pigs[END_REF][START_REF] Rehfeldt | Consequences of birth weight for postnatal growth performance and carcass quality in pigs as related to myogenesis[END_REF]. It may also be because leaner sows would have more difficulties mobilizing their body reserves to support the needs of embryos/piglets during gestation and lactation, as well as because of an unfavorable effect of a nutritional unbalance in late lactation on the subsequent ovulation or implantation rates.

Estimates of genetic correlations between OR and the different ultimate pH traits varied from 0.11 to 0.25 and from -0.05 to 0.37 between SP and ultimate pH traits, but also had very large SE (from 0.17 to 0.21 and 0.23 to 0.25, respectively), making it difficult to draw a conclusion regarding the genetic relationships between OR and meat quality. No other estimate of genetic correlations between OR or PS and meat quality traits is, to our knowledge, available in the literature to compare with results reported herein. More generally, few estimates of genetic relationships between reproduction and meat quality traits are available in the literature. Hermesch et al. (2000b) reported a negative genetic Table 6. Least squares (LS) and mixed model (MM) estimates (±SE) of line differences for carcass and meat quality traits after 6 generations of selection correlation (-0.26) between ultimate pH and number of piglets born alive, indicating that an increased litter size is associated with a smaller ultimate pH. Genetic correlations between reflectance, lightness, water-holding capacity, and OR or PS traits were inconsistent and mostly of small magnitude, indicating no clear genetic relationships between these traits. However, genetic correlation between OR or PS and water-holding capacity traits were weak and negative. In contrast, Hermesch et al. (2000b) presented slightly positive genetic correlations between water-holding capacity and number of piglets born alive in first and second parity, but this correlation was not confirmed for the number of piglets born alive in third parity litters.

Line differences in the last generations of the experiment indicate that correlated responses to selection on carcass traits are limited. It has to be noted that the tendency toward an increase in loin weight in the PS line is not consistent with the weakly unfavorable genetic correlations between PS and carcass leanness traits. No clear explanation for this discrepancy was found. Conversely, the tendency toward an increase in ultimate pH is consistent with genetic parameter estimates. Larzul et al. (1999b) estimated correlated responses to selection for a reduced muscle glycolytic potential, which is strongly correlated with ultimate pH, in a Large White population issued from the same base population and obtained little correlated response for litter traits.

In conclusion, this study provides first estimates of genetic correlations of components of litter size (i.e., OR and PS) with carcass and meat quality traits. They globally confirm that both groups of traits are weakly correlated, although there might be a small antagonism between OR or PS and carcass leanness and a weak positive (i.e., favorable) relationship with meat ultimate pH.

Table 2 .

 2 Descriptive statistics for carcass and meat quality traits

	Carcass trait	Abbreviation	n 1	Mean	SD
	Dressing percentage, %	DP	584	76.6	2.5
	Carcass length, mm	LGTH	593	956	28
	Average backfat thickness, mm	ABT	593	28.1	4.6
	Weight of joint, kg
	Backfat	BFWT	593	4.02	0.75
	Ham	HAWT	593	8.83	0.57
	Loin	LOWT	593	10.99	0.85
	Shoulder	SWT	593	4.99	0.37
	Feet	FWT	593	1.08	0.1
	Leaf fat	LFWT	471	0.54	0.17
	Belly	BEWT	593	4.64	0.53
	Head	HEWT	588	6.76	0.58
	Estimated lean meat content, %	LMC	588	51.6	3.3
	Meat quality trait 2
	pH24 of LM	pH24_LM	480	5.58	0.23
	pH24 of gluteus superficialis muscle	pH24_GS	480	5.90	0.29
	pH24 of biceps femoris muscle	pH24_BF	480	5.59	0.23
	pH24 of adductor femoris muscle	pH24_AF	480	5.69	0.26
	Reflectance of:
	Gluteus superficialis, a* value	a*GS	484	8.20	2.44
	Biceps femoris, a* value	a*BF	484	7.63	2.37
	Gluteus superficialis, b* value	b*GS	484	7.20	2.42
	Biceps femoris, b* value	b*BF	484	6.39	2.30
	Gluteus superficialis, L* value	L*GS	484	46.1	4.2
	Biceps femoris, L* value	L*BF	484	47.0	4.0
	Water-holding capacity, s
	Gluteus superficialis	WHC_GS	484	170.8	26.6
	Biceps femoris	WHC_BF	484	112.2	63.0
	1 n = number of records.

2 pH24 = pH recorded 24 h after slaughter.

Table 3 .

 3 Models of analyses, significance, and total variation explained by fixed-effect terms (R 2 )

					Fixed effect		Covariate	Random effect 2
					Parity	Contemporary		BW at
	Trait 1	R 2	Sex	Line	of dam	group 3	Inbreeding	slaughter	Animal	Litter
	Dressing percentage	0.33	***	ns	ns	***	ns	ns	X	-
	Carcass length	0.24	***	ns	*	**	ns	***	X	X
	Average backfat thickness	0.27	***	ns	***	***	ns	***	X	X
	Backfat weight	0.34	***	ns	***	**	***	***	X	X
	Ham weight	0.50	***	ns	**	***	ns	***	X	X
	Loin weight	0.35	***	+	*	***	*	***	X	X
	Shoulder weight	0.30	*	+	*	***	ns	***	X	X
	Feet weight	0.41	**	ns	ns	***	**	***	X	-
	Leaf fat weight	0.19	***	ns	ns	**	**	***	X	-
	Belly weight	0.33	ns	ns	ns	***	**	***	X	-
	Head weight	0.37	ns	ns	ns	***	ns	***	X	-
	Estimated lean meat content	0.29	***	ns	***	**	ns	ns	X	X
	pH24 of LM muscle	0.38	**	+	ns	***	ns	ns	X	X
	pH24 of GS muscle	0.35	**	*	ns	***	ns	*	X	X
	pH24 of BF muscle	0.33	**	*	ns	***	ns	*	X	X
	pH24 of AF muscle	0.38	***	+	ns	***	ns	ns	X	X
	a* value, GS muscle	0.48	ns	ns	ns	***	ns	*	X	-
	a* value, BF muscle	0.42	ns	ns	ns	***	ns	*	X	X
	b* value, GS muscle	0.57	ns	+	ns	***	ns	*	X	X
	b* value, BF muscle	0.55	ns	ns	ns	***	ns	ns	X	X
	L* value, GS muscle	0.20	ns	+	ns	**	ns	ns	X	X
	L* value, BF muscle	0.13	ns	ns	ns	**	ns	ns	X	X
	WHC, GS muscle	0.18	ns	*	ns	**	ns	*	X	-
	WHC, BF muscle	0.21	ns	ns	ns	**	ns	ns	X	-

Table 4 .

 4 Estimates of heritability (±SE) and common litter effect for carcass and meat quality traits

			Common	Phenotypic
	Trait 1	Heritability ± SE	litter effect 2	variance
	Dressing percentage, %	0.08 ± 0.04	-	6.8
	Carcass length, mm	0.58 ± 0.10	0.09	747
	Average backfat thickness, mm	0.70 ± 0.10	0.04	19.9
	Backfat weight, kg	0.68 ± 0.08	0.04	0.45
	Ham weight, kg	0.48 ± 0.11	0.08	0.24
	Loin weight, kg	0.55 ± 0.09	0.07	0.56
	Shoulder weight, kg	0.28 ± 0.06	0.03	0.12
	Feet weight, kg	0.66 ± 0.09	-	0.01
	Leaf fat weight, kg	0.72 ± 0.08	-	0.03
	Belly weight, kg	0.69 ± 0.10	-	0.23
	Head weight, kg	0.38 ± 0.08	-	0.74
	Estimated lean meat content, %	0.74 ± 0.10	0.04	8.9
	pH24 of LM muscle	0.35 ± 0.12	0.08	0.05
	pH24 of GS muscle	0.34 ± 0.13	0.07	0.05
	pH24 of BF muscle	0.29 ± 0.11	0.07	0.06
	pH24 of AF muscle	0.45 ± 0.12	0.08	0.11
	a* value, GS muscle	0.34 ± 0.12	-	3.78
	a* value, BF muscle	0.45 ± 0.13	0.09	4.26
	b* value, GS muscle	0.28 ± 0.11	0.02	3.05
	b* value, BF muscle	0.26 ± 0.11	0.02	3.28
	L* value, GS muscle	0.30 ± 0.11	0.04	16.4
	L* value, BF muscle	0.13 ± 0.11	0.09	15.3
	WHC, GS muscle, s	0.11 ± 0.06	-	645
	WHC, BF muscle, s	0.03 ± 0.06	-	651
	1

The thesis work of A. Rosendo is funded by the Consejo Nacional de Ciencia y Tecnologia, the Secretaria de Educacion Publica, Programa de Becas Complemantarias, México, and Colegio de Postgraduados,

Mexico. 2 Present address: Colegio de Postgraduados, Campus Veracruz km 26.5 carr. Fed. Veracruz-Xalapa, Predio Tepetates, Mpio. De Manlio Fabio Altamirano CP 91700, Apdo. Postal 421, Veracruz,

Mexico.3 Corresponding author: jean-pierre.bidanel@jouy.inra.fr Received July 21, 2009. Accepted October 26, 2009.

at INRA Institut National de la Recherche Agronomique on August 21, 2013 www.journalofanimalscience.org Downloaded from

Selection for litter size components in pigs

ne = not estimated because of convergence problems.

GS, BF, and AF = gluteus superficialis, biceps femoris, and adductor femoris muscles, respectively. pH24 = pH recorded 24 h after slaughter. a*, b*, and L* = reflectance measurements: redness/greenness, yellowness/blueness, and lightness, respectively. WHC = water-holding capacity. 2 Difference between selected and control lines; †P < 0.10; *P < 0.05. Selection for litter size components in pigs