Supreme glutathione-dependent ketosteroid isomerase in the yellow-fever transmitting mosquito Aedes aegypti.

The steroid hormone ecdysone is essential for the reproduction and survival of insects. The hormone is synthesized from dietary sterols such as cholesterol, yielding ecdysone in a series of consecutive enzymatic reactions. In the insect orders Lepidoptera and Diptera a glutathione transferase called Noppera-bo (Nobo) plays an essential, but biochemically uncharacterized, role in ecdysteroid biosynthesis. The Nobo enzyme is consequently a possible target in harmful dipterans, such as disease-carrying mosquitoes. Flavonoid compounds inhibit Nobo and have larvicidal effects in the yellow-fever transmitting mosquito Aedes aegypti, but the enzyme is functionally incompletely characterized. We here report that within a set of glutathione transferase substrates the double-bond isomerase activity with 5-androsten-3,17-dione stands out with an extraordinary specific activity of 4000 µmol min-1 mg-1. We suggest that the authentic function of Nobo is catalysis of a chemically analogous ketosteroid isomerization in ecdysone biosynthesis.

Importance of 3ß-hydroxysteroid dehydrogenases and their clinical use in prostate cancer.

Androgen receptor signaling is crucial for the development of treatment resistance in prostate cancer. Among steroidogenic enzymes, 3ß-hydroxysteroid dehydrogenases (3ßHSDs) play critical roles in extragonadal androgen synthesis, especially 3ßHSD1. Increased expression of 3ßHSDs is observed in castration-resistant prostate cancer tumors compared with primary prostate tumors, indicating their involvement in castration resistance. Recent studies link 3ßHSD1 to resistance to androgen receptor signaling inhibitors. The regulation of 3ßHSD1 expression involves various factors, including transcription factors, microenvironmental influences, and posttranscriptional modifications. Additionally, the clinical significance of HSD3B1 genotypes, particularly the rs1047303 variant, has been extensively studied. The impact of HSD3B1 genotypes on treatment outcomes varies according to the therapy administered, suggesting the potential of HSD3B1 genotyping for personalized medicine. Targeting 3ßHSDs may be a promising strategy for prostate cancer management. Overall, understanding the roles of 3ßHSDs and their genetic variations may enable the development and optimization of novel treatments for prostate cancer.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis.

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Association Between Adrenal-Restrictive HSD3B1 Inheritance and Hormone-Independent Subtypes of Endometrial and Breast Cancer.

BACKGROUND: The germline variant rs1047303 (HSD3B1[1245A/C]), restricting or enabling production of potent androgens and estrogens from adrenal precursors, affects outcomes of castration-resistant prostate cancer and is associated with estrogen receptor positivity in postmenopausal breast cancer. Like breast cancer, endometrial cancer is another malignancy with hormone-dependent and hormone-independent subtypes. We hypothesized that adrenal-restrictive HSD3B1 genotype would associate with hormone-independent cancer subtypes. METHODS: We employed a previously described classification of tumors in The Cancer Genome Atlas into genomic clusters. We determined HSD3B1 genotype frequencies by endometrial cancer genomic cluster and calculated the odds per adrenal-restrictive A allele for the largely hormone-independent copy-number (CN) high subtype vs other subtypes. An equivalent analysis was performed for the genomically similar, hormone-independent basal breast cancer subtype. Last, we performed survival analyses for UK Biobank participants with endometrial cancer by HSD3B1 genotype. All statistical tests were 2-sided. RESULTS: The adrenal-restrictive HSD3B1(1245A) allele was associated with the CN-high endometrial cancer subtype (odds ratio [OR] = 1.63, 95% confidence interval [CI] = 1.14 to 2.32; P = .007). Similarly, HSD3B1(1245A) was associated with the basal breast cancer subtype (OR = 1.54, 95% CI = 1.13 to 2.08; P = .006). In the UK Biobank, endometrial cancer patients homozygous for HSD3B1(1245A) had worse overall (hazard ratio [HR] = 1.39, 95% CI = 1.16 to 1.68; P < .001) and cancer-specific (HR = 1.39, 95% CI = 1.14 to 1.70; P = .001) survival, consistent with the A allele being enriched in the more aggressive CN-high subtype. CONCLUSIONS: These findings suggest roles for adrenal-restrictive vs adrenal-permissive steroidogenesis, by way of rs1047303 genotype, in the development of and/or outcomes from at least 3 commonly hormone-associated types of cancer: prostate, breast, and endometrial.

Role of single nucleotide polymorphisms of the HSD3B1 gene (rs6203 and rs33937873) in the prediction of prostate cancer risk.

3­ß­hydroxysteroid dehydrogenase 1 (HSD3B1) is shown to affect dihydrotestosterone level in prostatic tissue which is a risk factor for prostate cancer (PC). The present study aimed to determine whether rs33937873 (G313A) and rs6203 (C338T) single nucleotide polymorphisms (SNP) in HSD3B1 gene was a potential risk factor for PC susceptibility and can predict the recurrence of PC in Egyptian patients. A total of 186 Egyptian patients were selected with incident primary PC and compared with 180 age healthy controls. The frequencies and the main effect of rs33937873 and rs6203 in HSD3B1 were compared and investigated between the patients and control using genotyping technique and statistical analysis. The mutant GA genotype of G313A in rs33937873 SNP was considered as an independent risk for PC in the multivariate regression analysis [odds ratio (OR)=2.7, 95% confidence intervals (CI): 1.2­5.5, P=0.01] together with positive history of hypertension (HTN) (OR=6.2, 95% CI: 3.2­12.1, P=0.0001) and begin prostatic hyperplasia (BPH; OR=8.9, 95% CI: 4.5­17.5, P=0.0001). Conversely, in rs6203 (C338T), C allele is considered as major risk allele in the development of PC (OR=1.8, 95% CI: 1.3­2.4, P=0.0003). The univariate logistic regression analyses indicated that CC genotype of rs6203 was a PC risk factor (OR=1.9, 95% CI: 1.3­2.9, P=0.002). In addition, the frequency of the A­C haplotype established by rs33937873­rs6203 was also significantly higher for PC (P=0.013). The predication of PC recurrence was associated only with positive family history (OR=7.7, 95% CI: 2.3­25.9, P=0.001) and not for The G313A and C338T SNPs. These results suggested that the two HSD3B1 polymorphisms rs33937873 and rs6203 may modify the risk of PC, particularly among patients with HTN and history of BPH, suggesting them as prominent future markers for prediction of PC risk.

Association of HSD3B1 Genotype and Clinical Outcomes in Postmenopausal Estrogen-Receptor-Positive Breast Cancer.

BACKGROUND: Homozygous inheritance of a single-nucleotide polymorphism (1245A > C) in HSD3B1 results in an adrenal permissive phenotype of increased adrenal steroid precursor conversion to potent androgens. This is associated with poor outcomes in prostate cancer. We hypothesized that inheritance of the HSD3B1 adrenal permissive genotype would similarly negatively impact breast cancer outcomes. PATIENTS AND METHODS: Germline HSD3B1 was sequenced in 644 postmenopausal women diagnosed between 2004 and 2015 with stage I-III estrogen receptor-positive (ER+), HER2/neu-negative (HER2-) breast cancer enrolled in a population-based study in western Washington. Primary endpoint was distant metastatic recurrence according to genotype. Secondary endpoint was breast cancer-specific survival. Hazard ratios (HR) were calculated using cause-specific Cox regression accounting for competing risks. RESULTS: Adrenal restrictive genotype (homozygous wild type) was most prevalent (47%), followed by heterozygous (44%) and adrenal permissive (9%). There were no significant differences comparing demographic, tumor, or treatment characteristics apart from higher frequency of adrenal permissive genotype among non-Hispanic white participants (p = 0.04). After accounting for competing risks, the cumulative incidence of distant metastatic recurrence (15 events) was significantly higher among participants with adrenal permissive compared with the adrenal restrictive genotype (HR 4.9, 95% CI 1.32-18.4, p = 0.02). The adrenal permissive genotype was also predictive of breast cancer-specific mortality (HR 3.5, 95% CI 1.27-9.59, p = 0.02). CONCLUSIONS: Inheritance of the HSD3B1 adrenal permissive genotype is associated with increased incidence of distant metastasis and higher cause-specific mortality in postmenopausal ER+/HER2- breast cancer. Further research is necessary to understand the effect of excess adrenal androgen metabolism in promoting breast cancer growth and progression.

Letrozole protects against cadmium-induced inhibition of spermatogenesis via LHCGR and Hsd3b6 to activate testosterone synthesis in mice.

The heavy metal cadmium is proposed to be one of the environmental endocrine disruptors of spermatogenesis. Cadmium-induced inhibition of spermatogenesis is associated with a hormone secretion disorder. Letrozole is an aromatase inhibitor that increases peripheral androgen levels and stimulates spermatogenesis. However, the potential protective effects of letrozole on cadmium-induced reproductive toxicity remain to be elucidated. In this study, male mice were administered CdCl2 (4 mg/kg BW) orally by gavage alone or in combination with letrozole (0.25 mg/kg BW) for 30 days. Cd exposure caused a significant decreases in body weight, sperm count, motility, vitality, and plasma testosterone levels. Histopathological changes revealed extensive vacuolization and decreased spermatozoa in the lumen. However, in the Cd + letrozole group, letrozole treatment compensated for deficits in sperm parameters (count, motility, and vitality) induced by Cd. Letrozole treatment significantly increased serum testosterone levels, which were reduced by Cd. Histopathological studies revealed a systematic array of all germ cells, a preserved basement membrane and relatively less vacuolization. For a mechanistic examination, RNA-seq was used to profile alterations in gene expression in response to letrozole. Compared with that in the Cd-treated group, RNA-Seq analysis showed that 214 genes were differentially expressed in the presence of letrozole. Gene ontology (GO) enrichment analysis and KEGG signaling pathway analysis showed that steroid biosynthetic processes were the processes most affected by letrozole treatment. Furthermore, we found that the expression of the testosterone synthesis-related genes LHCGR (luteinizing hormone/choriogonadotropin receptor) and Hsd3b6 (3 beta- and steroid delta-isomerase 6) was significantly downregulated in Cd-treated testes, but these genes maintained similar expression levels in letrozole-treated testes as those in the control group. However, the transcription levels of inflammatory cytokines, such as IL-1ß and IL-6, and oxidative stress-related genes (Nrf2, Nqo1, and Ho-1) showed no changes. The present study suggests that the potential protective effect of letrozole on Cd-induced reproductive toxicity might be mediated by the upregulation of LHCGR and Hsd3b6, which would beneficially increase testosterone synthesis to achieve optimum protection of sperm quality and spermatogenesis.

Retrotransposition disrupting EBP in a girl and her mother with X-linked dominant chondrodysplasia punctata.

X-linked dominant chondrodysplasia punctata (CDPX2) is a rare congenital disorder caused by pathogenic variants in EBP on Xp11.23. We encountered a girl and her mother with CDPX2-compatible phenotypes including punctiform calcification in the neonatal period of the girl, and asymmetric limb shortening and ichthyosis following the Blaschko lines in both subjects. Although Sanger direct sequencing failed to reveal a disease-causing variant in EBP, whole genome sequencing (WGS) followed by Manta analysis identified a ~ 4.5 kb insertion at EBP exon 2 of both subjects. The insertion was associated with the hallmarks of retrotransposition such as an antisense poly(A) tail, a target site duplication, and a consensus endonuclease cleavage site, and the inserted sequence harbored full-length SVA_F1 element with 5'- and 3'-transductions containing the Alu sequence. The results imply the relevance of retrotransposition to the human genetic diseases and the usefulness of WGS in the identification of retrotransposition.

A Testis-Specific Long Noncoding RNA, Start, Is a Regulator of Steroidogenesis in Mouse Leydig Cells.

The testis expresses many long noncoding RNAs (lncRNAs), but their functions and overview of lncRNA variety are not well understood. The mouse Prss/Tessp locus contains six serine protease genes and two lncRNAs that have been suggested to play important roles in spermatogenesis. Here, we found a novel testis-specific lncRNA, Start (Steroidogenesis activating lncRNA in testis), in this locus. Start is 1822 nucleotides in length and was found to be localized mostly in the cytosol of germ cells and Leydig cells, although nuclear localization was also observed. Start-knockout (KO) mice generated by the CRISPR/Cas9 system were fertile and showed no morphological abnormality in adults. However, in adult Start-KO testes, RNA-seq and qRT-PCR analyses revealed an increase in the expression of steroidogenic genes such as Star and Hsd3b1, while ELISA analysis revealed that the testosterone levels in serum and testis were significantly low. Interestingly, at 8 days postpartum, both steroidogenic gene expression and testosterone level were decreased in Start-KO mice. Since overexpression of Start in two Leydig-derived cell lines resulted in elevation of the expression of steroidogenic genes including Star and Hsd3b1, Start is likely to be involved in their upregulation. The increase in expression of steroidogenic genes in adult Start-KO testes might be caused by a secondary effect via the androgen receptor autocrine pathway or the hypothalamus-pituitary-gonadal axis. Additionally, we observed a reduced number of Leydig cells at 8 days postpartum. Collectively, our results strongly suggest that Start is a regulator of steroidogenesis in Leydig cells. The current study provides an insight into the overall picture of the function of testis lncRNAs.

Parallel molecular mechanisms for enzyme temperature adaptation.

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

The association between missense polymorphisms in SRD5A2 and HSD3B1 and treatment failure with abiraterone for castration-resistant prostate cancer.

Missense polymorphism in HSD3B1, encoding 3ß-hydroxysteroid dehydrogenase-1, was associated with outcome after abiraterone treatment. Other androgen-metabolizing enzymes may be involved in therapeutic effect in abiraterone. In this study, we investigated the significance of polymorphisms in genes involved in androgen and abiraterone metabolisms in prostate cancer patients treated with abiraterone. A total of 99 Japanese male castration-resistant prostate cancer patients treated with abiraterone between 2014 and 2018 were included. Genomic DNA was obtained from whole blood samples, and genotyping on SRD5A2 (rs523349), CYP17A1 (rs743572), CYP17A1 (rs2486758), and AKR1C3 (rs12529) was performed by PCR-based technique. Among the 99 patients, 32 (32.3%), 49 (49.5%), and 18 patients (18.2%) carried GG, GC, and CC alleles in SRD5A2, respectively. CC allele was associated with lower risk of treatment failure (hazard ratio, 0.43; 95% confidence interval, 0.20-0.87; P = 0.017) on multivariate analyses, compared with GG/GC alleles. In the combination model using HSD3B1 and SRD5A2 polymorphisms, compared with the combination of AA in HSD3B1 and GG/GC in SRD5A2, other combinations were associated with lower risk of treatment failure (hazard ratio, 0.34; 95% confidence interval, 0.17-0.62; P = 0.0003) on multivariate analyses. This study showed that SRD5A2 genetic variation was associated with the risk of treatment failure in abiraterone. Combinational use of genetic variation in HSD3B1 with SRD5A2 genetic variation augmented the ability of prognostic stratification.

Effects of electromagnetic field (EMF) radiation on androgen synthesis and release from the pig endometrium during the fetal peri-implantation period.

An electromagnetic field (EMF) may have effects on female reproduction. This study was conducted to determine whether EMF [50 and 120 Hz, 2 and 4 h of incubation in the presence or absence of progesterone (P4, 10-5 M)] affects androgen synthesis and release from the pig endometrium. Endometrial slices were collected from pigs (n = 5) during the fetal peri-implantation period (i.e., days 15-16 of gestation) and treated in vitro with EMF. The selected endometrial slices were treated with P4 to determine whether this hormone has effects on protection of the tissue from EMF radiation. The CYP17A1 and HSD3B1 mRNA transcript abundance, steroid 17αhydroxylase/17, 20-lyase (cytochrome P450c17) and hydroxyΔ5steroid dehydrogenase/3ß and steroidΔisomerase (3ßHSD) protein abundance were examined using Real-Time PCR and Western Blot procedures, respectively. In media collected after incubation, the concentrations of androstenedione (A4) and testosterone (T) were quantified used a RIA. When P4 was added to the culture medium, EMF radiation had suppressive effects on endometrial T release after 2 and 4 h of incubation when the EMF treatment was occurring and increased A4 release after 4 h of incubation with EMF at 120 Hz. When there was no inclusion of P4, release of A4 was decreased after 2 h of EMF treatment at 120 Hz and after 4 h of EMF treatment at 50 and 120 Hz. Progesterone did not have functions that protected the pig endometrium against EMF radiation during the fetal peri-implantation period.

Progesterone dose during synchronization treatment alters luteinizing hormone receptor and steroidogenic enzyme mRNA abundances in granulosa cells of Nellore heifers.

The objective was to investigate effects of progesterone (P4) dose on abundance of luteinizing hormone receptor (LHCGR), aromatase (CYP19A1), 3ß-hydroxysteroid dehydrogenase (HSD3B1), and other steroidogenic mRNA transcripts in granulosa cells from dominant follicles. Nellore heifers were assigned to one of six groups: new, first-use controlled internal drug release device (CIDR1) inserted for 5 days (Large-P4-dose-D5; n = 7) or 6 days (Large-P4-dose-D6; n = 8), prostaglandin (PG)F2α administered on D0 and 1 previously-used CIDR (CIDR3) inserted for 5 days (Small- P4-dose-D5; n = 8) or 6 days (Small-P4-dose-D6; n = 8), CIDR1 inserted on D0 and removed plus PGF2α on D5 (Large-P4-dose-proestrus (PE); n = 7), and CIDR3 and PGF2α on D0 and 1, CIDR3 removed plus PGF2α on D5 (Small-P4-dose-PE; n = 7). Duration of P4 treatment (D5 compared to D6) affected abundances of CYP19A1 mRNA transcripts, with there being greater abundances on D6 than D5 (P ≤ 0.05). Heifers treated with the large dose of P4 had a smaller dominant follicle, less serum and intra-follicular estradiol (E2) concentrations (P ≤ 0.05) and lesser LHCGR, CYP19A1, and HSD3B1 transcript abundances (P ≤ 0.05). Heifers treated to induce PE had a larger follicle diameter (P = 0.09), greater intra-follicular E2 concentrations and larger abundances of CYP19A1 mRNA transcript (P ≤ 0.05) than heifers of the D6 group. Overall, treatment with larger doses of P4 resulted in lesser abundances of LHCGR, HSD3B1, and CYP19A1 mRNA transcripts; thus, potentially leading to development of smaller dominant follicles and lesser E2 concentrations.

Prostate Cancer Foundation Hormone-Sensitive Prostate Cancer Biomarker Working Group Meeting Summary.

Androgen deprivation therapy remains the backbone therapy for the treatment of metastatic hormone-sensitive prostate cancer (mHSPC). In recent years, several treatments, including docetaxel, abiraterone + prednisone, enzalutamide, and apalutamide, have each been shown to demonstrate survival benefit when used upfront along with androgen deprivation therapy. However, treatment selection for an individual patient remains a challenge. There is no high level clinical evidence for treatment selection among these choices based on biological drivers of clinical disease. In August 2020, the Prostate Cancer Foundation convened a working group to meet and discuss biomarkers for hormone-sensitive prostate cancer, the proceedings of which are summarized here. This meeting covered the state of clinical and biological evidence for systemic therapies in the mHSPC space, with emphasis on charting a course for the generation, interrogation, and clinical implementation of biomarkers for treatment selection.

Mutations in 3ß-hydroxysteroid-δ8, δ7-isomerase paradoxically benefit epidermal permeability barrier homeostasis in mice.

Inherited or acquired blockade of distal steps in the cholesterol synthetic pathway results in ichthyosis, due to reduced cholesterol production and/or the accumulation of toxic metabolic precursors, while inhibition of epidermal cholesterol synthesis compromises epidermal permeability barrier homeostasis. We showed here that 3ß-hydroxysteroid-δ8, δ7-isomerase-deficient mice (TD), an analog for CHILD syndrome in humans, exhibited not only lower basal transepidermal water loss rates, but also accelerated permeability barrier recovery despite the lower expression levels of mRNA for epidermal differentiation marker-related proteins and lipid synthetic enzymes. Moreover, TD mice displayed low skin surface pH, paralleled by increased expression levels of mRNA for sodium/hydrogen exchanger 1 (NHE1) and increased antimicrobial peptide expression, compared with wild-type (WT) mice, which may compensate for the decreased differentiation and lipid synthesis. Additionally, in comparison with WT controls, TD mice showed a significant reduction in ear thickness following challenges with either phorbol ester or oxazolone. However, TD mice exhibited growth retardation. Together, these results demonstrate that 3ß-hydroxysteroid-δ8, δ7-isomerase deficiency does not compromise epidermal permeability barrier in mice, suggesting that alterations in epidermal function depend on which step of the cholesterol synthetic pathway is interrupted. But whether these findings in mice could be mirrored in humans remains to be determined.

HSD3B1 variant and androgen-deprivation therapy outcome in prostate cancer.

OBJECTIVE: Rate-limiting enzyme 3b-hydroxysteroid dehydrogenase type 1 (3ßHSD1) encoded by HSD3B1 catalyzes the transition of dehydroepiandrosterone (DHEA) to dihydrotestosterone (DHT). The HSD3B1 (1245C) variant renders 3bHSD1 of resistant to ubiquitination and degradation, leading to a large amount of protein accumulation in the cell. Multiple clinical studies have shown that this mutation was correlated with resistance to androgen-deprivation therapy in prostate cancer. However, the results were not consistent depending on different treatment strategy and in some researches, the number of observed cases was relatively small. METHODS: To determine the effects of HSD3B1 (1245C) variant on resistance to androgen-deprivation therapy in prostate cancer, we performed a meta-analysis of the available literature. Electronic database searches identified appropriately designed studies that detected HSD3B1 in prostate cancer. We conducted a systematic search of studies in the following databases: PubMed, and EMBASE published until August 10, 2020 using the following search terms: (HSD3B1 AND ((((prostate cancer) OR prostatic neoplasm) OR prostatic carcinoma) OR prostatic cancer). RESULTS: Eight researches were included in this research. The result validated that the HSD3B1 (1245C) variant allele was associated with a shorter PFS (HR, 1.97; 95% CI, 1.39-2.79; P = 0.0001) (homozygous wild-type group) in men with prostate cancer when treated with ADT, however, a higher PFS (HR, 0.68; 95% CI, 0.48-0.96; P = 0.03) when treated with ADT and CYP17A1 inhibitor. CONCLUSION: The HSD3B1 (1245C) variant is a predictor of ADT plus CYP17A1 inhibitor response in prostate cancer.

Transcriptome analysis of sheep follicular development during prerecruitment, dominant, and mature stages after FSH superstimulation.

Sheep is usually a monovular animal; superovulation technology is used to increase the number of offspring per individual and shorten generation intervals. To date, mature FSH superstimulatory treatments have been successfully used in sheep breeding, but much remains unknown about genes, pathways, and biological functions involved in follicular development. Therefore, in this study, we performed transcriptome profiling of small follicles (SFs; 2-2.5 mm), medium follicles (MFs; 3.5-4.5 mm), and large follicles (LFs; > 6 mm) in Mongolian ewes after FSH superstimulation. Furthermore, we identified differentially expressed genes and performed Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology enrichment analyses in 3 separate pairwise comparisons. We found that ovarian steroidogenesis was significantly enriched in the SFs versus MFs analysis; the associated genes, cytochrome P450 family 19 (CYP19) and Hydroxy-delta-5-steroid dehydrogenase 3 beta- and steroid delta-isomerase 1 (HSD3B1), were significantly upregulated. Moreover, proline metabolism, glutathione metabolism, and PPAR signaling pathways were significantly enriched in the LFs versus SFs analysis; the associated genes, glutamate-cysteine ligase modifier subunit (GCLM) and cystathionine gamma-lyase (CTH), were significantly upregulated, whereas peroxisome proliferator-activated receptor gamma (PPARγ) was significantly downregulated. In summary, our study provides basic data and possible biological direction to further explore the molecular mechanism of sheep follicular development after FSH superstimulation.