In Vivo Gene Expression in Granulosa Cells During Pig Terminal Follicular Development. A
Total Page:16
File Type:pdf, Size:1020Kb
In vivo gene expression in granulosa cells during pig terminal follicular development. A. Bonnet, Kim-Anh Lê Cao, M. Sancristobal, F. Benne, C. Robert-Granié, G. Law-So, S. Fabre, Philippe Besse, E. de Billy, Hélène Quesnel, et al. To cite this version: A. Bonnet, Kim-Anh Lê Cao, M. Sancristobal, F. Benne, C. Robert-Granié, et al.. In vivo gene expression in granulosa cells during pig terminal follicular development.. Reproduction, BioScientifica, 2008, 136 (2), pp.211-24. 10.1530/REP-07-0312. hal-00960020 HAL Id: hal-00960020 https://hal.archives-ouvertes.fr/hal-00960020 Submitted on 30 May 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. REPRODUCTIONRESEARCH In vivo gene expression in granulosa cells during pig terminal follicular development A Bonnet, K A Leˆ Cao1,3, M SanCristobal, F Benne, C Robert-Granie´1, G Law-So, S Fabre2, P Besse3, E De Billy, H Quesnel4, F Hatey and G Tosser-Klopp INRA, UMR 444, Ge´ne´tique Cellulaire, F-31326 Castanet-Tolosan Cedex, France, 1INRA, UR, Station d’Ame´lioration Ge´ne´tique des Animaux, F-31326 Castanet-Tolosan Cedex, France, 2INRA, UMR6175, Physiologie de la Reproduction et des Comportements, CNRS, Universite´ de Tours, Haras Nationaux, F-37380 Nouzilly, France, 3Universite´ Paul Sabatier, UMR 5219, Institut de Mathe´matiques, F-31062 Toulouse Cedex 9, France and 4INRA, UMR INRA/AgroCampus Rennes: Syste`mes d’Elevage, Nutrition Animale et Humaine, F-35590 Saint Gilles, France Correspondence should be addressed to A Bonnet; Email: [email protected] Abstract Ovarian antral follicular development is clearly dependent on pituitary gonadotrophins FSH and LH. Although the endocrine mechanism that controls ovarian folliculogenesis leading to ovulation is quite well understood, the detailed mechanisms and molecular determinants in the different follicular compartments remain to be clarified. The aim of this study was to identify the genes differentially expressed in pig granulosa cells along the terminal ovarian follicle growth, to gain a comprehensive view of these molecular mechanisms. First, we developed a specific micro-array using cDNAs from suppression subtractive hybridization libraries (345 contigs) obtained by comparison of three follicle size classes: small, medium and large antral healthy follicles. In a second step, a transcriptomic analysis using cDNA probes from these three follicle classes identified 79 differentially expressed transcripts along the terminal follicular growth and 26 predictive genes of size classes. The differential expression of 18 genes has been controlled using real-time PCR experiments validating the micro-array analysis. Finally, the integration of the data using Ingenuity Pathways Analysis identified five gene networks providing descriptive elements of the terminal follicular development. Specifically, we observed: (1) the down-expression of ribosomal protein genes, (2) the genes involved in lipid metabolism and (3) the down-expression of cell morphology and ion-binding genes. In conclusion, this study gives new insight into the gene expression during pig terminal follicular growth in vivo and suggested, in particular, a morphological change in pig granulosa cells accompanying terminal follicular growth. Reproduction (2008) 136 211–224 Introduction Those factors control follicular growth either directly or The growth and development of ovarian follicles leading indirectly: for example, bone morphogenetic protein or to ovulation require a series of coordinated events that the insulin-like growth factor (IGF) systems (Monget et al. lead to follicular somatic cell differentiation and oocyte 2002, Mazerbourg et al. 2003, Shimasaki et al. 2004) development. It includes two successive periods, modify the sensitivity of follicular cells to FSH, and preantral (primordial, primary and secondary follicles) epidermal growth factor signalling network potentializes and antral (early antral, antral and preovulatory follicles) LH action (Hattori et al. 1995). Follicular development is follicular developments. Follicles undergo morpho- thus a complex process and requires the coordinate logical and functional changes as they progress towards expression of a large number of genes. The mechanisms ovulation. Among all growing follicles, only a small underlying this development have been intensively proportion of them (less than 1%) will ovulate. Most of studied mainly in granulosa cells because they constitute them undergo a degenerative process called atresia that an important compartment in the mammalian ovarian occurs at the different developmental stages. The antral follicle. They actively participate in the endocrine follicular development depends on a complex regulatory function of the ovaries by secreting oestradiol or network with endocrine regulation by pituitary gonado- progesterone under FSH or LH stimulation (Duda 1997). trophins follicle-stimulating hormone (FSH) and luteiniz- In pigs, in vitro experiments or in situ hybridization ing hormone (LH) but also with autocrine and paracrine have allowed the individual description of gene pathways (Hsueh 1986) including the action of steroids expression in granulosa cells during the growth of antral and peptides (Hillier & Miro 1993, Drummond 2006). follicles. They referred mainly to the role of FSH, the IGF q 2008 Society for Reproduction and Fertility DOI: 10.1530/REP-07-0312 ISSN 1470–1626 (paper) 1741–7899 (online) Online version via www.reproduction-online.org 212 A Bonnet and others system and growth factors on gene expression sequences (35%) corresponded to the whole insert regulation, and the expression of genes involved in (presence of vector flanking sequences). The average steroidogenesis such as CYP19A (Chan & Tan 1987) and insert sequence length was 640 bp. STAR (Balasubramanian et al. 1997). Despite continuing All sequences were deposited in EMBL public progress in the area of ovarian biology, many of the database (accession numbers: CR939025–CR940296; specific mechanisms involved in follicular development, CT971504–CT971571; CT990474–CT990531) and sub- including the initiation of primordial follicle growth, mitted to an assembly process (Sigenae contig assembly; antrum formation, follicular growth/selection and Pig V3 p.sc.3, 30/01/2006). Our 1378 sequences were follicular atresia, remain to be elucidated in greater part of 345 contigs, 63 of which were new sequences. detail. Then, global approaches have been undertaken to The contig analysis revealed a global redundancy of identify new genes involved in antral follicle maturation. 75%. This redundancy was explained by the high Our previous data obtained on porcine granulosa cells proportion of sequences coming from only three genes, in vitro by suppression subtractive hybridization (SSH) or CYP19A, GSTA and PGFS1 representing 15, 13.5 and differential display PCR suggested a role of FSH in 10% of the sequences respectively. There was no contig extracellular matrix synthesis, chromatin remodelling, overlap between the forward and reverse libraries for SF regulation of transcription activity and protection against versus LF comparison and around 10–18% for SF versus atresia (Clouscard-Martinato et al. 1998, Bonnet et al. MF comparison. 2006b). Moreover, different cDNA libraries obtained The PCR products of the 1697 clones (345 genes/con- from whole follicles of different size classes were tigs) resulting from SSH experiments were spotted onto generated to create a catalogue of differentially nylon membrane along with 1056 already sequenced expressed genes along antral follicle development PCR products (954 genes/contigs) to generate our micro- using sequence frequency in each library (Jiang et al. array platform (GEO accession number GPL3978). 2004). However, despite the development of DNA micro-array techniques in the pig (Tuggle et al. 2007), only few transcriptomic data are available concerning Micro-array analysis the ovarian function. In order to identify genes whose expression differs along Thus, using the pig as a model and a transcriptomic the terminal growth and maturation of porcine ovarian approach, the aim of this study was to identify follicles, 14 complex cDNA radiolabelled probes syn- differentially expressed genes along terminal ovarian thesized from the three different size class follicles were follicular development, before the onset of preovulatory hybridized onto our specific cDNA micro-array (1275 LH surge. The experiment has focused especially on genes/contigs). Data (GEO accession number GSE5798 the time of expression of functional LH receptors in dataset) were pre-processed resulting in a list of 1564 granulosa cells that occur in 4–5 mm follicle size expressed cDNAs, corresponding to 515 different genes and classically considered as a maturation marker in or contigs of sequences assembled by SIGENAE. follicle development (May & Schomberg 1984). Then, A mixed linear model was applied to the 1564 cDNAs we developed first a micro-array using cDNAs from to quantify the hybridization signal intensity measured SSH libraries obtained