- Research article
- Open Access
Identification of transcripts involved in meiosis and follicle formation during ovine ovary development
BMC Genomicsvolume 9, Article number: 436 (2008)
The key steps in germ cell survival during ovarian development are the entry into meiosis of oogonies and the formation of primordial follicles, which then determine the reproductive lifespan of the ovary. In sheep, these steps occur during fetal life, between 55 and 80 days of gestation, respectively. The aim of this study was to identify differentially expressed ovarian genes during prophase I meiosis and early folliculogenesis in sheep.
In order to elucidate the molecular events associated with early ovarian differentiation, we generated two ovary stage-specific subtracted cDNA libraries using SSH. Large-scale sequencing of these SSH libraries identified 6,080 ESTs representing 2,535 contigs. Clustering and assembly of these ESTs resulted in a total of 2,101 unique sequences depicted in 1,305 singleton (62.11%) and 796 contigs (37.9%) ESTs (clusters). BLASTX evaluation indicated that 99% of the ESTs were homologous to various known genes/proteins in a broad range of organisms, especially ovine, bovine and human species. The remaining 1% which exhibited any homology to known gene sequences was considered as novel. Detailed study of the expression patterns of some of these genes using RT-PCR revealed new promising candidates for ovary differentiation genes in sheep.
We showed that the SSH approach was relevant to determining new mammalian genes which might be involved in oogenesis and early follicle development, and enabled the discovery of new potential oocyte and granulosa cell markers for future studies. These genes may have significant implications regarding our understanding of ovarian function in molecular terms, and for the development of innovative strategies to both promote and control fertility.
Differentiation of the mammalian gonad is a dynamic process which occurs during embryogenesis and continues throughout an individual's life. During this process, both germ and somatic cells acquire sex-specific characteristics. Meiosis is initiated at different time points in males and females; female germ cells begin meiosis during embryogenesis while male germ cells do not enter meiosis until puberty. Recent findings have indicated that the entry of female germ cells into meiosis is mediated by retinoid responsive genes [1, 2]. A second major event in female gonad differentiation occurs when somatic pre-granulosa cells surround individual germ cells to form primordial follicles which constitute a reservoir of oocytes , available for subsequent growth and differentiation. This event occurs during fetal life or after birth, depending on the species.
During the last decade, substantial progress has been made in identifying the genes that control the assembly and initial growth of ovarian follicles. Different approaches, such as transcriptome analysis of fetal ovaries at different developmental stages and gene inactivation in the mouse, have revealed pathways leading to the formation and differentiation of follicles [4–10]. The transcription of several oocyte-specific genes, such as growth and differentiation factor 9 (Gdf9) and bone morphogenetic protein 15 (Bmp15) commences during early folliculogenesis [11–13]. Regulation of these genes is in part due to the expression of two oocyte-specific transcription factors, Figla  and Nobox . It has been also shown recently that Sohlh1 disruption perturbs follicular formation, partly by causing the down-regulation of Nobox and Figla .
In order to identify novel genes involved in follicle differentiation, microarray-generated expression profiles of the mature ovary have been compared with those of newborns, and several novel genes in ovarian development and follicle formation have thus been discovered [14–16].
By contrast, few genes are known to be related to meiotic prophase I in female mammals. In most cases, meiotic genes were initially discovered in the context of meiotic recombination in yeast. Mammalian counterparts were then discovered, predominantly by means of approaches such as gene-targeting in mice and mutational and chromosomal analysis in human patients with fertility disorders [17–24]. The application of microarray technology has enabled the comparison of transcriptional activity in developing male and female gonads at different stages, some of which contain female germ cells in meiotic prophase I [25–27]. More recently, a microarray study compared genome-wide expression in embryonic mouse ovaries, the aim being to identify novel genes associated with early female meiosis, and reported that around 200 genes were differentially expressed in E13.5 mouse ovaries compared with E11.5 ovaries .
Another means of isolating differentially expressed transcripts consists in using subtraction techniques. A highly effective method, termed suppression subtractive hybridization (SSH), has been applied to generating subtracted cDNA libraries of various systems. This method has also been used to isolate ovary-specific genes in several species [29–32].
In sheep, differences between the somatic components of ovaries and testes are microscopically evident by embryonic day 35 (D35). Germ cells, however, remain indistinguishable between the sexes until D55, when ovarian germ cells initiate prophase of meiosis (between days 49 and 56) and testicular germ cells arrest in the G0/G1 stage of the mitotic cell cycle . The structural integrity of ovigerous cords is maintained until D75, at which time primordial follicles first begin to emerge from the interface of the cortex and medulla . The maximum number of germ cells present within the fetal ovary of sheep is found at D75 (i.e., 800,000). Type 2 follicles (an oocyte with at least one but not two complete layers of cuboidal granulosa cells) are first observed at around D90, and in some cases, type 3 (an oocyte with 2 to <4 layers of cuboidal granulosa cells) and type 4 (an oocyte with 4 to >6 layers of cuboidal granulosa cells with no evidence of antrum formation) follicles are present by D120. Antral follicles are observed in the fetal sheep ovary at around D135, which is about 12 days before parturition . Few studies have described the expression profiles of genes controlling meiosis and follicular development in sheep [36–40].
In this study, we report on the generation of two sheep fetal ovary cDNA libraries, obtained by applying SSH using mRNAs isolated from D55 and D82 female gonads In the first subtraction D55 was used as tester and D82 as driver and the forward subtractive library was constructed (M, Meiosis). In the second one the tester and driver were interchanged and the reverse subtractive library was constructed (F, Folliculogenesis). These two stages correspond to the onset of meiotic prophase I in female germ cells and primordial follicle formation, respectively. 7,296 cDNA clones were selected and 6,080 were sequenced from both libraries, Then 5,025 valid sequences were subjected to database searches to annotate the putative functions of the genes thus represented. Differentially expressed genes were identified by in silico analysis and RT-PCR. Expression analysis confirmed that some of these genes were specifically expressed in fetal ovaries. We further demonstrated that two of the selected ESTs [EMBL:CU651916 and EMBL:CU654655] were expressed in the female gonad during meoitic prophase I period. Furthermore, about 1% (21/2101) of these ovary-expressing transcripts were novel. They were analyzed for their expression in gonads at different developmental stages. The expression patterns of these 21 cDNAs were validated by RT-PCR, which enabled the identification of two transcripts with ovary-specific expression.
Generation of two ovary stage-specific cDNA libraries
Messenger RNAs isolated from D55 and D82 fetal ovaries were used to construct two subtractive cDNA libraries using SSH. The M (meiosis) library, supposed to be enriched in cDNAs representing genes preferentially expressed during prophase I meiosis, was obtained using cDNAs from D55 mRNAs as the 'tester' and cDNAs from D82 mRNAs as the 'driver'. The F (folliculogenesis) library was obtained by reversing the 'tester' and 'driver' mRNA. Sequencing was carried out for 6,080 cDNAs clones selected at random from both libraries: 3,203 from the M library and 2,877 from the F library. After removing poor quality sequences, a total of 2,474 ESTs from the M library and 2,551 from the F library were conserved for further analyses (Table 1). The insert size of the clones ranged from 99 to 1362 base pairs, and the average insert size was approximately 450 base pairs.
A database named M/F db was generated using all ESTs from the M and F libraries . To generate contigs regrouping identical ESTs, two successive criteria were applied. The first step built clusters of all 5,025 sequences from the M/F db and publicly available ovine ESTs, sharing at least 75 bp at an identity rate of 96% using MegaBlast. The second step constructed coherent contigs from the previous clusters using the CAP3 assembly program . This analysis revealed a total of 2,101 independent contigs (82.8% of the total number of contigs) for both libraries, 1,361 contigs only being found in the M library and 740 in the F library (Table 1). These comparisons resulted in the assembly of 2,101 independent sequences including 796 clusters (contigs) and 1,305 singleton ESTs (Fig. 1). The percentage of unique sequences in each library ranged from 42% (F) to 73% (M). This high proportion of sequences represented by single ESTs seemed to indicate an efficient normalization of the libraries; only 3.7% of the genes were represented by more than 5 ESTs. Furthermore, only 8.5% of the sequences were represented by ESTs from both libraries (Table 1), as expected for libraries that were indeed normalized.
Each sequenced clone was analyzed using the BLAST program and classified as known, homologous, or novel after its evaluation against entries within the publicly accessible non-redundant NCBI database (Fig. 2). Known genes were classified as such if homology was high (identity of 90% or more) between the subtracted SSH-generated clones and the Ovis aries database: UniGene Build #11. They accounted for more than 60% in both libraries (62.2% and 65.3% in the M and F libraries, respectively).
Contigs were classified as homologous sequences if they possessed homology with bovine or human gene sequences, relative to entries in the Bos taurus: UniGene Build #84, Homo sapiens: UniGene Build #202 NCBI database. Homologous contigs accounted for 36.7% in the M library and 33.9% in the F library. Novel ESTs were designated as such if no significant match was obtained in the different databases with known genes (Fig. 2). Fifteen were found without gene annotation in the M library (15/1361, 1.1%) and six in the F library (6/740, 0.8%). The two most represented genes in the M library were Cytochrome C oxidase subunit 3 (EC 184.108.40.206) and Bos taurus hypothetical protein, LOC783526, with 13 and 12 ESTs respectively. Eleven genes were found more than 10 times in the F library [see Additional file 1].
Functional annotation of ESTs
Comparison of the 2,101 sequences from both SSH libraries with the non-redundant protein UniProt/SwissProt database revealed that 1,173 contigs (55.8%) displayed significant similarity (<1E-5) to existing protein sequences. Furthermore, among these 2,101 contigs, only 543 had one GO term (25.8%). These 543 sequences fell into the three principal GO categories: cellular component, biological process and molecular function. In both libraries, the three most widely represented biological processes were intracellular signaling cascade (40–45%), cell surface receptor (20–24%) and DNA repair (13–21%) (Fig. 3). However, the small number of GO annotated sequences restricted functional annotation of our libraries.
Identification of differentially expressed ovary-specific transcripts
In order to evaluate the specificity of the transcripts involved in both the developmental processes studied here, i.e. meiosis and follicle formation, individual library specific ESTs were studied. A bioinformatics comparison of our SSH library sequences with those published in the literature showed that at least 10 genes previously described as being expressed in mice fetal ovaries or adult male meiotic germ cells were present in our libraries. The detection of several genes already known to be expressed during meiosis, such as Stk31 [EMBL:CU655270], Smc3 [EMBL:CU653776], Dmrt7 [EMBL:CU652558], Zp148 [EMBL:CU655301], Mutl1 [EMBL:CU654415], Dazl [EMBL:CU654830], Stathmin [EMBL:CU652110], Msh2 [EMBL:CU654850] and Msh5 [EMBL:CU654104] in the M library increased confidence in the quality of the subtraction procedure. In a similar manner, sequences corresponding to the Foxo3A [EMBL:CU638379], Lhx8 [EMBL:CU637755] and BMP4 [EMBL:CU653220] genes involved in folliculogenesis were present in the F library.
RT-PCR screening was used to validate in our animal model the ovary-specific expression of several genes previously described in the context of mouse gonad transcriptome analysis. This screening contained (i) a mixture of cDNAs from four different fetal somatic tissues (kidney, liver, heart and brain) and (ii) cDNA from male and female gonads at two different developmental stages (D55 and D82). We examined seven genes from the M library using this mini-profiling procedure, with Dmc1 as the positive control of the meiotic gene (Fig. 4). Five genes (TEX11, TEX14, MOV10L1, Maelstrom, STK31) produced a band-pattern consistent with the profile expected for a gonad-specific transcript or a predominantly gonad-expressed transcript. Three of them (TEX11, MOV10L1, STK31) were more strongly expressed in female gonads than in males (Fig. 4). In order to confirm their specific expression pattern, TEX11 was subjected to further real-time quantitative RT-PCR expression analysis on a panel of several fetal stages for each sex. Its expression profile appeared to be specific to meiotic prophase I and was similar to typical meiotic genes such as DMC1 (Fig. 5).
Interestingly, 21 novel transcripts were also detected in the two libraries. Primers were designed for each of them [see Additional file 2] and RT-PCR analysis was performed using mini profiles (Fig. 6). The expression patterns of these 21 transcripts were classified into three categories: the first corresponded to transcripts expressed in both male and female gonads and somatic tissues, the second to transcripts specifically expressed in gonads of both sexes, and the third to those specifically expressed in the ovary (Fig. 6). Transcripts presenting a differential expression pattern between sexes or stages were subjected to real-time quantitative RT-PCR expression analysis (Fig. 7 and data not shown). This analysis showed that three transcripts were stage and sex -specifically expressed.
Use of the SSH method to study meiosis and early folliculogenesis
Of the different strategies available to study transcriptome variations between two biological situations, Suppression Subtractive Hybridization (SSH)  is an efficient method to identify genes that are differentially regulated. Because it includes a normalization step, it enables the isolation of low abundance, differentially expressed transcripts. It does not require previous sequence knowledge, which is an advantage when the genome of the studied species has not been sequenced, and it can be started using PCR-amplified cDNAs. It is thus particularly well suited to biological situations such as embryonic or fetal development where specific genes are expressed and tiny amounts of RNA are available .
The present study aimed to analyze changes in the gene expression profile of sheep ovary during development using the SSH technique. The choice of SSH was based on the assumption that we would be able to discover factors expressed at lower levels than the genes classically detected by microarray techniques, particularly since commercial microarrays are only available for a small number of species that does not include sheep.
M and F cDNA-libraries were generated by SSH using mRNAs isolated from fetal ovaries at two developmental stages, and a total of 5,025 ESTs from both libraries were found to represent a maximum of 2,535 genes. The level of redundancy in the libraries was approximately 50%.
We used an approach based on the systematic sequencing of SSH products and gene analysis by functional categories. Out of the 2,101 uncommon ESTs that were cloned during this study, about 82% corresponded to contigs and singletons found to be unique, demonstrating the effectiveness of the experimental design utilized for these SSH libraries. Of the two libraries, the one dedicated to meiosis was of better quality than the F library, i.e. it contained more independent contigs (1578 versus 957), more non redundant contigs (1361/740) and more singletons (994/311).
In total, 2,101 non-redundant contigs were sequenced and analyzed. Nearly all of them presented a UniGene nucleotidic annotation (99%) and around two-thirds of the sequences were present in sheep genome databases, while one third presented homology with human or bovine genome databases. These results should be placed in context insofar as in the 99% of sequences with nucleic annotation, a large proportion displayed no protein annotation or possessed a hypothetical or unknown function. Major efforts regarding annotation are now necessary in order to improve our understanding of the functional role of these clones.
The cDNAs isolated from the M library included several transcripts that had previously been reported to be involved in the meiosis process . Examples included STK31 , SMC3 , DMRT7 , ZFP148 MSH5 , MSH2 , Stathmin , DAZL A , and MUTL1 . Similar data had been observed regarding the F library, with the identification of FOXO3A , SOHLH1  and BMP4 . The detection of previously identified genes supported the validity of our experimental model and the use of SSH as the analytical method. Although the utilization of SSH in the current study was successful in identifying previously characterized genes, several expected genes were not present within the target cDNA library. Among these, Stra8 [1, 24], DMC1  for meiosis genes and BMP15  and Figla , for folliculogenesis genes were not found within our subtracted libraries. The absence of these genes from the libraries may have been due to the use of PCR, and some of the cDNAs may not have been amplified as efficiently as others and thus have been lost from the SSH starting material. A similar inability to identify all expected known genes after SSH was recently reported by others and ascribed to incomplete representation of the total cDNA repertoire . We can also suppose that the plating and sequencing of a larger number of clones could allow isolation of weakly expressed transcripts.
In addition, 1% of the transcripts thus isolated were unknown. They may correspond to additional exons of previously known genes, resulting from differential splicing specific to ovine species or to new genes not yet described (e.g. corresponding to non coding RNAs). Further 5' and 3'RACE analyses will be necessary to obtain full-length transcripts and to perform other sequence comparisons and functional analysis.
Analysis of functional gene ontology
In order to understand the pathways involved in the initiation of prophase I meiosis and follicle formation, GO was used to cluster genes according to their biological function. The resulting classification mainly highlighted three biological processes, each representing more than 10% of the contigs taken into account: intracellular signaling cascade, cell surface receptors and DNA repair. These results should be modulated by the fact that only 25% of the contigs analyzed displayed GO annotation. However, this partial analysis was consistent with the fact that signaling molecules and DNA repair proteins are known to be important to regulating germ cell development and meiosis [56, 57].
RT-PCR analysis of three (STK31, MOV10L1, TEX11) of seven selected ESTs from the M library confirmed their differential expression between males and females. These three transcripts exhibited regulation in the sex, organ and stage from which they were cloned; i.e. female, ovary and meiosis prophase I stage. Real-time quantitative RT-PCR analyses of two of them (MOV10L1 [EMBL:CU654655] and TEX11 [EMBL:CU651916]) revealed that they were over-expressed in ovaries from days 55 to 82 and presented similar patterns to DMC1, a typical meiotic gene. This suggested that the M library contained ESTs representing genes differentially expressed in meiosis prophase 1 period. TEX11, initially described as being expressed during adult spermatogenesis, has recently been demonstrated on female meiotic chromosomes in mice , findings in agreement with our results. Moreover, it is noteworthy that MOV10L1 transcripts have never been described in female fetal gonads; they were previously described as being expressed in adult testis . It is also interesting to note that the expression profile of the ovine Pecanex 1 gene differed from that described in the mouse . It had previously been shown that the mouse Pecanex 1 transcript was only present in the germ line and not in the somatic cells of the testis, reaching its peak at the pachytene stage of meiotic prophase. However, in sheep, we detected Pecanex 1 in the fetal testis when male germ cells were not in meiosis.
The expression of unknown transcripts revealed that 20% of them were specific to gonads and 14% to ovary. Two of these unknown transcripts from the M library were more strongly expressed in the ovary and not expressed in the testis during fetal life. Their expression profiles were similar to those described in the goat for aromatase, activin β genes and 3β HSD , with a high level of transcripts at around D55 followed by a decline and then a further rise at around D75–82. This similarity of profiles suggests that these unknown genes may be involved in the regulation of steroid production. By refining our methods for sequence homology analysis, we found nucleotidic sequence homology between both CU652878 [EMBL:CU652878] and CU638235 [EMBL:CU638235] and two distinct bovine chromosome 6 genomic contigs (NW_001495163 and NW_001495159, respectively). No genes have been localized near these homologous regions.
This report also presents a medium scale sequencing of cDNAs representing genes in sheep. Indeed, of the 2,535 non-redundant ESTs selected from the M and F libraries, 771 were not reported in the Ovis aries EST database collection and therefore constituted an unique contribution to the Ovis aries EST pool. Though modest by comparison with databases for some other ruminants, such as cattle, our sequence collection for sheep may serve as a source for comparative studies in the bovids, an animal family of considerable importance to agronomic research.
Using normalized cDNA subtraction as a transcript profiling tool to identify differentially expressed transcripts, we isolated several genes that may contribute toward understanding the mechanisms involved in meiosis and ovarian follicle differentiation in a mono-ovulatory species. Defects in these two processes are a leading cause of both infertility and birth defects. Further analyses of some selected genes in terms of their spatiotemporal expression patterns at the mRNA and protein levels in the developing ovarian follicle will be necessary to gain a clearer understanding of their functional roles.
In addition, 60 mer oligonucleotides from each of the 2101 unique ESTs will be designed and thus constitute a custom-made macroarray dedicated to ovine ovary differentiation. This tool will be used to evaluate modifications to gene expression under various nutritional and environmental conditions and to analyze their impacts on ovarian function.
Animals and tissue samples
Procedures for handling animals during this study were carried out in accordance with the INRA guidelines for the Care and Use of Agricultural Animals. The animals came from the INRA experimental farm at Brouessy, UCEA . Cyclic Préalpes ewes were synchronized using an intravaginal progestogen sponge. Animals were inseminated at Day 0. Pregnant female tracts at different developmental stages were collected at slaughter and rapidly dissected to extract fetuses. Fetuses were collected on Days 44, 55, 60, 65, 75, 82, 90 and 114 of gestation. The gonads of adult animals were collected at slaughter. After puberty, gonads from both sexes were detached, then cut into small pieces and frozen immediately in liquid nitrogen (about 1/4 for the ovary and 1/20 for the testis). Other organs, such as the liver, kidney, brain and heart, were also collected, frozen immediately and then stored at -80°C.
Total RNA was extracted from each sample using TRIzol LS reagent (Invitrogen Life Technologies, Cergy-Pontoise, France), according to the method described by Chomczynski and Sacchi . The quantitative yield of purification and the quality of RNAs were checked by spectrophotometrical scanning (Nanodrop ND1000 LABTECH) and by examination following 0.8% agarose gel electrophoresis to ensure equal purity and an absence of degradation in all samples.
One hundred micrograms of total RNA from each developmental stage (D55 and D82) were extracted, representing respectively the "meiosis" and "follicle formation" stages. PolyA+ mRNA were then fractionated from total RNA using the NucleoTrap mRNA Kit (BD Biosciences Clontech, Palo Alto, CA, USA) which combines a spin-column filter with oligo (dT) latex bead technology. Poly A+ RNA was treated with 40 U of RNase-free DNaseI (Roche, Saclay, France) for 2 h at 37°C to prevent genomic DNA contamination. After phenol/chloroform extraction and ethanol precipitation, the DNase-treated poly A+ RNA was resuspended in RNase-free water.
cDNAs were synthesized from 2 μg of polyA+ using 20 U AMV Reverse Transcriptase (BD Biosciences Clontech, Palo Alto, CA, USA) to generate the first cDNA-strand. Second-cDNA stands were produced with the DNA polymerase I (6 U/μl) (BD Biosciences Clontech, Palo Alto, CA, USA) provided in the kit. RsaI was added to these synthesized DS cDNA which were then incubated at 37°C. Tester and driver double strain cDNAs were prepared from each pool of gonads. According to the manufacturer's instructions, SSH experiments  were performed with the following steps: (i) Rsa I digestion of tester and driver cDNA; (ii) adaptor ligation (8 h at 68°C) to the testers: two tester populations were created with different adaptors (1 and 2R), driver cDNA had no adaptors; (iii) first hybridization: each tester cDNA was hybridized at 68°C overnight with an excess of the counterpart driver cDNA; (iv) second hybridization: the 2 samples from the first hybridization were newly hybridized with fresh driver cDNA; (v) two consecutive PCR amplification steps with primers designed according to the adaptor sequence were performed in order to enrich the differentially expressed sequences and reduce background. The first comprised 27 cycles (94°C for 30 sec, 66°C for 30 sec, 72°C for 1 min 30 sec) using the end of the adaptors as primers. The second PCR comprised 10 or 12 cycles (F and M libraries, respectively) using the nested primers (included in adaptors 1 and 2R) (94°C for 30 sec, 68°C for 30 sec, 72°C for 1 min 30 sec). All the primers used during SSH experiments are listed in additional file 3 [see Additional file 3].
The efficiency of subtraction was analyzed using PCR by comparison with the abundance of cDNAs before and after subtraction on the constitutively expressed gene G3PDH (94°C for 30 sec, 60°C for 30 sec, 72°C for 1 min 30 sec).
Five microliters aliquots were collected from 50-μl PCR reactions after 18, 23, 28, 33 and 38 cycles. Amplification was quantified by densitometry. The efficiency of subtraction was confirmed by comparing the abundance of cDNAs before and after subtraction using PCR with bovine GAPDH [see Additional file 4].
cDNA cloning, sequencing, data analysis and dbEST submission
The subtracted PCR products generated by SSH were purified using the GFX kit (GE Healthcare, Orsay, France), cloned into the pGEMT easy vector (Promega France, Charbonnière les Bains, France) and then transformed into Escherichia coli DH5-alpha T1 Max Efficiency Cells (Invitrogen, Cergy-Pontoise, France). Positive clones were picked out according to blue/white selection and 7296 colonies were transferred robotically into 384 well plates at the INRA National Biological Resources Centre for Animal Genomics , where they were grown in 2xYT medium overnight at 37°C. All individual colonies were spotted on nylon membrane and then fixed with NaOH, neutralized by Tris-HCl and stored at 4°C prior to hybridization.
Clones were randomly selected and single-pass sequenced with a vector primer from sequencing companies (Genome Express,  and MWG, ). On average, 500 clones were sequenced and then hybridized on the macro-membrane containing the libraries in order to identify any redundant clones. After hybridization, another batch of 500 clones was sequenced, hybridized and so on.
Using SIGENAE bioinformatics facilities , EST quality was computed using a home-made workflow system and the ESTs were then cleaned of vector and primer sequences and sequences containing contaminants such as E. Coli, Yeast, Mitochondria, Ribosome. SSH cDNA libraries usually contained insertion sizes within a range of 200–800 bp. A sequence was only selected for further analysis if it had a valid status, i.e. containing at least 100 bp outside repeats with a phred score over 20, without any vectors and adaptors so as to prevent chimerical contig linking sequences, for instance because of unclipped vectors.
Contigs and EST singletons of SIGENAE Sheep clustering 5  were annotated for sequence homologies by running BLASTX or BLASTN against the following databases for different mammalian species (human, bovine, ovine): UniProtKB/Swiss-Prot Release 53.0 of 29-May-2007, UniProtKB/TrEMBL Release 36.0 of 29-May-2007, ProDom Release CG67, Sigenae Sheep Contigs V4, TIGR Cattle BtGI 12, UniGene Cattle Build #84, UniGene Human Build #202, UniGene Sheep Build #11, Ensembl Cattle Transcripts Btau_3.1 44 and Ensembl Human Transcripts NCBI36 44.
Only sequences with a PHRED score over 20 on at least 100 bp were released in the EST division of the EMBL-EBI (European Molecular Biology Laboratory – European Bioinformatics Institute) Nucleotide Sequence Database with accession numbers.
Accession numbers in the EMBL Nucleotide Sequence Database ran from CU637754 to CU638663 and from CU651679 to CU655811, except for CU652171, CU652212, CU652215, CU652218, CU652484, CU652528, CU652531, CU652642, CU652698, CU652701, CU652704, CU652776, CU652780, CU652834, CU652837, CU652840, CU652843 and CU652893.
Contig similarity search and functional assignments
Similarity searches of each contig sequence were performed using BLASTX or BLASTN programs with different cut-off e-value criteria, depending on the databases used. These e-values were 1e-02 with the ProDom, TIGR and UniGene databases; 1e-05 against UniProt databases; 1e-10 against Ensembl transcripts and 1e-30 against other SIGENAE contigs. Functional annotation was carried out using Gene Ontology . Next, the GO categorization scheme of classification by biological process, cellular component and molecular function was used to generate representative pie diagrams .
Semi-quantitative reverse transcription-PCR (Semi-quantitative RT-PCR)
Total RNA was isolated from the ovaries, testis, liver, heart, brain and kidney of fetal or adult origin using the TRIzol LS reagent (Invitrogen Life Technologies, Cergy-Pontoise, France) according to the manufacturer's instructions. RNA extraction, DNase treatment, and cDNA synthesis were carried out as previously described . Briefly, 1 μl (corresponding to 0.25 μg reverse transcribed total RNA) of each RT mix was amplified in 25 μl PCR using 0.5 U Taq polymerase (TaKaRa, Lonza, Verviers, Belgium), 200 μM of each dNTP, and 150 nM of each primer. The sequences of the primers used and the PCR conditions utilized are shown in additional files 2 and 5 [see additional file 2], [see additional file 5].
The GAPDH gene was used as a cDNA quantity control. PCR conditions were as follows: 25 cycles of 94°C for 30 sec, 58°C for 30 s and 72°C for 30 sec, followed by a single cycle at 72°C for 10 min.
Real-time quantitative RT-PCR of differentially expressed transcripts
cDNA (50 ng) from ovary and testis was obtained by the reverse transcription of 5 μg RNA for each sample using Superscript II (Invitrogen, Cergy-Pontoise, France) and mixed with 12.5 μl (1×) of ABsolute blue QPCR SYBR Green ROX mix (Abgene, Les Ulis, France), 300 nM from each primer [see Additional file 6], in a total volume of 25 μl. The reaction mixture was finally transferred into a 96-well optical reaction plate, sealed with appropriate optical caps, and run on the ABI Prism 7700 HT apparatus (Applied Biosystems). Standard controls of both the specificity and efficiency of the qPCR assays were performed. All expression data were normalized to the HPRT1 expression level. The concentration of the target genes were divided by the concentration of the HPRT1 gene thus showing the relative expression of the mRNA in the tissue. This calculation was performed for all stages analyzed. The sequences of the primers used, and the qPCR conditions prevailing, are shown in additional file 6 [see Additional file 6]. All qPCR experiments are made in duplicate and for some stages (D60, 82, 90 and adult), different animals, RNA extractions, and RT were used. SEM were calculated when two independent qPCR from 2 different biological materials were done.
Koubova J, Menke DB, Zhou Q, Capel B, Griswold MD, Page DC: Retinoic acid regulates sex-specific timing of meiotic initiation in mice. Proc Natl Acad Sci USA. 2006, 103 (8): 2474-2479. 10.1073/pnas.0510813103.
Bowles J, Koopman P: Retinoic acid, meiosis and germ cell fate in mammals. Development. 2007, 134 (19): 3401-3411. 10.1242/dev.001107.
Anderson LD, Hirshfield AN: An overview of follicular development in the ovary: from embryo to the fertilized ovum in vitro. Md Med J. 1992, 41 (7): 614-620.
Joshi S, Davies H, Sims LP, Levy SE, Dean J: Ovarian gene expression in the absence of FIGLA, an oocyte-specific transcription factor. BMC Dev Biol. 2007, 7: 67-10.1186/1471-213X-7-67.
Soyal SM, Amleh A, Dean J: FIGalpha, a germ cell-specific transcription factor required for ovarian follicle formation. Development. 2000, 127 (21): 4645-4654.
Pangas SA, Choi Y, Ballow DJ, Zhao Y, Westphal H, Matzuk MM, Rajkovic A: Oogenesis requires germ cell-specific transcriptional regulators Sohlh1 and Lhx8. Proc Natl Acad Sci USA. 2006, 103 (21): 8090-8095. 10.1073/pnas.0601083103.
Lin YN, Matzuk MM: High-throughput discovery of germ-cell-specific genes. Semin Reprod Med. 2005, 23 (3): 201-212. 10.1055/s-2005-872448.
Rajkovic A, Pangas SA, Ballow D, Suzumori N, Matzuk MM: NOBOX deficiency disrupts early folliculogenesis and oocyte-specific gene expression. Science. 2004, 305 (5687): 1157-1159. 10.1126/science.1099755.
Uda M, Ottolenghi C, Crisponi L, Garcia JE, Deiana M, Kimber W, Forabosco A, Cao A, Schlessinger D, Pilia G: Foxl2 disruption causes mouse ovarian failure by pervasive blockage of follicle development. Hum Mol Genet. 2004, 13 (11): 1171-1181. 10.1093/hmg/ddh124.
Durlinger AL, Gruijters MJ, Kramer P, Karels B, Ingraham HA, Nachtigal MW, Uilenbroek JT, Grootegoed JA, Themmen AP: Anti-Mullerian hormone inhibits initiation of primordial follicle growth in the mouse ovary. Endocrinology. 2002, 143 (3): 1076-1084. 10.1210/en.143.3.1076.
Dong J, Albertini DF, Nishimori K, Kumar TR, Lu N, Matzuk MM: Growth differentiation factor-9 is required during early ovarian folliculogenesis. Nature. 1996, 383 (6600): 531-535. 10.1038/383531a0.
Dean J: Oocyte-specific genes regulate follicle formation, fertility and early mouse development. J Reprod Immunol. 2002, 53 (1–2): 171-180. 10.1016/S0165-0378(01)00100-0.
Elvin JA, Yan C, Matzuk MM: Oocyte-expressed TGF-beta superfamily members in female fertility. Mol Cell Endocrinol. 2000, 159 (1–2): 1-5. 10.1016/S0303-7207(99)00185-9.
Herrera L, Ottolenghi C, Garcia-Ortiz JE, Pellegrini M, Manini F, Ko MS, Nagaraja R, Forabosco A, Schlessinger D: Mouse ovary developmental RNA and protein markers from gene expression profiling. Dev Biol. 2005, 279 (2): 271-290. 10.1016/j.ydbio.2004.11.029.
Fayad T, Levesque V, Sirois J, Silversides DW, Lussier JG: Gene expression profiling of differentially expressed genes in granulosa cells of bovine dominant follicles using suppression subtractive hybridization. Biol Reprod. 2004, 70 (2): 523-533. 10.1095/biolreprod.103.021709.
Kocabas AM, Crosby J, Ross PJ, Otu HH, Beyhan Z, Can H, Tam WL, Rosa GJ, Halgren RG, Lim B: The transcriptome of human oocytes. Proc Natl Acad Sci USA. 2006, 103 (38): 14027-14032. 10.1073/pnas.0603227103.
Critchlow HM, Payne A, Griffin DK: Genes and proteins involved in the control of meiosis. Cytogenet Genome Res. 2004, 105 (1): 4-10. 10.1159/000078002.
Pittman DL, Cobb J, Schimenti KJ, Wilson LA, Cooper DM, Brignull E, Handel MA, Schimenti JC: Meiotic prophase arrest with failure of chromosome synapsis in mice deficient for Dmc1, a germline-specific RecA homolog. Mol Cell. 1998, 1 (5): 697-705. 10.1016/S1097-2765(00)80069-6.
Yoshida K, Kondoh G, Matsuda Y, Habu T, Nishimune Y, Morita T: The mouse RecA-like gene Dmc1 is required for homologous chromosome synapsis during meiosis. Mol Cell. 1998, 1 (5): 707-718. 10.1016/S1097-2765(00)80070-2.
Edelmann W, Cohen PE, Kneitz B, Winand N, Lia M, Heyer J, Kolodner R, Pollard JW, Kucherlapati R: Mammalian MutS homologue 5 is required for chromosome pairing in meiosis. Nat Genet. 1999, 21 (1): 123-127. 10.1038/5075.
Kneitz B, Cohen PE, Avdievich E, Zhu L, Kane MF, Hou H, Kolodner RD, Kucherlapati R, Pollard JW, Edelmann W: MutS homolog 4 localization to meiotic chromosomes is required for chromosome pairing during meiosis in male and female mice. Genes Dev. 2000, 14 (9): 1085-1097.
Baudat F, Manova K, Yuen JP, Jasin M, Keeney S: Chromosome synapsis defects and sexually dimorphic meiotic progression in mice lacking Spo11. Mol Cell. 2000, 6 (5): 989-998. 10.1016/S1097-2765(00)00098-8.
Mandon-Pepin B, Touraine P, Kuttenn F, Derbois C, Rouxel A, Matsuda F, Nicolas A, Cotinot C, Fellous M: Genetic investigation of four meiotic genes in women with premature ovarian failure. Eur J Endocrinol. 2008, 158 (1): 107-115. 10.1530/EJE-07-0400.
Baltus AE, Menke DB, Hu YC, Goodheart ML, Carpenter AE, de Rooij DG, Page DC: In germ cells of mouse embryonic ovaries, the decision to enter meiosis precedes premeiotic DNA replication. Nat Genet. 2006, 38 (12): 1430-1434. 10.1038/ng1919.
Grimmond S, Van Hateren N, Siggers P, Arkell R, Larder R, Soares MB, de Fatima Bonaldo M, Smith L, Tymowska-Lalanne Z, Wells C: Sexually dimorphic expression of protease nexin-1 and vanin-1 in the developing mouse gonad prior to overt differentiation suggests a role in mammalian sexual development. Hum Mol Genet. 2000, 9 (10): 1553-1560. 10.1093/hmg/9.10.1553.
Nef S, Schaad O, Stallings NR, Cederroth CR, Pitetti JL, Schaer G, Malki S, Dubois-Dauphin M, Boizet-Bonhoure B, Descombes P: Gene expression during sex determination reveals a robust female genetic program at the onset of ovarian development. Dev Biol. 2005, 287 (2): 361-377.
Small CL, Shima JE, Uzumcu M, Skinner MK, Griswold MD: Profiling gene expression during the differentiation and development of the murine embryonic gonad. Biol Reprod. 2005, 72 (2): 492-501. 10.1095/biolreprod.104.033696.
Olesen C, Nyeng P, Kalisz M, Jensen TH, Moller M, Tommerup N, Byskov AG: Global gene expression analysis in fetal mouse ovaries with and without meiosis and comparison of selected genes with meiosis in the testis. Cell Tissue Res. 2007, 328 (1): 207-221. 10.1007/s00441-006-0205-5.
Miyakoshi K, Murphy MJ, Yeoman RR, Mitra S, Dubay CJ, Hennebold JD: The identification of novel ovarian proteases through the use of genomic and bioinformatic methodologies. Biol Reprod. 2006, 75 (6): 823-835. 10.1095/biolreprod.106.052290.
Bonnet A, Frappart PO, Dehais P, Tosser-Klopp G, Hatey F: Identification of differential gene expression in in vitro FSH treated pig granulosa cells using suppression subtractive hybridization. Reprod Biol Endocrinol. 2006, 4: 35-10.1186/1477-7827-4-35.
Hourvitz A, Gershon E, Hennebold JD, Elizur S, Maman E, Brendle C, Adashi EY, Dekel N: Ovulation-selective genes: the generation and characterization of an ovulatory-selective cDNA library. J Endocrinol. 2006, 188 (3): 531-548. 10.1677/joe.1.06231.
Ndiaye K, Fayad T, Silversides DW, Sirois J, Lussier JG: Identification of downregulated messenger RNAs in bovine granulosa cells of dominant follicles following stimulation with human chorionic gonadotropin. Biol Reprod. 2005, 73 (2): 324-333. 10.1095/biolreprod.104.038026.
Sawyer HR, Smith P, Heath DA, Juengel JL, Wakefield SJ, McNatty KP: Formation of ovarian follicles during fetal development in sheep. Biol Reprod. 2002, 66 (4): 1134-1150. 10.1095/biolreprod66.4.1134.
Smith P, O WS, Hudson NL, Shaw L, Heath DA, Condell L, Phillips DJ, McNatty KP: Effects of the Booroola gene (FecB) on body weight, ovarian development and hormone concentrations during fetal life. J Reprod Fertil. 1993, 98 (1): 41-54.
McNatty KP, Quirke QL, Fidler AE, Smith P, Heath DA, Tisdall D: Ovarian development: fetus to puberty. Hormones and Women's Health. Edited by: Salamonsen LA. 2000, Gordon and Breach Science Publishers, 9-22.
Juengel JL, Heath DA, Quirke LD, McNatty KP: Oestrogen receptor alpha and beta, androgen receptor and progesterone receptor mRNA and protein localisation within the developing ovary and in small growing follicles of sheep. Reproduction. 2006, 131 (1): 81-92. 10.1530/rep.1.00704.
Mandon-Pepin B, Oustry-Vaiman A, Vigier B, Piumi F, Cribiu E, Cotinot C: Expression profiles and chromosomal localization of genes controlling meiosis and follicular development in the sheep ovary. Biol Reprod. 2003, 68 (3): 985-995. 10.1095/biolreprod.102.008557.
Logan KA, McNatty KP, Juengel JL: Expression of wilms' tumor gene and protein localization during ovarian formation and follicular development in sheep. Biol Reprod. 2003, 68 (2): 635-643. 10.1095/biolreprod.102.008946.
Quirke LD, Juengel JL, Tisdall DJ, Lun S, Heath DA, McNatty KP: Ontogeny of steroidogenesis in the fetal sheep gonad. Biol Reprod. 2001, 65 (1): 216-228. 10.1095/biolreprod65.1.216.
Tisdall DJ, Fidler AE, Smith P, Quirke LD, Stent VC, Heath DA, McNatty KP: Stem cell factor and c-kit gene expression and protein localization in the sheep ovary during fetal development. J Reprod Fertil. 1999, 116 (2): 277-291.
Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.
Diatchenko L, Lau YF, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci USA. 1996, 93 (12): 6025-6030. 10.1073/pnas.93.12.6025.
Bui LC, Leandri RD, Renard JP, Duranthon V: SSH adequacy to preimplantation mammalian development: scarce specific transcripts cloning despite irregular normalisation. BMC Genomics. 2005, 6: 155-10.1186/1471-2164-6-155.
Wang PJ, McCarrey JR, Yang F, Page DC: An abundance of X-linked genes expressed in spermatogonia. Nat Genet. 2001, 27 (4): 422-426. 10.1038/86927.
Prieto I, Tease C, Pezzi N, Buesa JM, Ortega S, Kremer L, Martinez A, Martinez AC, Hulten MA, Barbero JL: Cohesin component dynamics during meiotic prophase I in mammalian oocytes. Chromosome Res. 2004, 12 (3): 197-213. 10.1023/B:CHRO.0000021945.83198.0e.
Kim S, Namekawa SH, Niswander LM, Ward JO, Lee JT, Bardwell VJ, Zarkower D: A mammal-specific Doublesex homolog associates with male sex chromatin and is required for male meiosis. PLoS Genet. 2007, 3 (4): e62-10.1371/journal.pgen.0030062.
Takeuchi A, Mishina Y, Miyaishi O, Kojima E, Hasegawa T, Isobe K: Heterozygosity with respect to Zfp148 causes complete loss of fetal germ cells during mouse embryogenesis. Nat Genet. 2003, 33 (2): 172-176. 10.1038/ng1072.
Qin J, Baker S, Te Riele H, Liskay RM, Arnheim N: Evidence for the lack of mismatch-repair directed antirecombination during mouse meiosis. J Hered. 2002, 93 (3): 201-205. 10.1093/jhered/93.3.201.
Amat JA, Fields KL, Schubart UK: Stage-specific expression of phosphoprotein p19 during spermatogenesis in the rat. Mol Reprod Dev. 1990, 26 (4): 383-390. 10.1002/mrd.1080260414.
Saunders PT, Turner JM, Ruggiu M, Taggart M, Burgoyne PS, Elliott D, Cooke HJ: Absence of mDazl produces a final block on germ cell development at meiosis. Reproduction. 2003, 126 (5): 589-597. 10.1530/rep.0.1260589.
Eaker S, Cobb J, Pyle A, Handel MA: Meiotic prophase abnormalities and metaphase cell death in MLH1-deficient mouse spermatocytes: insights into regulation of spermatogenic progress. Dev Biol. 2002, 249 (1): 85-95. 10.1006/dbio.2002.0708.
Hosaka T, Biggs WH, Tieu D, Boyer AD, Varki NM, Cavenee WK, Arden KC: Disruption of forkhead transcription factor (FOXO) family members in mice reveals their functional diversification. Proc Natl Acad Sci USA. 2004, 101 (9): 2975-2980. 10.1073/pnas.0400093101.
Juengel JL, Reader KL, Bibby AH, Lun S, Ross I, Haydon LJ, McNatty KP: The role of bone morphogenetic proteins 2, 4, 6 and 7 during ovarian follicular development in sheep: contrast to rat. Reproduction. 2006, 131 (3): 501-513. 10.1530/rep.1.00958.
McNatty KP, Juengel JL, Reader KL, Lun S, Myllymaa S, Lawrence SB, Western A, Meerasahib MF, Mottershead DG, Groome NP: Bone morphogenetic protein 15 and growth differentiation factor 9 co-operate to regulate granulosa cell function in ruminants. Reproduction. 2005, 129 (4): 481-487. 10.1530/rep.1.00517.
Li GM: Mechanisms and functions of DNA mismatch repair. Cell Res. 2008, 18 (1): 85-98. 10.1038/cr.2007.115.
Best D, Sahlender DA, Walther N, Peden AA, Adams IR: Sdmg1 is a conserved transmembrane protein associated with germ cell sex determination and germline-soma interactions in mice. Development. 2008, 135 (8): 1415-1425. 10.1242/dev.019497.
Yang F, Gell K, Heijden van der GW, Eckardt S, Leu NA, Page DC, Benavente R, Her C, Hoog C, McLaughlin KJ: Meiotic failure in male mice lacking an X-linked factor. Genes Dev. 2008, 22 (5): 682-691. 10.1101/gad.1613608.
Geisinger A, Alsheimer M, Baier A, Benavente R, Wettstein R: The mammalian gene pecanex 1 is differentially expressed during spermatogenesis. Biochim Biophys Acta. 2005, 1728 (1–2): 34-43.
Pannetier M, Fabre S, Batista F, Kocer A, Renault L, Jolivet G, Mandon-Pepin B, Cotinot C, Veitia R, Pailhoux E: FOXL2 activates P450 aromatase gene transcription: towards a better characterization of the early steps of mammalian ovarian development. J Mol Endocrinol. 2006, 36 (3): 399-413. 10.1677/jme.1.01947.
Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162 (1): 156-159. 10.1016/0003-2697(87)90021-2.
CRB GADIE. [http://w3.jouy.inra.fr/unites/lreg/CRB/index.html]
Genome Express. [http://www.genome-express.com/]
Gene Ontology. [http://www.geneontology.org/GO.downloads.shtm]
Gene Ontology Annotation Plotting. [http://wego.genomics.org.cn/]
We would like to thank Céline Ducroix-Crépy (CRB-GADIE Jouy-en-Josas, France) for library arraying and bacteria membrane preparations and Véronique Duranthon (UMR-BDR, Jouy-en-Josas, France) for helpful discussions. We also are indebted to the UCEA members responsible for our sheep breed. All sequence analyses were carried out in fruitful collaboration with the AGENAE bioinformatics team (SIGENAE). Part of this work was financed by an INRA grant: AIP sequençage n° P00258.
AB is supported by a doctoral fellowship from the Ministère délégué à l'Enseignement Supérieur et à la Recherche, France.
The authors declare that they have no competing interests.
AB carried out screening procedures, conducted all expression validations, performed gene ontology annotations and helped to draft the manuscript. BMP participated in design of the study, carried out gonad dissections, mRNA extractions, constructed the SSH and helped to draft the manuscript. CCa carried out sequence and bioinformatics analyses. EPo participated in data analysis and real-time quantitative RT-PCR experiments; EPa participated in data analysis and interpretation of the results. CCo conceived the study, supervised the experiments and helped to write the manuscript. All authors read and approved the final manuscript.
Adrienne Baillet, Béatrice Mandon-Pépin contributed equally to this work.
Electronic supplementary material
Additional file 4: Evaluation of subtraction efficiency. SSH subtraction efficiency was determined by analyzing the amount of GAPDH present in both the unsubtracted starting cDNA and subtracted target cDNA through the use of increasing numbers of PCR cycles. In both libraries, the GAPDH cDNA fragment was absent of the subtracted samples. In the unsubtracted samples the amplified band was observed only following 18 PCR cycles. (PDF 25 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
- Suppression Subtractive Hybridization
- Primordial Follicle
- Meiotic Prophase
- Female Gonad
- Suppression Subtractive Hybridization Library