Granulosa cells provide transcriptomic information on ovarian follicle dynamics in the southern white rhinoceros

Granulosa cells provide transcriptomic information on ovarian follicle dynamics in the southern white rhinoceros | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Granulosa cells provide transcriptomic information on ovarian follicle dynamics in the southern white rhinoceros Kristin Klohonatz, Barbara Durrant, Marc-André Sirard, Elena Ruggeri This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4765147/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Much remains unknown about the reproductive physiology of southern white rhinoceros (SWR) and the effect of ovarian stimulation prior to ovum pickup (OPU) have not been fully elucidated. Granulosa cells (GC) provide valuable insight into follicle growth and oocyte maturation status. The goals of this study were to evaluate transcriptomic changes in GC from three stages of follicle development and to identify biomarkers possibly associated with follicular growth and maturation because of ovarian stimulation. GC collected from SWRs following OPU were assigned stages based upon follicle size. Total RNA was isolated, and cDNA libraries were prepared and sequenced on a NovaSeq 6000. All bioinformatics analyses were performed utilizing the Galaxy web platform. Reads were aligned to CerSimCot1.0, and the manual curation was performed with EquCab3.0. Overall, 39,455 transcripts (21,612 genes) were identified across follicle stages, and manual curation yielded a 61% increase in gene identification from the original annotation. Granulosa cells from preovulatory follicles expressed the highest number of unique transcripts. The following seven biomarkers were determined based upon cluster analysis and patterns of expression: COL1A1 , JMY, FBXW11 , NRG1, TMPO , MACIR and COL4A1 . These data can be used to potentially evaluate the effects of different ovarian stimulation protocols on follicle dynamics, improve OPU results, and support conservation efforts in this species. Animal Science Wildlife Biology Bioinformatics Molecular Genetics Animal Physiology southern white rhinoceros OPU ARTs oocytes granulosa cells reproduction endangered species Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Wild rhinoceroses are facing rapid population decline worldwide due to poaching and habitat loss. To preserve these keystone species, assisted reproductive technologies (ARTs) are indispensable. ARTs have been developed and applied to rhinoceros species with variable success rates due to the unique reproductive characteristics of each species ( 1 – 9 ). In general, the success rate of ART for rhinoceros species is lower than that for domestic species ( 10 – 14 ), mainly due to the lack of fundamental reproductive knowledge, challenges in adapting protocols for domestic species, and limited opportunities to collect gametes. In the past decade, considerable conservation effort has been directed toward the rescue of the near-extinct northern white rhinoceros (NWR), with successful application of ARTs and the production of blastocysts in vitro from the last two female NWRs ( 1 , 5 ). The closely related subspecies, the southern white rhinoceros (SWR), has been key to the success of developing these techniques for the NWR ( 2 , 6 , 15 ). Moreover, recent studies have shown that the NWR and SWR populations are distinct but genetically compatible for future efforts at genetic rescue; therefore, the less threatened SWR population is an ideal model for developing new reproductive strategies to help this species ( 15 ). Recent work has explored similarities and differences in the anatomy and reproductive physiology of the rhinoceros and the domestic horse due to their taxonomic proximity ( 16 ). The mare and the rhinoceros both have bicornuate uteri and similar ovarian structures, such as ovarian and follicle sizes. In addition, the mare and rhinoceros share similar hormonal estrous cycle features (estradiol, progesterone and LH)( 16 ). Therefore, based upon their reproductive similarities, the mare is a valuable candidate model for rhinoceros reproductive questions. With the recent success of ovum pick-up (OPU) in the SWR ( 1 , 2 , 5 , 17 ), it is imperative to fully understand the estrous cycle to develop optimal ovarian stimulation protocols. Ovarian ultrasonography has enabled the identification of cycle lengths in the SWR ( 3 , 4 , 18 ), but the mechanisms of follicle growth, oocyte maturation, ovulation, and embryo development are only partly understood ( 2 ). Like in horses and cows, follicle growth begins with increasing FSH levels in response to pulsatile GnRH stimulation ( 16 , 19 ). Once selection of a dominant follicle occurs, FSH levels decrease, and all nondominant follicles begin to undergo atresia. The dominant follicle becomes luteinizing hormone (LH)-dependent and will ultimately ovulate or plateau in the presence, or absence, of progesterone (i.e., the presence or absence of a corpus luteum). The developmental competence of an oocyte in vitro is strongly correlated with follicle development ( 19 ). Therefore, proper ovarian stimulation prior to OPU is crucial for aspirating the oocytes with the highest developmental potential. Ovarian stimulation is routinely performed in horses and cattle prior to OPU ( 20 , 21 ). In horses, although not widely successful, stimulation to retrieve mature oocytes during OPU is performed ( 22 – 24 ). This method of ovarian stimulation in horses utilizes GnRH analogs but results in only one or two preovulatory follicles. An alternative method of OPU for horses is to aspirate immature oocytes for in vitro maturation ( 25 ). Due to the limited opportunities to perform OPU on SWRs, aspirating a single preovulatory follicle is not ideal; therefore, stimulation methods for producing multiple large follicles are necessary. In cattle, superstimulation protocols utilize both FSH and LH, depending on the desired outcome [multiple large follicles (OPU) or timed artificial insemination (ovulation)] ( 26 ). For both methods, a CIDR that secretes progesterone may be used to synchronize estrus prior to stimulation ( 27 ). After CIDR removal, for OPU, multiple FSH injections are given, resulting in an increased number of large follicles. For timed artificial insemination, after CIDR removal, injections of prostaglandin and GnRH are provided to synchronize ovulation ( 26 ). In the SWR, a single GnRH injection successfully induces ovulation when the follicle length reaches 30–35 mm ( 18 ). In previous studies, the administration of the GnRH analog Histrelin increased follicle numbers with alternate day injections prior to OPU ( 1 , 5 ). However, the effects of ovarian stimulation on follicle gene expression in rhinoceros are unknown. Given the influence of the follicle stage on oocyte meiotic competence, it is crucial to understand the gene expression patterns associated with these stages and how they are influenced by ovarian stimulation prior to OPU. With limited access to rhinoceros oocytes, granulosa cells (GC) can provide insight into the physiologic status of follicles and oocytes as they undergo dynamic changes during exogenous ovarian stimulation ( 28 – 30 ). Our previous study identified the transcriptomic signature of a pool of granulosa cells from the SWR, but to date, no studies have evaluated GC across follicle stages ( 2 ). Little has been done in horses to understand gene expression patterns in GC across different follicle stages; hence, cattle are another valuable domestic model for providing comparative transcriptomic data generated from the SWR. A thorough analysis of bovine GC across different follicle stages and correlating gene expression patterns to oocyte and meiotic competence has been completed, and biomarkers for the physiologic status of the oocytes have been identified ( 19 , 28 , 31 , 32 ). It is especially important to determine biomarkers for these stages to develop appropriate ovarian stimulation protocols for ovum pickup (OPU). The aim of this study was to develop a stimulation protocol that results in the largest number of follicles at the appropriate stage for collecting oocytes capable of in vitro meiotic resumption after OPU. The goals of this study were to identify transcriptional changes across follicle development and develop a panel of biomarkers for use in future studies to evaluate transcriptional changes due to different ovarian stimulation protocols. First, we determined the transcriptional signatures of follicles in different stages of development in the SWR. This allowed us to identify gaps in knowledge of the reproductive physiology and ovarian function of the SWR, determine key drivers of follicle development and oocyte maturation in this species, and ultimately develop an in vitro system that best emulates in vivo conditions. Results Part 1. Additions to the NWR Genome The northern white rhinoceros (NWR) genome is newly assembled and annotated ( 15 ) but still contains many gaps, primarily in the annotation. Manual curation of both the sequences identified in the annotation and those identified via de novo sequencing with Stringtie was performed to improve the understanding of the data generated with this study. Manual curation involved searching for the specific sequences identified in the dataset through NCBI BLAST ( 33 ) and comparing the homology to the SWR and the domestic horse (as the genome is highly homologous and well annotated). If a sequence was at least 80% homologous to a gene identified via NCBI, that gene was assigned to that sequence. In this study, improvements were made to provide proper gene nomenclature, assignment of genes to scaffolds and sequences, and identification of genes that were not included in the original annotation. The original annotation included 14,274 transcripts that were expressed in granulosa cells according to this analysis. Throughout the annotation, the protein nomenclature was corrected and changed to gene nomenclature. Through manual curation, 6,763 transcripts were corrected from protein to gene nomenclature, with an overall improvement in nomenclature of 47%. In addition, 70% of the transcript sequences were reassigned to different genes after being searched for gene homology. Through de novo sequencing, an additional 7,429 genes were added to the original annotation, with an increase of 61% more genes from the original annotation (Fig. 1 ). These improvements in the annotation will be incorporated in a follow-up annotation to be made publicly available with the genome. Supplemental Table 1 contains all the sequences identified as a gene as well as their corresponding gene reference. Part 2. Analysis of granulosa cells Prior to analysis, principal component analysis (PCA) was performed to determine the variability of the samples within each follicle stage and between the different stages. Figure 2 shows that within each stage, the replicates had very little variability between the two samples at each follicle stage. In addition, there was clear clustering of each stage, independent of the other stages. A heatmap showed the same trend, with granulosa cells (GC) from preovulatory (P) follicles being in their own expression group and GC from growing (G) and dominant (D) follicles clustering closer together (Fig. 2 ). Across all the samples (n = 6), 39,455 different transcripts were detected, which corresponded to 21,621 genes, with multiple genes expressing multiple transcripts (Supplemental Table 1). The first analysis was to determine the genes that were specific to different follicle stages (growing, dominant, or preovulatory). The growing follicles contained 427 unique transcripts (corresponding to 385 genes), with only 8% of the transcripts identified being unannotated. A total of 812 transcripts were specifically identified in the dominant follicle stage, resulting in a total of 677 genes, and 104 transcripts (13%) were unannotated. Preovulatory follicles contained the largest number of follicle stage-specific transcripts, 945, corresponding to 433 genes. Of the 945 transcripts identified, 484 (51%) were unannotated. Importantly, the most highly expressed genes specific to each follicle type were associated with different transcription factors. The canonical pathways most enriched in growing follicles were involved in the regulation of TP53 activity through phosphorylation, the regulation of TP53 expression and degradation, the signaling of VEGF, the cell cycle checkpoint, and HDR through homologous recombination or single-strand annealing. The genes in preovulatory follicles targeted the following canonical pathways: signaling via the insulin receptor, iron uptake and transport, the iron homeostasis signaling pathway, amino acid regulation of mTORC1, and the adrenergic receptor signaling pathway. The next analyses evaluated the differentially expressed genes (DEG) between each follicle category. The comparison with the least number of differentially expressed transcripts was growing versus dominant, with 855 transcripts (resulting in 802 different genes). The genes annotated with the greatest significant differences were RNPC3 and ZDHHC5 , both more highly expressed in growing follicles. Figure 3 shows the representative volcano plot for this analysis, highlighting the genes with the largest significant differences. Table 1 shows the top 10 genes differentially expressed between growing and dominant follicles. Table 1 Top 10 DEGs when comparing growing (G) vs. dominant (D) follicles. A positive fold change indicates that the gene was more highly expressed in G, and a negative fold change indicates that the gene was more highly expressed in D. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated. Gene Symbol Annotation Assignment Fold Change P-Value RNPC3 Annotated 5.030277375 1.03E-22 ZDHHC5 Annotated 6.697711575 3.29E-19 REV3L Annotated -3.122418916 4.86E-18 MSTRG.16139.2 Unannotated -14.3228372 2.34E-17 CP2DE_g9976.t1 Unannotated -5.49818783 3.22E-17 RYR2 Annotated 4.597912368 7.13E-17 MSTRG.5093.1 Unannotated 3.681205492 1.77E-16 HUWE1 Annotated -5.088883652 2.11E-16 AKAP5 Annotated -14.03465133 2.11E-16 FANCG Annotated -14.39311232 5.43E-16 When comparing growing and preovulatory follicles, the largest number of transcripts (4,581) were significantly differentially expressed, corresponding to 3,980 genes. The genes annotated with the greatest significant differences were CDH11 and TMEM51 (both of which were more highly expressed in growing follicles, as shown in Fig. 3 ). Table 2 lists the top 10 significantly differentially expressed transcripts. Table 2 Top 10 DEGs when comparing growing (G) vs. preovulatory (P) follicles. A positive fold change indicates that the gene was more highly expressed in P, and a negative fold change indicates that the gene was more highly expressed in G. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated. Gene Symbol Annotation Assignment Fold Change P-Value CDH11 Annotated -8.022803609 1.05E-16 TMEM51 Annotated -16.3725893 1.19E-14 B3GNT4 Annotated -16.06868082 5.49E-14 MSTRG.2860.2 Unannotated -11.97253106 1.04E-13 GLCE Annotated 6.57298424 3.46E-13 WDR11 Annotated -5.837110353 4.85E-13 PHC3 Annotated 5.981697786 1.66E-12 KHSRP Annotated -15.14411634 1.90E-12 RPL6 Annotated -15.2562075 4.13E-12 LRP1B Annotated -14.99039555 8.71E-12 The final comparison was dominant versus preovulatory follicles. There were 4,315 transcripts differentially expressed within this comparison (3,704 genes). The annotated genes most differentially expressed were HSPG2 and CFAP52 , both of which are more highly expressed in dominant follicles. Table 3 lists the top 10 differentially expressed transcripts, and Fig. 4 shows the genes with the most differential expression. Supplemental Table 2 details all the differentially expressed transcripts identified in each comparison. Table 3 Top 10 DEGs when comparing dominant (D) vs. preovulatory (P) follicles. A positive fold change indicates that the gene was more highly expressed in P, and a negative fold change indicates that the gene was more highly expressed in D. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated. Gene Symbol Annotation Assignment Fold Change P-Value HSPG2 Annotated -20.063746 9.50707E-23 CFAP52 Annotated -7.7004087 1.68824E-19 MSTRG.2860.2 Unannotated -12.934659 8.26415E-18 KHSRP Annotated -15.24234 3.979E-15 RPL6 Annotated -15.156424 7.45097E-15 TMEM51 Annotated -15.438493 7.68241E-15 SUMO1 Annotated -15.366785 9.58062E-15 MSTRG.9751.1 Unannotated 7.62408891 1.98344E-14 RPL15 Annotated -14.491793 1.02625E-13 PTPRM Annotated -14.556105 9.71413E-14 Although not specific to any follicle type, MACIR was the most highly expressed gene in each follicle category. Table 4 contains the most highly expressed genes in each category, and interestingly, the same genes are present in multiple follicle types. Table 4 The most highly expressed genes in each stage. The presence of multiple genes indicates that different transcripts were among the most highly expressed for that stage. Growing Dominant Pre-Ovulatory MACIR MACIR MACIR MSTRG.15261.1 MSTRG.15261.1 MSTRG.15261.1 COL1A1 COL1A1 COL1A1 SERPINE2 SERPINE2 MSTRG.13580.1 HSPG2 HSPG2 SERPINE2 HSPG2 SLC23A1 KATNIP KATNIP VCAN SLC23A1 MSTRG.13580.1 HSPG2 SERPINE2 SLC23A1 MSTRG.13580.1 VCAN VCAN HSPG2 MSTRG.7397.1 Further evaluation of the DEG in each comparison led to the identification of the most enhanced canonical pathways for each comparison. Table 5 shows the canonical pathways affected by the most transcripts in each category, and the canonical pathways with the greatest fold changes and z scores are shown in Fig. 4 . Table 5 Canonical pathways affected by the most transcripts in each follicle stage. Growing vs. Dominant Growing vs. Preovulatory Dominant vs. Preovulatory RHO GTPase cycle Rho GTPase cycle Selenoamino acid metabolism Ribonucleotide reductase signaling pathway Selenoamino acid metabolism Response of EIF2AK4 (GCN2) to amino acid deficiency Role of tissue factor in cancer Response of EIF2AK4 (GCN2) to amino acid deficiency Rho GTPase cycle Angiopoietin signaling SRP-dependent cotranslational protein targeting to membrane Nonsense-mediated decay (NMD) Production of nitric oxide and reactive oxygen species in macrophages Cell cycle checkpoints Eukaryotic translation termination To understand the expression patterns across follicle types, a cluster analysis was performed, revealing seven unique patterns of expression (Fig. 5 ). Each cluster contained varying numbers of transcripts and differentially expressed transcripts within different comparisons, ranging from 173 transcripts (cluster four) to 1620 transcripts (cluster five). The expression of 173 transcripts in cluster four increased from G to D to P, with the largest number of differentially expressed transcripts occurring between the G and P follicle stages. Pathways associated with this cluster included focal adhesion, the T-cell receptor signaling pathway, the Pi3K-Akt signaling pathway, and the MAPK signaling pathway. This cluster closely resembled the gene expression patterns associated with follicle differentiation. Cluster five contained the largest number of transcripts (1620) whose expression largely decreased from the D to the P follicle stage. The pathways associated with these transcripts include cellular senescence, EGFR tyrosine kinase inhibitor resistance, and endocytosis. Discussion This is the first study evaluating the transcriptomic changes in growing (G), dominant (D) and preovulatory (P) follicles in southern white rhinoceros (SWR). This information is needed to advance assisted reproductive technologies (ARTs) by determining appropriate methods for manipulating ovarian function prior to oocyte collection to obtain oocytes with the highest developmental potential. Through this analysis, we generated candidate gene markers for the physiological status of the follicle at the time of OPU. Although this study was performed on SWR granulosa cells (GC), the SWR genome is incomplete and contains large gaps in annotation. Therefore, the northern white rhinoceros (NWR) genome and currently available annotations are the only resources for generating transcriptomic data on white rhinoceroses, but extensive manual curation is still necessary ( 15 ). In addition, the NWR and SWR are genetically similar at the chromosomal and mitochondrial genome levels, enabling the NWR genome to be used for this analysis ( 15 ). The original annotation for the NWR contained 14,274 different transcripts (12,252 genes) that were expressed in the GC from this dataset. After adding the transcripts from de novo sequencing, a total of 39,455 different transcripts (21,621 genes) were expressed. This was an overall improvement of 119% in transcripts and 61% in genes compared to the original annotation. In addition, it is important to note that through manual curation, multiple transcripts were reassigned to different genes, and transcripts with unknown associations were identified. This approach resulted in a substantial improvement to the available annotation, and these new annotations will be available publicly in the future. Due to their scarcity, reproductive cells are often underrepresented when generating annotations for genomes. Many new genes were annotated in the granulosa cells due to the massive transcriptional changes that occurred within the cells over time and their dynamic behavior. This additional information to the annotation will benefit other researchers in reproductive physiology in this species, and it emphasizes the importance of using reproductive cells, such as granulosa cells, as a resource for conservation research and genome improvement and accuracy. One of our objectives was to identify the genes most highly expressed across the follicle stages. Interestingly, all follicle stages had the same top three most highly expressed genes. These three genes were involved in macrophage immunometabolism regulator ( MACIR ), MSTRG15261. 1, and collagen type I alpha 1 chain ( COL1A1 ). The most highly expressed gene, MACIR , formally known as C5ORF30 , has a well-documented role in modulating the immune response and regulating macrophage function ( 34 , 35 ). In addition, blocking MACIR enhances the expression of genes involved in regulating cell migration, adhesion, and angiogenesis ( 36 ). In cattle GC, MACIR has been shown to be upregulated in follicles undergoing late-stage atresia ( 37 ). This is a potential indicator that the follicles in this study were under stress and already progressed through advanced, irreversible atresia. The second most highly expressed gene, MSTRG15261.1 , was found to be unannotated in the SWR and in the domestic horse. The third most highly expressed gene was COL1A1 , which is the most abundant collagen in the mammalian ovary and is regulated by transforming growth factor beta 1 ( TGF-β1 ) ( 38 ). In addition, COL1A1 promotes the secretion of estrogen while also inhibiting the increase in progesterone, resulting in follicle growth ( 39 ) ( 40 ). Most importantly, COL1A1 plays an important role in oocyte maturation and embryo development through its association with the PI3K/Akt signaling pathway in mice and pigs ( 41 – 44 ). Interestingly, the most highly expressed gene was involved in follicle death, and the third highest was involved in follicle growth. This could indicate that prior to OPU, follicles received mixed signals during ovarian treatment. Cluster analysis was performed to generate gene expression patterns across the follicle stages to identify genes that follow conventional literature for these stages and those that are not. Due to the limited opportunities to perform OPU in SWRs, the current preferred approach is to stimulate the growth of multiple follicles to increase the oocyte recovery rate. Therefore, unstimulated animals were not considered for this elective procedure in our study. Due to the lack of an unstimulated rhinoceros for comparison, the gene expression patterns were compared to published literature generated from unstimulated cows ( 28 ). The two main cluster types we focused on were patterns associated with follicle growth (clusters four and seven) and differentiation (clusters three, five, and six). For cluster three, the gene expression levels were consistent between G and D and increased from D to P; therefore, these genes are associated with follicle differentiation. Of the top 20 genes from this cluster, 13 have known expression patterns across follicle stages ( 28 ); interestingly, five of those genes should be downregulated in P, yet their expression is increased in our dataset. COL4A1 (collagen type IV alpha 1 chain) was found to be differentially expressed in this cluster compared to cattle data ( 28 ). In cattle, COL4A1 is a gene that stays in relatively low abundance throughout the follicle stages, with a decrease in P and an increase only observed in atretic follicles ( 28 ). COL4A1 is directly regulated by forkhead box L2 (FOXL2), a transcription factor specifically expressed in granulosa cells ( 45 ). High levels of FOXL2 decrease COL4A1 ( 45 ), which is required for the proliferation of GC, allowing follicle remodeling, and should remain high during folliculogenesis to allow for proliferation and differentiation ( 46 ). In our data, we observed a significant decrease in FOXL2 expression between D and P, which led to a significant increase in COL4A1 expression between D and P, a sign of late-stage atresia associated with follicular demise ( 31 ). In addition, high levels of collagen in the follicle basement membrane have been linked to an imbalance in gonadotropin stimulation ( 47 ). It is hypothesized that this delay in maturation due to increased collagen may lead to asynchrony during maturation and negative effects on oocyte developmental competence ( 47 ). Therefore, we can hypothesize that high levels of COL4A1 in our SWRs could be caused by incomplete or suboptimal stimulation prior to OPU. Junction mediating and regulatory protein p53 cofactor ( JMY ) is a transcriptional cofactor and an actin nucleation factor ( 48 ) ( 49 ). Actin nucleation occurs during oocyte polarization and affects the microtubule and microfilament cytoskeleton by activating Arp 2/3 (for actin remodeling) and the assembly of filaments in the cytoskeletal spindle ( 48 , 50 , 51 ). If JMY is not present, the spindle fails to migrate to the cortex, and the oocyte arrests with a centrally located spindle, highlighting the crucial role of JMY during peripheral spindle migration during oocyte maturation ( 49 ). In both mouse and porcine oocytes, JMY decreases during oocyte maturation ( 49 , 52 ). This same pattern of expression was observed in cluster five of our dataset, with JMY being highly significant within that cluster. Cluster four included the genes whose expression increased from G to D to P. A gene of interest in this pattern was F-box and WD repeat domain containing 11 ( FBXW11 ) due to its regulation by estradiol, its role in meiotic resumption, and its normal pattern of expression compared to our data. FBXW11 is regulated by many different factors, one of which is helicase for meiosis 1 ( HFM1 ) ( 53 ). Low levels of HFM1 result in high levels of FBXW11 , which corresponds to the pattern of expression we observed for both of these genes; however, according to the literature in mice, HFM1 expression should be high for oocytes to resume meiosis ( 53 ). In addition, FBXW11 is regulated by estradiol, as indicated by high levels of estradiol (as found in late stage developing follicles), which results in low levels of FBXW11 ( 28 , 54 ). The primary function of FBXW11 , formally known as β-TRCP2 , is to degrade cell division cycle 25A ( CDC25A ), which is crucial for the resumption of meiosis (at the appropriate time) ( 55 , 56 ). CDC25A must be downregulated for oocytes to be released from MII arrest ( 57 ), which we observed in the rhinoceros data. Therefore, from the literature, it appears that the timing of the increase in FBXW11 is important in relation to the resumption of oocyte meiosis. NRG1 (neuregulin 1) is a member of the EGF-like factor family that is induced by gonadotropins ( 58 , 59 ). There is a dramatic increase in NRG1 expression due to LH induction, which regulates proper oocyte maturation and developmental competence ( 60 ). Recent literature has shown that NRG1 supplementation during IVM can promote oocyte nuclear maturation and blastocyst development and hatching in cows ( 60 ). During development, amphiregulin ( AREG ) and NRG1 work synergistically ( AREG is an accelerator, and NRG1 is a brake) to ensure that the required temporal changes within the oocyte occur for proper meiotic resumption and competence ( 59 ). It is important that these genes are balanced to ensure that oocytes have high developmental competence. According to our data, NRG1 is part of cluster five and remains consistent from G to D but decreases from D to P, which is the opposite of what the literature states regarding high levels of NRG1 being needed for oocyte competence acquisition ( 28 , 60 ). NRG1 is a novel factor that may impact oocyte maturation; therefore, it should be considered a biomarker in GC for oocyte maturation competence and should be added to in vitro maturation medium. Cluster six was classified as a cluster related to differentiation because there was a decrease in gene expression from G to D and a large increase from D to P. A gene that followed this cluster pattern but, based upon the literature should not be, was thymopoeitin ( TMPO ), also known as lamina-associated polypeptide 2 (LAP2) ( 61 ). According to both bovine and porcine data, TMPO decreases as follicle progression progresses ( 28 , 62 ), which we previously reported from G to D but subsequently deviates with an increase from D to P. A specific region in TMPO [regions 32–36 known as thymopentin ( TP5 )] has been studied extensively due to its role in the activation of the pluripotency marker LIN28a in oocytes ( 63 , 64 ); therefore, high levels of TMPO result in high levels of LIN28a . LIN28a expression decreases until the 8-cell stage of embryo development, when a dramatic increase is then noted ( 63 ). In contrast, LIN28a levels should be low until after resumption of meiosis, and high LIN28a levels in oocytes are associated with MII-arrested oocytes in mice ( 63 , 65 ). High levels of LAP2 have also been associated with DNA damage in oocytes and programmed cell death in granulosa cells in mice ( 61 ). Overall, these data indicate that the P follicles in this dataset may have DNA damage, potentially leading to oocyte arrest at MII. Based upon all the cluster analyses, the largest digression from the literature appears to be in the D to P transition, corresponding to differentiation. The five main drivers of differentiation in cattle are TGF-β1 , tumor protein 53 ( TP53 ), estradiol, hepatocyte nuclear factor 4 alpha ( HNF4A ) and retinoic acid ( RA ) ( 28 ). Our dataset included all but RA , and the expression of each gene significantly increased from D to P. RA has a positive regulatory effect on granulosa cell proliferation and oocyte maturation ( 66 ). RA levels in growing follicles are required for the follicle to reach appropriate final maturation (differentiation) in response to gonadotropins and affect gonadotropin receptor expression in granulosa cells in cows ( 67 ). It can therefore be hypothesized that with this ovarian treatment protocol, follicles are not receiving the proper stimuli to fully differentiate and produce oocytes with the highest developmental competence potential. Based upon these results, it appears that the ovarian treatment protocol used on these animals did not result in GC or oocytes with the highest developmental competence potential. It is important to note that these comparisons are based primarily on data from unstimulated cows. We hypothesize that the gene expression changes we observed could be an outcome of the simulation protocol, but it is unknown whether some of these deviations are due to the unique expression pattern of the rhinoceros from other species. Because we performed oocyte group culture in vitro , we cannot directly correlate the in vivo GC analyzed with the corresponding oocyte maturation outcome; hence, in the future, we will match GC to the in vitro maturation outcome of oocytes by performing individual culture. Overall, for the animals treated with this stimulation protocol, we achieved an in vitro maturation rate of 33%, but none of those oocytes fertilized or developed into embryos. Although a percentage of retrieved oocytes were capable of extruding the first polar body and reaching metaphase II (MII) after in vitro maturation, we hypothesize that unpaired cytoplasmic maturation, erroneous chromosomal segregation, or morphological characteristics of atresia ( 68 ) left the oocyte incapable of fertilization and becoming a zygote. Incomplete cytoplasmic maturation can lead to failed sperm decondensation, erroneous male pronuclear formation, and ultimately unsuccessful fertilization ( 69 ). However, additional studies are needed to characterize the effects of ovarian stimulation on oocyte developmental competence acquisition/failure and cytoplasmic maturation in this species. A balance between the quantity and quality of oocytes obtained during OPU is suggested, and low, frequent doses of GnRH might be a valuable treatment. With this dataset, we generated a list of potential biomarkers of the GC for the physiologic status of the follicle to determine which oocytes had the highest potential for developmental competence. These markers include COL1A1 (follicle growth) ( 41 – 44 ), JMY and FBXW11 (oocyte preparedness to resume meiosis) ( 49 , 53 , 57 ), NRG1 (oocyte maturation potential) ( 60 ), TMPO (DNA damage) ( 61 ), and MACIR and COL4A1 (markers for cell death and atresia) ( 31 , 46 ). However, further testing will be needed to correlate these markers with the status of the follicles we collected during OPU and ultimately with oocyte maturation and cleavage. It is important to note that one limitation in this study is grouping the follicle stages were grouped based upon follicle size at OPU. Nevertheless, for this study, we were not able to obtain follicles at all stages from the same animal due to the unknown status of the estrous cycle prior to ovarian stimulation and OPU. We attempted to overcome this limitation by obtaining technical replicates from different granulosa cell aliquots of the preovulatory follicle. In addition, the animals in this study were not ultrasounded regularly prior to ovarian treatment and OPU; therefore, the trajectory of the follicle prior to OPU was unknown (i.e., was it increasing in size or was it already regressing and undergoing atresia). We attempted to overcome this by evaluating multiple follicles per stage, and in the future, we will attempt to monitor these animals prior to OPU. In conclusion, this was the first study to evaluate the transcriptomic differences in GC between the different stages of follicle development in the SWR, an important chapter in the reproductive knowledge of this species. In this threatened species, much remains unknown regarding ovarian physiology and follicle activity; nevertheless, it is important to continue addressing these reproductive questions to support assisted reproductive technologies for rhinoceros species. Finally, these data provide multiple biomarkers to potentially improve OPU results. This body of work contains novel information on the SWR and enriches the reproductive knowledge of the species, which will support conservation efforts for this species. Materials and Methods Animal management and ovum pickup (OPU) All procedures, experiments, and methods were reviewed and approved by San Diego Zoo Wildlife Alliance’s Institutional Animal Care and Use Committee (IACUC; protocol number 18–018). In addition, guidelines set forth by ARRIVE (Animal Research: Reporting of In Vivo Experiments, https://arriveguidelines.org/arrive-guidelines ) were followed. All methods for this study were performed in accordance with all relevant guidelines and regulations of the IACUC and ARRIVE. Three southern white rhinoceros females (aged 7, 11, and 24 at the time of OPU) were utilized for this study. While under the care of the SDZWA, all rhinoceroses have undergone health assessments to ensure that no reproductive pathologies, tumors, or cysts were present and were continually monitored. Complete information regarding each individual rhinoceros is included in supplemental table 3. All females were stimulated prior to transrectal OPU as previously described ( 2 , 17 ), facilitating the in vivo collection of follicular fluid aspirate from follicles of distinct sizes. Prior to OPU, females received synthetic chlormadinone acetate (CMA) at 30 mg/day for 35 days (days 0–34). The day following CMA withdrawal (Day 35), all animals received 3 mg of deslorelin [gonadotropin-releasing hormone (GnRH) analog] followed by 2 mg of deslorelin every 48 h prior to OPU (Day 37 and 39; 7 mg total) via intramuscular injection. OPU was performed on day 41. All rhinoceroses were anesthetized using a combination of etorphine (2 mcg/kg), butorphanol (20–24 mcg/kg), medetomidine (23–25 mcg/kg), and azaperone (14–16 mcg/kg) administered intramuscularly using a remote drug delivery system. During initial positioning to facilitate intubation, propofol was administered intravenously (0.5 mg/kg). Following OPU procedures, anesthesia reversal was achieved using atipamezole delivered intramuscularly at a 5:1 ratio to medetomidine (116–125 mcg/kg) and naltrexone delivered intramuscularly at a target 50:1 ratio to etorphine (53–61 mcg/kg). Following fecal removal, rinsing, and disinfection, OPU was achieved using a customized, ultrasound-guided probe accompanied by three double-lumen needles. The follicles were ablated, the contents were aspirated, and the samples were rinsed with a warm (37°C) flushing solution (Vigro) containing 12.5 I.U./mL of heparin. During OPU, follicle diameters were measured via trans-rectal ultrasonography. Follicular fluid and granulosa cells were separated based upon follicle stage/size. The follicle stages were growing (11–17 mm), dominant (25–29 mm), and preovulatory (30–34 mm). Granulosa cell collection and RNA isolation/quantification After OPU and the oocyte was collected from the collection fluid, free-floating mural granulosa cells were collected at the follicle stage, pipetted directly into RNAlater (Thermo Fisher Scientific, Waltham, MA) for 24 hours at 4°C and subsequently transferred to -80°C until RNA isolation. Prior to isolation, the granulosa cells were thawed, mixed 1:1 with cold PBS (Sigma Aldrich, St. Louis, MO) and centrifuged at 3,000 × g for 10 min. The supernatant was discarded, and the pellet was resuspended in cold PBS. The samples were centrifuged again at 3,000 × g for 5 min to remove all remnants of RNAlater. Total RNA was isolated from granulosa cells using an Arcturus PicoPure RNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA) per the manufacturer’s instructions. Two biological replicates (one sample from two separate individuals) were used for both growing and dominant follicles. Neither of these two animals had a preovulatory follicle present at the time of OPU, and a third animal provided two technical replicates (the GC were stored in numerous separate aliquots) for the preovulatory stage. In total, six samples were used for this study (n = 2 per follicle type). Granulosa cells were incubated with extraction buffer for 30 min before centrifugation to remove debris and extracellular material. The cell extract was incubated with ethanol, bound to, and washed on preconditioned purification columns. Total RNA was recovered into elution buffer and quantified using a Qubit 4 Fluorometer (Thermo Fischer Scientific, Waltham, MA. After quantification, the samples were evaluated on a 4150 TapeStation System (Agilent Technologies, Santa Clara, CA) to determine RNA integrity number (RIN) values. Only samples with a RIN greater than 6.0 were used for RNA-Seq analysis. RNA sequencing Library preparation and RNA sequencing were performed at the University of California San Diego Institute for Genomic Medicine Center. Following the manufacturer’s protocol, RNA-sequencing (cDNA) libraries were prepared using the Illumina TruSeq Stranded Total RNA Prep with Ribo-Zero Plus. A library was prepared for each sample using 50 ng of total RNA. Briefly, RNA was rRNA depleted and fragmented, first- and second-strand cDNA synthesis was performed, and adapters were ligated. Each sample was amplified with a specific barcoded PCR primer for sample identification purposes. The prepared libraries were sequenced as 100 base pair paired end reads on an Illumina NovaSeq 6000 (Illumina, San Diego, CA). The raw data files were uploaded to the Gene Expression Omnibus under accession number GSE261038. Bioinformatic and statistical analysis Bioinformatic analysis was performed on the Galaxy web platform using the public server usegalaxy.org ( 70 ). Initial assessment of sequence quality was performed by FastQC, and the sequences were aggregated for comparison by MultiQC ( 71 , 72 ). The adapter sequences and reads with quality scores below the quality threshold of 25 and with fewer than 20 base pairs were removed from analysis via Trimmomatic ( 73 ). HISAT2 was used to align reads to the northern white rhinoceros (NWR) genome CerSimCot1.0 (GCA_021442165.1)( 15 , 74 ). The southern white rhinoceros genome (CerSimSim1.0; GCA_000283155.1) was not used for this analysis due to large gaps in the sequencing and extremely limited annotation. Transcript assembly and quantification for both annotated and unannotated transcripts were initially performed with Stringtie, and each animal sample was analyzed individually ( 75 ). The resulting transcripts for individual samples were merged to create a final annotation file representing the intersection and union of all samples ( 75 ). Stringtie was subsequently applied to the complete merged annotation to generate transcript read counts. Transcript counts were exported from Galaxy and imported into RStudio for statistical analysis ( 76 ). Transcript read counts were analyzed between each follicle type utilizing DESeq2 within R ( 77 ). The average expression value for replicates within a follicle stage was ≥ 3, except in the cases where the expression value was 0 for the replicates within a follicle type. When not expressed in a follicle type, the expression value for both replicates was 0. The data were normalized internally using the median ratio method of DESeq2. Benjamini Hochberg false discovery rate adjustment was used. Significance was assessed at P ≤ 0.05. Differentially expressed transcripts were grouped into different cluster patterns according to their expression profiles among the different developmental stages using the Mfuzz Bioconductor package ( 78 ). Pathway analysis and visualization were completed using Qiagen Ingenuity Pathway Analysis ( 79 ). Manual curation was required to assign sequences identified to corresponding genes and to ensure that the HGNC/VGNC naming conventions were maintained. Declarations Competing Interests The authors have no competing interests to declare. Funding Sources This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Author Contributions E.R. and B.D. collected the samples. E.R. processed and prepared the samples for RNA sequencing. E.R. and K.K. performed the bioinformatic analyses and wrote the manuscript. M.A.S. aided in comparative data analysis and interpretation. B.D. aided in the writing and revision of the manuscript. Acknowledgments We thank Dr. Candace Williams at the San Diego Zoo Wildlife Alliance and Dr. Ahmed Gad and Nico Menjivar at Colorado State University for aiding in the data analysis and visualization. In addition, we would like to acknowledge the Conservation Science Wildlife Health team, the Veterinary Services, and the Rhino Rescue Center team at San Diego Zoo Wildlife Alliance for aiding in animal management, care, and sample collection. Data Availability Statement The datasets generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO) under accession number GSE261038. References Hildebrandt TB, Holtze S, Colleoni S, Hermes R, Stejskal J, Lekolool I et al (2023) In vitro fertilization (IVF) program in white rhinoceros. Reproduction (Cambridge, England). :REP-23-0087 Ruggeri E, Klohonatz K, Sirard M-A, Durrant B, Coleman S (2023) Genomic insights into southern white rhinoceros (Ceratotherium simum simum) reproduction: revealing granulosa cell gene expression. Theriogenology Wild. :100055 Pennington PM, Marshall KL, Capiro JM, Howard L, Durrant BS (2020) Pregnancies following long luteal phases in southern white rhinoceros (Ceratotherium simum simum). Zoo Biol 39(2):141–144 Pennington PM, Durrant BS (2019) Assisted reproductive technologies in captive rhinoceroses. Mammal Rev 49(1):1–15 Hildebrandt TB, Hermes R, Colleoni S, Diecke S, Holtze S, Renfree MB et al (2018) Embryos and embryonic stem cells from the white rhinoceros. Nat Commun 9(1):2589 Saragusty J, Diecke S, Drukker M, Durrant B, Friedrich Ben-Nun I, Galli C et al (2016) Rewinding the process of mammalian extinction. Zoo Biol 35(4):280–292 Roth TL (2006) A review of the reproductive physiology of rhinoceros species in captivity. Int Zoo Yearbook 40(1):130–143 Patton ML, Swaisgood RR, Czekala NM, White AM, Fetter GA, Montagne JP et al (1999) Reproductive cycle length and pregnancy in the southern white rhinoceros (Ceratotherium simum simum) as determined by fecal pregnane analysis and observations of mating behavior. Zoo Biol 18(2):111–127 Roth TL (2023) That was then, this is now–over two decades of progress in rhinoceros reproductive science and technology. Theriogenology Wild. :100065 Ferre LB, Kjelland ME, Strobech LB, Hyttel P, Mermillod P, Ross PJ (2020) Review: Recent advances in bovine in vitro embryo production: reproductive biotechnology history and methods. Animal 14(5):991–1004 Chang EM, Song HS, Lee DR, Lee WS, Yoon TK (2014) In vitro maturation of human oocytes: Its role in infertility treatment and new possibilities. Clin Exp Reprod Med 41(2):41–46 Thomas MR, Sparks AE, Ryan GL, Van Voorhis BJ (2010) Clinical predictors of human blastocyst formation and pregnancy after extended embryo culture and transfer. Fertil Steril 94(2):543–548 Hinrichs K (2010) The equine oocyte: factors affecting meiotic and developmental competence. Mol Reprod Dev 77(8):651–661 Galli C, Colleoni S, Duchi R, Lagutina I, Lazzari G (2007) Developmental competence of equine oocytes and embryos obtained by in vitro procedures ranging from in vitro maturation and ICSI to embryo culture, cryopreservation and somatic cell nuclear transfer. Anim Reprod Sci 98(1–2):39–55 Wang G, Brändl B, Rohrandt C, Hong K, Pang A, Lee J et al Chromosome-level genome assembly of the functionally extinct northern white rhinoceros (Ceratotherium simum cottoni). bioRxiv. 2021:2021.12.11.472206. Meuffels-Barkas J, Wilsher S, Allen WRT, Ververs C, Lueders I (2023) Comparative reproduction of the female horse, elephant and rhinoceros: implications for advancing Assisted Reproductive Technologies (ART). Reprod Fertil. ;4(3) Ruggeri E, Young C, Ravida N, Sirard MA, Krisher R, de la Rey M et al (2022) Glucose consumption and gene expression in granulosa cells collected before and after in vitro oocyte maturation in the southern white rhinoceros (Ceratotherium simum simum). Reprod Fertil Dev Pennington PM, Marshall KL, Capiro JM, Felton RG, Durrant BS (2019) Ovulation induction in anovulatory southern white rhinoceros (Ceratotherium simum simum) without altrenogest. Conserv Physiol 7(1):coz033 Girard A, Dufort I, Douville G, Sirard MA (2015) Global gene expression in granulosa cells of growing, plateau and atretic dominant follicles in cattle. Reprod Biol Endocrinol 13:17 Squires E (2020) Current Reproductive Technologies Impacting Equine Embryo Production. J Equine Vet Sci 89:102981 De Roover R, Feugang JM, Bols PE, Genicot G, Hanzen C (2008) Effects of ovum pick-up frequency and FSH stimulation: a retrospective study on seven years of beef cattle in vitro embryo production. Reprod Domest Anim 43(2):239–245 Carnevale EM, Sessions DR (2012) In vitro production of equine embryos. J Equine Vet Sci 32(7):367–371 Carnevale EM, Stokes J, Squires EL, Campos-Chillon LF, Altermatt J, Suh TK (2007) Clinical use of intracytoplasmic sperm injection in horses. Proc Am Assoc Equine Pract 53:560 Carnevale EM, Maclellan LJ (2006) Collection, evaluation, and use of oocytes in equine assisted reproduction. Vet Clin North Am Equine Pract 22(3):843–856 Benammar A, Derisoud E, Vialard F, Palmer E, Ayoubi JM, Poulain M et al (2021) The Mare: A Pertinent Model for Human Assisted Reproductive Technologies? Animals (Basel). ;11(8) Bo GA, Mapletoft RJ (2020) Superstimulation of ovarian follicles in cattle: Gonadotropin treatment protocols and FSH profiles. Theriogenology 150:353–359 Mantovani AP, Reis E, Gacek F, Bó G, Binelli M, Baruselli PS (2018) Prolonged use of a progesterone-releasing intravaginal device (CIDR®) for induction of persistent follicles in bovine embryo recipients. Anim Reprod (AR) 2(4):272–277 Khan DR, Fournier E, Dufort I, Richard FJ, Singh J, Sirard MA (2016) Meta-analysis of gene expression profiles in granulosa cells during folliculogenesis. Reproduction 151(6):R103–R110 Cheng J, Huang J, Yuan S, Zhou S, Yan W, Shen W et al (2017) Circular RNA expression profiling of human granulosa cells during maternal aging reveals novel transcripts associated with assisted reproductive technology outcomes. PLoS ONE 12(6):e0177888 Li Z, Wang J, Zhao Y, Ma D, Zhao M, Li N et al (2021) scRNA-seq of ovarian follicle granulosa cells from different fertility goats reveals distinct expression patterns. Reprod Domest Anim 56(5):801–811 Douville G, Sirard MA (2014) Changes in granulosa cells gene expression associated with growth, plateau and atretic phases in medium bovine follicles. J Ovarian Res 7:50 Nivet AL, Vigneault C, Blondin P, Sirard MA (2013) Changes in granulosa cells' gene expression associated with increased oocyte competence in bovine. Reproduction 145(6):555–565 Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K et al (2009) BLAST+: architecture and applications. BMC Bioinformatics 10:421 Dorris ER, Tazzyman SJ, Moylett J, Ramamoorthi N, Hackney J, Townsend M et al (2019) The Autoimmune Susceptibility Gene C5orf30 Regulates Macrophage-Mediated Resolution of Inflammation. J Immunol 202(4):1069–1078 Muthana M, Hawtree S, Wilshaw A, Linehan E, Roberts H, Khetan S et al (2015) C5orf30 is a negative regulator of tissue damage in rheumatoid arthritis. Proceedings of the National Academy of Sciences. ;112(37):11618-23 Serwin K, Kozak M, Mazurek-Mochol M, Piotrowska K, Czerewaty M, Pawlik A (2023) Human macrophage immunometabolism regulator (MACIR) in patients with periodontitis. Immunobiology 228(6):152760 Hatzirodos N, Hummitzsch K, Irving-Rodgers HF, Harland ML, Morris SE, Rodgers RJ (2014) Transcriptome profiling of granulosa cells from bovine ovarian follicles during atresia. BMC Genomics 15:40 Li H, Chang HM, Shi Z, Leung PCK (2020) The p38 signaling pathway mediates the TGF-beta1-induced increase in type I collagen deposition in human granulosa cells. FASEB J 34(11):15591–15604 Huet C, Monget P, Pisselet C, Monniaux D (1997) Changes in extracellular matrix components and steroidogenic enzymes during growth and atresia of antral ovarian follicles in the sheep. Biol Reprod 56(4):1025–1034 Huet C, Pisselet C, Mandon-Pepin B, Monget P, Monniaux D (2001) Extracellular matrix regulates ovine granulosa cell survival, proliferation and steroidogenesis: relationships between cell shape and function. J Endocrinol 169(2):347–360 Shimada M, Ito J, Yamashita Y, Okazaki T, Isobe N (2003) Phosphatidylinositol 3-kinase in cumulus cells is responsible for both suppression of spontaneous maturation and induction of gonadotropin-stimulated maturation of porcine oocytes. J Endocrinol 179(1):25–34 Hoshino Y, Yokoo M, Yoshida N, Sasada H, Matsumoto H, Sato E (2004) Phosphatidylinositol 3-kinase and Akt participate in the FSH-induced meiotic maturation of mouse oocytes. Mol Reprod Dev 69(1):77–86 Hall SE, Upton RM, McLaughlin EA, Sutherland JM (2018) Phosphoinositide 3-kinase/protein kinase B (PI3K/AKT) and Janus kinase/signal transducer and activator of transcription (JAK/STAT) follicular signalling is conserved in the mare ovary. Reprod Fertility Dev 30(4):624–633 Riley JK, Carayannopoulos MO, Wyman AH, Chi M, Ratajczak CK, Moley KH (2005) The PI3K/Akt pathway is present and functional in the preimplantation mouse embryo. Dev Biol 284(2):377–386 Du H, Guo Y, Wu X, Gong Y (2022) FOXL2 regulates the expression of the Col4a1 collagen gene in chicken granulosa cells. Mol Reprod Dev 89(2):95–103 Marongiu M, Deiana M, Marcia L, Sbardellati A, Asunis I, Meloni A et al (2016) Novel action of FOXL2 as mediator of Col1a2 gene autoregulation. Dev Biol 416(1):200–211 Adriaenssens T, Mazoyer C, Segers I, Wathlet S, Smitz J (2009) Differences in collagen expression in cumulus cells after exposure to highly purified menotropin or recombinant follicle-stimulating hormone in a mouse follicle culture model. Biol Reprod 80(5):1015–1025 Zuchero JB, Coutts AS, Quinlan ME, Thangue NBL, Mullins RD (2009) p53-cofactor JMY is a multifunctional actin nucleation factor. Nat Cell Biol 11(4):451–459 Sun S-C, Sun Q-Y, Kim N-H (2011) JMY is required for asymmetric division and cytokinesis in mouse oocytes. MHR: Basic Sci reproductive Med 17(5):296–304 Brunet S, Verlhac MH (2011) Positioning to get out of meiosis: the asymmetry of division. Hum Reprod Update 17(1):68–75 Namgoong S, Kim NH (2016) Roles of actin binding proteins in mammalian oocyte maturation and beyond. Cell Cycle 15(14):1830–1843 Lin Z, Xu Y-N, Namgoong S, Kim N-H (2014) JMY functions as actin nucleation-promoting factor and mediator for p53-mediated DNA damage in porcine oocytes. PLoS ONE 9(10):e109385 Pu D, Zhong C, Wang H, Yuan X, He Y, Cong J et al (2023) Regulation of FUS ubiquitination and localization by HFM1 is essential for oocyte meiosis prophase I progression in mice Xie J, Jin Y, Wang G (2019) The role of SCF ubiquitin-ligase complex at the beginning of life. Reproductive Biology Endocrinol 17(1):1–9 Silverman JS, Skaar JR, Pagano M (2012) SCF ubiquitin ligases in the maintenance of genome stability. Trends Biochem Sci 37(2):66–73 Isoda M, Kanemori Y, Nakajo N, Uchida S, Yamashita K, Ueno H et al (2009) The extracellular signal-regulated kinase-mitogen-activated protein kinase pathway phosphorylates and targets Cdc25A for SCF beta-TrCP-dependent degradation for cell cycle arrest. Mol Biol Cell 20(8):2186–2195 Oh JS, Susor A, Schindler K, Schultz RM, Conti M (2013) Cdc25A activity is required for the metaphase II arrest in mouse oocytes. J Cell Sci 126(5):1081–1085 Chowdhury I, Branch A, Mehrabi S, Ford BD, Thompson WE (2017) Gonadotropin-Dependent Neuregulin-1 Signaling Regulates Female Rat Ovarian Granulosa Cell Survival. Endocrinology 158(10):3647–3660 Noma N, Kawashima I, Fan HY, Fujita Y, Kawai T, Tomoda Y et al (2011) LH-induced neuregulin 1 (NRG1) type III transcripts control granulosa cell differentiation and oocyte maturation. Mol Endocrinol 25(1):104–116 Dellaqua TT, Vigaro RA, Janini LCZ, Dal Canto M, Renzini MM, Lodde V et al (2023) Neuregulin 1 (NRG1) modulates oocyte nuclear maturation during IVM and improves post-IVF embryo development. Theriogenology 195:209–216 Miao X, Guo R, Williams A, Lee C, Ma J, Wang PJ et al (2023) Replication Protein A1 is essential for DNA damage repair during mammalian oogenesis. bioRxiv Kulus J, Kulus M, Kranc W, Jopek K, Zdun M, Jozkowiak M et al (2021) Transcriptomic Profile of New Gene Markers Encoding Proteins Responsible for Structure of Porcine Ovarian Granulosa Cells. Biology (Basel). ;10(11) Vogt EJ, Meglicki M, Hartung KI, Borsuk E, Behr R (2012) Importance of the pluripotency factor LIN28 in the mammalian nucleolus during early embryonic development. Development 139(24):4514–4523 Liu T, Jing F, Huang P, Geng Z, Xu J, Li J et al (2021) Thymopentin alleviates premature ovarian failure in mice by activating YY2/Lin28A and inhibiting the expression of let-7 family microRNAs. Cell Prolif 54(8):e13089 Flemr M, Moravec M, Libova V, Sedlacek R, Svoboda P (2014) Lin28a is dormant, functional, and dispensable during mouse oocyte-to-embryo transition. Biol Reprod 90(6):131 Kawai T, Yanaka N, Richards JS, Shimada M (2016) De Novo-Synthesized Retinoic Acid in Ovarian Antral Follicles Enhances FSH-Mediated Ovarian Follicular Cell Differentiation and Female Fertility. Endocrinology 157(5):2160–2172 Damdimopoulou P, Chiang C, Flaws JA (2019) Retinoic acid signaling in ovarian folliculogenesis and steroidogenesis. Reprod Toxicol 87:32–41 Hinrichs K (1991) The relationship of follicle atresia to follicle size, oocyte recovery rate on aspiration, and oocyte morphology in the mare. Theriogenology 36(2):157–168 Fulka J Jr., First NL, Moor RM (1998) Nuclear and cytoplasmic determinants involved in the regulation of mammalian oocyte maturation. Mol Hum Reprod 4(1):41–49 Afgan E, Baker D, Batut B, van den Beek M, Bouvier D, Cech M et al (2018) The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res 46(W1):W537–W44 Andrews S (2010) FastQC: a quality control tool for high throughput sequence data Ewels P, Magnusson M, Lundin S, Kaller M (2016) MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32(19):3047–3048 Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30(15):2114–2120 Kim D, Paggi JM, Park C, Bennett C, Salzberg SL (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37(8):907–915 Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol 33(3):290–295 R Core Team R (2013) R: A language and environment for statistical computing Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15(12):550 Kumar L (2007) Mfuzz: a software package for soft clustering of microarray data. Bioinformation 2(1):5–7 Kramer A, Green J, Pollard J Jr., Tugendreich S (2014) Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 30(4):523–530 Additional Declarations The authors declare no competing interests. Supplementary Files SupplementalTable1.pdf All the sequences identified as a gene as well as their corresponding gene reference. SupplementalTable2.pdf All the differentially expressed transcripts identified in each comparison SupplementalTable3.pdf All information on the animal is based upon the time of OPU Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4765147","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":328915008,"identity":"dd8ef4ce-e0a4-49aa-b4fa-028b7044d1d0","order_by":0,"name":"Kristin Klohonatz","email":"","orcid":"https://orcid.org/0000-0002-9790-1950","institution":"University of Pennsylvania","correspondingAuthor":false,"prefix":"","firstName":"Kristin","middleName":"","lastName":"Klohonatz","suffix":""},{"id":328915009,"identity":"7d9b4efa-a0dd-4908-aaed-03cab9931b68","order_by":1,"name":"Barbara Durrant","email":"","orcid":"https://orcid.org/0000-0001-6019-9865","institution":"San Diego Zoo Wildlife Alliance","correspondingAuthor":false,"prefix":"","firstName":"Barbara","middleName":"","lastName":"Durrant","suffix":""},{"id":328915010,"identity":"44c3cae8-fd64-4744-9984-5cc196d59029","order_by":2,"name":"Marc-André Sirard","email":"","orcid":"https://orcid.org/0000-0001-8667-6682","institution":"Université Laval","correspondingAuthor":false,"prefix":"","firstName":"Marc-André","middleName":"","lastName":"Sirard","suffix":""},{"id":328915011,"identity":"da08be30-d817-44fd-8603-89a94ef80f89","order_by":3,"name":"Elena Ruggeri","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA2ElEQVRIiWNgGAWjYBADHn4GBgMStUg2kKqFweAAsVr4pc8Yv/xScU/G+PzhzZ95ahjkzPsX4Nci2ZdjZi1zppjH7EZagTHPMQZjmRsPCLjnDI+ZsWRbAlALj0FyDhtD4gyJA/i12IO1/EvgMe4/Y3A45x8RWgx4eIwffmxI4DFgyDFszm0DauFvwK9F4gxbGTPDsQQeiRtpxcx/+ySMJSTw62Dg72He/PFHTYI9f//hzR9nfLORk+An4DAgYJPmQbIViBIIamH++APVYsK2jIJRMApGwcgCALeNPkLyqoimAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0003-3543-1903","institution":"San Diego Zoo Wildlife Alliance","correspondingAuthor":true,"prefix":"","firstName":"Elena","middleName":"","lastName":"Ruggeri","suffix":""}],"badges":[],"createdAt":"2024-07-18 22:07:52","currentVersionCode":1,"declarations":{"humanSubjects":false,"vertebrateSubjects":true,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":true},"doi":"10.21203/rs.3.rs-4765147/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4765147/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":60833007,"identity":"ec7317d7-9329-40e0-8d60-725ac1c61701","added_by":"auto","created_at":"2024-07-22 15:16:34","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":5713,"visible":true,"origin":"","legend":"\u003cp\u003eImprovement of the SWR annotation before and after annotation.\u003c/p\u003e","description":"","filename":"placeholderimage.png","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/4bb361120f13a2c523c70503.png"},{"id":60833009,"identity":"6741fb91-2509-420c-a87c-fe1cdc02203a","added_by":"auto","created_at":"2024-07-22 15:16:34","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":156754,"visible":true,"origin":"","legend":"\u003cp\u003eAnalysis of data grouping by follicle type and variability using principal component analysis (PCA) and heatmapping.\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/84d4bc62a882cad9478e7fe5.png"},{"id":60833965,"identity":"33c1f2e0-18b6-4a7f-b83f-13a335742025","added_by":"auto","created_at":"2024-07-22 15:32:34","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":236701,"visible":true,"origin":"","legend":"\u003cp\u003eVolcano plots representing the three primary comparisons: growing (G) vs. dominant (D), D vs. preovulatory (P), and G vs. P. The genes marked have the greatest significant differences between stages. In the comparison, the red dots are more highly expressed in the first stage, and the green dots are more highly expressed in the second stage.\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/68abc5824f9d31509d17556d.png"},{"id":60833010,"identity":"220a1196-ee41-49dc-9be2-e90f8ea9a2b9","added_by":"auto","created_at":"2024-07-22 15:16:34","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":707725,"visible":true,"origin":"","legend":"\u003cp\u003eTop canonical pathways for each comparison. The threshold was set to -log(p-value)\u0026gt;1.3, but the z score thresholds were varied due to the large number of canonical pathways. A.) Represents growing vs. dominant follicles with an absolute z score ³ 2. B.) Represents growing vs. preovulatory follicles with an absolute z score ³ 5. C.) Represents dominant vs. preovulatory follicles with an absolute z score ³ 5.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/5816cff6a5f11e371434d939.png"},{"id":60833438,"identity":"97205c20-6407-4f4c-89f4-1d22dd3f9673","added_by":"auto","created_at":"2024-07-22 15:24:34","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":394547,"visible":true,"origin":"","legend":"\u003cp\u003eOverall expression patterns of the seven different clusters. These clusters were determined by the significantly differentially expressed genes between each follicle stage. Each stage is represented as growing (G), dominant (D), or preovulatory (P).\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/25a55ee696f2fe1c016e7374.png"},{"id":60834476,"identity":"e41fc3dd-5c60-49cb-9804-5a3f810a7469","added_by":"auto","created_at":"2024-07-22 15:40:36","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2012502,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/783ff653-2809-4380-b8d0-2be8635cc566.pdf"},{"id":60833013,"identity":"1e9db1a0-bd79-4e48-9541-0877cdb6e5b2","added_by":"auto","created_at":"2024-07-22 15:16:34","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":3761790,"visible":true,"origin":"","legend":"\u003cp\u003eAll the sequences identified as a gene as well as their corresponding gene reference.\u003c/p\u003e","description":"","filename":"SupplementalTable1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/f09f44d6ae747810f4e212cd.pdf"},{"id":60833015,"identity":"c90ce985-ad78-4de5-b953-ef1efc589de2","added_by":"auto","created_at":"2024-07-22 15:16:35","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":3879158,"visible":true,"origin":"","legend":"\u003cp\u003eAll the differentially expressed transcripts identified in each comparison\u003c/p\u003e","description":"","filename":"SupplementalTable2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/adc292958256fae17072cc91.pdf"},{"id":60833008,"identity":"d13b2ada-0a7e-48da-b26a-746999f5177f","added_by":"auto","created_at":"2024-07-22 15:16:34","extension":"pdf","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":15970,"visible":true,"origin":"","legend":"\u003cp\u003eAll information on the animal is based upon the time of OPU\u003c/p\u003e","description":"","filename":"SupplementalTable3.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4765147/v1/74b374cf76d0c0e22c99b0e6.pdf"}],"financialInterests":"The authors declare no competing interests.","formattedTitle":"\u003cp\u003eGranulosa cells provide transcriptomic information on ovarian follicle dynamics in the southern white rhinoceros\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003eWild rhinoceroses are facing rapid population decline worldwide due to poaching and habitat loss. To preserve these keystone species, assisted reproductive technologies (ARTs) are indispensable. ARTs have been developed and applied to rhinoceros species with variable success rates due to the unique reproductive characteristics of each species (\u003cspan additionalcitationids=\"CR2 CR3 CR4 CR5 CR6 CR7 CR8\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e). In general, the success rate of ART for rhinoceros species is lower than that for domestic species (\u003cspan additionalcitationids=\"CR11 CR12 CR13\" citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e), mainly due to the lack of fundamental reproductive knowledge, challenges in adapting protocols for domestic species, and limited opportunities to collect gametes. In the past decade, considerable conservation effort has been directed toward the rescue of the near-extinct northern white rhinoceros (NWR), with successful application of ARTs and the production of blastocysts \u003cem\u003ein vitro\u003c/em\u003e from the last two female NWRs (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e). The closely related subspecies, the southern white rhinoceros (SWR), has been key to the success of developing these techniques for the NWR (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e). Moreover, recent studies have shown that the NWR and SWR populations are distinct but genetically compatible for future efforts at genetic rescue; therefore, the less threatened SWR population is an ideal model for developing new reproductive strategies to help this species (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eRecent work has explored similarities and differences in the anatomy and reproductive physiology of the rhinoceros and the domestic horse due to their taxonomic proximity (\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e). The mare and the rhinoceros both have bicornuate uteri and similar ovarian structures, such as ovarian and follicle sizes. In addition, the mare and rhinoceros share similar hormonal estrous cycle features (estradiol, progesterone and LH)(\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e). Therefore, based upon their reproductive similarities, the mare is a valuable candidate model for rhinoceros reproductive questions.\u003c/p\u003e \u003cp\u003eWith the recent success of ovum pick-up (OPU) in the SWR (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e, \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e), it is imperative to fully understand the estrous cycle to develop optimal ovarian stimulation protocols. Ovarian ultrasonography has enabled the identification of cycle lengths in the SWR (\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e), but the mechanisms of follicle growth, oocyte maturation, ovulation, and embryo development are only partly understood (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e). Like in horses and cows, follicle growth begins with increasing FSH levels in response to pulsatile GnRH stimulation (\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e, \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e). Once selection of a dominant follicle occurs, FSH levels decrease, and all nondominant follicles begin to undergo atresia. The dominant follicle becomes luteinizing hormone (LH)-dependent and will ultimately ovulate or plateau in the presence, or absence, of progesterone (i.e., the presence or absence of a corpus luteum). The developmental competence of an oocyte \u003cem\u003ein vitro\u003c/em\u003e is strongly correlated with follicle development (\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e). Therefore, proper ovarian stimulation prior to OPU is crucial for aspirating the oocytes with the highest developmental potential.\u003c/p\u003e \u003cp\u003eOvarian stimulation is routinely performed in horses and cattle prior to OPU (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e, \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e). In horses, although not widely successful, stimulation to retrieve mature oocytes during OPU is performed (\u003cspan additionalcitationids=\"CR23\" citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e). This method of ovarian stimulation in horses utilizes GnRH analogs but results in only one or two preovulatory follicles. An alternative method of OPU for horses is to aspirate immature oocytes for \u003cem\u003ein vitro\u003c/em\u003e maturation (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e). Due to the limited opportunities to perform OPU on SWRs, aspirating a single preovulatory follicle is not ideal; therefore, stimulation methods for producing multiple large follicles are necessary. In cattle, superstimulation protocols utilize both FSH and LH, depending on the desired outcome [multiple large follicles (OPU) or timed artificial insemination (ovulation)] (\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e). For both methods, a CIDR that secretes progesterone may be used to synchronize estrus prior to stimulation (\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e). After CIDR removal, for OPU, multiple FSH injections are given, resulting in an increased number of large follicles. For timed artificial insemination, after CIDR removal, injections of prostaglandin and GnRH are provided to synchronize ovulation (\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn the SWR, a single GnRH injection successfully induces ovulation when the follicle length reaches 30\u0026ndash;35 mm (\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e). In previous studies, the administration of the GnRH analog Histrelin increased follicle numbers with alternate day injections prior to OPU (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e). However, the effects of ovarian stimulation on follicle gene expression in rhinoceros are unknown. Given the influence of the follicle stage on oocyte meiotic competence, it is crucial to understand the gene expression patterns associated with these stages and how they are influenced by ovarian stimulation prior to OPU.\u003c/p\u003e \u003cp\u003eWith limited access to rhinoceros oocytes, granulosa cells (GC) can provide insight into the physiologic status of follicles and oocytes as they undergo dynamic changes during exogenous ovarian stimulation (\u003cspan additionalcitationids=\"CR29\" citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e). Our previous study identified the transcriptomic signature of a pool of granulosa cells from the SWR, but to date, no studies have evaluated GC across follicle stages (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e). Little has been done in horses to understand gene expression patterns in GC across different follicle stages; hence, cattle are another valuable domestic model for providing comparative transcriptomic data generated from the SWR. A thorough analysis of bovine GC across different follicle stages and correlating gene expression patterns to oocyte and meiotic competence has been completed, and biomarkers for the physiologic status of the oocytes have been identified (\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e, \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e, \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIt is especially important to determine biomarkers for these stages to develop appropriate ovarian stimulation protocols for ovum pickup (OPU). The aim of this study was to develop a stimulation protocol that results in the largest number of follicles at the appropriate stage for collecting oocytes capable of \u003cem\u003ein vitro\u003c/em\u003e meiotic resumption after OPU. The goals of this study were to identify transcriptional changes across follicle development and develop a panel of biomarkers for use in future studies to evaluate transcriptional changes due to different ovarian stimulation protocols. First, we determined the transcriptional signatures of follicles in different stages of development in the SWR. This allowed us to identify gaps in knowledge of the reproductive physiology and ovarian function of the SWR, determine key drivers of follicle development and oocyte maturation in this species, and ultimately develop an \u003cem\u003ein vitro\u003c/em\u003e system that best emulates \u003cem\u003ein vivo\u003c/em\u003e conditions.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePart 1. Additions to the NWR Genome\u003c/h2\u003e \u003cp\u003eThe northern white rhinoceros (NWR) genome is newly assembled and annotated (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e) but still contains many gaps, primarily in the annotation. Manual curation of both the sequences identified in the annotation and those identified via de novo sequencing with Stringtie was performed to improve the understanding of the data generated with this study. Manual curation involved searching for the specific sequences identified in the dataset through NCBI BLAST (\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e) and comparing the homology to the SWR and the domestic horse (as the genome is highly homologous and well annotated). If a sequence was at least 80% homologous to a gene identified via NCBI, that gene was assigned to that sequence. In this study, improvements were made to provide proper gene nomenclature, assignment of genes to scaffolds and sequences, and identification of genes that were not included in the original annotation. The original annotation included 14,274 transcripts that were expressed in granulosa cells according to this analysis. Throughout the annotation, the protein nomenclature was corrected and changed to gene nomenclature. Through manual curation, 6,763 transcripts were corrected from protein to gene nomenclature, with an overall improvement in nomenclature of 47%. In addition, 70% of the transcript sequences were reassigned to different genes after being searched for gene homology. Through de novo sequencing, an additional 7,429 genes were added to the original annotation, with an increase of 61% more genes from the original annotation (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). These improvements in the annotation will be incorporated in a follow-up annotation to be made publicly available with the genome. Supplemental Table\u0026nbsp;1 contains all the sequences identified as a gene as well as their corresponding gene reference.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003ePart 2. Analysis of granulosa cells\u003c/h2\u003e \u003cp\u003ePrior to analysis, principal component analysis (PCA) was performed to determine the variability of the samples within each follicle stage and between the different stages. Figure\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e shows that within each stage, the replicates had very little variability between the two samples at each follicle stage. In addition, there was clear clustering of each stage, independent of the other stages. A heatmap showed the same trend, with granulosa cells (GC) from preovulatory (P) follicles being in their own expression group and GC from growing (G) and dominant (D) follicles clustering closer together (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eAcross all the samples (n\u0026thinsp;=\u0026thinsp;6), 39,455 different transcripts were detected, which corresponded to 21,621 genes, with multiple genes expressing multiple transcripts (Supplemental Table\u0026nbsp;1). The first analysis was to determine the genes that were specific to different follicle stages (growing, dominant, or preovulatory). The growing follicles contained 427 unique transcripts (corresponding to 385 genes), with only 8% of the transcripts identified being unannotated. A total of 812 transcripts were specifically identified in the dominant follicle stage, resulting in a total of 677 genes, and 104 transcripts (13%) were unannotated. Preovulatory follicles contained the largest number of follicle stage-specific transcripts, 945, corresponding to 433 genes. Of the 945 transcripts identified, 484 (51%) were unannotated. Importantly, the most highly expressed genes specific to each follicle type were associated with different transcription factors. The canonical pathways most enriched in growing follicles were involved in the regulation of TP53 activity through phosphorylation, the regulation of TP53 expression and degradation, the signaling of VEGF, the cell cycle checkpoint, and HDR through homologous recombination or single-strand annealing. The genes in preovulatory follicles targeted the following canonical pathways: signaling via the insulin receptor, iron uptake and transport, the iron homeostasis signaling pathway, amino acid regulation of mTORC1, and the adrenergic receptor signaling pathway.\u003c/p\u003e \u003cp\u003eThe next analyses evaluated the differentially expressed genes (DEG) between each follicle category. The comparison with the least number of differentially expressed transcripts was growing versus dominant, with 855 transcripts (resulting in 802 different genes). The genes annotated with the greatest significant differences were \u003cem\u003eRNPC3\u003c/em\u003e and \u003cem\u003eZDHHC5\u003c/em\u003e, both more highly expressed in growing follicles. Figure\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e shows the representative volcano plot for this analysis, highlighting the genes with the largest significant differences. Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e shows the top 10 genes differentially expressed between growing and dominant follicles.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eTop 10 DEGs when comparing growing (G) vs. dominant (D) follicles. A positive fold change indicates that the gene was more highly expressed in G, and a negative fold change indicates that the gene was more highly expressed in D. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGene Symbol\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotation Assignment\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFold Change\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eP-Value\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRNPC3\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e5.030277375\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.03E-22\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eZDHHC5\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e6.697711575\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3.29E-19\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eREV3L\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-3.122418916\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e4.86E-18\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.16139.2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.3228372\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2.34E-17\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCP2DE_g9976.t1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-5.49818783\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3.22E-17\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRYR2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e4.597912368\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e7.13E-17\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.5093.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e3.681205492\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.77E-16\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eHUWE1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-5.088883652\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2.11E-16\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eAKAP5\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.03465133\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2.11E-16\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eFANCG\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.39311232\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e5.43E-16\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eWhen comparing growing and preovulatory follicles, the largest number of transcripts (4,581) were significantly differentially expressed, corresponding to 3,980 genes. The genes annotated with the greatest significant differences were \u003cem\u003eCDH11\u003c/em\u003e and \u003cem\u003eTMEM51\u003c/em\u003e (both of which were more highly expressed in growing follicles, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e lists the top 10 significantly differentially expressed transcripts.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eTop 10 DEGs when comparing growing (G) vs. preovulatory (P) follicles. A positive fold change indicates that the gene was more highly expressed in P, and a negative fold change indicates that the gene was more highly expressed in G. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGene Symbol\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotation Assignment\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFold Change\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eP-Value\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCDH11\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-8.022803609\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.05E-16\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eTMEM51\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-16.3725893\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.19E-14\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eB3GNT4\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-16.06868082\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e5.49E-14\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.2860.2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-11.97253106\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.04E-13\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eGLCE\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e6.57298424\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3.46E-13\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eWDR11\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-5.837110353\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e4.85E-13\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003ePHC3\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e5.981697786\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.66E-12\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eKHSRP\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.14411634\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.90E-12\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRPL6\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.2562075\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e4.13E-12\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eLRP1B\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.99039555\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e8.71E-12\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe final comparison was dominant versus preovulatory follicles. There were 4,315 transcripts differentially expressed within this comparison (3,704 genes). The annotated genes most differentially expressed were \u003cem\u003eHSPG2\u003c/em\u003e and \u003cem\u003eCFAP52\u003c/em\u003e, both of which are more highly expressed in dominant follicles. Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e lists the top 10 differentially expressed transcripts, and Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e shows the genes with the most differential expression. Supplemental Table\u0026nbsp;2 details all the differentially expressed transcripts identified in each comparison.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eTop 10 DEGs when comparing dominant (D) vs. preovulatory (P) follicles. A positive fold change indicates that the gene was more highly expressed in P, and a negative fold change indicates that the gene was more highly expressed in D. Annotation assignment indicates whether that sequence was assigned to a specific gene or was unannotated.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGene Symbol\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotation Assignment\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFold Change\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eP-Value\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-20.063746\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.50707E-23\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCFAP52\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-7.7004087\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.68824E-19\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.2860.2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-12.934659\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e8.26415E-18\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eKHSRP\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.24234\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3.979E-15\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRPL6\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.156424\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e7.45097E-15\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eTMEM51\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.438493\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e7.68241E-15\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSUMO1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-15.366785\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.58062E-15\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.9751.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnannotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e7.62408891\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.98344E-14\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eRPL15\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.491793\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.02625E-13\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003ePTPRM\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotated\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-14.556105\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.71413E-14\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eAlthough not specific to any follicle type, \u003cem\u003eMACIR\u003c/em\u003e was the most highly expressed gene in each follicle category. Table\u0026nbsp;\u003cspan refid=\"Tab4\" class=\"InternalRef\"\u003e4\u003c/span\u003e contains the most highly expressed genes in each category, and interestingly, the same genes are present in multiple follicle types.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab4\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 4\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe most highly expressed genes in each stage. The presence of multiple genes indicates that different transcripts were among the most highly expressed for that stage.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGrowing\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDominant\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePre-Ovulatory\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMACIR\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eMACIR\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eMACIR\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.15261.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.15261.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.15261.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCOL1A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eCOL1A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eCOL1A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSERPINE2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eSERPINE2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.13580.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eSERPINE2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eSLC23A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eKATNIP\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eKATNIP\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eVCAN\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eSLC23A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.13580.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eSERPINE2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSLC23A1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.13580.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eVCAN\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eVCAN\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eHSPG2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eMSTRG.7397.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eFurther evaluation of the DEG in each comparison led to the identification of the most enhanced canonical pathways for each comparison. Table\u0026nbsp;\u003cspan refid=\"Tab5\" class=\"InternalRef\"\u003e5\u003c/span\u003e shows the canonical pathways affected by the most transcripts in each category, and the canonical pathways with the greatest fold changes and z scores are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab5\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 5\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eCanonical pathways affected by the most transcripts in each follicle stage.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGrowing vs. Dominant\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGrowing vs. Preovulatory\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eDominant vs. Preovulatory\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRHO GTPase cycle\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eRho GTPase cycle\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSelenoamino acid metabolism\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRibonucleotide reductase signaling pathway\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSelenoamino acid metabolism\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eResponse of EIF2AK4 (GCN2) to amino acid deficiency\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRole of tissue factor in cancer\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eResponse of EIF2AK4 (GCN2) to amino acid deficiency\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eRho GTPase cycle\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAngiopoietin signaling\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSRP-dependent cotranslational protein targeting to membrane\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNonsense-mediated decay (NMD)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProduction of nitric oxide and reactive oxygen species in macrophages\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCell cycle checkpoints\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eEukaryotic translation termination\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eTo understand the expression patterns across follicle types, a cluster analysis was performed, revealing seven unique patterns of expression (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e). Each cluster contained varying numbers of transcripts and differentially expressed transcripts within different comparisons, ranging from 173 transcripts (cluster four) to 1620 transcripts (cluster five). The expression of 173 transcripts in cluster four increased from G to D to P, with the largest number of differentially expressed transcripts occurring between the G and P follicle stages. Pathways associated with this cluster included focal adhesion, the T-cell receptor signaling pathway, the Pi3K-Akt signaling pathway, and the MAPK signaling pathway. This cluster closely resembled the gene expression patterns associated with follicle differentiation. Cluster five contained the largest number of transcripts (1620) whose expression largely decreased from the D to the P follicle stage. The pathways associated with these transcripts include cellular senescence, EGFR tyrosine kinase inhibitor resistance, and endocytosis.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eThis is the first study evaluating the transcriptomic changes in growing (G), dominant (D) and preovulatory (P) follicles in southern white rhinoceros (SWR). This information is needed to advance assisted reproductive technologies (ARTs) by determining appropriate methods for manipulating ovarian function prior to oocyte collection to obtain oocytes with the highest developmental potential. Through this analysis, we generated candidate gene markers for the physiological status of the follicle at the time of OPU.\u003c/p\u003e \u003cp\u003eAlthough this study was performed on SWR granulosa cells (GC), the SWR genome is incomplete and contains large gaps in annotation. Therefore, the northern white rhinoceros (NWR) genome and currently available annotations are the only resources for generating transcriptomic data on white rhinoceroses, but extensive manual curation is still necessary (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e). In addition, the NWR and SWR are genetically similar at the chromosomal and mitochondrial genome levels, enabling the NWR genome to be used for this analysis (\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e). The original annotation for the NWR contained 14,274 different transcripts (12,252 genes) that were expressed in the GC from this dataset. After adding the transcripts from de novo sequencing, a total of 39,455 different transcripts (21,621 genes) were expressed. This was an overall improvement of 119% in transcripts and 61% in genes compared to the original annotation. In addition, it is important to note that through manual curation, multiple transcripts were reassigned to different genes, and transcripts with unknown associations were identified. This approach resulted in a substantial improvement to the available annotation, and these new annotations will be available publicly in the future. Due to their scarcity, reproductive cells are often underrepresented when generating annotations for genomes. Many new genes were annotated in the granulosa cells due to the massive transcriptional changes that occurred within the cells over time and their dynamic behavior. This additional information to the annotation will benefit other researchers in reproductive physiology in this species, and it emphasizes the importance of using reproductive cells, such as granulosa cells, as a resource for conservation research and genome improvement and accuracy. One of our objectives was to identify the genes most highly expressed across the follicle stages. Interestingly, all follicle stages had the same top three most highly expressed genes. These three genes were involved in macrophage immunometabolism regulator (\u003cem\u003eMACIR\u003c/em\u003e), \u003cem\u003eMSTRG15261.\u003c/em\u003e1, and collagen type I alpha 1 chain (\u003cem\u003eCOL1A1\u003c/em\u003e). The most highly expressed gene, \u003cem\u003eMACIR\u003c/em\u003e, formally known as \u003cem\u003eC5ORF30\u003c/em\u003e, has a well-documented role in modulating the immune response and regulating macrophage function (\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e, \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e). In addition, blocking \u003cem\u003eMACIR\u003c/em\u003e enhances the expression of genes involved in regulating cell migration, adhesion, and angiogenesis (\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e). In cattle GC, \u003cem\u003eMACIR\u003c/em\u003e has been shown to be upregulated in follicles undergoing late-stage atresia (\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e). This is a potential indicator that the follicles in this study were under stress and already progressed through advanced, irreversible atresia. The second most highly expressed gene, \u003cem\u003eMSTRG15261.1\u003c/em\u003e, was found to be unannotated in the SWR and in the domestic horse. The third most highly expressed gene was \u003cem\u003eCOL1A1\u003c/em\u003e, which is the most abundant collagen in the mammalian ovary and is regulated by transforming growth factor beta 1 (\u003cem\u003eTGF-β1\u003c/em\u003e) (\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e). In addition, \u003cem\u003eCOL1A1\u003c/em\u003e promotes the secretion of estrogen while also inhibiting the increase in progesterone, resulting in follicle growth (\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e) (\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e). Most importantly, \u003cem\u003eCOL1A1\u003c/em\u003e plays an important role in oocyte maturation and embryo development through its association with the PI3K/Akt signaling pathway in mice and pigs (\u003cspan additionalcitationids=\"CR42 CR43\" citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e). Interestingly, the most highly expressed gene was involved in follicle death, and the third highest was involved in follicle growth. This could indicate that prior to OPU, follicles received mixed signals during ovarian treatment.\u003c/p\u003e \u003cp\u003eCluster analysis was performed to generate gene expression patterns across the follicle stages to identify genes that follow conventional literature for these stages and those that are not. Due to the limited opportunities to perform OPU in SWRs, the current preferred approach is to stimulate the growth of multiple follicles to increase the oocyte recovery rate. Therefore, unstimulated animals were not considered for this elective procedure in our study. Due to the lack of an unstimulated rhinoceros for comparison, the gene expression patterns were compared to published literature generated from unstimulated cows (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). The two main cluster types we focused on were patterns associated with follicle growth (clusters four and seven) and differentiation (clusters three, five, and six). For cluster three, the gene expression levels were consistent between G and D and increased from D to P; therefore, these genes are associated with follicle differentiation. Of the top 20 genes from this cluster, 13 have known expression patterns across follicle stages (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e); interestingly, five of those genes should be downregulated in P, yet their expression is increased in our dataset. \u003cem\u003eCOL4A1\u003c/em\u003e (collagen type IV alpha 1 chain) was found to be differentially expressed in this cluster compared to cattle data (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). In cattle, \u003cem\u003eCOL4A1\u003c/em\u003e is a gene that stays in relatively low abundance throughout the follicle stages, with a decrease in P and an increase only observed in atretic follicles (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). \u003cem\u003eCOL4A1\u003c/em\u003e is directly regulated by forkhead box L2 (FOXL2), a transcription factor specifically expressed in granulosa cells (\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e). High levels of \u003cem\u003eFOXL2\u003c/em\u003e decrease \u003cem\u003eCOL4A1\u003c/em\u003e (\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e), which is required for the proliferation of GC, allowing follicle remodeling, and should remain high during folliculogenesis to allow for proliferation and differentiation (\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e). In our data, we observed a significant decrease in \u003cem\u003eFOXL2\u003c/em\u003e expression between D and P, which led to a significant increase in \u003cem\u003eCOL4A1\u003c/em\u003e expression between D and P, a sign of late-stage atresia associated with follicular demise (\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e). In addition, high levels of collagen in the follicle basement membrane have been linked to an imbalance in gonadotropin stimulation (\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e). It is hypothesized that this delay in maturation due to increased collagen may lead to asynchrony during maturation and negative effects on oocyte developmental competence (\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e). Therefore, we can hypothesize that high levels of \u003cem\u003eCOL4A1\u003c/em\u003e in our SWRs could be caused by incomplete or suboptimal stimulation prior to OPU.\u003c/p\u003e \u003cp\u003eJunction mediating and regulatory protein p53 cofactor (\u003cem\u003eJMY\u003c/em\u003e) is a transcriptional cofactor and an actin nucleation factor (\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e) (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e). Actin nucleation occurs during oocyte polarization and affects the microtubule and microfilament cytoskeleton by activating Arp 2/3 (for actin remodeling) and the assembly of filaments in the cytoskeletal spindle (\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e, \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e, \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e). If \u003cem\u003eJMY\u003c/em\u003e is not present, the spindle fails to migrate to the cortex, and the oocyte arrests with a centrally located spindle, highlighting the crucial role of JMY during peripheral spindle migration during oocyte maturation (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e). In both mouse and porcine oocytes, JMY decreases during oocyte maturation (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e, \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e). This same pattern of expression was observed in cluster five of our dataset, with JMY being highly significant within that cluster.\u003c/p\u003e \u003cp\u003eCluster four included the genes whose expression increased from G to D to P. A gene of interest in this pattern was F-box and WD repeat domain containing 11 (\u003cem\u003eFBXW11\u003c/em\u003e) due to its regulation by estradiol, its role in meiotic resumption, and its normal pattern of expression compared to our data. \u003cem\u003eFBXW11\u003c/em\u003e is regulated by many different factors, one of which is helicase for meiosis 1 (\u003cem\u003eHFM1\u003c/em\u003e) (\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e). Low levels of \u003cem\u003eHFM1\u003c/em\u003e result in high levels of \u003cem\u003eFBXW11\u003c/em\u003e, which corresponds to the pattern of expression we observed for both of these genes; however, according to the literature in mice, \u003cem\u003eHFM1\u003c/em\u003e expression should be high for oocytes to resume meiosis (\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e). In addition, \u003cem\u003eFBXW11\u003c/em\u003e is regulated by estradiol, as indicated by high levels of estradiol (as found in late stage developing follicles), which results in low levels of \u003cem\u003eFBXW11\u003c/em\u003e (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e). The primary function of \u003cem\u003eFBXW11\u003c/em\u003e, formally known as \u003cem\u003eβ-TRCP2\u003c/em\u003e, is to degrade cell division cycle 25A (\u003cem\u003eCDC25A\u003c/em\u003e), which is crucial for the resumption of meiosis (at the appropriate time) (\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e, \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e). \u003cem\u003eCDC25A\u003c/em\u003e must be downregulated for oocytes to be released from MII arrest (\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e), which we observed in the rhinoceros data. Therefore, from the literature, it appears that the timing of the increase in \u003cem\u003eFBXW11\u003c/em\u003e is important in relation to the resumption of oocyte meiosis.\u003c/p\u003e \u003cp\u003e \u003cem\u003eNRG1\u003c/em\u003e (neuregulin 1) is a member of the EGF-like factor family that is induced by gonadotropins (\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e, \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e). There is a dramatic increase in NRG1 expression due to LH induction, which regulates proper oocyte maturation and developmental competence (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e). Recent literature has shown that \u003cem\u003eNRG1\u003c/em\u003e supplementation during IVM can promote oocyte nuclear maturation and blastocyst development and hatching in cows (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e). During development, amphiregulin (\u003cem\u003eAREG\u003c/em\u003e) and \u003cem\u003eNRG1\u003c/em\u003e work synergistically (\u003cem\u003eAREG\u003c/em\u003e is an accelerator, and \u003cem\u003eNRG1\u003c/em\u003e is a brake) to ensure that the required temporal changes within the oocyte occur for proper meiotic resumption and competence (\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e). It is important that these genes are balanced to ensure that oocytes have high developmental competence. According to our data, \u003cem\u003eNRG1\u003c/em\u003e is part of cluster five and remains consistent from G to D but decreases from D to P, which is the opposite of what the literature states regarding high levels of \u003cem\u003eNRG1\u003c/em\u003e being needed for oocyte competence acquisition (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e). \u003cem\u003eNRG1\u003c/em\u003e is a novel factor that may impact oocyte maturation; therefore, it should be considered a biomarker in GC for oocyte maturation competence and should be added to \u003cem\u003ein vitro\u003c/em\u003e maturation medium.\u003c/p\u003e \u003cp\u003eCluster six was classified as a cluster related to differentiation because there was a decrease in gene expression from G to D and a large increase from D to P. A gene that followed this cluster pattern but, based upon the literature should not be, was thymopoeitin (\u003cem\u003eTMPO\u003c/em\u003e), also known as lamina-associated polypeptide 2 (LAP2) (\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e). According to both bovine and porcine data, \u003cem\u003eTMPO\u003c/em\u003e decreases as follicle progression progresses (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e), which we previously reported from G to D but subsequently deviates with an increase from D to P. A specific region in \u003cem\u003eTMPO\u003c/em\u003e [regions 32\u0026ndash;36 known as thymopentin (\u003cem\u003eTP5\u003c/em\u003e)] has been studied extensively due to its role in the activation of the pluripotency marker \u003cem\u003eLIN28a\u003c/em\u003e in oocytes (\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e, \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e); therefore, high levels of \u003cem\u003eTMPO\u003c/em\u003e result in high levels of \u003cem\u003eLIN28a\u003c/em\u003e. \u003cem\u003eLIN28a\u003c/em\u003e expression decreases until the 8-cell stage of embryo development, when a dramatic increase is then noted (\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e). In contrast, \u003cem\u003eLIN28a\u003c/em\u003e levels should be low until after resumption of meiosis, and high \u003cem\u003eLIN28a\u003c/em\u003e levels in oocytes are associated with MII-arrested oocytes in mice (\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e, \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e). High levels of \u003cem\u003eLAP2\u003c/em\u003e have also been associated with DNA damage in oocytes and programmed cell death in granulosa cells in mice (\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e). Overall, these data indicate that the P follicles in this dataset may have DNA damage, potentially leading to oocyte arrest at MII.\u003c/p\u003e \u003cp\u003eBased upon all the cluster analyses, the largest digression from the literature appears to be in the D to P transition, corresponding to differentiation. The five main drivers of differentiation in cattle are \u003cem\u003eTGF-β1\u003c/em\u003e, tumor protein 53 (\u003cem\u003eTP53\u003c/em\u003e), estradiol, hepatocyte nuclear factor 4 alpha (\u003cem\u003eHNF4A\u003c/em\u003e) and retinoic acid (\u003cem\u003eRA\u003c/em\u003e) (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e). Our dataset included all but \u003cem\u003eRA\u003c/em\u003e, and the expression of each gene significantly increased from D to P. \u003cem\u003eRA\u003c/em\u003e has a positive regulatory effect on granulosa cell proliferation and oocyte maturation (\u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e). \u003cem\u003eRA\u003c/em\u003e levels in growing follicles are required for the follicle to reach appropriate final maturation (differentiation) in response to gonadotropins and affect gonadotropin receptor expression in granulosa cells in cows (\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e). It can therefore be hypothesized that with this ovarian treatment protocol, follicles are not receiving the proper stimuli to fully differentiate and produce oocytes with the highest developmental competence potential.\u003c/p\u003e \u003cp\u003eBased upon these results, it appears that the ovarian treatment protocol used on these animals did not result in GC or oocytes with the highest developmental competence potential. It is important to note that these comparisons are based primarily on data from unstimulated cows. We hypothesize that the gene expression changes we observed could be an outcome of the simulation protocol, but it is unknown whether some of these deviations are due to the unique expression pattern of the rhinoceros from other species. Because we performed oocyte group culture \u003cem\u003ein vitro\u003c/em\u003e, we cannot directly correlate the \u003cem\u003ein vivo\u003c/em\u003e GC analyzed with the corresponding oocyte maturation outcome; hence, in the future, we will match GC to the \u003cem\u003ein vitro\u003c/em\u003e maturation outcome of oocytes by performing individual culture. Overall, for the animals treated with this stimulation protocol, we achieved an \u003cem\u003ein vitro\u003c/em\u003e maturation rate of 33%, but none of those oocytes fertilized or developed into embryos. Although a percentage of retrieved oocytes were capable of extruding the first polar body and reaching metaphase II (MII) after \u003cem\u003ein vitro\u003c/em\u003e maturation, we hypothesize that unpaired cytoplasmic maturation, erroneous chromosomal segregation, or morphological characteristics of atresia (\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e) left the oocyte incapable of fertilization and becoming a zygote. Incomplete cytoplasmic maturation can lead to failed sperm decondensation, erroneous male pronuclear formation, and ultimately unsuccessful fertilization (\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e). However, additional studies are needed to characterize the effects of ovarian stimulation on oocyte developmental competence acquisition/failure and cytoplasmic maturation in this species. A balance between the quantity and quality of oocytes obtained during OPU is suggested, and low, frequent doses of GnRH might be a valuable treatment.\u003c/p\u003e \u003cp\u003eWith this dataset, we generated a list of potential biomarkers of the GC for the physiologic status of the follicle to determine which oocytes had the highest potential for developmental competence. These markers include \u003cem\u003eCOL1A1\u003c/em\u003e (follicle growth) (\u003cspan additionalcitationids=\"CR42 CR43\" citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e), \u003cem\u003eJMY\u003c/em\u003e and \u003cem\u003eFBXW11\u003c/em\u003e (oocyte preparedness to resume meiosis) (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e, \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e), \u003cem\u003eNRG1\u003c/em\u003e (oocyte maturation potential) (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e), \u003cem\u003eTMPO\u003c/em\u003e (DNA damage) (\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e), and \u003cem\u003eMACIR\u003c/em\u003e and \u003cem\u003eCOL4A1\u003c/em\u003e (markers for cell death and atresia) (\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e, \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e). However, further testing will be needed to correlate these markers with the status of the follicles we collected during OPU and ultimately with oocyte maturation and cleavage.\u003c/p\u003e \u003cp\u003eIt is important to note that one limitation in this study is grouping the follicle stages were grouped based upon follicle size at OPU. Nevertheless, for this study, we were not able to obtain follicles at all stages from the same animal due to the unknown status of the estrous cycle prior to ovarian stimulation and OPU. We attempted to overcome this limitation by obtaining technical replicates from different granulosa cell aliquots of the preovulatory follicle. In addition, the animals in this study were not ultrasounded regularly prior to ovarian treatment and OPU; therefore, the trajectory of the follicle prior to OPU was unknown (i.e., was it increasing in size or was it already regressing and undergoing atresia). We attempted to overcome this by evaluating multiple follicles per stage, and in the future, we will attempt to monitor these animals prior to OPU.\u003c/p\u003e \u003cp\u003eIn conclusion, this was the first study to evaluate the transcriptomic differences in GC between the different stages of follicle development in the SWR, an important chapter in the reproductive knowledge of this species. In this threatened species, much remains unknown regarding ovarian physiology and follicle activity; nevertheless, it is important to continue addressing these reproductive questions to support assisted reproductive technologies for rhinoceros species. Finally, these data provide multiple biomarkers to potentially improve OPU results. This body of work contains novel information on the SWR and enriches the reproductive knowledge of the species, which will support conservation efforts for this species.\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eAnimal management and ovum pickup (OPU)\u003c/h2\u003e \u003cp\u003e All procedures, experiments, and methods were reviewed and approved by San Diego Zoo Wildlife Alliance\u0026rsquo;s Institutional Animal Care and Use Committee (IACUC; protocol number 18\u0026ndash;018). In addition, guidelines set forth by ARRIVE (Animal Research: Reporting of In Vivo Experiments, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://arriveguidelines.org/arrive-guidelines\u003c/span\u003e\u003cspan address=\"https://arriveguidelines.org/arrive-guidelines\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) were followed. All methods for this study were performed in accordance with all relevant guidelines and regulations of the IACUC and ARRIVE. Three southern white rhinoceros females (aged 7, 11, and 24 at the time of OPU) were utilized for this study. While under the care of the SDZWA, all rhinoceroses have undergone health assessments to ensure that no reproductive pathologies, tumors, or cysts were present and were continually monitored. Complete information regarding each individual rhinoceros is included in supplemental table 3. All females were stimulated prior to transrectal OPU as previously described (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e, \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e), facilitating the \u003cem\u003ein vivo\u003c/em\u003e collection of follicular fluid aspirate from follicles of distinct sizes. Prior to OPU, females received synthetic chlormadinone acetate (CMA) at 30 mg/day for 35 days (days 0\u0026ndash;34). The day following CMA withdrawal (Day 35), all animals received 3 mg of deslorelin [gonadotropin-releasing hormone (GnRH) analog] followed by 2 mg of deslorelin every 48 h prior to OPU (Day 37 and 39; 7 mg total) via intramuscular injection. OPU was performed on day 41. All rhinoceroses were anesthetized using a combination of etorphine (2 mcg/kg), butorphanol (20\u0026ndash;24 mcg/kg), medetomidine (23\u0026ndash;25 mcg/kg), and azaperone (14\u0026ndash;16 mcg/kg) administered intramuscularly using a remote drug delivery system. During initial positioning to facilitate intubation, propofol was administered intravenously (0.5 mg/kg). Following OPU procedures, anesthesia reversal was achieved using atipamezole delivered intramuscularly at a 5:1 ratio to medetomidine (116\u0026ndash;125 mcg/kg) and naltrexone delivered intramuscularly at a target 50:1 ratio to etorphine (53\u0026ndash;61 mcg/kg).\u003c/p\u003e \u003cp\u003eFollowing fecal removal, rinsing, and disinfection, OPU was achieved using a customized, ultrasound-guided probe accompanied by three double-lumen needles. The follicles were ablated, the contents were aspirated, and the samples were rinsed with a warm (37\u0026deg;C) flushing solution (Vigro) containing 12.5 I.U./mL of heparin. During OPU, follicle diameters were measured via trans-rectal ultrasonography. Follicular fluid and granulosa cells were separated based upon follicle stage/size. The follicle stages were growing (11\u0026ndash;17 mm), dominant (25\u0026ndash;29 mm), and preovulatory (30\u0026ndash;34 mm).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eGranulosa cell collection and RNA isolation/quantification\u003c/h2\u003e \u003cp\u003eAfter OPU and the oocyte was collected from the collection fluid, free-floating mural granulosa cells were collected at the follicle stage, pipetted directly into RNAlater (Thermo Fisher Scientific, Waltham, MA) for 24 hours at 4\u0026deg;C and subsequently transferred to -80\u0026deg;C until RNA isolation. Prior to isolation, the granulosa cells were thawed, mixed 1:1 with cold PBS (Sigma Aldrich, St. Louis, MO) and centrifuged at 3,000 \u0026times; g for 10 min. The supernatant was discarded, and the pellet was resuspended in cold PBS. The samples were centrifuged again at 3,000 \u0026times; g for 5 min to remove all remnants of RNAlater. Total RNA was isolated from granulosa cells using an Arcturus PicoPure RNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA) per the manufacturer\u0026rsquo;s instructions. Two biological replicates (one sample from two separate individuals) were used for both growing and dominant follicles. Neither of these two animals had a preovulatory follicle present at the time of OPU, and a third animal provided two technical replicates (the GC were stored in numerous separate aliquots) for the preovulatory stage. In total, six samples were used for this study (n\u0026thinsp;=\u0026thinsp;2 per follicle type). Granulosa cells were incubated with extraction buffer for 30 min before centrifugation to remove debris and extracellular material. The cell extract was incubated with ethanol, bound to, and washed on preconditioned purification columns. Total RNA was recovered into elution buffer and quantified using a Qubit 4 Fluorometer (Thermo Fischer Scientific, Waltham, MA. After quantification, the samples were evaluated on a 4150 TapeStation System (Agilent Technologies, Santa Clara, CA) to determine RNA integrity number (RIN) values. Only samples with a RIN greater than 6.0 were used for RNA-Seq analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eRNA sequencing\u003c/h2\u003e \u003cp\u003eLibrary preparation and RNA sequencing were performed at the University of California San Diego Institute for Genomic Medicine Center. Following the manufacturer\u0026rsquo;s protocol, RNA-sequencing (cDNA) libraries were prepared using the Illumina TruSeq Stranded Total RNA Prep with Ribo-Zero Plus. A library was prepared for each sample using 50 ng of total RNA. Briefly, RNA was rRNA depleted and fragmented, first- and second-strand cDNA synthesis was performed, and adapters were ligated. Each sample was amplified with a specific barcoded PCR primer for sample identification purposes. The prepared libraries were sequenced as 100 base pair paired end reads on an Illumina NovaSeq 6000 (Illumina, San Diego, CA). The raw data files were uploaded to the Gene Expression Omnibus under accession number GSE261038.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eBioinformatic and statistical analysis\u003c/h2\u003e \u003cp\u003eBioinformatic analysis was performed on the Galaxy web platform using the public server usegalaxy.org (\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e). Initial assessment of sequence quality was performed by FastQC, and the sequences were aggregated for comparison by MultiQC (\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e, \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e72\u003c/span\u003e). The adapter sequences and reads with quality scores below the quality threshold of 25 and with fewer than 20 base pairs were removed from analysis via Trimmomatic (\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e). HISAT2 was used to align reads to the northern white rhinoceros (NWR) genome CerSimCot1.0 (GCA_021442165.1)(\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e74\u003c/span\u003e). The southern white rhinoceros genome (CerSimSim1.0; GCA_000283155.1) was not used for this analysis due to large gaps in the sequencing and extremely limited annotation. Transcript assembly and quantification for both annotated and unannotated transcripts were initially performed with Stringtie, and each animal sample was analyzed individually (\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e). The resulting transcripts for individual samples were merged to create a final annotation file representing the intersection and union of all samples (\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e). Stringtie was subsequently applied to the complete merged annotation to generate transcript read counts. Transcript counts were exported from Galaxy and imported into RStudio for statistical analysis (\u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e). Transcript read counts were analyzed between each follicle type utilizing DESeq2 within R (\u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e). The average expression value for replicates within a follicle stage was \u0026ge; 3, except in the cases where the expression value was 0 for the replicates within a follicle type. When not expressed in a follicle type, the expression value for both replicates was 0. The data were normalized internally using the median ratio method of DESeq2. Benjamini Hochberg false discovery rate adjustment was used. Significance was assessed at \u003cem\u003eP\u003c/em\u003e \u0026le; 0.05. Differentially expressed transcripts were grouped into different cluster patterns according to their expression profiles among the different developmental stages using the Mfuzz Bioconductor package (\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e78\u003c/span\u003e). Pathway analysis and visualization were completed using Qiagen Ingenuity Pathway Analysis (\u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e). Manual curation was required to assign sequences identified to corresponding genes and to ensure that the HGNC/VGNC naming conventions were maintained.\u003c/p\u003e \u003c/div\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eCompeting Interests\u003c/h2\u003e \u003cp\u003eThe authors have no competing interests to declare.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eFunding Sources\u003c/h2\u003e \u003cp\u003eThis research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.\u003c/p\u003e\u003ch2\u003eAuthor Contributions\u003c/h2\u003e \u003cp\u003eE.R. and B.D. collected the samples. E.R. processed and prepared the samples for RNA sequencing. E.R. and K.K. performed the bioinformatic analyses and wrote the manuscript. M.A.S. aided in comparative data analysis and interpretation. B.D. aided in the writing and revision of the manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgments\u003c/h2\u003e \u003cp\u003eWe thank Dr. Candace Williams at the San Diego Zoo Wildlife Alliance and Dr. Ahmed Gad and Nico Menjivar at Colorado State University for aiding in the data analysis and visualization. In addition, we would like to acknowledge the Conservation Science Wildlife Health team, the Veterinary Services, and the Rhino Rescue Center team at San Diego Zoo Wildlife Alliance for aiding in animal management, care, and sample collection.\u003c/p\u003e\u003ch2\u003eData Availability Statement\u003c/h2\u003e \u003cp\u003eThe datasets generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO) under accession number GSE261038.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eHildebrandt TB, Holtze S, Colleoni S, Hermes R, Stejskal J, Lekolool I et al (2023) In vitro fertilization (IVF) program in white rhinoceros. Reproduction (Cambridge, England). :REP-23-0087\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRuggeri E, Klohonatz K, Sirard M-A, Durrant B, Coleman S (2023) Genomic insights into southern white rhinoceros (Ceratotherium simum simum) reproduction: revealing granulosa cell gene expression. Theriogenology Wild. :100055\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePennington PM, Marshall KL, Capiro JM, Howard L, Durrant BS (2020) Pregnancies following long luteal phases in southern white rhinoceros (Ceratotherium simum simum). Zoo Biol 39(2):141\u0026ndash;144\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePennington PM, Durrant BS (2019) Assisted reproductive technologies in captive rhinoceroses. Mammal Rev 49(1):1\u0026ndash;15\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHildebrandt TB, Hermes R, Colleoni S, Diecke S, Holtze S, Renfree MB et al (2018) Embryos and embryonic stem cells from the white rhinoceros. Nat Commun 9(1):2589\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaragusty J, Diecke S, Drukker M, Durrant B, Friedrich Ben-Nun I, Galli C et al (2016) Rewinding the process of mammalian extinction. Zoo Biol 35(4):280\u0026ndash;292\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRoth TL (2006) A review of the reproductive physiology of rhinoceros species in captivity. Int Zoo Yearbook 40(1):130\u0026ndash;143\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePatton ML, Swaisgood RR, Czekala NM, White AM, Fetter GA, Montagne JP et al (1999) Reproductive cycle length and pregnancy in the southern white rhinoceros (Ceratotherium simum simum) as determined by fecal pregnane analysis and observations of mating behavior. Zoo Biol 18(2):111\u0026ndash;127\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRoth TL (2023) That was then, this is now\u0026ndash;over two decades of progress in rhinoceros reproductive science and technology. Theriogenology Wild. :100065\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFerre LB, Kjelland ME, Strobech LB, Hyttel P, Mermillod P, Ross PJ (2020) Review: Recent advances in bovine in vitro embryo production: reproductive biotechnology history and methods. Animal 14(5):991\u0026ndash;1004\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChang EM, Song HS, Lee DR, Lee WS, Yoon TK (2014) In vitro maturation of human oocytes: Its role in infertility treatment and new possibilities. Clin Exp Reprod Med 41(2):41\u0026ndash;46\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eThomas MR, Sparks AE, Ryan GL, Van Voorhis BJ (2010) Clinical predictors of human blastocyst formation and pregnancy after extended embryo culture and transfer. Fertil Steril 94(2):543\u0026ndash;548\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHinrichs K (2010) The equine oocyte: factors affecting meiotic and developmental competence. Mol Reprod Dev 77(8):651\u0026ndash;661\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGalli C, Colleoni S, Duchi R, Lagutina I, Lazzari G (2007) Developmental competence of equine oocytes and embryos obtained by in vitro procedures ranging from in vitro maturation and ICSI to embryo culture, cryopreservation and somatic cell nuclear transfer. Anim Reprod Sci 98(1\u0026ndash;2):39\u0026ndash;55\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang G, Br\u0026auml;ndl B, Rohrandt C, Hong K, Pang A, Lee J et al Chromosome-level genome assembly of the functionally extinct northern white rhinoceros (Ceratotherium simum cottoni). bioRxiv. 2021:2021.12.11.472206.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeuffels-Barkas J, Wilsher S, Allen WRT, Ververs C, Lueders I (2023) Comparative reproduction of the female horse, elephant and rhinoceros: implications for advancing Assisted Reproductive Technologies (ART). Reprod Fertil. ;4(3)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRuggeri E, Young C, Ravida N, Sirard MA, Krisher R, de la Rey M et al (2022) Glucose consumption and gene expression in granulosa cells collected before and after in vitro oocyte maturation in the southern white rhinoceros (Ceratotherium simum simum). Reprod Fertil Dev\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePennington PM, Marshall KL, Capiro JM, Felton RG, Durrant BS (2019) Ovulation induction in anovulatory southern white rhinoceros (Ceratotherium simum simum) without altrenogest. Conserv Physiol 7(1):coz033\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGirard A, Dufort I, Douville G, Sirard MA (2015) Global gene expression in granulosa cells of growing, plateau and atretic dominant follicles in cattle. Reprod Biol Endocrinol 13:17\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSquires E (2020) Current Reproductive Technologies Impacting Equine Embryo Production. J Equine Vet Sci 89:102981\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDe Roover R, Feugang JM, Bols PE, Genicot G, Hanzen C (2008) Effects of ovum pick-up frequency and FSH stimulation: a retrospective study on seven years of beef cattle in vitro embryo production. Reprod Domest Anim 43(2):239\u0026ndash;245\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarnevale EM, Sessions DR (2012) In vitro production of equine embryos. J Equine Vet Sci 32(7):367\u0026ndash;371\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarnevale EM, Stokes J, Squires EL, Campos-Chillon LF, Altermatt J, Suh TK (2007) Clinical use of intracytoplasmic sperm injection in horses. Proc Am Assoc Equine Pract 53:560\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarnevale EM, Maclellan LJ (2006) Collection, evaluation, and use of oocytes in equine assisted reproduction. Vet Clin North Am Equine Pract 22(3):843\u0026ndash;856\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBenammar A, Derisoud E, Vialard F, Palmer E, Ayoubi JM, Poulain M et al (2021) The Mare: A Pertinent Model for Human Assisted Reproductive Technologies? Animals (Basel). ;11(8)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBo GA, Mapletoft RJ (2020) Superstimulation of ovarian follicles in cattle: Gonadotropin treatment protocols and FSH profiles. Theriogenology 150:353\u0026ndash;359\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMantovani AP, Reis E, Gacek F, B\u0026oacute; G, Binelli M, Baruselli PS (2018) Prolonged use of a progesterone-releasing intravaginal device (CIDR\u0026reg;) for induction of persistent follicles in bovine embryo recipients. Anim Reprod (AR) 2(4):272\u0026ndash;277\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKhan DR, Fournier E, Dufort I, Richard FJ, Singh J, Sirard MA (2016) Meta-analysis of gene expression profiles in granulosa cells during folliculogenesis. Reproduction 151(6):R103\u0026ndash;R110\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCheng J, Huang J, Yuan S, Zhou S, Yan W, Shen W et al (2017) Circular RNA expression profiling of human granulosa cells during maternal aging reveals novel transcripts associated with assisted reproductive technology outcomes. PLoS ONE 12(6):e0177888\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Z, Wang J, Zhao Y, Ma D, Zhao M, Li N et al (2021) scRNA-seq of ovarian follicle granulosa cells from different fertility goats reveals distinct expression patterns. Reprod Domest Anim 56(5):801\u0026ndash;811\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDouville G, Sirard MA (2014) Changes in granulosa cells gene expression associated with growth, plateau and atretic phases in medium bovine follicles. J Ovarian Res 7:50\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNivet AL, Vigneault C, Blondin P, Sirard MA (2013) Changes in granulosa cells' gene expression associated with increased oocyte competence in bovine. Reproduction 145(6):555\u0026ndash;565\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCamacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K et al (2009) BLAST+: architecture and applications. BMC Bioinformatics 10:421\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDorris ER, Tazzyman SJ, Moylett J, Ramamoorthi N, Hackney J, Townsend M et al (2019) The Autoimmune Susceptibility Gene C5orf30 Regulates Macrophage-Mediated Resolution of Inflammation. J Immunol 202(4):1069\u0026ndash;1078\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMuthana M, Hawtree S, Wilshaw A, Linehan E, Roberts H, Khetan S et al (2015) C5orf30 is a negative regulator of tissue damage in rheumatoid arthritis. Proceedings of the National Academy of Sciences. ;112(37):11618-23\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSerwin K, Kozak M, Mazurek-Mochol M, Piotrowska K, Czerewaty M, Pawlik A (2023) Human macrophage immunometabolism regulator (MACIR) in patients with periodontitis. Immunobiology 228(6):152760\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHatzirodos N, Hummitzsch K, Irving-Rodgers HF, Harland ML, Morris SE, Rodgers RJ (2014) Transcriptome profiling of granulosa cells from bovine ovarian follicles during atresia. BMC Genomics 15:40\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi H, Chang HM, Shi Z, Leung PCK (2020) The p38 signaling pathway mediates the TGF-beta1-induced increase in type I collagen deposition in human granulosa cells. FASEB J 34(11):15591\u0026ndash;15604\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuet C, Monget P, Pisselet C, Monniaux D (1997) Changes in extracellular matrix components and steroidogenic enzymes during growth and atresia of antral ovarian follicles in the sheep. Biol Reprod 56(4):1025\u0026ndash;1034\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuet C, Pisselet C, Mandon-Pepin B, Monget P, Monniaux D (2001) Extracellular matrix regulates ovine granulosa cell survival, proliferation and steroidogenesis: relationships between cell shape and function. J Endocrinol 169(2):347\u0026ndash;360\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShimada M, Ito J, Yamashita Y, Okazaki T, Isobe N (2003) Phosphatidylinositol 3-kinase in cumulus cells is responsible for both suppression of spontaneous maturation and induction of gonadotropin-stimulated maturation of porcine oocytes. J Endocrinol 179(1):25\u0026ndash;34\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHoshino Y, Yokoo M, Yoshida N, Sasada H, Matsumoto H, Sato E (2004) Phosphatidylinositol 3-kinase and Akt participate in the FSH-induced meiotic maturation of mouse oocytes. Mol Reprod Dev 69(1):77\u0026ndash;86\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHall SE, Upton RM, McLaughlin EA, Sutherland JM (2018) Phosphoinositide 3-kinase/protein kinase B (PI3K/AKT) and Janus kinase/signal transducer and activator of transcription (JAK/STAT) follicular signalling is conserved in the mare ovary. Reprod Fertility Dev 30(4):624\u0026ndash;633\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRiley JK, Carayannopoulos MO, Wyman AH, Chi M, Ratajczak CK, Moley KH (2005) The PI3K/Akt pathway is present and functional in the preimplantation mouse embryo. Dev Biol 284(2):377\u0026ndash;386\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDu H, Guo Y, Wu X, Gong Y (2022) FOXL2 regulates the expression of the Col4a1 collagen gene in chicken granulosa cells. Mol Reprod Dev 89(2):95\u0026ndash;103\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMarongiu M, Deiana M, Marcia L, Sbardellati A, Asunis I, Meloni A et al (2016) Novel action of FOXL2 as mediator of Col1a2 gene autoregulation. Dev Biol 416(1):200\u0026ndash;211\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAdriaenssens T, Mazoyer C, Segers I, Wathlet S, Smitz J (2009) Differences in collagen expression in cumulus cells after exposure to highly purified menotropin or recombinant follicle-stimulating hormone in a mouse follicle culture model. Biol Reprod 80(5):1015\u0026ndash;1025\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZuchero JB, Coutts AS, Quinlan ME, Thangue NBL, Mullins RD (2009) p53-cofactor JMY is a multifunctional actin nucleation factor. Nat Cell Biol 11(4):451\u0026ndash;459\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSun S-C, Sun Q-Y, Kim N-H (2011) JMY is required for asymmetric division and cytokinesis in mouse oocytes. MHR: Basic Sci reproductive Med 17(5):296\u0026ndash;304\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBrunet S, Verlhac MH (2011) Positioning to get out of meiosis: the asymmetry of division. Hum Reprod Update 17(1):68\u0026ndash;75\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNamgoong S, Kim NH (2016) Roles of actin binding proteins in mammalian oocyte maturation and beyond. Cell Cycle 15(14):1830\u0026ndash;1843\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLin Z, Xu Y-N, Namgoong S, Kim N-H (2014) JMY functions as actin nucleation-promoting factor and mediator for p53-mediated DNA damage in porcine oocytes. PLoS ONE 9(10):e109385\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePu D, Zhong C, Wang H, Yuan X, He Y, Cong J et al (2023) Regulation of FUS ubiquitination and localization by HFM1 is essential for oocyte meiosis prophase I progression in mice\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXie J, Jin Y, Wang G (2019) The role of SCF ubiquitin-ligase complex at the beginning of life. Reproductive Biology Endocrinol 17(1):1\u0026ndash;9\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSilverman JS, Skaar JR, Pagano M (2012) SCF ubiquitin ligases in the maintenance of genome stability. Trends Biochem Sci 37(2):66\u0026ndash;73\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIsoda M, Kanemori Y, Nakajo N, Uchida S, Yamashita K, Ueno H et al (2009) The extracellular signal-regulated kinase-mitogen-activated protein kinase pathway phosphorylates and targets Cdc25A for SCF beta-TrCP-dependent degradation for cell cycle arrest. Mol Biol Cell 20(8):2186\u0026ndash;2195\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOh JS, Susor A, Schindler K, Schultz RM, Conti M (2013) Cdc25A activity is required for the metaphase II arrest in mouse oocytes. J Cell Sci 126(5):1081\u0026ndash;1085\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChowdhury I, Branch A, Mehrabi S, Ford BD, Thompson WE (2017) Gonadotropin-Dependent Neuregulin-1 Signaling Regulates Female Rat Ovarian Granulosa Cell Survival. Endocrinology 158(10):3647\u0026ndash;3660\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNoma N, Kawashima I, Fan HY, Fujita Y, Kawai T, Tomoda Y et al (2011) LH-induced neuregulin 1 (NRG1) type III transcripts control granulosa cell differentiation and oocyte maturation. Mol Endocrinol 25(1):104\u0026ndash;116\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDellaqua TT, Vigaro RA, Janini LCZ, Dal Canto M, Renzini MM, Lodde V et al (2023) Neuregulin 1 (NRG1) modulates oocyte nuclear maturation during IVM and improves post-IVF embryo development. Theriogenology 195:209\u0026ndash;216\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMiao X, Guo R, Williams A, Lee C, Ma J, Wang PJ et al (2023) Replication Protein A1 is essential for DNA damage repair during mammalian oogenesis. bioRxiv\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKulus J, Kulus M, Kranc W, Jopek K, Zdun M, Jozkowiak M et al (2021) Transcriptomic Profile of New Gene Markers Encoding Proteins Responsible for Structure of Porcine Ovarian Granulosa Cells. Biology (Basel). ;10(11)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVogt EJ, Meglicki M, Hartung KI, Borsuk E, Behr R (2012) Importance of the pluripotency factor LIN28 in the mammalian nucleolus during early embryonic development. Development 139(24):4514\u0026ndash;4523\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu T, Jing F, Huang P, Geng Z, Xu J, Li J et al (2021) Thymopentin alleviates premature ovarian failure in mice by activating YY2/Lin28A and inhibiting the expression of let-7 family microRNAs. Cell Prolif 54(8):e13089\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFlemr M, Moravec M, Libova V, Sedlacek R, Svoboda P (2014) Lin28a is dormant, functional, and dispensable during mouse oocyte-to-embryo transition. Biol Reprod 90(6):131\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKawai T, Yanaka N, Richards JS, Shimada M (2016) De Novo-Synthesized Retinoic Acid in Ovarian Antral Follicles Enhances FSH-Mediated Ovarian Follicular Cell Differentiation and Female Fertility. Endocrinology 157(5):2160\u0026ndash;2172\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDamdimopoulou P, Chiang C, Flaws JA (2019) Retinoic acid signaling in ovarian folliculogenesis and steroidogenesis. Reprod Toxicol 87:32\u0026ndash;41\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHinrichs K (1991) The relationship of follicle atresia to follicle size, oocyte recovery rate on aspiration, and oocyte morphology in the mare. Theriogenology 36(2):157\u0026ndash;168\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFulka J Jr., First NL, Moor RM (1998) Nuclear and cytoplasmic determinants involved in the regulation of mammalian oocyte maturation. Mol Hum Reprod 4(1):41\u0026ndash;49\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAfgan E, Baker D, Batut B, van den Beek M, Bouvier D, Cech M et al (2018) The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res 46(W1):W537\u0026ndash;W44\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAndrews S (2010) FastQC: a quality control tool for high throughput sequence data\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEwels P, Magnusson M, Lundin S, Kaller M (2016) MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32(19):3047\u0026ndash;3048\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30(15):2114\u0026ndash;2120\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKim D, Paggi JM, Park C, Bennett C, Salzberg SL (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37(8):907\u0026ndash;915\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol 33(3):290\u0026ndash;295\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eR Core Team R (2013) R: A language and environment for statistical computing\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLove MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15(12):550\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKumar L (2007) Mfuzz: a software package for soft clustering of microarray data. Bioinformation 2(1):5\u0026ndash;7\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKramer A, Green J, Pollard J Jr., Tugendreich S (2014) Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 30(4):523\u0026ndash;530\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"San Diego Zoo Wildlife Alliance","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"southern white rhinoceros, OPU, ARTs, oocytes, granulosa cells, reproduction, endangered species","lastPublishedDoi":"10.21203/rs.3.rs-4765147/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4765147/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eMuch remains unknown about the reproductive physiology of southern white rhinoceros (SWR) and the effect of ovarian stimulation prior to ovum pickup (OPU) have not been fully elucidated. Granulosa cells (GC) provide valuable insight into follicle growth and oocyte maturation status. The goals of this study were to evaluate transcriptomic changes in GC from three stages of follicle development and to identify biomarkers possibly associated with follicular growth and maturation because of ovarian stimulation. GC collected from SWRs following OPU were assigned stages based upon follicle size. Total RNA was isolated, and cDNA libraries were prepared and sequenced on a NovaSeq 6000. All bioinformatics analyses were performed utilizing the Galaxy web platform. Reads were aligned to CerSimCot1.0, and the manual curation was performed with EquCab3.0. Overall, 39,455 transcripts (21,612 genes) were identified across follicle stages, and manual curation yielded a 61% increase in gene identification from the original annotation. Granulosa cells from preovulatory follicles expressed the highest number of unique transcripts. The following seven biomarkers were determined based upon cluster analysis and patterns of expression: \u003cem\u003eCOL1A1\u003c/em\u003e, \u003cem\u003eJMY, FBXW11\u003c/em\u003e, \u003cem\u003eNRG1, TMPO\u003c/em\u003e, \u003cem\u003eMACIR\u003c/em\u003e and \u003cem\u003eCOL4A1\u003c/em\u003e. These data can be used to potentially evaluate the effects of different ovarian stimulation protocols on follicle dynamics, improve OPU results, and support conservation efforts in this species.\u003c/p\u003e","manuscriptTitle":"Granulosa cells provide transcriptomic information on ovarian follicle dynamics in the southern white rhinoceros","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-07-22 15:16:30","doi":"10.21203/rs.3.rs-4765147/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"61c1502f-83ef-4890-a0bd-5dcf72b4698e","owner":[],"postedDate":"July 22nd, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":34817269,"name":"Animal Science"},{"id":34817270,"name":"Wildlife Biology"},{"id":34817271,"name":"Bioinformatics"},{"id":34817272,"name":"Molecular Genetics"},{"id":34817273,"name":"Animal Physiology"}],"tags":[],"updatedAt":"2024-07-22T15:16:30+00:00","versionOfRecord":[],"versionCreatedAt":"2024-07-22 15:16:30","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4765147","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4765147","identity":"rs-4765147","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}