Identification and modelling of human breast cancer subtypes

Résumé Le cancer du sein est le cancer le plus commun chez les femmes et est responsable de presque 30% de tous les nouveaux cas de cancer en Europe. On estime le nombre de décès liés au cancer du sein en Europe est à plus de 130.000 par an. Ces chiffres expliquent l'impact social considérable de cette maladie. Les objectifs de cette thèse étaient: (1) d'identifier les prédispositions et les mécanismes biologiques responsables de l'établissement des sous-types spécifiques de cancer du sein; (2) les valider dans un modèle ín vivo "humain-dans-souris"; et (3) de développer des traitements spécifiques à chaque sous-type de cancer du sein identifiés. Le premier objectif a été atteint par l'intermédiaire de l'analyse des données d'expression de gènes des tumeurs, produite dans notre laboratoire. Les données obtenues par puces à ADN ont été produites à partir de 49 biopsies des tumeurs du sein provenant des patientes participant dans l'essai clinique EORTC 10994/BIG00-01. Les données étaient très riches en information et m'ont permis de valider des données précédentes des autres études d'expression des gènes dans des tumeurs du sein. De plus, cette analyse m'a permis d'identifier un nouveau sous-type biologique de cancer du sein. Dans la première partie de la thèse, je décris I identification des tumeurs apocrines du sein par l'analyse des puces à ADN et les implications potentielles de cette découverte pour les applications cliniques. Le deuxième objectif a été atteint par l'établissement d'un modèle de cancer du sein humain, basé sur des cellules épithéliales mammaires humaines primaires (HMECs) dérivées de réductions mammaires. J'ai choisi d'adapter un système de culture des cellules en suspension basé sur des mammosphères précédemment décrit et pat décidé d'exprimer des gènes en utilisant des lentivirus. Dans la deuxième partie de ma thèse je décris l'établissement d'un système de culture cellulaire qui permet la transformation quantitative des HMECs. Par la suite, j'ai établi un modèle de xénogreffe dans les souris immunodéficientes NOD/SCID, qui permet de modéliser la maladie humaine chez la souris. Dans la troisième partie de ma thèse je décris et je discute les résultats que j'ai obtenus en établissant un modèle estrogène-dépendant de cancer du sein par transformation quantitative des HMECs avec des gènes définis, identifiés par analyse de données d'expression des gènes dans le cancer du sein. Les cellules transformées dans notre modèle étaient estrogène-dépendantes pour la croissance, diploïdes et génétiquement normales même après la culture cellulaire in vitro prolongée. Les cellules formaient des tumeurs dans notre modèle de xénogreffe et constituaient des métastases péritonéales disséminées et du foie. Afin d'atteindre le troisième objectif de ma thèse, j'ai défini et examiné des stratégies de traitement qui permettent réduire les tumeurs et les métastases. J'ai produit un modèle de cancer du sein génétiquement défini et positif pour le récepteur de l'estrogène qui permet de modéliser le cancer du sein estrogène-dépendant humain chez la souris. Ce modèle permet l'étude des mécanismes impliqués dans la formation des tumeurs et des métastases. Abstract Breast cancer is the most common cancer in women and accounts for nearly 30% of all new cancer cases in Europe. The number of deaths from breast cancer in Europe is estimated to be over 130,000 each year, implying the social impact of the disease. The goals of this thesis were first, to identify biological features and mechanisms --responsible for the establishment of specific breast cancer subtypes, second to validate them in a human-in-mouse in vivo model and third to develop specific treatments for identified breast cancer subtypes. The first objective was achieved via the analysis of tumour gene expression data produced in our lab. The microarray data were generated from 49 breast tumour biopsies that were collected from patients enrolled in the clinical trial EORTC 10994/BIG00-01. The data set was very rich in information and allowed me to validate data of previous breast cancer gene expression studies and to identify biological features of a novel breast cancer subtype. In the first part of the thesis I focus on the identification of molecular apacrine breast tumours by microarray analysis and the potential imptìcation of this finding for the clinics. The second objective was attained by the production of a human breast cancer model system based on primary human mammary epithelial cells {HMECs) derived from reduction mammoplasties. I have chosen to adopt a previously described suspension culture system based on mammospheres and expressed selected target genes using lentiviral expression constructs. In the second part of my thesis I mainly focus on the establishment of a cell culture system allowing for quantitative transformation of HMECs. I then established a xenograft model in immunodeficient NOD/SCID mice, allowing to model human disease in a mouse. In the third part of my thesis I describe and discuss the results that I obtained while establishing an oestrogen-dependent model of breast cancer by quantitative transformation of HMECs with defined genes identified after breast cancer gene expression data analysis. The transformed cells in our model are oestrogen-dependent for growth; remain diploid and genetically normal even after prolonged cell culture in vitro. The cells farm tumours and form disseminated peritoneal and liver metastases in our xenograft model. Along the lines of the third objective of my thesis I defined and tested treatment schemes allowing reducing tumours and metastases. I have generated a genetically defined model of oestrogen receptor alpha positive human breast cancer that allows to model human oestrogen-dependent breast cancer in a mouse and enables the study of mechanisms involved in tumorigenesis and metastasis.

SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

GENCODE: producing a reference annotation for ENCODE

Background: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manualannotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results.Results: The GENCODE gene features are divided into eight different categories of which onlythe first two (known and novel coding sequence) are confidently predicted to be protein-codinggenes. 5’ rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentallyverify the initial annotation. Of the 420 coding loci tested, 229 RACE products have beensequenced. They supported 5’ extensions of 30 loci and new splice variants in 50 loci. In addition,46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15putative transcripts. We assessed the comprehensiveness of the GENCODE annotation byattempting to validate all the predicted exon boundaries outside the GENCODE annotation. Outof 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only twoof them in intergenic regions.Conclusions: In total, 487 loci, of which 434 are coding, have been annotated as part of theGENCODE reference set available from the UCSC browser. Comparison of GENCODEannotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained withinthe two sets, which is a reflection of the high number of alternative splice forms with uniqueexons annotated. Over 50% of coding loci have been experimentally verified by 5’ RACE forEGASP and the GENCODE collaboration is continuing to refine its annotation of 1% humangenome with the aid of experimental validation.

Gene finding in the chicken genome

Background: Despite the continuous production of genome sequence for a number of organisms,reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularlytrue for genomes for which there is not a large collection of known gene sequences, such as therecently published chicken genome. We used the chicken sequence to test comparative andhomology-based gene-finding methods followed by experimental validation as an effective genomeannotation method.Results: We performed experimental evaluation by RT-PCR of three different computational genefinders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram wascomputed and each component of it was evaluated. The results showed that de novo comparativemethods can identify up to about 700 chicken genes with no previous evidence of expression, andcan correctly extend about 40% of homology-based predictions at the 5' end.Conclusions: De novo comparative gene prediction followed by experimental verification iseffective at enhancing the annotation of the newly sequenced genomes provided by standardhomology-based methods.

Conserved chromosomal clustering of genes governed by chromatin regulators in Drosophila

BACKGROUND: The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure. RESULTS: We have performed genome-wide expression studies of trx and ash2 mutants in Drosophila melanogaster. Using computational analysis of our microarray data, we have identified 25 clusters of genes potentially regulated by TRX. Most of these clusters consist of genes that encode structural proteins involved in cuticle formation. This organization appears to be a distinctive feature of the regulatory networks of TRX and other chromatin regulators, since we have observed the same arrangement in clusters after experiments performed with ASH2, as well as in experiments performed by others with NURF, dMyc, and ASH1. We have also found many of these clusters to be significantly conserved in D. simulans, D. yakuba, D. pseudoobscura and partially in Anopheles gambiae. CONCLUSION: The analysis of genes governed by chromatin regulators has led to the identification of clusters of functionally related genes conserved in other insect species, suggesting this chromosomal organization is biologically important. Moreover, our results indicate that TRX and other chromatin regulators may act globally on chromatin domains that contain transcriptionally co-regulated genes.

Use of a short hPDE6b promoter for rod gene transfer in a model of severe retinal cystrophy, the Rd10 mouse

Purpose: Gene therapy of severe retinal dystrophies directly affecting photoreceptor is still a challenge in terms of clinical application. One of the main hurdles is to generate high transgene expression specifically in rods or cones. In the present study, we are investigating the possibility to drive hPDE6b expression in the Rd10 mouse retina using a specific sequence of the human PDE6b promoter. Methods: Two 5' flanking fragments of the human PDE6b gene: (-93 to +53 (146 bp) and -297 to +53 (350 bp, see Di Polo and Farber, 1995) were cloned in different plasmids in order to check their expression in vitro and in vivo. These elements drove the activity of either luciferase (pGL3 plasmids) or EGFP (AAV2/8 backbone). Then, an AAV2/8 vector carrying the PDE6b cDNA was tested with subretinal injections at P9 in the Rd10 eyes. Eye fundus, OCT, ERG recordings and histological investigations were performed to assess the efficacy of the gene transfer. Results: The short PDE6b promoter containing 146bp (-93 to +53) showed the highest activity in the Y-79 cells, as described previously (Di Polo and Farber, 1995). Subretinal administrations of AAV2/8-PDE6bpromoter-EGFP allowed a rapid expression specifically in rods and not in cones. The expression is faster than a vector containing the CMV promoter. The AAV2/8-PDE6bpromoter-PDE6b and the control vector were injected at P9 in the Rd10 mouse retina and investigated 5 weeks post-injection. Out of 14 eyes, 6 presented an increased rod sensitivity of about 300 fold, and increased a- and b-wave responses in ERG recordings. Flicker stimulations revealed that cones are also functional. OCT images and histological analyses revealed an increased ONL size in the injected area. The retina treated with the therapeutic vector presented 4-6 rows of photoreceptors with outersegments containing PDE6b. In the control eyes, only 2-4 rows of photoreceptors with almost no OS were observed . Conclusions: The 146 bp promoter sequence (-93 to + 53) is the shortest regulatory element described to date which allows to obtain efficient rod-specific expression in the context of somatic gene transfer. This first result is of great interest for AAV vector design in general allowing more space for the accommodation of transgenes of interest and good expression in rods. Moreover we showed the proof of principle of the efficacy of AAV2/8-PDE6bp-PDE6b vector in the Rd10 mouse model of severe photoreceptor degeneration without using neither AAV mutated capsids, nor self-complementary vectors.

Characterization of human gene expression changes after olive oil ingestion: an exploratory approach

Olive oil consumption is protective against risk factors for cardiovascular and cancer diseases. A nutrigenomic approach was performed to assess whether changes in gene expression could occur in human peripheral blood mononuclear cells after oli ve oil ingestion at postprandial state. Six healthy male volunteers ingested, at fasting state, 50 ml of olive oil. Prior to intervention a 1-week washout period with a controlled diet and sunflower oil as the only source of fat was followed. During the 3 days before and on the intervention day, a very low-phenolic compound diet was followed. At baseline (0 h) and at post-ingestion (6 h), total RNA was isolated and gene expression (29,082 genes) was evaluated by microarray. From microarray data, nutrient-gene interactions were observed in genes related to metabolism, cellular processes, cancer, and atherosclerosis (e.g. USP48 by 2.16; OGT by 1.68-fold change) and associated processes such as inflammation (e.g. AKAP13 by 2.30; IL-10 by 1.66-fold change) and DNA damage (e.g. DCLRE1C by 1.47; POLK by 1.44- fold change). When results obtained by microarray were verified by qRT-PCR in nine genes, full concordance was achieved only in the case of up-regulated genes. Changes were observed at a real-life dose of olive oil, as it is daily consumed in some Mediterranean areas. Our results support the hypothesis that postprandial protective changes related to olive oil consumption could be mediated through gene expression changes.

Changes in microRNA expression contribute to pancreatic β-cell dysfunction in prediabetic NOD mice.

During the initial phases of type 1 diabetes, pancreatic islets are invaded by immune cells, exposing β-cells to proinflammatory cytokines. This unfavorable environment results in gene expression modifications leading to loss of β-cell functions. To study the contribution of microRNAs (miRNAs) in this process, we used microarray analysis to search for changes in miRNA expression in prediabetic NOD mice islets. We found that the levels of miR-29a/b/c increased in islets of NOD mice during the phases preceding diabetes manifestation and in isolated mouse and human islets exposed to proinflammatory cytokines. Overexpression of miR-29a/b/c in MIN6 and dissociated islet cells led to impairment in glucose-induced insulin secretion. Defective insulin release was associated with diminished expression of the transcription factor Onecut2, and a consequent rise of granuphilin, an inhibitor of β-cell exocytosis. Overexpression of miR-29a/b/c also promoted apoptosis by decreasing the level of the antiapoptotic protein Mcl1. Indeed, a decoy molecule selectively masking the miR-29 binding site on Mcl1 mRNA protected insulin-secreting cells from apoptosis triggered by miR-29 or cytokines. Taken together, our findings suggest that changes in the level of miR-29 family members contribute to cytokine-mediated β-cell dysfunction occurring during the initial phases of type 1 diabetes.

SV40-induced expression of calretinin protects mesothelial cells from asbestos cytotoxicity and may be a key factor contributing to mesothelioma pathogenesis.

The calcium-binding protein calretinin has emerged as a useful marker for the identification of mesotheliomas of the epithelioid and mixed types, but its putative role in tumor development has not been addressed previously. Although exposure to asbestos fibers is considered the main cause of mesothelioma, undoubtedly, not all mesothelioma patients have a history of asbestos exposure. The question as to whether the SV40 virus is involved as a possible co-factor is still highly debated. Here we show that increased expression of SV40 early gene products in the mesothelial cell line MeT-5A induces the expression of calretinin and that elevated calretinin levels strongly correlate with increased resistance to asbestos cytotoxicity. Calretinin alone mediates a significant part of this protective effect because cells stably transfected with calretinin cDNA were clearly more resistant to the toxic effects of crocidolite than mock-transfected control cells. Down-regulation of calretinin by antisense methods restored the sensitivity to asbestos toxicity to a large degree. The protective effect observed in clones with higher calretinin expression levels could be eliminated by phosphatidylinositol 3-kinase (PI3K) inhibitors, implying an important role for the PI3K/AKT signaling (survival) pathway in mediating the protective effect. Up-regulation of calretinin, resulting from either asbestos exposure or SV40 oncoproteins, may be a common denominator that leads to increased resistance to asbestos cytotoxicity and thereby contributes to mesothelioma carcinogenesis.

Repeated exposure to Lutzomyia intermedia sand fly saliva induces local expression of interferon-inducible genes both at the site of injection in mice and in human blood.

During a blood meal, Lutzomyia intermedia sand flies transmit Leishmania braziliensis, a parasite causing tegumentary leishmaniasis. In experimental leishmaniasis, pre-exposure to saliva of most blood-feeding sand flies results in parasite establishment in absence of any skin damages in mice challenged with dermotropic Leishmania species together with saliva. In contrast, pre-immunization with Lu. intermedia salivary gland sonicate (SGS) results in enhanced skin inflammatory exacerbation upon co-inoculation of Lu. intermedia SGS and L. braziliensis. These data highlight potential unique features of both L. braziliensis and Lu. intermedia. In this study, we investigated the genes modulated by Lu. intermedia SGS immunization to understand their potential impact on the subsequent cutaneous immune response following inoculation of both SGS and L. braziliensis. The cellular recruitment and global gene expression profile was analyzed in mice repeatedly inoculated or not with Lu. intermedia. Microarray gene analysis revealed the upregulation of a distinct set of IFN-inducible genes, an immune signature not seen to the same extent in control animals. Of note this INF-inducible gene set was not induced in SGS pre-immunized mice subsequently co-inoculated with SGS and L. braziliensis. These data suggest the parasite prevented the upregulation of this Lu. intermedia saliva-related immune signature. The presence of these IFN-inducible genes was further analyzed in peripheral blood mononuclear cells (PBMCs) sampled from uninfected human individuals living in a L. braziliensis-endemic region of Brazil thus regularly exposed to Lu. intermedia bites. PBMCs were cultured in presence or absence of Lu. intermedia SGS. Using qRT-PCR we established that the IFN-inducible genes induced in the skin of SGS pre-immunized mice, were also upregulated by SGS in PBMCs from human individuals regularly exposed to Lu. intermedia bites, but not in PBMCs of control subjects. These data demonstrate that repeated exposure to Lu. intermedia SGS induces the expression of potentially host-protective IFN-inducible genes.

Altered expression of the transcription factor Mef2c during retinal degeneration in Rpe65-/- mice.

Purpose. To investigate the role of the myocyte enhancer factor 2 (Mef2) transcription factor family in retinal diseases, Mef2c expression was assessed during retinal degeneration in the Rpe65(-/-) mouse model of Leber's congenital amaurosis (LCA). Mef2c-dependent expression of photoreceptor-specific genes was further addressed. Methods. Expression of Mef2 members was analyzed by oligonucleotide microarray, quantitative PCR (qPCR) and in situ hybridization. Mef2c-dependent transcriptional activity was assayed by luciferase assay in HEK293T cells. Results. Mef2c was the only Mef2 member markedly downregulated during retinal degeneration in Rpe65(-/-) mice. Mef2c mRNA level was decreased by more than 2 fold at 2 and 4 months and by 3.5 fold at 6 months in retinas of Rpe65(-/-) mice. Downregulation of Mef2c at the protein level was confirmed in Rpe65(-/-) retinas. The decrease in Mef2c mRNA levels in the developing Rpe65(-/-) retinas, from post-natal day (P)13 onward, was concomitant with the decreased expression of the rod-specific transcription factors Nrl and Nr2e3. Nrl was further shown to drive Mef2c transcriptional activity, supporting a physiological role for Mef2c in the retina. In addition, Mef2c appeared to act as a transcriptional repressor of its own expression, as well as those of the retina-specific retinal G-protein coupled receptor (Rgr), rhodopsin and M-opsin genes. Conclusions. These findings highlight the early altered regulation of the rod-specific transcriptional network in Rpe65-related disease. They further indicate that Mef2c may act as a novel transcription factor involved in the development and the maintenance of photoreceptor cells.

Dissection of the complex phenotype in cuticular mutants of Arabidopsis reveals a role of SERRATE as a mediator.

Mutations in LACERATA (LCR), FIDDLEHEAD (FDH), and BODYGUARD (BDG) cause a complex developmental syndrome that is consistent with an important role for these Arabidopsis genes in cuticle biogenesis. The genesis of their pleiotropic phenotypes is, however, poorly understood. We provide evidence that neither distorted depositions of cutin, nor deficiencies in the chemical composition of cuticular lipids, account for these features, instead suggesting that the mutants alleviate the functional disorder of the cuticle by reinforcing their defenses. To better understand how plants adapt to these mutations, we performed a genome-wide gene expression analysis. We found that apparent compensatory transcriptional responses in these mutants involve the induction of wax, cutin, cell wall, and defense genes. To gain greater insight into the mechanism by which cuticular mutations trigger this response in the plants, we performed an overlap meta-analysis, which is termed MASTA (MicroArray overlap Search Tool and Analysis), of differentially expressed genes. This suggested that different cell integrity pathways are recruited in cesA cellulose synthase and cuticular mutants. Using MASTA for an in silico suppressor/enhancer screen, we identified SERRATE (SE), which encodes a protein of RNA-processing multi-protein complexes, as a likely enhancer. In confirmation of this notion, the se lcr and se bdg double mutants eradicate severe leaf deformations as well as the organ fusions that are typical of lcr and bdg and other cuticular mutants. Also, lcr does not confer resistance to Botrytis cinerea in a se mutant background. We propose that there is a role for SERRATE-mediated RNA signaling in the cuticle integrity pathway.

Chicken BAFF--a highly conserved cytokine that mediates B cell survival.

Members of the tumor necrosis factor (TNF) family play key roles in the regulation of inflammation, immune responses and tissue homeostasis. Here we describe the identification of the chicken homologue of mammalian B cell activating factor of the TNF family (BAFF/BLyS). By searching a chicken EST database we identified two overlapping cDNA clones that code for the entire open reading frame of chicken BAFF (chBAFF), which contains a predicted transmembrane domain and a putative furin protease cleavage site like its mammalian counterparts. The amino acid identity between soluble chicken and human BAFF is 76%, considerably higher than for most other known cytokines. The chBAFF gene is most strongly expressed in the bursa of Fabricius. Soluble recombinant chBAFF produced by human 293T cells interacted with the mammalian cell-surface receptors TACI, BCMA and BAFF-R. It bound to chicken B cells, but not to other lymphocytes, and it promoted the survival of splenic chicken B cells in culture. Furthermore, bacterially expressed chBAFF induced the selective expansion of B cells in the spleen and cecal tonsils when administered to young chicks. Our results suggest that like its mammalian counterpart, chBAFF plays an important role in survival and/or proliferation of chicken B cells.

The peroxisome proliferator-activated receptor alpha regulates amino acid metabolism.

The peroxisome proliferator-activated receptor alpha is a ligand-activated transcription factor that plays an important role in the regulation of lipid homeostasis. PPARalpha mediates the effects of fibrates, which are potent hypolipidemic drugs, on gene expression. To better understand the biological effects of fibrates and PPARalpha, we searched for genes regulated by PPARalpha using oligonucleotide microarray and subtractive hybridization. By comparing liver RNA from wild-type and PPARalpha null mice, it was found that PPARalpha decreases the mRNA expression of enzymes involved in the metabolism of amino acids. Further analysis by Northern blot revealed that PPARalpha influences the expression of several genes involved in trans- and deamination of amino acids, and urea synthesis. Direct activation of PPARalpha using the synthetic PPARalpha ligand WY14643 decreased mRNA levels of these genes, suggesting that PPARalpha is directly implicated in the regulation of their expression. Consistent with these data, plasma urea concentrations are modulated by PPARalpha in vivo. It is concluded that in addition to oxidation of fatty acids, PPARalpha also regulates metabolism of amino acids in liver, indicating that PPARalpha is a key controller of intermediary metabolism during fasting.