Identification of downregulated messenger RNAs in bovine granulosa cells of dominant follicles following stimulation with human chorionic gonadotropin

The role of tribbles homolog 2 in cell proliferation

Tribbles homolog 2 (TRIB2), a pseudoserine/threonine kinase, is a member of the TRIB family. TRIB2 primarily regulates cell proliferation through its scaffold or adaptor effect on promoting the degradation of target proteins by E3 ligase-dependent ubiquitination and regulating mitogen-activated protein kinase (MAPK) and protein kinase B (AKT) signaling pathways. TRIB2 is not only involved in the physiological proliferation of cells (granulosa cells, myoblasts, naive T cells, and thymocytes) during normal development but also in the pathological proliferation of vascular smooth muscle cells and a variety of cancer cells (lung cancer cells, liver cancer cells, leukemia cells, pancreatic cancer cells, gastric cancer cells, prostate cancer cells, thyroid cancer cells, cervical cancer cells, melanoma cells, colorectal cancer cells, ovarian cancer cells and osteosarcoma cells) under disease conditions. Its expression level and functional role predominantly hinge on the specific tissue and cell type it targets. This review elucidates the specific mechanisms of TRIB2 in physiological and pathological cell proliferation from the perspective of different kinds of cells.

Genetic characterization and selection of litter size traits of Guizhou Black goat and Meigu goat

The aim of this study is to explore the genetic characteristics of Guizhou Black goats and Meigu goats and their relationship to reproductive performance through population structure analysis, genetic diversity assessment, and selection signal analysis. Blood samples of 19 Guizhou Black goats and 11 Meigu goats were collected for whole-genome high-throughput sequencing. Using PCA and ADMIXTURE analyses, their population structure and genetic relationships were revealed. Further genetic diversity analysis showed that although there is significant population differentiation, the levels of genetic diversity are similar. Subsequently, these goats were categorized into high-yield and low-yield groups based on their litter sizes, with 15 goats in each group. Then, a selection signal analysis was performed using FST and π ratios for 33,563 SNP loci. The results identified six candidate genes, including KCNIP4, GFRA2, and DGKH, which are significantly associated with high litter performance. These findings enhanced our understanding of the genetic characteristics and population structure of Guizhou Black goats and Meigu goats. Moreover, they provide an important theoretical foundation and scientific basis for further breeding improvements.

The molecular characteristic analysis of TRIB2 gene and its expressional patterns in Bos grunniens tissue and granulosa cells

Tribbles homolog 2 (TRIB2) plays an important role in the follicular development of female mammals. However, its expression and function in the yak (Bos grunniens) are still unclear. In this study, we predicted the molecular characteristics of TRIB2, and revealed its expression pattern in yak (Bos grunniens) tissues and ovarian granulosa cells. We cloned the full length of the yak TRIB2 gene obtained by RT-PCR was 1368 bp and the coding sequence (CDS) was 624 bp, encoding 207 amino acids (AA). Homology analysis showed that the yak TRIB2 is highly conserved among species. TRIB2 was detected to be extensively expressed in seven tissues of the yak liver, spleen, lung, kidney, ovary, oviduct and uterus by qPCR. The expression of TRIB2 mRNA in the ovary during gestation was significantly lower than that in the non-pregnant (p < 0.05). At each stage of follicle development, the TRIB2 mRNA in granulosa cells showed a significant upward trend with the development of follicles. The expression of TRIB2 gradually decreased with the increase of the culture time of the granulosa cells in vitro. In conclusion, these results suggest that TRIB2 may play an important role in the follicular development of yaks.

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.

Transcriptome Analysis Reveals Key miRNA-mRNA Pathways in Ovarian Tissues of Yunshang Black Goats With Different Kidding Numbers

The granulosa cell (GC) is the basic functional unit of follicles, and it is important for promoting follicle growth and sex hormones, as well as growth factor secretion in the process of reproduction. A variety of factors influence granulocyte proliferation, yet there are still many gaps to be filled in target and non-coding RNA regulation. In our study, the differentially expressed (DE) mRNAs and miRNAs were detected by using RNA-seq, and we constructed a mRNA-miRNA network related to goat prolificacy. Then, the goat primary GCs were isolated from the follicle for the function validation of candidate genes and their regulator miRNAs. A total of 2,968 DE mRNAs and 99 DE miRNAs were identified in the high- and low-prolificacy goat by RNA-seq, of which there were 1,553 upregulated and 1,415 downregulated mRNAs, and 80 upregulated and 19 downregulated miRNAs, respectively. JAK3 was identified as highly expressed in the low-prolificacy goats (3 times higher than high-prolificacy goats), and the integrated analysis showed that chi-miR-493-3p was a potential regulator of JAK3. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that JAK3 was involved in the PI3K-Akt signaling pathway, the Jak-STAT signaling pathway, and signaling pathways regulating pluripotency of stem cells. In particular, the PI3K-Akt signaling pathway was a typical pathway for cell proliferation, differentiation, apoptosis, and migration. We found that the chi-miR-493-3p targets JAK3 directly via RT-qPCR, dual fluorescence assays, and Western blot. Furthermore, the expression of JAK3 was significantly decreased by the chi-miR-493-3p mimic and increased by the chi-miR-493-3p inhibitor. The CCK-8 assay showed that overexpression of JAK3 promoted cell proliferation, while inhibiting JAK3 had the opposite effect. The expression of cell proliferation markers CDK4 and cyclin D2 also showed the same results. Moreover, the enzyme-linked immunosorbent assay showed that steroid hormones E₂ and PROG were increased by overexpressing JAK3 and decreased by inhibiting JAK3. Therefore, our results identified a chi-miR-439-3p-JAK3 regulatory pathway, which provided a new insight into the GC proliferation and prolificacy of goat.

Ankyrin-repeat and SOCS box-containing protein 9 (ASB9) regulates ovarian granulosa cells function and MAPK signaling

Ankyrin-repeat and SOCS box-containing proteins (ASB) interact with the elongin B-C adapter via their SOCS box domain and with the cullin and ring box proteins to form E3 ubiquitin ligase complexes within the protein ubiquitination pathway. ASB9 in particular is a differentially expressed gene in ovulatory follicles (OFs) induced by the luteinizing hormone (LH) surge or hCG injection in ovarian granulosa cells (GC) while downregulated in growing dominant follicles. Although ASB9 has been involved in biological processes such as protein modification, the signaling network associated with ASB9 in GC is yet to be fully defined. We previously identified and reported ASB9 interactions and binding partners in GC including PAR1, TAOK1, and TNFAIP6/TSG6. Here, we further investigate ASB9 effects on target binding partners regulation and signaling in GC. CRISPR/Cas9-induced inhibition of ASB9 revealed that ASB9 regulates PAR1, TAOK1, TNFAIP6 as well as genes associated with proliferation and cell cycle progression such as PCNA, CCND2, and CCNE2 while CCNA2 was not affected. Inhibition of ASB9 was also associated with increased GC number and decreased caspase3/7 activity, CASP3 expression, and BAX/BCL2 ratio. Furthermore, ASB9 induction in OF in vivo 24 h post-hCG is concomitant with a significant decrease in phosphorylation levels of MAPK3/1 while pMAPK3/1 levels increased following ASB9 inhibition in GC in vitro. Together, these results provide strong evidence for ASB9 as a regulator of GC activity and function by modulating MAPK signaling likely through specific binding partners such as PAR1, therefore controlling GC proliferation and contributing to GC differentiation into luteal cells.

Expression of steroidogenic enzymes and TGFβ superfamily members in follicular cells of prepubertal gilts with distinct endocrine profiles

Regulation of the transforming growth factor beta (TGFβ) superfamily by gonadotrophins in swine follicular cells is not fully understood. This study evaluated the expression of steroidogenic enzymes and members of the TGFβ superfamily in prepubertal gilts allocated to three treatments: 1200 IU eCG at D -3 (eCG); 1200 IU eCG at D -6 plus 500 IU hCG at D -3 (eCG + hCG); and the control, composed of untreated gilts. Blood samples and ovaries were collected at slaughter (D0) and follicular cells were recovered thereafter. Relative gene expression was determined by real-time PCR. Serum progesterone levels were greater in the eCG + hCG group compared with the other groups (P < 0.01). No differences were observed in the expression of BMP15, BMPR1A, BMPR2, FSHR, GDF9, LHCGR and TGFBR1 (P > 0.05). Gilts from the eCG group presented numerically greater mean expression of CYP11A1 mRNA than in the control group that approached statistical significance (P = 0.08) and greater expression of CYP19A1 than in both the eCG and the control groups (P < 0.05). Expression of BMPR1B was lower in the eCG + hCG treatment group compared with the control (P < 0.05). In conclusion, eCG treatment increased the relative expression of steroidogenic enzymes, whereas treatment with eCG + hCG increased serum progesterone levels. Although most of the evaluated TGFβ members were not regulated after gonadotrophin treatment, the downregulation of BMPR1B observed after treatment with eCG + hCG and suggests a role in luteinization regulation.

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.

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.

Interleukins' expression profile changes in granulosa cells of preovulatory follicles during the postpartum period in dairy cows

•

Interleukin (IL) 4 expression was significantly altered during the postpartum period.

•

IL8 and IL15 were the most significantly induced during the postpartum period.

•

Correlation between increased BHB levels and induction of proinflammatory cytokines.

•

Interleukins analyzed were differentially regulated during follicular development.

•

ILs induced during the post-partum period were also induced 24 h post-hCG injection.

The postpartum period in dairy cows is associated with a state of temporary negative energy balance and could induce functional changes into ovarian granulosa cells (GC) resulting in significant impact on the ovarian function and fertility. Yet, the regulation of interleukin receptors (ILRs) in GC as well as ILs expression profile during the postpartum period have not been fully investigated. We hypothesized that the postpartum period is associated with changes in ILs expression profile that could affect follicular development and ovulation rate. First, we aimed to investigate the expression and regulation of different IL and IL receptors in GC at different stages of follicular development and then analyse the changes in target ILs expression profile induced during the postpartum period. In the first objective, normal cycling cows were selected and GC were collected from small follicles (SF), dominant follicles at day 5 of the estrous cycle (DF), and ovulatory follicles, 24 h following hCG injection (OF). In the second objective, dairy cows between 50 and 70 days postpartum were randomly selected, and β-hydroxybutyrate (BHB) concentrations were measured in blood samples in order to assign cows to the BHB⁺ group (>1.4 mmol/L) or BHB⁻ group (<1.2 mmol/L). GC were collected from preovulatory follicles by transvaginal aspiration. Total RNA was extracted from GC of all groups for analysis of target ILs and ILRs expression. Steady-state mRNA levels of IL4R was strongest in the DF, while IL15R expression was greatest in the OF, and IL21R showed increased steady-state mRNA levels in the corpus luteum as compared to the different groups of follicles. Overall, expression of IL1A, IL1B, IL8, IL15, IL23 and TNFα was stronger in OF as compared to DF, while IL4 and IL10 expression was stronger in SF than in DF. Similarly, expression of IL1A, IL1B, IL8, IL15, IL23, and TNFα were significantly stronger in GC of BHB⁺ cows than in the control, while IL4 expression was significantly reduced in BHB⁺ as compared to control cows. We have established an IL expression profile, which suggest a correlation with BHB levels during the postpartum period. Additionally, we have demonstrated a differential regulation of target ILRs in GC at different stages of follicular development. Overall, these data provide a better understanding of the changes that could affect follicular development and ovulation during the postpartum period and lay the ground for further investigations.

Proteomic Analysis of Porcine Pre-ovulatory Follicle Differentiation Into Corpus Luteum

The luteinization of the follicular cells, following a LH surge, causes extensive molecular and structural changes in preovulatory follicles (POF) that lead to ovulation and ultimate formation of the corpus luteum (CL). The objective of this study was to identify proteins expressed in porcine POF before the LH surge and a new CL formed, 2-3 days after ovulation, and evaluate proteome changes associated with formation of the CL from a follicle. We used 2D-gel electrophoresis-based proteomics and tandem mass spectrometry followed by a functional analysis using Ingenuity Pathway analysis (IPA) to evaluate functional pathways associated with the luteinization process. Protein lysates were prepared from isolated POFs and from the newly formed CL. A total of 422 protein spots were identified in both structures. A total of 15 and 48 proteins or their proteoforms were detected only in the POFs and CL, respectively. An IPA analysis of a POF proteome showed that most of the follicular proteins were involved in cellular infiltration, endoplasmic stress responses, and the protein ubiquitination pathway. Most of the early luteal proteins were associated with steroid metabolism, cell death and survival, free radical scavenging, and the protein ubiquitination pathway. A comparison of a follicular proteome with that of an early luteal proteome revealed that 167 identified proteins or their proteoforms were differentially regulated between POFs and the newly formed CL (p < 0.05 and a fold change of >1.8). Proteins that were significantly more abundant in follicles included cAMP-dependent protein kinase, histone binding protein RBBP4, reticulocalbin, vimentin, and calumenin; more abundant luteal proteins included albumin, farnesyl diphosphate synthase, serine protease inhibitors, elongation factor-1, glutaredoxin, and selenium-binding protein. Proteins that were significantly altered with luteal formation were found to be associated with cholesterol biosynthesis, cell death and survival, and acute phase response. Moreover, upstream regulators of differentially abundant proteins in CL were identified that included insulin growth factor-1, sterol regulatory element-binding transcription factor-1, and nuclear factor erythroid-derived 2. We have identified novel proteins that advance our understanding of (1) processes associated with differentiation of POFs into the CL, (2) possible mechanisms of luteal cell survival, and (3) pathways regulating steroidogenesis in the newly formed CL.

Gonadotropin regulation of ankyrin-repeat and SOCS-box protein 9 (ASB9) in ovarian follicles and identification of binding partners

Ankyrin-repeat and SOCS-box protein 9 (ASB9) is a member of the large SOCS-box containing proteins family and acts as the specific substrate recognition component of E3 ubiquitin ligases in the process of ubiquitination and proteasomal degradation. We previously identified ASB9 as a differentially expressed gene in granulosa cells (GC) of bovine ovulatory follicles. This study aimed to further investigate ASB9 mRNA and protein regulation, identify binding partners in GC of bovine ovulatory follicles, and study its function. GC were obtained from small follicles (SF: 2–4 mm), dominant follicles at day 5 of the estrous cycle (DF), and ovulatory follicles, 24 hours following hCG injection (OF). Analyses by RT-PCR showed a 104-fold greater expression of ASB9 in GC of OF than in DF. Steady-state levels of ASB9 in follicular walls (granulosa and theca cells) analyzed at 0, 6, 12, 18 and 24 hours after hCG injection showed a significant induction of ASB9 expression at 12 and 18 hours, reaching a maximum induction of 10.2-fold at 24 hours post-hCG as compared to 0 hour. These results were confirmed in western blot analysis showing strongest ASB9 protein amounts in OF. Yeast two-hybrid screening of OF-cDNAs library resulted in the identification of 10 potential ASB9 binding partners in GC but no interaction was found between ASB9 and creatine kinase B (CKB) in these GC. Functional studies using CRISPR-Cas9 approach revealed that ASB9 inhibition led to increased GC proliferation and modulation of target genes expression. Overall, these results support a physiologically relevant role of ASB9 in the ovulatory follicle by targeting specific proteins likely for degradation, contributing to reduced GC proliferation, and could be involved in the final GC differentiation into luteal cells.

Characterization of the miRNA regulators of the human ovulatory cascade

Ovarian follicular development and ovulation are complex and tightly regulated processes that involve regulation by microRNAs (miRNAs). We previously identified differentially expressed mRNAs between human cumulus granulosa cells (CGCs) from immature early antral follicles (germinal vesicle - GV) and mature preovulatory follicles (metaphase II - M2). In this study, we performed an integrated analysis of the transcriptome and miRNome in CGCs obtained from the GV cumulus-oocyte complex (COC) obtained from IVM and M2 COC obtained from IVF. A total of 43 differentially expressed miRNAs were identified. Using Ingenuity IPA analysis, we identified 7288 potential miRNA-regulated target genes. Two hundred thirty-four of these target genes were also found in our previously generated ovulatory gene library while exhibiting anti-correlated expression to the identified miRNAs. IPA pathway analysis suggested that miR-21 and FOXM1 cooperatively inhibit CDC25A, TOP2A and PRC1. We identified a mechanism for the temporary inhibition of VEGF during ovulation by TGFB1, miR-16-5p and miR-34a-5p. The linkage bioinformatics analysis between the libraries of the coding genes from our preliminary study with the newly generated library of regulatory miRNAs provides us a comprehensive, integrated overview of the miRNA-mRNA co-regulatory networks that may play a key role in controlling post-transcriptomic regulation of the ovulatory process.

Transcriptome analysis of granulosa cells after conventional vs long FSH-induced superstimulation in cattle

BACKGROUND

Prolongation of superstimulatory treatment appears to be associated with a greater superovulatory response and with greater oocyte maturation in cattle. A genome-wide bovine oligo-microarray was used to compare the gene expression of granulosa cells collected from ovarian follicles after differing durations of the growing phase induced by exogenous FSH treatment. Cows were given a conventional (4-day) or long (7-day) superstimulatory treatment (25 mg FSH im at 12-h intervals; n = 6 per group), followed by prostaglandin treatment with last FSH and LH treatment 24 h later. Granulosa cells were harvested 24 h after LH treatment.

RESULTS

The expression of 416 genes was down-regulated and 615 genes was up-regulated in the long FSH group compared to the conventional FSH group. Quantification by RT-PCR of 7 genes (NTS, PTGS2, PTX3, RGS2, INHBA, CCND2 and LRP8) supported the microarrays data. Multigene bioinformatic analysis indicates that markers of fertility and follicle maturity were up-regulated in the long FSH group.

CONCLUSION

Using the large gene expression dataset generated by the genomic analysis and our previous associated with the growth phase and gene expression changes post LH, we can conclude that a prolonged FSH-induced growing phase is associated with transcriptomic characteristics of greater follicular maturity and may therefore be more appropriate for optimizing the superovulatory response and developmental competence of oocytes in cattle.

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.

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.

Copy number variation analysis detects novel candidate genes involved in follicular growth and oocyte maturation in a cohort of premature ovarian failure cases

STUDY QUESTION

Can spontaneous premature ovarian failure (POF) patients derived from population-based biobanks reveal the association between copy number variations (CNVs) and POF?

SUMMARY ANSWER

CNVs can hamper the functional capacity of ovaries by disrupting key genes and pathways essential for proper ovarian function.

WHAT IS KNOWN ALREADY

POF is defined as the cessation of ovarian function before the age of 40 years. POF is a major reason for female infertility, although its cause remains largely unknown.

STUDY DESIGN, SIZE, DURATION

The current retrospective CNV study included 301 spontaneous POF patients and 3188 control individuals registered between 2003 and 2014 at Estonian Genome Center at the University of Tartu (EGCUT) biobank.

PARTICIPANTS/MATERIALS, SETTING, METHODS

DNA samples from 301 spontaneous POF patients were genotyped by Illumina HumanCoreExome (258 samples) and HumanOmniExpress (43 samples) BeadChip arrays. Genotype and phenotype information was drawn from the EGCUT for the 3188 control population samples, previously genotyped with HumanCNV370 and HumanOmniExpress BeadChip arrays. All identified CNVs were subjected to functional enrichment studies for highlighting the POF pathogenesis. Real-time quantitative PCR was used to validate a subset of CNVs. Whole-exome sequencing was performed on six patients carrying hemizygous deletions that encompass genes essential for meiosis or folliculogenesis.

MAIN RESULTS AND THE ROLE OF CHANCE

Eleven novel microdeletions and microduplications that encompass genes relevant to POF were identified. For example, FMN2 (1q43) and SGOL2 (2q33.1) are essential for meiotic progression, while TBP (6q27), SCARB1 (12q24.31), BNC1 (15q25) and ARFGAP3 (22q13.2) are involved in follicular growth and oocyte maturation. The importance of recently discovered hemizygous microdeletions of meiotic genes SYCE1 (10q26.3) and CPEB1 (15q25.2) in POF patients was also corroborated.

LIMITATIONS, REASONS FOR CAUTION

This is a descriptive analysis and no functional studies were performed. Anamnestic data obtained from population-based biobank lacked clinical, biological (hormone levels) or ultrasonographical data, and spontaneous POF was predicted retrospectively by excluding known extraovarian causes for premature menopause.

WIDER IMPLICATIONS OF THE FINDINGS

The present study, with high number of spontaneous POF cases, provides novel data on associations between the genomic aberrations and premature menopause of ovarian cause and demonstrates that population-based biobanks are powerful source of biological samples and clinical data to reveal novel genetic lesions associated with human reproductive health and disease, including POF.

STUDY FUNDING/COMPETING INTEREST

This study was supported by the Estonian Ministry of Education and Research (IUT20-43, IUT20-60, IUT34-16, SF0180027s10 and 9205), Enterprise Estonia (EU30020 and EU48695), Eureka's EUROSTARS programme (NOTED, EU41564), grants from European Union's FP7 Marie Curie Industry-Academia Partnerships and Pathways (IAPP, SARM, |EU324509) and Horizon 2020 innovation programme (WIDENLIFE, 692065), Academy of Finland and the Sigrid Juselius Foundation.

Molecular determinants of a competent bovine corpus luteum: first- vs final-wave dominant follicles

Reproductive management in cattle requires the synchrony of follicle development and oestrus before insemination. However, ovulation of follicles that have not undergone normal physiological maturation can lead to suboptimal luteal function. Here, we investigated the expression of a targeted set of 47 genes in (a) a first-wave vs final-wave dominant follicle (DF; the latter destined to ovulate spontaneously) and (b) 6-day-old corpora lutea (CLs) following either spontaneous ovulation or induced ovulation of a first-wave DF to ascertain their functional significance for competent CL development. Both the mass and progesterone-synthesising capacity of a CL formed following induced ovulation of a first-wave DF were impaired. These impaired CLs had reduced expression of steroidogenic enzymes (e.g. STAR and HSD3B1), luteotrophic receptors (LHCGR) and angiogenic regulators (e.g. VEGFA) and increased expression of BMP2 (linked to luteolysis). Relative to final-wave DFs, characteristic features of first-wave DFs included reduced oestradiol concentrations and a reduced oestradiol:progesterone ratio in the face of increased expression of key steroidogenic enzymes (i.e. CYP11A1, HSD3B1 and CYP19A1) in granulosa cells and reduced expression of the HDL receptor SCARB1 in thecal cells. Transcripts for further components of the TGF and IGF systems (e.g. INHA, INHBA, IGF2R and IGFBP2) varied between the first- and final-wave DFs. These results highlight the importance of hormones such as progesterone interacting with local components of both the TGF and IGF systems to affect the maturation of the ovulatory follicle and functional competency of the subsequent CL.

Effects of intramuscular administration of folic acid and vitamin B12 on granulosa cells gene expression in postpartum dairy cows

The fertility of dairy cows is challenged during early lactation, and better nutritional strategies need to be developed to address this issue. Combined supplementation of folic acid and vitamin B12 improve energy metabolism in the dairy cow during early lactation. Therefore, the present study was undertaken to explore the effects of this supplement on gene expression in granulosa cells from the dominant follicle during the postpartum period. Multiparous Holstein cows received weekly intramuscular injection of 320 mg of folic acid and 10 mg of vitamin B12 (treated group) beginning 24 (standard deviation=4) d before calving until 56 d after calving, whereas the control group received saline. The urea plasma concentration was significantly decreased during the precalving period, and the concentration of both folate and vitamin B12 were increased in treated animals. Milk production and dry matter intake were not significantly different between the 2 groups. Plasma concentrations of folates and vitamin B12 were increased in treated animals. Daily dry matter intake was not significantly different between the 2 groups before [13.5 kg; standard error (SE)=0.5] and after (23.6 kg; SE=0.9) calving. Average energy-corrected milk tended to be greater in vitamin-treated cows, 39.7 (SE=1.4) and 38.1 (SE=1.3) kg/d for treated and control cows, respectively. After calving, average plasma concentration of β-hydroxybutyrate tended to be lower in cows injected with the vitamin supplement, 0.47 (SE=0.04) versus 0.55 (SE=0.03) for treated and control cows, respectively. The ovarian follicle ≥12 mm in diameter was collected by ovum pick-up after estrus synchronization. Recovered follicular fluid volumes were greater in the vitamin-treated group. A microarray platform was used to investigate the effect of treatment on gene expression of granulosa cells. Lower expression of genes involved in the cell cycle and higher expression of genes associated with granulosa cell differentiation before ovulation were observed. Selected candidate genes were analyzed by reverse transcription quantitative PCR. Although the effects of intramuscular injections of folic acid and vitamin B12 on lactational performance and metabolic status of animals were limited, ingenuity pathway analysis of gene expression in granulosa cells suggests a stimulation of cell differentiation in vitamin-treated cows, which may be the result of an increase in LH secretion.

Joint MiRNA/mRNA expression profiling reveals changes consistent with development of dysfunctional corpus luteum after weight gain

Obese women exhibit decreased fertility, high miscarriage rates and dysfunctional corpus luteum (CL), but molecular mechanisms are poorly defined. We hypothesized that weight gain induces alterations in CL gene expression. RNA sequencing was used to identify changes in the CL transcriptome in the vervet monkey (Chlorocebus aethiops) during weight gain. 10 months of high-fat, high-fructose diet (HFHF) resulted in a 20% weight gain for HFHF animals vs. 2% for controls (p = 0.03) and a 66% increase in percent fat mass for HFHF group. Ovulation was confirmed at baseline and after intervention in all animals. CL were collected on luteal day 7-9 based on follicular phase estradiol peak. 432 mRNAs and 9 miRNAs were differentially expressed in response to HFHF diet. Specifically, miR-28, miR-26, and let-7b previously shown to inhibit sex steroid production in human granulosa cells, were up-regulated. Using integrated miRNA and gene expression analysis, we demonstrated changes in 52 coordinately regulated mRNA targets corresponding to opposite changes in miRNA. Specifically, 2 targets of miR-28 and 10 targets of miR-26 were down-regulated, including genes linked to follicular development, steroidogenesis, granulosa cell proliferation and survival. To the best of our knowledge, this is the first report of dietary-induced responses of the ovulating ovary to developing adiposity. The observed HFHF diet-induced changes were consistent with development of a dysfunctional CL and provide new mechanistic insights for decreased sex steroid production characteristic of obese women. MiRNAs may represent novel biomarkers of obesity-related subfertility and potential new avenues for therapeutic intervention.

Effect in dedicator of cytokinesis 6 (DOCK6) on steroid production in theca cells of follicular cysts

Ovarian follicular cysts are one of the most common causes of reproductive failure in mammals. A comparative gene expression approach may aid in elucidating the causes of ovarian cyst disease. In the present study, the differential display technique was used to identify mRNA sequences that accumulate preferentially in theca cells of bovine cystic follicles. Dedicator of cytokinesis 6 (Dock6) expression was observed in the theca cells of cystic follicles. Small interfering RNA (siRNA) knockdown of Dock6 increased progesterone (P4) production and StAR expression in theca cells of high-estrogen follicular cysts, but did not affect androstenedione (A4) production. We propose that Dock6 may be a marker associated with the development of follicular cysts. Additionally, Dock6 may be involved in the development of cystic follicles by suppressing P4 production rather than increasing A4 production in theca cells.

Differential expression of lysosome-associated protein transmembrane-4 beta (LAPTM4B) in granulosa cells of ovarian follicles and in other bovine tissues

Background

LAPTM4B is a member of the lysosome-associated transmembrane protein superfamily that is differentially expressed in normal human tissues and upregulated in various types of carcinomas. These proteins are thought to be involved in the regulation of cell proliferation and survival. The objective of this study was to investigate the expression of bovine LAPTM4B during ovarian follicular development and in various bovine tissues.

Methods and results

Northern blot analysis revealed a 1.8 kb transcript, with highly variable steady state levels among tissues. RT-PCR analysis showed that LAPTM4B mRNA transcripts were low in granulosa cells of small antral follicles, increased in large dominant follicles, and decreased in ovulatory follicles following injection of human chorionic gonadotropin (hCG; P < 0.003). Ovulatory follicles collected at various times after hCG injection revealed a significant reduction of LAPTM4B mRNA starting at 18 h post-hCG (P < 0.029). Immunobloting analysis using antibodies generated against bovine LAPTM4B recognized proteins of 26.3 and 31.5 kDa in granulosa cells of developing follicles and corpus luteum. Further analyses of affinity-purified His-tag LAPTM4B overexpressed in HEK cells showed that the 31.5 kDa protein represented the ubiquinated isoform of the 26.3 kDa native protein. The 26.3 kDa protein was differentially expressed showing highest amounts in dominant follicles and lowest amounts in ovulatory follicles 24 h post-hCG. Immunohistochemical analyses of LAPTM4B showed marked heterogeneity of labeling signal among tissues, with LAPTM4B mainly localized to perinuclear vesicles, in keeping with its putative lysosomal membrane localization.

Conclusion

This study reports for the first time that bovine LAPTM4B in granulosa cells is present in both unubiquinated and ubiquinated forms, and is differentially expressed in developing ovarian follicles, suggesting a possible role in terminal follicular growth.

Global gene expression in granulosa cells of growing, plateau and atretic dominant follicles in cattle

BACKGROUND

The physiological state of the dominant follicle is important as it may be linked to the competence of the oocyte within. The objective of this study was to analyze, by transcriptomic analysis, the changes occurring in granulosa cells from dominant follicles at different phases of follicular growth.

METHODS

Granulosa cells were collected from slaughterhouse dairy cattle follicles with a diameter greater than 9 mm, and were classified at different phases of follicle growth based on flow cytometry profiles of DNA content after staining with propidium iodide. Three phases were identified based on the proportion of cells in -G1 (less than 2n DNA), G0-G1 (2n DNA) or S-M (more than 2n DNA) and follicles were thus allocated to the growing, plateau or atresia group. Between group analysis (BGA) showed clear segregation of the three groups, and the groups were contrasted against each other in a loop design to identify differently expressed genes. Ingenuity Pathway Analysis (IPA) was used to identify the functions and upstream regulators associated with the observed differently expressed genes.

RESULTS

Major differences were observed between the growth phases. Granulosa cells from follicles in the plateau phase had increased expression of TYRO3 and downregulation of JAM2 compared to growing follicles, supporting the idea of a shift from proliferation to differentiation. On the other hand, genes regulating the response to oxidative stress (VNN1) and angiogenesis (ANGPT2) were upregulated in granulosa cells from atretic follicles. While the predicted activated functions in cells at the plateau stage compared to cells at the growing stage included synthesis and transport of molecules, the predictions for atretic follicles relative to plateau ones included an increase in apoptosis and cell death.

CONCLUSION

Consistent with previous studies, these observations allowed us to match the presence of specific gene transcripts to a particular physiological status and consequently to classify follicles. The results also demonstrated that the plateau phase is not a simple 'in between' status between growth and atresia, as several characteristics are unique to this stage.

Tribbles expression in cumulus cells is related to oocyte maturation and fatty acid metabolism

BACKGROUND

In mammals, the Tribbles family includes widely expressed serine-threonine kinase-like proteins (TRIB1, TRIB2 and TRIB3) that are involved in multiple biological processes including cell proliferation and fatty acid (FA) metabolism. Our recent studies highlighted the importance of FA metabolism in cumulus cells (CC) during oocyte maturation in vertebrates and reported a higher TRIB1 expression in CC surrounding in vivo mature oocytes as compared to immature ooocytes in mice and cows. The objective was to investigate Tribbles expression patterns in bovine CC during in vitro maturation (IVM) and to examine their roles in the cumulus-oocyte complex.

METHODS

Tribbles gene expression was analyzed in bovine and murine CC using quantitative RT-PCR. Proteins were detected using Western blot and intracellular localization was assessed by immunofluorescence. Bovine COCs were treated with etomoxir, an inhibitor of FA oxidation (FAO) which blocks CPT1 activity, during 6 h and 18 h IVM. Oocyte meiotic stage was assessed and expression of Tribbles and lipid metabolism genes was quantified in CC.

RESULTS AND DISCUSSION

TRIB1 and TRIB3 were more strongly expressed whereas TRIB2 was less expressed in CC surrounding the oocytes from preovulatory follicles than in CC of immature ones. In CC, Tribbles were located in the cytoplasm and nucleus; in mitotic cells TRIB2 and TRIB3 were detected in the spindle. In the oocyte, Tribbles proteins were detected in the ooplasm; also TRIB2 and TRIB3 were more accumulated in the germinal vesicle. In bovine CC, expression of TRIB1 and TRIB3 was transiently increased at a time preceding oocyte meiosis resumption in vitro. Treatment with etomoxir (150 μM) during IVM resulted in a significant reduction of oocyte maturation rate and decreased MAPK3/1 phosphorylation in the oocytes. In CC, 18 h IVM of etomoxir treatment significantly increased expression of TRIB1, TRIB3, CPTA1 (enzyme regulating FA entry in mitochondria for FAO) and CD36 (thrombospondin receptor involved in FA transport). Under the same conditions, expression of TRIB2 and ACACA (Acetyl coenzyme A carboxylase involved in FA synthesis) decreased in CC. All considered, Tribbles family members may be involved in cell proliferation and in FAO signaling in CC and participate in oocyte meiotic resumption regulation.

Follicle-stimulating hormone receptor (FSHR) alternative skipping of exon 2 or 3 affects ovarian response to FSH

Genes critical for fertility are highly conserved in mammals. Interspecies DNA sequence variation, resulting in amino acid substitutions and post-transcriptional modifications, including alternative splicing, are a result of evolution and speciation. The mammalian follicle-stimulating hormone receptor (FSHR) gene encodes distinct species-specific forms by alternative splicing. Skipping of exon 2 of the human FSHR was reported in women of North American origin and correlated with low response to ovarian stimulation with exogenous follicle-stimulating hormone (FSH). To determine whether this variant correlated with low response in women of different genetic backgrounds, we performed a blinded retrospective observational study in a Turkish cohort. Ovarian response was determined as low, intermediate or high according to retrieved oocyte numbers after classifying patients in four age groups (<35, 35-37, 38-40, >40). Cumulus cells collected from 96 women undergoing IVF/ICSI following controlled ovarian hyperstimulation revealed four alternatively spliced FSHR products in seven patients (8%): exon 2 deletion in four patients; exon 3 and exons 2 + 3 deletion in one patient each, and a retention of an intron 1 fragment in one patient. In all others (92%) splicing was intact. Alternative skipping of exons 2, 3 or 2 + 3 were exclusive to low responders and was independent of the use of agonist or antagonist. Interestingly, skipping of exon 3 occurs naturally in the ovaries of domestic cats--a good comparative model for human fertility. We tested the signaling potential of human and cat variants after transfection in HEK293 cells and FSH stimulation. None of the splicing variants initiated cAMP signaling despite high FSH doses, unlike full-length proteins. These data substantiate the occurrence of FSHR exon skipping in a subgroup of low responders and suggest that species-specific regulation of FSHR splicing plays diverse roles in mammalian ovarian function.

Possible role of IGF2 receptors in regulating selection of 2 dominant follicles in cattle selected for twin ovulations and births

Abundance of IGF-2 receptor (IGF2R), FSH receptor (FSHR), and LH receptor (LHCGR) mRNA in granulosa cells (GCs) or theca cells (TCs) or both cells as well as estradiol (E2), progesterone (P4), and androstenedione concentrations in follicular fluid were compared in cows genetically selected (Twinner) or not selected (control) for multiple ovulations and twin births. Cows were slaughtered at day 3 to 4 (day 3) and day 5 to 6 (day 5) of an estrous cycle, and ovaries, follicular fluid, GCs, and TCs were collected. The two largest (F1 and F2) E2-active (EA) and E2-inactive (EI) follicles were selected according to their E2-to-P4 ratio and diameter. Androstenedione levels in EA F1 and F2 follicles were 5-fold greater (P < 0.05) in Twinner cows than in control cows on day 3 but did not differ on day 5. Twinner cows also had greater (P < 0.05) E2 and P4 concentrations, whereas steroid levels in EI follicles did not differ (P > 0.10) between genotypes. In EA F2 follicles, IGF2R levels in GCs were greater (P < 0.05) in control cows than in Twinner cows on day 3 and day 5, whereas IGF2R mRNA in TCs did not differ (P > 0.10). On day 3, FSHR mRNA levels were greater (P < 0.05) in GCs of EA F1 and EI F2 follicles of control cows than of Twinner cows. LH receptor mRNA expression was less in GCs and greater in TCs of EA F2 follicles in control cows than in Twinner cows (P < 0.05). We hypothesize that reduced GC IGF2R expression in F2 follicles of Twinner cows may play a role in the development of 2 or more dominant follicles.

Research resource: preovulatory LH surge effects on follicular theca and granulosa transcriptomes

The molecular mechanisms that regulate the pivotal transformation processes observed in the follicular wall following the preovulatory LH surge, are still not established, particularly for cells of the thecal layer. To elucidate thecal cell (TC) and granulosa cell (GC) type-specific biologic functions and signaling pathways, large dominant bovine follicles were collected before and 21 hours after an exogenous GnRH-induced LH surge. Antral GCs (aGCs; aspirated by follicular puncture) and membrane-associated GCs (mGCs; scraped from the follicular wall) were compared with TC expression profiles determined by mRNA microarrays. Of the approximately 11 000 total genes expressed in the periovulatory follicle, only 2% of thecal vs 25% of the granulosa genes changed in response to the LH surge. The majority of the 203 LH-regulated thecal genes were also LH regulated in GCs, leaving a total of 57 genes as LH-regulated TC-specific genes. Of the 57 thecal-specific LH-regulated genes, 74% were down-regulated including CYP17A1 and NR5A1, whereas most other genes are being identified for the first time within theca. Many of the newly identified up-regulated thecal genes (eg, PTX3, RND3, PPP4R4) were also up-regulated in granulosa. Minimal expression differences were observed between aGCs and mGCs; however, transcripts encoding extracellular proteins (NID2) and matrix modulators (ADAMTS1, SASH1) dominated these differences. We also identified large numbers of unknown LH-regulated GC genes and discuss their putative roles in ovarian function. This Research Resource provides an easy-to-access global evaluation of LH regulation in TCs and GCs that implicates numerous molecular pathways heretofore unknown within the follicle.

Follicular-phase concentrations of progesterone, estradiol-17β, LH, FSH, and a PGF2α metabolite and daily clustering of prolactin pulses, based on hourly blood sampling and hourly detection of ovulation in heifers

Circulating concentrations of hormones were determined each hour in 13 heifers from the end of the luteolytic period to ovulation (follicular phase, 3.5 days). Diameter of the preovulatory follicle was determined every 8 hours, and the time of ovulation was determined hourly. The diameter of the preovulatory follicle decreased 0.8 ± 0.1 mm/h in heifers when there was 1 to 3 hours between the last two diameter measurements before ovulation. The concentration of progesterone (P4) after the end of the luteolytic period (P4 < 1 ng/mL) changed (P < 0.0001), as shown by a continued decrease until Hour -57 (Hour 0 = ovulation), then was maintained at approximately 0.2 ng/mL until 2 hours before the peak of the LH surge at Hour -26, and then a decrease to 0.1 ng/mL along with a decrease in estradiol-17β. Concentrations of LH gradually increased (P < 0.007) and concentrations of FSH gradually decreased (P < 0.0001) after the end of luteolysis until the beginning nadirs of the respective preovulatory surges. A cluster of prolactin (PRL) pulses occurred (P < 0.0001) each day with approximately 24 hours between the maximum value of successive clusters. Hourly concentrations of a PGF2α metabolite decreased (P < 0.007) until Hour -40, but did not differ among hours thereafter. Novel observations included the gradual increase in LH and decrease in FSH until the beginning of the preovulatory surges and follicle diameter decrease a few hours before ovulation. Results supported the following hypotheses: (1) change in the low circulating P4 concentrations during the follicular phase are temporally associated with change in LH concentrations; and (2) PRL pulses occur in a cluster each day during the follicular phase of the estrous cycle.

Gene expression profiling of granulosa cells from PCOS patients following varying doses of human chorionic gonadotropin

PURPOSE

Human chorionic gonadotrophin (hCG) has been used to induce ovulation and oocyte maturation. Although the most common dose of hCG used in IVF is 10,000 IU, there are reports that suggest 5,000 IU is sufficient to yield similar results. The objective of this study is to evaluate the dose dependent differences in gene expression of granulosa cells following various doses of hCG treatment.

METHODS

Patients with polycystic ovarian syndrome (PCOS) were stimulated for IVF treatment. The hCG injection was either withheld or given at 5,000 or 10,000 IU. Granulosa cells from the follicular fluids have been collected for RNA isolation and analyzed using Affymetrix genechip arrays.

RESULTS

Unsupervised hierarchical clustering based on whole gene expression revealed two distinct groups of patients in this experiment. All untreated patients were clustered together whereas hCG-treated patients separated to a different group regardless of the dose. A large number of the transcripts were similarly up- or down-regulated across both hCG doses (2229 and 1945 transcripts, respectively). However, we observed dose-dependent statistically significant differences in gene expression in only 15 transcripts.

CONCLUSIONS

Although hCG injection caused a major change in the gene expression profile of granulosa cells, 10,000 IU hCG resulted in minimal changes in the gene expression profiles of granulosa cells as compared with 5,000 IU. Thus, based on our results, we suggest the use of 10,000 IU hCG should be reconsidered in PCOS patients.

The cloning, characterization, and expression profiling of the LRP8 gene in duck (Anas platyrhynchos)

Low-density lipoprotein receptor-related protein 8 (LRP8) is a member of the low-density lipoprotein receptor gene family that functions in body lipoprotein homeostasis. In this study, reverse transcription-polymerase chain reaction, rapid amplification of cDNA ends, and real-time PCR were performed to characterize the duck LRP8 gene. The cDNA of duck LRP8 contained a 14-bp 5' UTR, a 2754-bp open reading frame, and a 189-bp 3' UTR. The duck LRP8 encoded a protein of 917 amino acid residues composed of five functional domains and resembling other members of the LDLR family, and it displayed high nucleotide and amino acid homology to the LRP8 sequences in other avian species. The mRNA expression level of LRP8 was greater in duck extra-hepatic adipose tissue than in the liver. The peak expression values of LRP8 in both liver and adipose tissues occurred at week 1 and were significantly higher than the values observed during any other week (p < 0.05). Differences in the expression patterns of LRP8 mRNA from weeks 2 to 8 of growth were observed in different organs. A consistent low expression was observed in the liver, and fluctuating expression was observed in the subcutaneous adipose tissue (up- and then down-regulated) and abdominal adipose tissue (down-, then up-, then down-regulated). These findings suggest that LRP8 might play more important roles in regulating lipid metabolism in extra-hepatic adipose tissues than in the liver during early growth after hatching in the duck.

Vascular endothelial growth factor-induced expression of its receptors and activation of the MAPK signaling pathway during ovine oocyte maturation in vitro

The vascular endothelial growth factor (VEGF) has beneficial effects on ovine oocytes during in vitro maturation and their subsequent early embryonic development, but the biochemical pathway underlying this effect has not been elucidated. Therefore, the focus of the present study was to investigate the activation of the mitogen-activated protein kinase (MAPK) pathway in response to the addition of VEGF to the maturation medium, and to study the subcellular localization of VEGF and its receptors during ovine oocyte maturation. We concluded that: (1) VEGF mainly localized in the cytoplasm, whereas its receptors, fms-tyrosine kinase-1 and kinase domain region (KDR), were localized on the plasma membrane of oocytes; (2) the addition of 5 ng/mL VEGF increased the percentage of oocytes with extruded first polar bodies (50.9 ± 2.2% vs. 34.6 ± 2.9%; treatment vs. control, respectively; P < 0.01) and the rate of oocytes competent to undergo nuclear maturation (70.6 ± 0.9% vs. 62.9 ± 1.9%, P < 0.01); and (3) as the expression of VEGF, fms-tyrosine kinase-1, and KDR increased after supplementation with 5 ng/mL, expression of VEGF, mitogen-activated protein kinase kinase (MEK), and MAPK mRNA, as well as MAPK phosphorylation, were stimulated in a time-dependent manner. We inferred that, in a paracrine manner, exogenous VEGF bound to KDR, its main receptor, and then activated the MAPK signaling pathway, which promoted maturation of ovine oocytes. However, the VEGF system also had an autocrine regulatory loop that contributed to creating an environment optimal for oocyte maturation.

Impact of the LH surge on granulosa cell transcript levels as markers of oocyte developmental competence in cattle

In the case of in vitro embryonic production, it is known that not all oocytes detain the developmental capacity to form an embryo. This capacity appears to be acquired through completion of folliculogenesis, during which the oocyte and follicular cells influence their respective destinies. The differentiation status of granulosa cells (GCs) could therefore offer an indicator of oocyte quality. The aim of this study was to compare mRNA transcript abundance in GCs associated with oocytes that subsequently reach or not the blastocyst stage. GCs were collected from cattle following an ovarian stimulation protocol that did or did not include the administration of LH. GCs were classified according to the developmental stage achieved by the associated oocytes. Transcript abundance was measured by microarray. Follicles (n=189) obtained from cows before and after the LH surge were essentially similar and the rates of oocytes reaching the blastocyst stage were not significantly different (52 vs 41%), but blastocyst quality was significantly better in the post-LH-surge group. In GCs from the pre-LH-surge group and associated with developmentally competent oocytes, 18 overexpressed and 22 underexpressed transcripts were found, including novel uncharacterized transcripts, whereas no differentially expressed transcripts were associated with developmentally different oocytes in the post-LH-surge group. The novel transcriptomic response associated with LH appeared to mask the difference. Based on oocyte developmental competence, the period prior to the LH surge appears best suited for studying competence-associated mRNA transcripts in bovine follicle cells.

Extensive effects of in vitro oocyte maturation on rhesus monkey cumulus cell transcriptome

The elaboration of a quality oocyte is integrally linked to the correct developmental progression of cumulus cell phenotype. In humans and nonhuman primates, oocyte quality is diminished with in vitro maturation. To determine the changes in gene expression in rhesus monkey cumulus cells (CC) that occur during the final day prior to oocyte maturation and how these changes differ between in vitro (IVM) and in vivo maturation (VVM), we completed a detailed comparison of transcriptomes using the Affymetrix gene array. We observed a large number of genes differing in expression when comparing IVM-CC and VVM-CC directly but a much larger number of differences when comparing the transitions from the prematuration to the post-IVM and post-VVM states. We observed a truncation or delay in the normal pattern of gene regulation but also remarkable compensatory changes in gene expression during IVM. Among the genes affected by IVM are those that contribute to productive cell-cell interactions between cumulus cell and oocyte and between cumulus cells. Numerous genes involved in lipid metabolism are incorrectly regulated during IVM, and the synthesis of sex hormones appears not to be suppressed during IVM. We identified a panel of 24 marker genes, the expression of which should provide the foundation for understanding how IVM can be improved for monitoring IVM conditions and for diagnosing oocyte quality.

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.

Differential gene expression of serine protease inhibitors in bovine ovarian follicle: possible involvement in follicular growth and atresia

BACKGROUND

SERPINs (serine protease inhibitors) regulate proteases involving fibrinolysis, coagulation, inflammation, cell mobility, cellular differentiation and apoptosis. This study aimed to investigate differentially expressed genes of members of the SERPIN superfamily between healthy and atretic follicles using a combination of microarray and quantitative real-time PCR (QPCR) analysis. In addition, we further determined mRNA and protein localization of identified SERPINs in estradiol (E2)-active and E2-inactive follicles by in situ hybridization and immunohistochemistry.

METHODS

We performed microarray analysis of healthy (10.7 +/- 0.7 mm) and atretic (7.8 +/- 0.2 mm) follicles using a custom-made bovine oligonucleotide microarray to screen differentially expressed genes encoding SERPIN superfamily members between groups. The expression profiles of six identified SERPIN genes were further confirmed by QPCR analysis. In addition, mRNA and protein localization of four SERPINs was investigated in E2-active and E2-inactive follicles using in situ hybridization and immunohistochemistry.

RESULTS

We have identified 11 SERPIN genes expressed in healthy and atretic follicles by microarray analysis. QPCR analysis confirmed that mRNA expression of four SERPINs (SERPINA5, SERPINB6, SERPINE2 and SERPINF2) was greater in healthy than in atretic follicles, while two SERPINs (SERPINE1 and SERPING1) had greater expression in atretic than in healthy follicles. In situ hybridization showed that SERPINA5, SERPINB6 and SERPINF2 mRNA were localized in GCs of E2-active follicles and weakly expressed in GCs of E2-inactive follicles. SERPING1 mRNA was localized in both GCs and the theca layer (TL) of E2-inactive follicles and a weak hybridization signal was also detected in both GCs and TL of E2-active follicles. Immunohistochemistry showed that SERPINA5, SERPINB6 and SERPINF2 were detected in GCs of E2-active and E2-inactive follicles. SERPING1 protein was localized in both GCs and the TL of E2-active and E2-inactive follicles.

CONCLUSIONS

Our results demonstrate a characteristic expression of SERPIN superfamily member genes in bovine healthy and atretic follicles. The cell-type-and stage-specific expression of SERPINs may be associated with bovine follicular growth and atresia.

New insights into the pathogenesis of cystic follicles in cattle: microarray analysis of gene expression in granulosa cells

Ovarian follicular growth and development are regulated by extraovarian and intraovarian factors, which influence granulosa cell proliferation and differentiation. However, the molecular mechanisms that drive follicular growth are not completely understood. Ovarian follicular cysts are one of the most common causes of reproductive failure in dairy cattle. Nevertheless, the primary cause of cyst formation has not been clearly established. A gene expression comparison may aid in elucidating the causes of ovarian cyst disease. Our objective was to identify differentially expressed genes in ovarian granulosa cells between normal dominant and cystic follicles of cattle. Granulosa cells and follicular fluid were isolated from dominant and cystic follicles collected via either ultrasound-guided aspiration from dairy cows (n = 24) or slaughterhouse ovaries from beef cows (n = 23). Hormonal analysis for progesterone, estradiol, and androstenedione in follicular fluid was performed by RIA. Total RNA was extracted and hybridized to 6 Affymetrix GeneChip Bovine Genome Arrays (Affymetrix, Santa Clara, CA). Abundance of mRNA for differentially expressed selected genes was determined through quantitative real-time reverse-transcription PCR. Follicular cysts showed greater (P < 0.05) progesterone, lesser (P < 0.05) estradiol, and no differences (P > 0.10) in androstenedione concentrations compared with noncystic follicles. A total of 163 gene sequences were differentially expressed (P < 0.01), with 19 upregulated and 144 downregulated. From selected target genes, quantitative real-time reverse-transcription PCR confirmed angiogenin, PGE(2) receptor 4, and G-protein coupled receptor 34 genes as upregulated in cystic follicles, and Indian hedgehog protein precursor and secreted frizzled-related protein 4 genes as downregulated in cystic follicles. Further research is required to elucidate the role of these factors in follicular development and cyst formation.

Molecular characterization and expression of DERL1 in bovine ovarian follicles and corpora lutea

The endoplasmic reticulum (ER) is a major site of protein synthesis and facilitates the folding and assembly of newly synthesized proteins. Misfolded proteins are retrotranslocated across the ER membrane and destroyed at the proteasome. DERL1 is an important protein involved in the retrotranslocation and degradation of a subset of misfolded proteins from the ER. We characterized a 2617 bp cDNA from bovine granulosa cells that corresponded to bovine DERL1. Two transcripts of 3 and 2.6 kb were detected by Northern blot analysis, and showed variations in expression among tissues. During follicular development, DERL1 expression was greater in day 5 dominant follicles compared to small follicles, ovulatory follicles, or corpus luteum (CL). Within the CL, DERL1 mRNA expression was intermediate in midcycle, and lowest in late cycle as compared to early in the estrous cycle. Western blot analyses demonstrated the presence of DERL1 in the bovine CL at days 5, 11, and 18 of the estrous cycle. Co-immunoprecipitation using luteal tissues showed that DERL1 interacts with class I MHC but not with VIMP or p97 ATPase. The interaction between DERL1 and MHC I suggests that, in the CL, DERL1 may regulate the integrity of MHC I molecules that are transported to the ER membrane. Furthermore, the greater expression of DERL1 mRNA is associated with the active follicular development and early luteal stages, suggesting a role of DERL1 in tissue remodeling events and maintenance of function in reproductive tissues.

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.

Temporal relationships of the LH surge and ovulation to echotexture and power Doppler signals of blood flow in the wall of the preovulatory follicle in heifers

Changes in echotexture and blood flow in the wall of preovulatory follicles in heifers were studied in relation to the LH surge and ovulation in gonadotrophin-releasing hormone-induced (n = 7; Experiment 1) and spontaneous (n = 8; Experiment 2) ovulators. Ultrasonographic examinations and blood sampling were performed either every hour (Experiment 1) or every 6 h (Experiment 2). The interval from LH peak to ovulation in induced and spontaneous ovulators was 27.1 ± 0.3 and 34.5 ± 1.5 h, respectively. Follicle diameter did not increase between the LH peak and ovulation. In the induced ovulators, serration of the stratum granulosum was detected in one (14%), two (29%), three (43%) and four (57%) heifers at 4, 3, 2 and 1 h before ovulation, respectively. An initial increase in blood flow (P < 0.001) encompassed the LH peak in both experiments. In the induced ovulators, blood flow increased (P < 0.02) to maximum 3 h after the LH peak, maintained a plateau for 5 h, decreased (P < 0.05) between 9 and 14 h, increased (P < 0.05) again between 19 and 21 h and then decreased (P < 0.01) between 25 and 26 h (1 h before ovulation). The biphasic increase and decrease in blood flow and serration of the granulosum in the wall of the preovulatory follicle in cattle are novel findings.

Equine thrombospondin II and secreted protein acidic and cysteine-rich in a model of normal and pathological wound repair

Wound healing in horses is complicated, particularly when wounds are on the limb. The objectives of this study were to clone equine thrombospondin II (THBS2) and secreted protein acidic and cysteine-rich (SPARC) cDNAs and to compare the spatiotemporal expression of mRNAs and proteins during repair of body and limb wounds. These molecules were targeted in view of their potential biological contribution to angiogenesis, which is exacerbated during the repair of limb wounds in horses. Cloning was achieved by screening size-selected cDNA libraries previously derived from 7-day-old wounds. Expression was studied in unwounded skin and in samples from 1, 2, 3, 4, and 6 wk old wounds of the body and limb. Temporal gene expression was determined by semiquantitative RT-PCR, while protein expression was mapped immunohistochemically. The temporal pattern of expression for both genes was similar; wounding caused immediate upregulation of mRNA, which did not return to baseline by the end of the study, and overexpression was noted in body relative to limb wounds. Immunostaining for THBS2 and SPARC was induced by wounding, though no differences in stain location or intensity were detected between body and limb wounds. This study is the first to characterize equine cDNA for THBS2 and SPARC and to document mRNA expression over the different phases of repair. THBS2 and SPARC might modulate angiogenesis during wound healing in the horse, which could protect against the disproportionate fibroplasia commonly afflicting limb wounds and leading to the development of exuberant granulation tissue.

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.

Role of transforming growth factor-β1 in gene expression and activity of estradiol and progesterone-generating enzymes in FSH-stimulated bovine granulosa cells

Survival and inhibitory factors regulate steroidogenesis and determine the fate of developing follicles. The objective of this study was to determine the role of transforming growth factor-β1 (TGFB1) in the regulation of estradiol-17β (E2) and progesterone (P4) secretion in FSH-stimulated bovine granulosa cells. Granulosa cells were obtained from 2 to 5 mm follicles and cultured in serum-free medium. FSH dose (1 and 10 ng/ml for 6 days) and time in culture (2, 4, and 6 days with 1 ng/ml FSH) increased E2secretion and mRNA expression of E2-related enzymes cytochrome P450 aromatase (CYP19A1) and 17β-hydroxysteroid dehydrogenase type 1 (HSD17B1), but notHSD17B7. TGFB1 in the presence of FSH (1 ng/ml) inhibited E2secretion, and decreased mRNA expression of FSH receptor(FSHR),CYP19A1, andHSD17B1, but notHSD17B7. FSH dose did not affect P4secretion and mRNA expression of 3β-hydroxysteroid dehydrogenase (HSD3B) and α-glutathioneS-transferase (GSTA), but inhibited the amount of steroidogenic acute regulatory protein(STAR)mRNA. Conversely, P4and mRNA expression ofSTAR, cytochrome P450 side-chain cleavage(CYP11A1),HSD3B, andGSTAincreased with time in culture. TGFB1 inhibited P4secretion and decreased mRNA expression ofSTAR,CYP11A1,HSD3B, andGSTA. TGFB1 modified the formation of granulosa cell clumps and reduced total cell protein. Finally, TGFB1 decreased conversion of androgens to E2, but did not decrease the conversion of estrone (E1) to E2and pregnenolone to P4. Overall, these results indicate that TGFB1 counteracts stimulation of E2and P4synthesis in granulosa cells by inhibiting key enzymes involved in the conversion of androgens to E2and cholesterol to P4without shutting down HSD17B reducing activity and HSD3B activity.

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.