Gene expression profiling of differentially expressed genes in granulosa cells of bovine dominant follicles using suppression subtractive hybridization

Janus Kinase 3 phosphorylation and the JAK/STAT pathway are positively modulated by follicle-stimulating hormone (FSH) in bovine granulosa cells

Janus kinase 3 (JAK3) is a member of the JAK family of tyrosine kinase proteins involved in cytokine receptor-mediated intracellular signal transduction through the JAK/STAT signaling pathway. JAK3 was previously shown as differentially expressed in granulosa cells (GC) of bovine pre-ovulatory follicles suggesting that JAK3 could modulate GC function and activation/inhibition of downstream targets. We used JANEX-1, a JAK3 inhibitor, and FSH treatments and analyzed proliferation markers, steroidogenic enzymes and phosphorylation of target proteins including STAT3, CDKN1B/p27^Kip1 and MAPK8IP3/JIP3. Cultured GC were treated with or without FSH in the presence or not of JANEX-1. Expression of steroidogenic enzyme CYP11A1, but not CYP19A1, was upregulated in GC treated with FSH and both were significantly decreased when JAK3 was inhibited. Proliferation markers CCND2 and PCNA were reduced in JANEX-1-treated GC and upregulated by FSH. Western blots analyses showed that JANEX-1 treatment reduced pSTAT3 amounts while JAK3 overexpression increased pSTAT3. Similarly, FSH treatment increased pSTAT3 even in JANEX-1-treated GC. UHPLC-MS/MS analyses revealed phosphorylation of specific amino acid residues within JAK3 as well as CDKN1B and MAPK8IP3 suggesting possible activation or inhibition post-FSH or JANEX-1 treatments. We show that FSH activates JAK3 in GC, which could phosphorylate target proteins and likely modulate other signaling pathways involving CDKN1B and MAPK8IP3, therefore controlling GC proliferation and steroidogenic activity.

Bone morphogenetic protein 2 upregulates SERPINE2 expression through noncanonical SMAD2/3 and p38 MAPK signaling pathways in human granulosa-lutein cells

Serine protease inhibitor-E2 (SERPINE2) is highly expressed in the granulosa cells of growing follicles and the dynamic changes in SERPINE2 expression are correlated with follicular development and ovulation in several mammals, including mice, cattle, sheep, and humans. Bone morphogenetic proteins (BMPs) and their functional receptors are extensively expressed in the ovary and play critical roles in the regulation of ovarian folliculogenesis and luteal function. To date, whether BMPs regulate the expression of SERPINE2 during human follicular development remains to be elucidated. The aim of this study was to investigate the effects of BMPs on the regulation of SERPINE2 expression (a major regulator of plasminogen activators [PA]) and the underlying mechanisms using primary and immortalized human granulosa-lutein (hGL) cells. Our results demonstrated that these BMPs (BMP2, BMP4, BMP6, BMP7, and BMP15) induced differential upregulation of SERPINE2 expression. In this regard, BMP2 is the major modulator that has the best cellular activity, which further decreased the production of urokinase PA and tissue PA in hGL cells. In addition to canonical SMAD1/5/8 signaling, BMP2 also activates noncanonical SMAD2/3 and p38 mitogen-activated protein kinase (MAPK) signaling. Using two inhibition approaches (kinase receptor inhibitors and siRNA-mediated knockdown), we found that SMAD2/3-SMAD4 and p38 MAPK, but not SMAD1/5/8 signaling, was involved in the BMP2-induced upregulation of SERPINE2 expression via activin receptor-like kinase 3. These findings deepen our understanding of the differential effect of BMPs in regulating follicular function and provide new insights of the molecular mechanisms by which BMP2 regulates the expression of SERPINE2 in human granulosa cells.

Tribbles Pseudokinase 2 (TRIB2) Regulates Expression of Binding Partners in Bovine Granulosa Cells

Members of the Tribbles (TRIB) family of pseudokinases are critical components of intracellular signal transduction pathways in physiological and pathological processes. TRIBs, including TRIB2, have been previously shown as signaling mediators and scaffolding proteins regulating numerous cellular events such as proliferation, differentiation and cell death through protein stability and activity. However, the signaling network associated with TRIB2 and its binding partners in granulosa cells during ovarian follicular development is not fully defined. We previously reported that TRIB2 is differentially expressed in growing dominant follicles while downregulated in ovulatory follicles following the luteinizing hormone (LH) surge or human chorionic gonadotropin (hCG) injection. In the present study, we used the yeast two-hybrid screening system and in vitro coimmunoprecipitation assays to identify and confirm TRIB2 interactions in granulosa cells (GCs) of dominant ovarian follicles (DFs), which yielded individual candidate binding partners including calmodulin 1 (CALM1), inhibin subunit beta A (INHBA), inositol polyphosphate phosphatase-like 1 (INPPL1), 5'-nucleotidase ecto (NT5E), stearoyl-CoA desaturase (SCD), succinate dehydrogenase complex iron sulfur subunit B (SDHB) and Ras-associated protein 14 (RAB14). Further analyses showed that all TRIB2 binding partners are expressed in GCs of dominant follicles but are differentially regulated throughout the different stages of follicular development. CRISPR/Cas9-driven inhibition along with pQE-driven overexpression of TRIB2 showed that TRIB2 differently regulates expression of binding partners, which reveals the importance of TRIB2 in the control of gene expression linked to various biological processes such as proliferation, differentiation, cell migration, apoptosis, calcium signaling and metabolism. These data provide a larger view of potential TRIB2-regulated signal transduction pathways in GCs and provide strong evidence that TRIB2 may act as a regulator of target genes during ovarian follicular development.

Decreased ascorbic acid biosynthesis in response to PMSG in the pre-pubertal female rat ovary

Ascorbic acid (AA) is known to be an important antioxidant serving as a cofactor in collagen synthesis, and thus facilitates follicular growth in the ovary. Many studies have shown that AA is synthesized in the liver and transported to other organs including ovary, however, there is no direct evidence of ascorbic acid synthesis in the ovary. Hence, we examined the expression pattern of different proteins (SMP30/GNL and GULO) involved in the AA synthesis in pre-pubertal rat, which showed significant expression of these proteins, suggesting the synthesis of AA in the ovary. Accumulation of AA in the ovary during follicular growth has been well demonstrated. However, the effect of Pregnant Mare Serum Gonadotropin (PMSG) on the AA synthesis in the ovary has not been studied in detail. Hence to decipher the effect, different doses of PMSG were injected subcutaneously into the pre-pubertal female rats, and ovarian AA level was measured after 48 h. A significant increase in AA content was observed in PMSG treated animal groups. Further, to understand the mechanism underlying ovarian AA accumulation, the expression levels of SMP30/GNL and GULO genes were measured. Expression of both the genes was significantly suppressed, which suggested a lowered AA synthesis in the PMSG treated rat ovary. For further understanding, mRNA expression of AA transporters SVCT1 and SVCT2 encoded by SLC23A1 and SLC23A2 genes respectively were measured, which showed increased level of SVCT1 expression. These observations suggested that the increased AA content might not be due to increased synthesis of AA within the ovary but possibly due to increased uptake from blood during the stimulation of follicular growth.

Functional effects of Tribbles homolog 2 in bovine ovarian granulosa cells†

Tribbles homologs (TRIB) 1, 2, and 3 represent atypical members of the serine/threonine kinase superfamily. We previously identified TRIB2 as a differentially expressed gene in granulosa cells (GCs) of bovine preovulatory follicles. The current study aimed to further investigate TRIB2 regulation and study its function in the ovary. GCs were collected from follicles at different developmental stages: small antral follicles (SF), dominant follicles (DF) at day 5 of the estrous cycle, and hCG-induced ovulatory follicles (OFs). RT-qPCR analyses showed greater expression of TRIB2 in GC of DF as compared to OF and a significant downregulation of TRIB2 steady-state mRNA amounts by hCG/LH, starting at 6 h through 24 h post-hCG as compared to 0 h. Specific anti-TRIB2 polyclonal antibodies were generated and western blot analysis confirmed TRIB2 downregulation by hCG at the protein level. In vitro studies showed that FSH stimulates TRIB2 expression in GC. Inhibition of TRIB2 using CRISPR/Cas9 resulted in a significant increase in PCNA expression and an increase in steroidogenic enzyme CYP19A1 expression, while TRIB2 overexpression tended to decrease GC proliferation. TRIB2 inhibition also resulted in a decrease in transcription factors connective tissue growth factor (CTGF) and ankyrin repeat domain-containing protein 1 (ANKRD1) expression, while TRIB2 overexpression increased CTGF and ANKRD1. Additionally, western blot analyses showed reduction in ERK1/2 (MAPK3/1) and p38MAPK (MAPK14) phosphorylation levels following TRIB2 inhibition, while TRIB2 overexpression increased p-ERK1/2 and p-p38MAPK. These results provide evidence that TRIB2 modulates MAPK signaling in GC and that TRIB2 could act as a regulator of GC proliferation and function, which could affect steroidogenesis during follicular development.

The Role of MicroRNAs in Mammalian Fertility: From Gametogenesis to Embryo Implantation

The genetic codes inscribed during two key developmental processes, namely gametogenesis and embryogenesis, are believed to determine subsequent development and survival of adult life. Once the embryo is formed, its further development mainly depends on its intrinsic characteristics, maternal environment (the endometrial receptivity), and the embryo–maternal interactions established during each phase of development. These developmental processes are under strict genetic regulation that could be manifested temporally and spatially depending on the physiological and developmental status of the cell. MicroRNAs (miRNAs), one of the small non-coding classes of RNAs, approximately 19–22 nucleotides in length, are one of the candidates for post-transcriptional developmental regulators. These tiny non-coding RNAs are expressed in ovarian tissue, granulosa cells, testis, oocytes, follicular fluid, and embryos and are implicated in diverse biological processes such as cell-to-cell communication. Moreover, accumulated evidences have also highlighted that miRNAs can be released into the extracellular environment through different mechanisms facilitating intercellular communication. Therefore, understanding miRNAs mediated regulatory mechanisms during gametogenesis and embryogenesis provides further insights about the molecular mechanisms underlying oocyte/sperm formation, early embryo development, and implantation. Thus, this review highlights the role of miRNAs in mammalian gametogenesis and embryogenesis and summarizes recent findings about miRNA-mediated post-transcriptional regulatory mechanisms occurring during early mammalian development.

Non-coding RNA regulation in reproduction: Their potential use as biomarkers

Non-coding RNAs (ncRNAs) are crucial regulatory elements in most biological processes and reproduction is also controlled by them. The different types of ncRNAs, as well as the high complexity of these regulatory pathways, present a complex scenario; however, recent studies have shed some light on these questions, discovering the regulatory function of specific ncRNAs on concrete reproductive biology processes. This mini review will focus on the role of ncRNAs in spermatogenesis and oogenesis, and their potential use as biomarkers for reproductive diseases or for reproduction success.

The Orphan Nuclear Receptors Steroidogenic Factor-1 and Liver Receptor Homolog-1: Structure, Regulation, and Essential Roles in Mammalian Reproduction

Nuclear receptors are intracellular proteins that act as transcription factors. Proteins with classic nuclear receptor domain structure lacking identified signaling ligands are designated orphan nuclear receptors. Two of these, steroidogenic factor-1 (NR5A1, also known as SF-1) and liver receptor homolog-1 (NR5A2, also known as LRH-1), bind to the same DNA sequences, with different and nonoverlapping effects on targets. Endogenous regulation of both is achieved predominantly by cofactor interactions. SF-1 is expressed primarily in steroidogenic tissues, LRH-1 in tissues of endodermal origin and the gonads. Both receptors modulate cholesterol homeostasis, steroidogenesis, tissue-specific cell proliferation, and stem cell pluripotency. LRH-1 is essential for development beyond gastrulation and SF-1 for genesis of the adrenal, sexual differentiation, and Leydig cell function. Ovary-specific depletion of SF-1 disrupts follicle development, while LRH-1 depletion prevents ovulation, cumulus expansion, and luteinization. Uterine depletion of LRH-1 compromises decidualization and pregnancy. In humans, SF-1 is present in endometriotic tissue, where it regulates estrogen synthesis. SF-1 is underexpressed in ovarian cancer cells and overexpressed in Leydig cell tumors. In breast cancer cells, proliferation, migration and invasion, and chemotherapy resistance are regulated by LRH-1. In conclusion, the NR5A orphan nuclear receptors are nonredundant factors that are crucial regulators of a panoply of biological processes, across multiple reproductive tissues.

Reproductive failure in mice expressing transgenic follicle-stimulating hormone is not caused by loss of oocyte quality

Human female reproductive aging features declining ovarian follicle reserve and oocyte quality, and rising levels of circulating follicle-stimulating hormone (FSH). We determined the effects of elevated FSH on oocyte-embryo development in mature mice exhibiting premature infertility caused by progressively rising transgenic human FSH (TgFSH) levels. Oocyte-embryo developmental competence and quality were examined using oocyte maturation and aneuploidy rates, biomarkers of oocyte quality, and reciprocal embryo transfers assessed for implantation and pregnancy. In vitro maturation suggested that TgFSH exposure only hindered oocyte developmental competence in old females, as significantly more oocytes from ≥12-month-old TgFSH females remained at germinal vesicle stage compared with age-matched control oocytes. Aneuploidy rates were equivalent in oocytes from aging TgFSH compared with wildtype females. Cumulus cell expression levels of candidate biomarker Inhba, Egfr, and Rgs2 transcripts were elevated in associated aneuploid vs euploid oocytes from both TgFSH and wildtype females. In vivo, embryos transferred from subfertile 6-month-old TgFSH females to wildtype recipients yielded normal implantation rates and more pups born compared with controls. Transfer of wildtype embryos rescued the fertility of 6-month-old TgFSH-recipient females, although pup birth weight was reduced in TgFSH vs wildtype recipients. Our current findings show that elevated FSH had minimal disruption of either embryo developmental capacity or uterine function when examined in isolation, and the subfertility of TgFSH female mice was not caused by altered oocyte aneuploidy or quality.

Gene expression profiling of upregulated mRNAs in granulosa cells of bovine ovulatory follicles following stimulation with hCG

BACKGROUND

Ovulation and luteinization of follicles are complex biological processes initiated by the preovulatory luteinizing hormone surge. The objective of this study was to identify genes that are differentially expressed in bovine granulosa cells (GC) of ovulatory follicles.

METHODS

Granulosa cells were collected during the first follicular wave of the bovine estrous cycle from dominant follicles (DF) and from ovulatory follicles (OF) obtained 24 h following injection of human chorionic gonadotropin (hCG). A granulosa cell subtracted cDNA library (OF-DF) was generated using suppression subtractive hybridization and screened.

RESULTS

Detection of genes known to be upregulated in bovine GC during ovulation, such as ADAMTS1, CAV1, EGR1, MMP1, PLAT, PLA2G4A, PTGES, PTGS2, RGS2, TIMP1, TNFAIP6 and VNN2 validated the physiological model and analytical techniques used. For a subset of genes that were identified for the first time, gene expression profiles were further compared by semiquantitative RT-PCR in follicles obtained at different developmental stages. Results confirmed an induction or upregulation of the respective mRNAs in GC of OF 24 h after hCG-injection compared with those of DF for the following genes: ADAMTS9, ARAF, CAPN2, CRISPLD2, FKBP5, GFPT2, KIT, KITLG, L3MBLT3, MRO, NUDT10, NUDT11, P4HA3, POSTN, PSAP, RBP1, SAT1, SDC4, TIMP2, TNC and USP53. In bovine GC, CRISPLD2 and POSTN mRNA were found as full-length transcript whereas L3MBLT3 mRNA was alternatively spliced resulting in a truncated protein missing the carboxy-terminal end amino acids, ⁷⁷⁴KNSHNEL⁷⁸⁰. Conversely, L3MBLT3 is expressed as a full-length mRNA in a bovine endometrial cell line. The ⁷⁷⁴KNSHNEL⁷⁸⁰ sequence is well conserved in all mammalian species and follows a SAM domain known to confer protein/protein interactions, which suggest a key function for these amino acids in the epigenetic control of gene expression.

CONCLUSIONS

We conclude that we have identified novel genes that are upregulated by hCG in bovine GC of OF, thereby providing novel insight into peri-ovulatory regulation of genes that contribute to ovulation and/or luteinization processes.

Stable reference genes in granulosa cells of bovine dominant follicles during follicular growth, FSH stimulation and maternal aging

The aim of the present study was to determine a set of reference genes in granulosa cells of dominant follicles that are suitable for relative gene expression analyses during maternal and follicular aging. Granulosa cells of growing and preovulatory dominant follicles were collected from aged and young cows (maternal aging study) and from FSH-stimulated follicles developing under different durations of FSH treatment (follicular aging study). The mRNA levels of the two commonly used reference genes (GAPDH, ACTB) and four novel genes (UBE2D2, EIF2B2, SF3A1, RNF20) were analysed using cycle threshold values. Results revealed that mRNA levels of GAPDH, ACTB, EIF2B2, RNF20, SF3A1 and UBE2D2 were similar (P>0.05) between dominant follicle type, age and among follicles obtained after FSH-stimulation, but differed (P=0.005) due to mRNA processing (i.e. with versus without amplification). The stability of reference genes was analysed using GeNorm, DeltaCT and NormFinder programs and comprehensive ranking order was determined using RefFinder. The mRNA levels of GAPDH and ACTB were less stable than those of UBE2D2 and EIF2B2. The geometric mean of multiple genes (UBE2D2, EIF2B2, GAPDH and SF3A1) is a more appropriate reference control than the use of a single reference gene to compare relative gene expression among dominant and FSH-stimulated follicles during maternal and/or follicular aging studies.

MicroRNA Mediating Networks in Granulosa Cells Associated with Ovarian Follicular Development

Ovaries, which provide a place for follicular development and oocyte maturation, are important organs in female mammals. Follicular development is complicated physiological progress mediated by various regulatory factors including microRNAs (miRNAs). To demonstrate the role of miRNAs in follicular development, this study analyzed the expression patterns of miRNAs in granulosa cells through investigating three previous datasets generated by Illumina miRNA deep sequencing. Furthermore, via bioinformatic analyses, we dissected the associated functional networks of the observed significant miRNAs, in terms of interacting with signal pathways and transcription factors. During the growth and selection of dominant follicles, 15 dysregulated miRNAs and 139 associated pathways were screened out. In comparison of different styles of follicles, 7 commonly abundant miRNAs and 195 pathways, as well as 10 differentially expressed miRNAs and 117 pathways in dominant follicles in comparison with subordinate follicles, were collected. Furthermore, SMAD2 was identified as a hub factor in regulating follicular development. The regulation of miR-26a/b on smad2 messenger RNA has been further testified by real time PCR. In conclusion, we established functional networks which play critical roles in follicular development including pivotal miRNAs, pathways, and transcription factors, which contributed to the further investigation about miRNAs associated with mammalian follicular development.

Differential regulation of Janus kinase 3 (JAK3) in bovine preovulatory follicles and identification of JAK3 interacting proteins in granulosa cells

Background

Janus kinase 3 (JAK3) is a member of the membrane-associated non-receptor tyrosine kinase protein family and is considered predominantly expressed in hematopoietic cells. We previously identified JAK3 as a differentially expressed gene in granulosa cells (GC) of bovine preovulatory follicles. The present study aimed to further investigate JAK3 regulation, to identify protein binding partners and better understand its mode of action in bovine reproductive cells.

Results

GC were obtained from small follicles (SF), dominant follicles at day 5 of the estrous cycle (DF), and ovulatory follicles, 24 h following hCG injection (OF). RT-PCR analyses showed greatest expression of JAK3 in GC of DF, while JAK3 expression was downregulated in OF (P < 0.0001). In addition, there was a 5- and 20-fold reduction of JAK3 steady-state mRNA levels in follicular walls, respectively at 12 and 24 hours post-hCG as compared to 0 h (P < 0.05). Similarly, JAK3 expression was downregulated by the endogenous LH surge. These results were confirmed in western blot analysis showing weakest JAK3 protein amounts in OF as compared to DF. Yeast two-hybrid screening of a DF-cDNA library resulted in the identification of JAK3 partners in GC that were confirmed by co-immunoprecipitation and included leptin receptor overlapping transcript-like 1 (LEPROTL1), inhibin beta A (INHBA) and cyclin-dependent kinase inhibitor 1B (CDKN1B). In functional studies using bovine endometrial cells, JAK3 increased phosphorylation of STAT3 and cell viability, while the addition of JANEX-1 inhibited JAK3 actions.

Conclusion

These results support a physiologically relevant role of JAK3 in follicular development and provide insights into the mode of action and function of JAK3 in reproductive tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13048-016-0280-5) contains supplementary material, which is available to authorized users.

Transcriptome Analysis of Bovine Ovarian Follicles at Predeviation and Onset of Deviation Stages of a Follicular Wave

For two libraries (PDF1 and ODF1) using Illumina sequencing 44,082,301 and 43,708,132 clean reads were obtained, respectively. After being mapped to the bovine RefSeq database, 15,533 genes were identified to be expressed in both types of follicles (cut-off RPKM > 0.5), of which 719 were highly expressed in bovine follicles (cut-off RPKM > 100). Furthermore, 83 genes were identified as being differentially expressed in ODF1 versus PDF1, where 42 genes were upregulated and 41 genes were downregulated. KEGG pathway analysis revealed two upregulated genes in ODF1 versus PDF1, CYP11A1, and CYP19A1, which are important genes in the steroid hormone biosynthesis pathway. This study represents the first investigation of transcriptome of bovine follicles at predeviation and onset of deviation stages and provides a foundation for future investigation of the regulatory mechanisms involved in follicular development in cattle.

Expression and regulation of regulator of G-protein signaling protein-2 (RGS2) in equine and bovine follicles prior to ovulation: molecular characterization of RGS2 transactivation in bovine granulosa cells

The luteinizing hormone preovulatory surge stimulates several signal pathways essential for ovulation, and the regulator of G-protein signaling protein-2 (RGS2) is thought to be involved in this process. The objectives of this study were to characterize the regulation of RGS2 transcripts in equine and bovine follicles prior to ovulation and to determine its transcriptional control in bovine granulosa cells. To assess the regulation of equine RGS2 prior to ovulation, RT-PCR was performed using total RNA extracted from equine follicles collected at various times after human chorionic gonadotropin (hCG) injection. Results showed that RGS2 mRNA levels were very low at 0 h but markedly increased 12-39 h post-hCG (P < 0.05). In the bovine species, results revealed that RGS2 mRNA levels were low in small and dominant follicles and in ovulatory follicles obtained at 0 h, but markedly increased in ovulatory follicles 6-24 h post-hCG (P < 0.05). To study the molecular control of RGS2 expression, primary cultures of bovine granulosa cells were used. Stimulation with forskolin induced an up-regulation of RGS2 mRNA in vitro. Studies using 5'-deletion mutants identified a minimal region containing full-length basal and forskolin-inducible RGS2 promoter activities. Site-directed mutagenesis indicated that these activities were dependent on CRE and ETS1 cis-elements. Electrophoretic mobility shift assays confirmed the involvement of these elements and revealed their interactions with CREB1 and ETS1 proteins. Chromatin immunoprecipitation assays confirmed endogenous interactions of these proteins with the RGS2 promoter in granulosa cells. Forskolin-inducible RGS2 promoter activity and mRNA expression were markedly decreased by PKA and ERK1/2 inhibitors, and treatment with an antagonist of PGR (RU486) and inhibitors of PTGS2 (NS398) and EGFR (PD153035) blocked the forskolin-dependent RGS2 transcript expression, suggesting the importance of RGS2 in ovulation. Collectively, this study reports for the first time the gonadotropin-dependent up-regulation of RGS2 in equine and bovine preovulatory follicles and presents some of the regulatory controls involved in RGS2 gene expression in granulosa cells.

Expression of microRNAs and their target genes and pathways associated with ovarian follicle development in cattle

Development of ovarian follicles is controlled at the molecular level by several gene products whose precise expression leads to regression or ovulation of follicles. MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression through sequence-specific base pairing with target messenger RNAs (mRNAs) causing translation repression or mRNA degradation. The aim of this study was to identify miRNAs expressed in theca and/or granulosa layers and their putative target genes/pathways that are involved in bovine ovarian follicle development. By using miRCURY microarray (Exiqon) we identified 14 and 49 differentially expressed miRNAs (P < 0.01) between dominant and subordinate follicles in theca and granulosa cells, respectively. The expression levels of four selected miRNAs were confirmed by qRT-PCR. To identify target prediction and pathways of differentially expressed miRNAs we used Union of Genes option in DIANA miRPath v.2.0 software. The predicted targets for these miRNAs were enriched for pathways involving oocyte meiosis, Wnt, TGF-beta, ErbB, insulin, P13K-Akt, and MAPK signaling pathways. This study identified differentially expressed miRNAs in the theca and granulosa cells of dominant and subordinate follicles and implicates them in having important roles in regulating known molecular pathways that determine the fate of ovarian follicle development.

Granulosa cell function and oocyte competence: Super-follicles, super-moms and super-stimulation in cattle

The review presents an overview of studies that examined the effects of follicular aging and maternal aging in the bovine model. The first of three main sections is a discussion of the developmental competence of oocytes from (1) the ovulatory follicle of 2-wave and 3-wave estrous cycles, (2) dominant follicles that develop under high or low LH pulse frequency, and (3) natural versus FSH-stimulated ovulatory follicles. The second section highlights the effects of maternal aging. Maternal aging in cattle is associated with (1) elevated circulating FSH concentrations, (2) reduced response to superstimulatory treatment, and (3) markedly decreased early embryonic development in cows >12 year of age. The third and final section on superstimulation protocols addresses the effects of the duration of FSH stimulation and withdrawal (i.e., FSH "starvation" or "coasting") on oocyte competence. Ovarian superstimulation for 4 days altered the expression of genes related to angiogenesis, and activated oxidative stress-response genes. Extending the duration of FSH stimulation from 4 to 7 days resulted in a greater and more synchronous ovulatory response and optimal oocyte maturation. The highest rates of blastocyst development in vitro were obtained when FSH support was discontinued for 44 to 68h and granulosa cell SMAD7 mRNA was predictive of this period. Longer periods of FSH starvation resulted in a loss of oocyte competence or ovulatory capability. By extending the bovine model to the transcriptome level, new approaches and treatments may be devised to resolve subfertility in women and animals.

Signaling by the extracellular matrix protein Reelin promotes granulosa cell proliferation in the chicken follicle

Chicken oocytes develop in follicles and reach an enormous size because of a massive uptake of yolk precursors such as very low density lipoprotein and vitellogenin. Oocyte growth is supported by theca cells and granulosa cells, which establish dynamic and highly organized cell layers surrounding the oocyte. The signaling processes orchestrating the development of these layered structures are largely unknown. Here we demonstrate that the Reelin pathway, which determines the development of layered neuronal structures in the brain, is also active in chicken follicles. Reelin, which is expressed in theca cells, triggers a signal in granulosa cells via apolipoprotein E receptor 2 and the very low density lipoprotein receptor, resulting in the phosphorylation of disabled-1 and consecutive activation of the phosphatidylinositol 3-kinase/Akt pathway. This signaling pathway supports the proliferation of differentiated granulosa cells to keep up with the demand of cells to cover the rapidly increasing surface of the giant germ cell.

Transcriptome profiling of granulosa cells of bovine ovarian follicles during growth from small to large antral sizes

Background

At later stages of folliculogenesis, the mammalian ovarian follicle contains layers of epithelial granulosa cells surrounding an antral cavity. During follicle development granulosa cells replicate, secrete hormones and support the growth of the oocyte. In cattle, the follicle needs to grow > 10 mm in diameter to allow an oocyte to ovulate, following which the granulosa cells cease dividing and differentiate into the specialised cells of the corpus luteum. To better understand the molecular basis of follicular growth and granulosa cell maturation, we undertook transcriptome profiling of granulosa cells from small (< 5 mm; n = 10) and large (> 10 mm, n = 4) healthy bovine follicles using Affymetrix microarrays (24,128 probe sets).

Results

Principal component analysis for the first two components and hierarchical clustering showed clustering into two groups, small and large, with the former being more heterogeneous. Size-frequency distributions of the coefficient of variation of the signal intensities of each probe set also revealed that small follicles were more heterogeneous than the large. IPA and GO enrichment analyses revealed that processes of axonal guidance, immune signalling and cell rearrangement were most affected in large follicles. The most important networks were associated with: (A) Notch, SLIT/ROBO and PI3K signalling, and (B) ITGB5 and extracellular matrix signalling through extracellular signal related kinases (ERKs). Upstream regulator genes which were predicted to be active in large follicles included STAT and XBP1. By comparison, developmental processes such as those stimulated by KIT, IHH and MEST were most active in small follicles. MGEA5 was identified as an upstream regulator in small follicles. It encodes an enzyme that modifies the activity of many target proteins, including those involved in energy sensing, by removal of N-acetylglucosamine from serine and threonine residues.

Conclusions

Our data suggest that as follicles enlarge more genes and/or pathways are activated than are inactivated, and gene expression becomes more uniform. These findings could be interpreted that either the cells in large follicles are more uniform in their gene expression, or that follicles are more uniform or a combination of both and that additional factors, such as LH, are additionally controlling the granulosa cells.

Expression, regulation, and promoter activation of vanin-2 (VNN2) in bovine follicles prior to ovulation

Vanin-2 (VNN2) is known to be involved in inflammation and leukocyte migration, but its regulation in follicles remains unknown. The objectives of this work were to study the regulation of VNN2 transcripts in bovine follicles prior to ovulation and to characterize the control of its expression in bovine granulosa cells. VNN2 expression was studied using total RNA extracted from granulosa cells of small follicles (2-4 mm in diameter), dominant follicles obtained on Day 5 of the estrous cycle, ovulatory follicles obtained 0-24 h after human chorionic gonadotropin (hCG), and corpora lutea on Day 5 of the cycle. The results from RT-PCR analyses showed that levels of VNN2 mRNA were high in ovulatory follicles 24 h post-hCG but low in the other tissues. In ovulatory follicles, levels of VNN2 mRNA were low at 0 h but significantly up-regulated 12-24 h post-hCG. To determine factors controlling VNN2 gene expression, established primary cultures of granulosa cells isolated from bovine dominant follicles were used. Treatment with forskolin elevated VNN2 mRNA expression as observed in vivo. Mutation studies identified the minimal region conferring basal and forskolin-stimulated VNN2 promoter activities, which were dependent on chicken ovalbumin upstream promoter-transcription factor (COUP-TF), GATA, and Ebox cis-elements. Electrophoretic mobility shift assays identified COUP-TF, GATA4, and upstream stimulating factor proteins as key factors interacting with these elements. Chromatin immunoprecipitation assays confirmed basal and forskolin-induced interactions between these proteins and the VNN2 promoter in bovine granulosa cell cultures. VNN2 promoter activity and mRNA expression were markedly stimulated by forskolin and overexpression of the catalytic subunit of PKA, but inhibited by PKA and ERK1/2 inhibitors. Collectively, the findings from this study describe for the first time the gonadotropin/forskolin-dependent up-regulation of VNN2 transcripts in granulosa cells of preovulatory follicles and provide insights into some of the molecular bases of VNN2 gene expression in follicular cells.

FSH in vitro versus LH in vivo: similar genomic effects on the cumulus

The use of gonadotropins to trigger oocyte maturation both in vivo and in vitro has provided precious and powerful knowledge that has significantly increased our understanding of the ovarian function. Moreover, the efficacy of most assisted reproductive technologies (ART) used in both humans and livestock species relies on gonadotropin input, mainly FSH and LH. Despite the significant progress achieved and the huge impact of gonadotropins, the exact molecular pathways of the two pituitary hormones, FSH and LH, still remain poorly understood. Moreover, these pathways may not be the same when moving from the in vivo to the in vitro context. This misunderstanding of the intricate synergy between these two hormones leads to a lack of consensus about their use mainly in vitro or in ovulation induction schedules in vivo. In order to optimize their use, additional work is thus required with a special focus on comparing the in vitro versus the in vivo effects. In this context, this overview will briefly summarize the downstream gene expression pathways induced by both FSH in vitro and LH in vivo in the cumulus compartment. Based on recent microarray comparative analysis, we are reporting that in vitro FSH stimulation on cumulus cells appears to achieve at least part of the gene expression activity after in vivo LH stimulation. We are then proposing that the in vitro FSH-response of cumulus cells have similitudes with the in vivo LH-response.

Molecular characterization and transcriptional regulation of a disintegrin and metalloproteinase with thrombospondin motif 1 (ADAMTS1) in bovine preovulatory follicles

The ovulatory process involves a complex remodeling of the extracellular matrix during which a desintegrin and metalloproteinase with thrombospondin motif 1 (ADAMTS1) is thought to play a key role, but its transcriptional regulation in bovine follicles remains largely unknown. The objectives of this study were to characterize the regulation of ADAMTS1 in bovine follicles before ovulation and to determine its transcriptional control in bovine granulosa cells. Regulation of ADAMTS1 was assessed using total RNA isolated from bovine preovulatory follicles obtained at various times after human chorionic gonadotropin treatment. Results from RT-PCR analyses showed that levels of ADAMTS1 mRNA were very low at 0 hours but increased at 6 to 24 hours after human chorionic gonadotropin in granulosa cells. To determine the regulatory mechanisms controlling ADAMTS1 gene expression in vitro, primary cultures of bovine granulosa cells were established, and treatment with forskolin up-regulated ADAMTS1 mRNA levels. Promoter activity assays, 5'-deletion, and site-directed mutagenesis identified a minimal region conferring full-length basal and forskolin-stimulated ADAMTS1 promoter activities, with both being dependent on Ebox cis-acting elements. EMSAs revealed upstream stimulating factor (USF) proteins as key trans-activating factors interacting with Ebox. Chromatin immunoprecipitation assays confirmed such interactions between USF and Ebox in vivo, and USF binding to Ebox elements was increased by forskolin treatment. ADAMTS1 promoter activity and mRNA expression were increased by forskolin and overexpression of the catalytic subunit of protein kinase A, but not by cotreatment with inhibitors of protein kinase A, ERK1/2, and epidermal growth factor receptor signaling pathways. Furthermore, treatment with a soluble epidermal growth factor induced ADAMTS1 mRNA expression in granulosa cells. Collectively, results from this study describe the gonadotropin/forskolin-dependent up-regulation of ADAMTS1 mRNA in granulosa cells of bovine preovulatory follicles in vivo and in vitro and identify for the first time some of the molecular mechanisms responsible for ADAMTS1 promoter activation in follicular cells of a large monoovulatory species.

Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle

Background

Identification of single nucleotide polymorphisms (SNPs) for specific genes involved in reproduction might improve reliability of genomic estimates for these low-heritability traits. Semen from 550 Holstein bulls of high (≥ 1.7; n = 288) or low (≤ −2; n = 262) daughter pregnancy rate (DPR) was genotyped for 434 candidate SNPs using the Sequenom MassARRAY® system. Three types of SNPs were evaluated: SNPs previously reported to be associated with reproductive traits or physically close to genetic markers for reproduction, SNPs in genes that are well known to be involved in reproductive processes, and SNPs in genes that are differentially expressed between physiological conditions in a variety of tissues associated in reproductive function. Eleven reproduction and production traits were analyzed.

Results

A total of 40 SNPs were associated (P < 0.05) with DPR. Among these were genes involved in the endocrine system, cell signaling, immune function and inhibition of apoptosis. A total of 10 genes were regulated by estradiol. In addition, 22 SNPs were associated with heifer conception rate, 33 with cow conception rate, 36 with productive life, 34 with net merit, 23 with milk yield, 19 with fat yield, 13 with fat percent, 19 with protein yield, 22 with protein percent, and 13 with somatic cell score. The allele substitution effect for SNPs associated with heifer conception rate, cow conception rate, productive life and net merit were in the same direction as for DPR. Allele substitution effects for several SNPs associated with production traits were in the opposite direction as DPR. Nonetheless, there were 29 SNPs associated with DPR that were not negatively associated with production traits.

Conclusion

SNPs in a total of 40 genes associated with DPR were identified as well as SNPs for other traits. It might be feasible to include these SNPs into genomic tests of reproduction and other traits. The genes associated with DPR are likely to be important for understanding the physiology of reproduction. Given the large number of SNPs associated with DPR that were not negatively associated with production traits, it should be possible to select for DPR without compromising production.

Identification of ovarian genes regulated by follicle-stimulating hormone (Fsh) in vitro during early secondary oocyte growth in coho salmon

Follicle-stimulating hormone (Fsh) function in fishes is poorly understood. This study aimed to reveal Fsh-regulated genes in coho salmon previtellogenic ovarian follicles in vitro. Four suppression subtractive hybridization libraries were generated with RNA isolated from Fsh-treated and control follicles or follicle cell-enriched tissue fractions. Fsh induced steroidogenesis and dynamically upregulated several genes predominantly expressed in follicle cells, including WAP domain-containing protease, connexin 34.3, clusterin (clu1, clu2), fibronectin, wilms tumor 2-like, and influenza virus NS1A-binding protein a. Genes downregulated by Fsh included connective tissue growth factor, alcohol dehydrogenase 8-like, and serine/threonine-protein kinase pim-1. This study demonstrates for the first time in fishes that Fsh influences the expression of a unique suite of ovarian genes involved in processes like cell communication, survival and differentiation, and extracellular matrix remodeling. Collectively, these findings suggest that Fsh and/or steroids induce differentiation of granulosa cells and remodeling of the follicle in preparation for onset of vitellogenesis.

Profiling of differentially expressed genes using suppression subtractive hybridization in an equine model of chronic asthma

Background

Gene expression analyses are used to investigate signaling pathways involved in diseases. In asthma, they have been primarily derived from the analysis of bronchial biopsies harvested from mild to moderate asthmatic subjects and controls. Due to ethical considerations, there is currently limited information on the transcriptome profile of the peripheral lung tissues in asthma.

Objective

To identify genes contributing to chronic inflammation and remodeling in the peripheral lung tissue of horses with heaves, a naturally occurring asthma-like condition.

Methods

Eleven adult horses (6 heaves-affected and 5 controls) were studied while horses with heaves were in clinical remission (Pasture), and during disease exacerbation induced by a 30-day natural antigen challenge during stabling (Challenge). Large peripheral lung biopsies were obtained by thoracoscopy at both time points. Using suppression subtractive hybridization (SSH), lung cDNAs of controls (Pasture and Challenge) and asymptomatic heaves-affected horses (Pasture) were subtracted from cDNAs of horses with heaves in clinical exacerbation (Challenge). The differential expression of selected genes of interest was confirmed using quantitative PCR assay.

Results

Horses with heaves, but not controls, developed airway obstruction when challenged. Nine hundred and fifty cDNA clones isolated from the subtracted library were screened by dot blot array and 224 of those showing the most marked expression differences were sequenced. The gene expression pattern was confirmed by quantitative PCR in 15 of 22 selected genes. Novel genes and genes with an already defined function in asthma were identified in the subtracted cDNA library. Genes of particular interest associated with asthmatic airway inflammation and remodeling included those related to PPP3CB/NFAT, RhoA, and LTB4/GPR44 signaling pathways.

Conclusions

Pathways representing new possible targets for anti-inflammatory and anti-remodeling therapies for asthma were identified. The findings of genes previously associated with asthma validate this equine model for gene expression studies.

Gene expression profiling of preovulatory follicle in the buffalo cow: effects of increased IGF-I concentration on periovulatory events

The preovulatory follicle in response to gonadotropin surge undergoes dramatic biochemical, and morphological changes orchestrated by expression changes in hundreds of genes. Employing well characterized bovine preovulatory follicle model, granulosa cells (GCs) and follicle wall were collected from the preovulatory follicle before, 1, 10 and 22 h post peak LH surge. Microarray analysis performed on GCs revealed that 450 and 111 genes were differentially expressed at 1 and 22 h post peak LH surge, respectively. For validation, qPCR and immunocytochemistry analyses were carried out for some of the differentially expressed genes. Expression analysis of many of these genes showed distinct expression patterns in GCs and the follicle wall. To study molecular functions and genetic networks, microarray data was analyzed using Ingenuity Pathway Analysis which revealed majority of the differentially expressed genes to cluster within processes like steroidogenesis, cell survival and cell differentiation. In the ovarian follicle, IGF-I is established to be an important regulator of the above mentioned molecular functions. Thus, further experiments were conducted to verify the effects of increased intrafollicular IGF-I levels on the expression of genes associated with the above mentioned processes. For this purpose, buffalo cows were administered with exogenous bGH to transiently increase circulating and intrafollicular concentrations of IGF-I. The results indicated that increased intrafollicular concentrations of IGF-I caused changes in expression of genes associated with steroidogenesis (StAR, SRF) and apoptosis (BCL-2, FKHR, PAWR). These results taken together suggest that onset of gonadotropin surge triggers activation of various biological pathways and that the effects of growth factors and peptides on gonadotropin actions could be examined during preovulatory follicle development.

Transitional gene expression profiling in ovarian follicle during ovulation in normal-cycle rats

Evaluation of ovarian toxicity requires an understanding of the physiological changes related to the estrous cycle in the ovary. The authors investigated the transitional gene expression profile of ovulatory follicles in rats that show normal estrous cyclicity. Ovaries were collected at 10:00 and 22:00 on the proestrus day and at 10:00 on the estrus day. Ovarian follicles or early corpora lutea were isolated using laser microdissection, and extracted total RNA was analyzed using microarray technology. Clustering analysis revealed four different expression patterns: transient up- or down-regulation only at 22:00 on the proestrus day (pattern 1), up- or down-regulation only at 10:00 on the estrus day (pattern 2), continuous increase at 22:00 on the proestrus day and at 10:00 on the estrus day (pattern 3), and up- or down-regulation at 22:00 on the proestrus day and level maintenance at 10:00 on the estrus day (pattern 4). In addition, these probe sets were functionally categorized in each pattern using the Ingenuity Pathways Analysis database. These data will aid in understanding the physiology of ovulation and may be useful in assessing ovarian toxicity and its mechanism, such as in investigations of chemical-induced ovulatory impairment.

SERPINE2, an inhibitor of plasminogen activators, is highly expressed in the human endometrium during the secretory phase

BACKGROUND

SERPINE2, also known as protease nexin-1, belongs to the serine protease inhibitor (SERPIN) superfamily. It is one of the potent SERPINs that modulates the activity of plasminogen activators (PAs). PAs and their SERPIN inhibitors, such as SERPINB2 and SERPINE1, were expressed in the human endometrium and were implicated in implantation. However, expression data about SERPINE2 in the human endometrium is still unknown. Thus, we conducted an investigation to reveal the spatiotemporal and cellular expression of SERPINE2 in the human uterus during the menstrual cycle.

METHODS

Seven patients who underwent a hysterectomy and samples of 120 archived patients' endometrial curettage or parts of the uterus that were formalin-fixed and embedded in paraffin. Western blotting was performed to evaluate the specificity and sensitivity of the antibody. Immunohistochemistry was conducted to localize the SERPINE2 expression site. Quantitative analysis was conducted to evaluate expression levels of SERPINE2 in various sub-phases of the menstrual cycle.

RESULTS

The SERPINE2 protein was primarily detected in the uterine fluid during the mid- and late-secretory phases of the menstrual cycle. It was predominantly expressed in the luminal and glandular epithelium, less in the myometrium, and only dispersedly in certain stromal cells throughout the menstrual cycle. A quantitative analysis of expression levels of SERPINE2 in the glandular epithelium revealed that it was highly expressed in the endometrium during the secretory phase compared to the proliferative phase.

CONCLUSIONS

The SERPINE2 protein is highly expressed in the endometrium during the secretory phase, indicating that it may participate in tissue remodeling involved in implantation.

Spatiotemporal expression of the serine protease inhibitor, SERPINE2, in the mouse placenta and uterus during the estrous cycle, pregnancy, and lactation

BACKGROUND

SERPINE2, also known as glia-derived nexin or protease nexin-1, belongs to the serine protease inhibitor (SERPIN) superfamily. It is one of the potent serpins that modulates the activity of the plasminogen activator (PA) and was implicated in tissue remodeling. In this study, we investigated the expression patterns of SERPINE2 in the mouse placenta and uterus during the estrous cycle, pregnancy, and lactation.

METHODS

SERPINE2 was purified from mouse seminal vesicle secretion using liquid chromatography (LC) and identified by LC/tandem mass spectrometry. The antiserum against the SERPINE2 protein was raised in rabbits. To reveal the uterine and placental expression of SERPINE2, tissues at various stages were collected for real-time PCR quantification, Western blotting, and immunohistochemical staining.

RESULTS

Serpine2 mRNA was the major PA inhibitor in the placenta and uterus during the estrous cycle, pregnancy, and lactation, although Serpine1 mRNA had higher expression levels than Serpine2 mRNA in the placenta. Plat seemed to be the major PA in the mouse uterus and placenta. Antiserum against the SERPINE2 protein specifically recognized two forms of SERPINE2 and an extra 75-kDa protein, which was probably a complex of SERPINE2 with a certain protease, from among thousands of protein components in the tissue extract as demonstrated by Western blotting. In the uterus, SERPINE2 was primarily localized in luminal and glandular epithelial cells but it also was detected in circular and longitudinal smooth muscle cells during the estrous cycle and lactation. It was prominently expressed in decidual stroma cells, the metrial gland, and endometrial epithelium of the pregnant uterus. In the placenta, SERPINE2 was expressed in trophoblasts of the labyrinth and spongiotrophoblasts. However, its expression was remarkably reduced in giant cells which existed in the giant cell-decidual junction zone. In contrast, prominent expression of SERPINE2 seemed to be detected on clusters of glycogen cells near the junction zone. In addition, yolk sac membranes also showed high expression of SERPINE2.

CONCLUSIONS

These findings indicate that SERPINE2 is a major PA inhibitor in the placenta and uterus during the estrous cycle, pregnancy, and lactation. It may participate in the PA-modulated tissue remodeling process in the mouse placenta and uterus.

Cumulus cell gene expression following the LH surge in bovine preovulatory follicles: potential early markers of oocyte competence

Cumulus cells (CCs) are essential for oocytes to reach full development competency and become fertilized. Many major functional properties of CCs are triggered by gonadotropins and governed by the oocyte. Consequently, cumulus may reflect oocyte quality and is often used for oocyte selection. The most visible function of CCs is their ability for rapid extracellular matrix expansion after the LH surge. Although unexplained, LH induces the final maturation and improves oocyte quality. To study the LH signaling and gene expression cascade patterns close to the germinal vesicle breakdown, bovine CCs collected at 2 h before and 6 h after the LH surge were hybridized to a custom-made microarray to better understand the LH genomic action and find differentially expressed genes associated with the LH-induced oocyte final maturation. Functional genomic analysis of the 141 overexpressed and 161 underexpressed clones was performed according to their molecular functions, gene networks, and cell compartments. Following real-time PCR validation of our gene lists, some interesting pathways associated with the LH genomic action on CCs and their possible roles in oocyte final maturation, ovulation, and fertilization are discussed. A list of early potential markers of oocyte competency in vivo and in vitro is thereafter suggested. These early biomarkers are a preamble to understand the LH molecular pathways that trigger the final oocyte competence acquisition process in bovine.

Differential genome-wide gene expression profiling of bovine largest and second-largest follicles: identification of genes associated with growth of dominant follicles

BACKGROUND

Bovine follicular development is regulated by numerous molecular mechanisms and biological pathways. In this study, we tried to identify differentially expressed genes between largest (F1) and second-largest follicles (F2), and classify them by global gene expression profiling using a combination of microarray and quantitative real-time PCR (QPCR) analysis. The follicular status of F1 and F2 were further evaluated in terms of healthy and atretic conditions by investigating mRNA localization of identified genes.

METHODS

Global gene expression profiles of F1 (10.7 +/- 0.7 mm) and F2 (7.8 +/- 0.2 mm) were analyzed by hierarchical cluster analysis and expression profiles of 16 representative genes were confirmed by QPCR analysis. In addition, localization of six identified transcripts was investigated in healthy and atretic follicles using in situ hybridization. The healthy or atretic condition of examined follicles was classified by progesterone and estradiol concentrations in follicular fluid.

RESULTS

Hierarchical cluster analysis of microarray data classified the follicles into two clusters. Cluster A was composed of only F2 and was characterized by high expression of 31 genes including IGFBP5, whereas cluster B contained only F1 and predominantly expressed 45 genes including CYP19 and FSHR. QPCR analysis confirmed AMH, CYP19, FSHR, GPX3, PlGF, PLA2G1B, SCD and TRB2 were greater in F1 than F2, while CCL2, GADD45A, IGFBP5, PLAUR, SELP, SPP1, TIMP1 and TSP2 were greater in F2 than in F1. In situ hybridization showed that AMH and CYP19 were detected in granulosa cells (GC) of healthy as well as atretic follicles. PlGF was localized in GC and in the theca layer (TL) of healthy follicles. IGFBP5 was detected in both GC and TL of atretic follicles. GADD45A and TSP2 were localized in both GC and TL of atretic follicles, whereas healthy follicles expressed them only in GC.

CONCLUSION

We demonstrated that global gene expression profiling of F1 and F2 clearly reflected a difference in their follicular status. Expression of stage-specific genes in follicles may be closely associated with their growth or atresia. Several genes identified in this study will provide intriguing candidates for the determination of follicular growth.

Molecular mechanisms regulating bovine ovarian follicular selection

Transcription profiling of ovarian follicles. Understanding the mechanisms by which a single follicle is selected for further ovulation is important to control fertility in mammals. However, development of new treatments is limited by our poor understanding of molecular mechanisms regulating follicular selection. Our hypothesis is that genes involved in the control of cell proliferation and apoptosis are differentially regulated during follicular selection. Our objective was to identify these new genes. Bovine follicles were collected and gene expression levels were measured using microarrays. First, follicles were allocated to three groups, according to the time spent from the initiation of follicular wave to surgery (24 H, 36 H, and 48-60 H). Fifty-seven genes are differentially expressed at a false discovery rate of 5%. These genes are involved in the control of lipid metabolism (P-value = 0.0005), cell proliferation (0.007), cell death (0.003), cell morphology (0.003), and immune response (0.003). Follicles were also grouped into four categories, according to the expected time of deviation (early deviation; 8 mm, mid-deviation; 8.5 mm, late deviation; 9 mm, dominant follicles; >or=10 mm). One hundred and twenty eight genes are differentially expressed between these four groups, including genes involved in cell proliferation (0.00002), cell death (0.0006), cell-to-cell signaling (0.003), cell morphology (0.003), lipid metabolism (0.0004), and immune response (0.00007). The expression levels of 10 genes were confirmed using quantitative real time PCR. As expected, we identified new differentially regulated genes involved in the control of cell growth and apoptosis. We also discovered a potential role for immune cells, and in particular macrophages, in follicular selection.

A pilot study on transcriptome data analysis of folliculogenesis in pigs

Three different stages of pig antral follicles have been studied in a granulosa-cell transcriptome analysis on nylon microarrays (1152 clones). The data have been generated from seven RNA follicle pools and several technical replicates were made. The objective of this paper was to state the feasibility of a transcriptomic protocol for the study of folliculogenesis in the pig. A statistical analysis was chosen, relying on the linear mixed model (LMM) paradigm. Low variability within technical replicates was hence checked with a LMM. Relevant genes that might be involved in the studied process were then selected. For the most significant genes, statistical methods such as principal component analysis and unsupervised hierarchical clustering were applied to assess their relevance, and a random forest analysis proved their predictive value. The selection of genes was consistent with previous studies and also allowed the identification of new genes whose role in pig folliculogenesis will be further investigated.

Follicle development in mares

The mare provides a unique experimental model for studying follicle development in monovular species. Development of antral follicles in horses is characterized by the periodic growth of follicular waves which often involve the selection of a single dominant follicle. If properly stimulated, the dominant follicle will complete development and eventually ovulate a fertile oocyte. Regulation of follicular wave emergence and follicle selection involves an interplay between circulating gonadotropins and follicular factors that ensures that individual follicles are properly stimulated to grow (or to regress) at any given stage of follicular wave development. Periodic development of follicular waves continuously occurs during most of post-natal life in the mare and is influenced by factors such as stage of oestrous cycle, season, pregnancy, age, breed and individual so that different types of follicular waves (minor or major, ovulatory or anovulatory) and different levels of activity within waves may develop under different physiological conditions. Changes in gonadotropin levels and/or in the sensitivity of follicles to circulating gonadotropins seem to account largely for these physiological variations in follicle development.

Analysis of gene expression in granulosa cells of ovine antral growing follicles using suppressive subtractive hybridization

Follicular growth, development and ovulation are highly ordered processes that involve the expression of many genes under precise temporal and spatial regulation. However, information on stage-specific gene expression during the antral follicle phase in sheep is not well understood. In the present study, suppressive subtractive hybridization (SSH) was performed to screen genes that were differentially expressed in the granulosa cells between large follicles (LF, >5mm) and small follicles (SF, 3-5mm), and subtractive cDNA library was constructed. Furthermore, with dot-blot analysis, a total of 90 clones randomly selected from the library were proven to be differentially expressed in the granulosa cells. Among these, 38 exhibited high homology to known genes, 14 sequences were corresponding to novel expressed sequence tags (ESTs). Four ESTs, LAPTM4A, SERPINE2, GSTA1, and INHBA, were further examined the reproducibility of the SSH data by the real-time quantitative PCR. Results confirmed an increase expression of respective mRNA in granulosa cells of large follicles compared with that of small follicles. It is concluded that we have identified several genes (known or unknown) that may effect follicular growth, dominance or ovulation in ewes.

Identification of transcripts involved in meiosis and follicle formation during ovine ovary development

Background

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.

Results

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.

Conclusion

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.

Identification of potential markers of oocyte competence expressed in bovine cumulus cells matured with follicle-stimulating hormone and/or phorbol myristate acetate in vitro

Oocyte competence is the ability of the oocyte to complete maturation, undergo successful fertilization, and reach the blastocyst stage. Cumulus cells are indispensable for this process. Their removal significantly affects the blastocyst rates. Moreover, the properties and functions of cumulus cells are regulated by the oocyte. They also reflect the oocyte's degree of maturation. Our study was aimed at identifying markers of oocyte competence that are expressed in bovine cumulus cells. In a previous study in our laboratory, the blastocyst yield following FSH or phorbol myristate acetate (PMA) treatment was 45%%. Therefore, we tested four sets of conditions during the first 6 h of in vitro maturation (IVM): FSH (0.1 microg/ml), PMA (0.1 microM), FSH ++ PMA, and negative control. Extracts from each IVM treatment were hybridized against the same negative control on a microarray containing a partial library of differentially expressed transcripts in the cumulus of competent oocytes collected at 6 h after LH in vivo. Common positive clones between diffrentially treated cells were selected, and 15 candidates were validated by real-time PCR. Based on this, the main candidates expressed in cumulus cells and that could be valuable and indirect markers of oocyte competence are hyaluronan synthase 2 (HAS2), inhibin betaA (INHBA), epidermal growth factor receptor (EGFR), gremlin 1 (GREM1), betacellulin (BTC), CD44, tumor necrosis factor-induced protein 6 (TNFAIP6), and prostaglandin-endoperoxide synthase 2 (PTGS2). These biomarkers could be potential candidates to predict oocyte competence and to select higher-quality embryos for transfer. Additionally, these indirect predictors of oocyte competence and follicular health could improve our knowledge of gene expression patterns in the cumulus and yield insights into the molecular pathways controlling oocyte competence.

Differentiation of the bovine dominant follicle from the cohort upregulates mRNA expression for new tissue development genes

This study was designed to identify genes that regulate the transition from FSH- to LH-dependent development in the bovine dominant follicle (DF). Serial analysis of gene expression (SAGE) was used to compare the transcriptome of granulosa cells isolated from the most oestrogenic growing cohort follicle (COH), the newly selected DF and its largest subordinate follicle (SF) which is destined for atresia. Follicle diameter, follicular fluid oestradiol (E) and E:progesterone ratio confirmed follicle identity. Results show that there are 93 transcript species differentially expressed in DF granulosa cells, but only 8 of these encode proteins known to be involved in DF development. Most characterised transcripts upregulated in the DF are from tissue development genes that regulate cell differentiation, proliferation, apoptosis, signalling and tissue remodelling. Semiquantitative real-time PCR analysis confirmed seven genes with upregulated (P≤0.05) mRNA expression in DF compared with both COH and SF granulosa cells. Thus, the new genes identified by SAGE and real-time PCR, which show enhanced mRNA expression in the DF, may regulate proliferation (cyclin D2;CCND2), prevention of apoptosis or DNA damage (growth arrest and DNA damage-inducible, β;GADD45B), RNA synthesis (splicing factor, arginine/serine rich 9;SFRS9) and unknown processes associated with enhanced steroidogenesis (ovary-specific acidic protein; DQ004742) in granulosa cells of DF at the onset of LH-dependent development. Further studies are required to show whether the expression of identified genes is dysregulated when abnormalities occur during DF selection or subsequent development.

Isolation of genes associated with developmental competency of bovine oocytes

Eggs differ widely in their ability to develop into an embryo. To address this characteristic, the concept of developmental competency has been coined, defined as the ability or potential of an oocyte to undergo maturation, fertilization and development to blastocyst stages or live offspring. Developmental competency is acquired progressively during folliculogenesis and is linked to follicular size. In an effort to understand the molecular changes underlying differences in competency we compared oocytes derived from large follicles (>or=5mm) to those from small follicles (<or=2mm). We used an approach combining suppressive subtraction hybridization with a linear amplification step to identify genes upregulated in the more competent oocytes. Real-time RT-PCR quantification indicated highly significant upregulation for 10 genes. However, the observed changes did not exceed three-fold suggesting that the molecular causes for poor developmental capacity may be reliant on many small changes. In monovulatory species oocyte developmental competency is further modulated in a process termed follicular dominancy, whereby only one of a cohort of developing ovarian follicles continues to grow. In our second approach, we aimed to identify genes that may be involved in the choice of one follicle as becoming dominant and thus restricting the developmental competency to a single oocyte. This approach, focusing on granulosa cells, yielded a small set of five genes that could be verified to be reliable markers for dominant follicles. We have further analyzed one of these involving the activin/inhibin pathway. Lastly, in a third approach we are investigating the feasibility of using nuclear transfer (NT) to interrogate oocyte developmental competency.

Combining mouse mammary gland gene expression and comparative mapping for the identification of candidate genes for QTL of milk production traits in cattle

Background

Many studies have found segregating quantitative trait loci (QTL) for milk production traits in different dairy cattle populations. However, even for relatively large effects with a saturated marker map the confidence interval for QTL location by linkage analysis spans tens of map units, or hundreds of genes. Combining mapping and arraying has been suggested as an approach to identify candidate genes. Thus, gene expression analysis in the mammary gland of genes positioned in the confidence interval of the QTL can bridge the gap between fine mapping and quantitative trait nucleotide (QTN) determination.

Results

We hybridized Affymetrix microarray (MG-U74v2), containing 12,488 murine probes, with RNA derived from mammary gland of virgin, pregnant, lactating and involuting C57BL/6J mice in a total of nine biological replicates. We combined microarray data from two additional studies that used the same design in mice with a total of 75 biological replicates. The same filtering and normalization was applied to each microarray data using GeneSpring software. Analysis of variance identified 249 differentially expressed probe sets common to the three experiments along the four developmental stages of puberty, pregnancy, lactation and involution. 212 genes were assigned to their bovine map positions through comparative mapping, and thus form a list of candidate genes for previously identified QTLs for milk production traits. A total of 82 of the genes showed mammary gland-specific expression with at least 3-fold expression over the median representing all tissues tested in GeneAtlas.

Conclusion

This work presents a web tool for candidate genes for QTL (cgQTL) that allows navigation between the map of bovine milk production QTL, potential candidate genes and their level of expression in mammary gland arrays and in GeneAtlas. Three out of four confirmed genes that affect QTL in livestock (ABCG2, DGAT1, GDF8, IGF2) were over expressed in the target organ. Thus, cgQTL can be used to determine priority of candidate genes for QTN analysis based on differential expression in the target organ.

Follicle selection in cattle and horses: role of intrafollicular factors

The eminent event in follicle selection during a follicular wave in monovular species is diameter deviation, wherein one follicle continues to grow (developing dominant) and other follicles (subordinates) begin to regress. In cattle, the IGF system, oestradiol and LH receptors are involved in the intrafollicular events initiating deviation as indicated by the following: (1) concentrations of free IGF-I and oestradiol in the follicular fluid and number of LH receptors in the follicular wall increase more dramatically in the future dominant follicle than in the future subordinate follicles before the beginning of deviation and (2) ablation of the largest follicle (LF) or injection of recombinant human IGF (rhIGF)-I into the second LF at the expected beginning of deviation increases the concentrations of oestradiol in second LF before the expected beginning of deviation between second LF and third LF. In horses, an increase in free IGF-I, oestradiol, inhibin-A and activin-A is greater in the future dominant follicle than in other follicles before the beginning of deviation. However, free IGF-I is the only one of these four factors needed for the initiation of deviation in horses as indicated by the following: (1) ablation of LF at the expected beginning of deviation increases the concentrations of free IGF-I in second LF before the beginning of deviation between second LF and third LF but does not increase the other factors; (2) injection of rhIGF-I into second LF at the expected beginning of deviation causes second LF to continue to grow and become a codominant follicle and (3) injection of IGF-binding protein-3 into LF at the expected beginning of deviation causes LF to regress and second LF to become dominant. Thus, the dramatic changes in the IGF system in LF compared to other follicles before the beginning of deviation play a crucial role in the events that lead to the beginning of diameter deviation in both cattle and horses. Oestradiol and LH receptors also play a role in cattle. These intrafollicular events prepare the selected follicle for the decreasing availability of FSH and increasing availability of LH. The other follicles of the wave have the same future capability but do not have adequate time to attain a similar preparatory stage.

Steroidogenic changes and steady state amount of messenger RNA encoding steroidogenic enzymes, gonadotropin receptors and cell-death signalling in the dominant ovarian follicle during estradiol-induced atresia in cattle

Changes in steroidogenic function and associated gene expression were characterized in dominant ovarian follicles (DF) of cattle where follicles were induced to become atretic by systemic administration of estradiol benzoate (EB). In experiment 1, follicular fluid (FF) steroid concentrations in the DF were measured at 12-hourly time points for 48 h in heifers treated with 1 mg EB i.m./500 kg body weight (EB; n=20) as compared with untreated controls (C; n=19). Treatment with EB promoted a transient reduction in circulating FSH, a rapid (12 h) and sustained reduction in FF estradiol, a rapid (12 h) but transient reduction in FF progesterone and a delayed (36 h) increase in FF testosterone concentrations. In experiment 2, whole follicular wall tissue was collected from DF of mature non-lactating cows allocated to a 0 h control group (0 HC: n=7), a 24h control group (24 HC; n=7) or an EB-treated group where tissue was collected 24 h after administration of 1 mg EB i.m./500 kg body weight (EB; n=8). As for experiment 1, EB promoted a transient reduction in circulating FSH, a pronounced reduction in FF estradiol and a smaller but significant reduction in FF progesterone concentrations. Semi-quantitative RT-PCR on follicular wall tissue revealed that the loss in estrogen activity at 24 h after EB was associated with two-fold reduction in aromatase mRNA, with an apparent acceleration in loss of 17alpha-hydroxylase mRNA. Expression of genes for gonadotropin receptors (LHR and FSHR) and a cell-death signalling pathway (Fas antigen and Fas ligand) were unchanged during the initial 24h of EB-induced atresia. These results suggest that EB initiates atresia in dominant ovarian follicles through a rapid suppression of follicular estradiol synthesis, an effect associated with down-regulation of the aromatase gene. A transient suppression in circulating FSH following administration of EB appears to have initiated these events, and it is suggested that subsequent processes involved in atresia follow this loss in estrogenic function.