Advances in modeling of the bovine estrous cycle: Synchronization with PGF2α

Our model of the bovine estrous cycle is a set of ordinary differential equations which generates hormone profiles of successive estrous cycles with several follicular waves per cycle. It describes the growth and decay of the follicles and the corpus luteum, as well as the change of the key reproductive hormones, enzymes and processes over time. In this work we describe recent developments of this model towards the administration of prostaglandin F2α. We validate our model by showing that the simulations agree with observations from synchronization studies and with measured progesterone data after single dose administrations of synthetic prostaglandin F2α.     

A delay differential equation model of follicle waves in women

This article presents a mathematical model for hormonal regulation of the menstrual cycle which predicts the occurrence of follicle waves in normally cycling women. Several follicles of ovulatory size that develop sequentially during one menstrual cycle are referred to as follicle waves. The model consists of 13 nonlinear, delay differential equations with 51 parameters. Model simulations exhibit a unique stable periodic cycle and this menstrual cycle accurately approximates blood levels of ovarian and pituitary hormones found in the biological literature. Numerical experiments illustrate that the number of follicle waves corresponds to the number of rises in pituitary follicle stimulating hormone. Modifications of the model equations result in simulations which predict the possibility of two ovulations at different times during the same menstrual cycle and, hence, the occurrence of dizygotic twins via a phenomenon referred to as superfecundation. Sensitive parameters are identified and bifurcations in model behaviour with respect to parameter changes are discussed. Studying follicle waves may be helpful for improving female fertility and for understanding some aspects of female reproductive ageing.     

Validation of a mathematical model of the bovine estrous cycle for cows with different estrous cycle characteristics

A recently developed mechanistic mathematical model of the bovine estrous cycle was parameterized to fit empirical data sets collected during one estrous cycle of 31 individual cows, with the main objective to further validate the model. The a priori criteria for validation were (1) the resulting model can simulate the measured data correctly (i.e. goodness of fit), and (2) this is achieved without needing extreme, probably non-physiological parameter values. We used a least squares optimization procedure to identify parameter configurations for the mathematical model to fit the empirical in vivo measurements of follicle and corpus luteum sizes, and the plasma concentrations of progesterone, estradiol, FSH and LH for each cow. The model was capable of accommodating normal variation in estrous cycle characteristics of individual cows. With the parameter sets estimated for the individual cows, the model behavior changed for 21 cows, with improved fit of the simulated output curves for 18 of these 21 cows. Moreover, the number of follicular waves was predicted correctly for 18 of the 25 two-wave and three-wave cows, without extreme parameter value changes. Estimation of specific parameters confirmed results of previous model simulations indicating that parameters involved in luteolytic signaling are very important for regulation of general estrous cycle characteristics, and are likely responsible for differences in estrous cycle characteristics between cows.     

The Expression Pattern of microRNAs in Granulosa Cells of Subordinate and Dominant Follicles during the Early Luteal Phase of the Bovine Estrous Cycle

This study aimed to investigate the miRNA expression patterns in granulosa cells of subordinate (SF) and dominant follicle (DF) during the early luteal phase of the bovine estrous cycle. For this, miRNA enriched total RNA isolated from granulosa cells of SF and DF obtained from heifers slaughtered at day 3 and day 7 of the estrous cycle was used for miRNAs deep sequencing. The results revealed that including 17 candidate novel miRNAs, several known miRNAs (n = 291–318) were detected in SF and DF at days 3 and 7 of the estrous cycle of which 244 miRNAs were common to all follicle groups. The let-7 families, bta-miR-10b, bta-miR-26a, bta-miR-99b and bta-miR-27b were among abundantly expressed miRNAs in both SF and DF at both days of the estrous cycle. Further analysis revealed that the expression patterns of 16 miRNAs including bta-miR-449a, bta-miR-449c and bta-miR-222 were differentially expressed between the granulosa cells of SF and DF at day 3 of the estrous cycle. However, at day 7 of the estrous cycle, 108 miRNAs including bta-miR-409a, bta-miR-383 and bta-miR-184 were differentially expressed between the two groups of granulosa cell revealing the presence of distinct miRNA expression profile changes between the two follicular stages at day 7 than day 3 of the estrous cycle. In addition, unlike the SF, marked temporal miRNA expression dynamics was observed in DF groups between day 3 and 7 of the estrous cycle. Target gene prediction and pathway analysis revealed that major signaling associated with follicular development including Wnt signaling, TGF-beta signaling, oocyte meiosis and GnRH signaling were affected by differentially expressed miRNAs. Thus, this study highlights the miRNA expression patterns of granulosa cells in subordinate and dominant follicles that could be associated with follicular recruitment, selection and dominance during the early luteal phase of the bovine estrous cycle.     

Estrous behavior in dairy cows: identification of underlying mechanisms and gene functions

Selection in dairy cattle for a higher milk yield has coincided with declined fertility. One of the factors is reduced expression of estrous behavior. Changes in systems that regulate the estrous behavior could be manifested by altered gene expression. This literature review describes the current knowledge on mechanisms and genes involved in the regulation of estrous behavior. The endocrinological regulation of the estrous cycle in dairy cows is well described. Estradiol (E2) is assumed to be the key regulator that synchronizes endocrine and behavioral events. Other pivotal hormones are, for example, progesterone, gonadotropin releasing hormone and insulin-like growth factor-1. Interactions between the latter and E2 may play a role in the unfavorable effects of milk yield-related metabolic stress on fertility in high milk-producing dairy cows. However, a clear understanding of how endocrine mechanisms are tied to estrous behavior in cows is only starting to emerge. Recent studies on gene expression and signaling pathways in rodents and other animals contribute to our understanding of genes and mechanisms involved in estrous behavior. Studies in rodents, for example, show that estrogen-induced gene expression in specific brain areas such as the hypothalamus play an important role. Through these estrogen-induced gene expressions, E2 alters the functioning of neuronal networks that underlie estrous behavior, by affecting dendritic connections between cells, receptor populations and neurotransmitter releases. To improve the understanding of complex biological networks, like estrus regulation, and to deal with the increasing amount of genomic information that becomes available, mathematical models can be helpful. Systems biology combines physiological and genomic data with mathematical modeling. Possible applications of systems biology approaches in the field of female fertility and estrous behavior are discussed.     

Reproductive cycles in pigs

The oestrous cycle in pigs spans a period of 18-24 days. It consists of a follicular phase of 5-7 days and a luteal phase of 13-15 days. During the follicular phase, small antral follicles develop into large, pre-ovulatory follicles. Being a polytocous species, the pig may ovulate from 15 to 30 follicles, depending on age, nutritional status and other factors. During the luteal phase, follicle development is less pronounced, although there is probably a considerable turnover of primordial to early antral follicles that fail to further develop due to progesterone inhibition of gonadotrophic hormones. Nevertheless, formation of the early antral follicle pool during this stage probably has a major impact on follicle dynamics in the follicular phase in terms of number and quality of follicles. Generally, gilts are mated at their second or third estrous cycle after puberty. After farrowing, pigs experience a lactational anoestrus period, until they are weaned and the follicular phase is initiated, resulting in oestrus and ovulation 4-7 days after weaning. This paper describes the major endocrine processes during the follicular and luteal phases that precede and follow ovulation. The role of nutrition and metabolic status on these processes are briefly discussed.     

Changes in ovarian follicles following acute infection with bovine viral diarrhea virus

Bovine viral diarrhea virus (BVDV) has been associated with several reproductive problems in cattle, including poor fertility, early embryonic deaths, abortion and congenital anomalies. Little is known about the cause of poor fertility in cows acutely infected with BVDV. The purpose of this study was to identify changes in ovarian function following acute infection with noncytopathic BVDV. The ovaries of 5 BVDV sero-negative and virus-negative pubertal heifers were monitored daily for 4 consecutive estrous cycles. The position and diameter of all follicles (> 5 mm) and luteal structures were recorded. Daily plasma samples were collected to measure peripheral progesterone and estradiol levels. Each heifer was infected intranasally with noncytopathic BVDV following ovulation of the second estrous cycle. The maximum diameter and growth rate of dominant anovulatory and ovulatory follicles were significantly reduced following acute BVDV infection. Similarly, the number of subordinate follicles associated with both the anovulatory and ovulatory follicle was reduced following infection. There were no significant differences in other follicle or luteal dynamic parameters or in peripheral progesterone or estradiol levels. Ovarian follicular growth was different during the first 2 estrous cycles following acute infection with BVDV when compared with the 2 estrous cycles preceding infection. These differences may be important in explaining reduced fertility in herds with acute BVDV infection.     

Alterations in peripheral concentrations of inhibin A in cattle studied using a time-resolved immunofluorometric assay: relationship with estradiol and follicle-stimulating hormone in various reproductive conditions

The aims of this study were to develop a sensitive and specific assay for bovine inhibin A using europium and to investigate the endocrine role of inhibin A in various reproductive conditions by characterizing the relationship between profiles of inhibin A, FSH, and estradiol and follicle growth during the postpartum period, during the intact estrous cycle, and in cows with follicular cysts. The time-resolved immunofluorometric assay (Tr-IFMA) for bovine inhibin A, using purified polyclonal antibodies to alpha and beta(A) subunits, was specific for bovine inhibin A and did not cross-react with bovine activin A, activin AB, activin B, pro-alphaC or human recombinant inhibin B. The detection limit of the IFMA was 3.3 pg/ml expressed in terms of bovine 32-kDa inhibin A. Dose-response curves of plasma samples obtained from intact and FSH-stimulated cows and cystic cows were parallel to the standard without any preassay processing of samples. Plasma inhibin A levels increased (P < 0.01) concomitant with emergence of nonovulatory or ovulatory follicular waves during the postpartum period. In cystic cows, plasma inhibin A was sustained at high levels for a longer period, associated with growth of persistent dominant follicles. The highest levels of inhibin A were noted during the growth phase of normal and persistent dominant follicles; however, inhibin A levels declined (P < 0.01) as these dominant follicles ceased to grow or ovulated. An inverse relationship between patterns of plasma inhibin A and FSH existed during each follicular wave in the three physiologic conditions. Increases in plasma inhibin A levels were associated with increases in plasma estradiol levels during most follicular waves; however, there was no increase in plasma estradiol level and no relationship between patterns of estradiol and FSH during follicular waves observed during the early postpartum period or midluteal phase of the estrous cycle. In conclusion, the Tr-IFMA does not require pretreatment of samples and can be used for precise measurement of bovine inhibin A without interference with free inhibin alpha subunits. Inhibin A, produced primarily during growth of the dominant follicle, functions as a negative feedback regulator for FSH secretion throughout the postpartum period and the estrous cycle, whereas estradiol appears to have a minor role in regulation of FSH compared with inhibin A, especially during the early postpartum period and midluteal phase of the estrous cycle. The results also indicate that a persistent dominant follicle sustains inhibin A production for a longer period than the dominant follicle emerging in the estrous cycle and establishes long-term dominance by suppressing emergence of a new follicular wave.     

Repeat breeding in dairy heifers: follicular dynamics and estrous cycle characteristics in relation to sexual hormone patterns

Repeat breeding occurs at an incidence of 10% in the Swedish dairy cow population.Evidence is available for a hormonal asynchrony around estrus in repeat-breeder heifers (RBH). This asynchrony seems to be the underlying cause for a series of dysfunctions such as prolonged standing estrus and delayed ovulation, leading to fertilization failure. For further determinations of repeat-breeder estrous cycle characteristics, seven strictly selected RBH and six virgin heifers (VH) were studied during 3-7 consecutive cycles, with particular attention paid to the estrous period. Follicular dynamics were studied by ultrasonography and related to estrous behavior and pattern of sexual hormones (progesterone, estradiol-17beta, and LH) in peripheral circulation. Mean group data were compared and a classification model was designed. The most prominent findings for RBH were prolonged duration of estrus, delayed LH peak, prolonged lifespan of the preovulatory follicle, and a late postovulatory rise in plasma progesterone. There was also a strong tendency for peri-ovulatory suprabasal progesterone levels in RBH. It is suspected that these deviations cause changes in the microenvironment of the preovulatory follicle, negatively affecting the final maturation of the oocyte. The heterogeneity of the RBH group underlines the multifactorial cause of the repeat-breeder syndrome. The VH formed a homogenous group with data varying within physiological limits. A classification model based on three characteristic variables managed to identify 81% of the VH and 79% of the RBH correctly. Results from this study propose that some heifers have general, consistent problems in synchronizing estrous events, displayed as varying symptoms in the course of consecutive estrous cycles. These subfertile animals could be classified as repeat-breeders.     

Sex hormone regulation of anti-bacterial activity in rat uterine secretions and apical release of anti-bacterial factor(s) by uterine epithelial cells in culture

In mature female rats, sex hormones regulate the reproductive (estrous) cycle to optimize mating and fertility. During the part of the estrous cycle when mating occurs, and when estrogen is the dominant sex hormone, the uterus is susceptible to infection with bacteria that can be deleterious for survival and fertility. The present study investigated whether sex hormones regulate innate immunity in the female reproductive tract by affecting the secretion of an anti-bacterial factor(s) in the rat uterus. Uterine fluids from intact rats at the proestrous stage of the estrous cycle significantly inhibited Staphylococcus aureus growth. When ovariectomized rats were treated with estradiol, anti-bacterial activity against both S. aureus and Escherichia coli increased in uterine secretions with hormone treatment. In contrast, rats injected with either progesterone and estradiol or progesterone alone displayed no bactericidal activity indicating that progesterone reversed the stimulatory effect of estradiol on anti-bacterial activity. In other studies, isolated uterine epithelial cells from intact animals were grown to confluence and high transepithelial resistance on cell inserts. Analysis of apical secretions indicated that a soluble factor(s) is released by polarized epithelial cells which inhibits bacterial growth. These results demonstrate that sex hormones influence the presence of a broad-spectrum bactericidal factor(s) in luminal secretions of the rat uterus. Further these studies suggest that epithelial cells which line the uterine lumen are a primary source of anti-bacterial activity.     

Modelling Hair Follicle Growth Dynamics as an Excitable Medium

The hair follicle system represents a tractable model for the study of stem cell behaviour in regenerative adult epithelial tissue. However, although there are numerous spatial scales of observation (molecular, cellular, follicle and multi follicle), it is not yet clear what mechanisms underpin the follicle growth cycle. In this study we seek to address this problem by describing how the growth dynamics of a large population of follicles can be treated as a classical excitable medium. Defining caricature interactions at the molecular scale and treating a single follicle as a functional unit, a minimal model is proposed in which the follicle growth cycle is an emergent phenomenon. Expressions are derived, in terms of parameters representing molecular regulation, for the time spent in the different functional phases of the cycle, a formalism that allows the model to be directly compared with a previous cellular automaton model and experimental measurements made at the single follicle scale. A multi follicle model is constructed and numerical simulations are used to demonstrate excellent qualitative agreement with a range of experimental observations. Notably, the excitable medium equations exhibit a wider family of solutions than the previous work and we demonstrate how parameter changes representing altered molecular regulation can explain perturbed patterns in Wnt over-expression and BMP down-regulation mouse models. Further experimental scenarios that could be used to test the fundamental premise of the model are suggested. The key conclusion from our work is that positive and negative regulatory interactions between activators and inhibitors can give rise to a range of experimentally observed phenomena at the follicle and multi follicle spatial scales and, as such, could represent a core mechanism underlying hair follicle growth.     

A complex mathematical model of the human menstrual cycle

Despite the fact that more than 100 million women worldwide use birth control pills and that half of the world's population is concerned, the menstrual cycle has so far received comparatively little attention in the field of mathematical modeling. The term menstrual cycle comprises the processes of the control system in the female body that, under healthy circumstances, lead to ovulation at regular intervals, thus making reproduction possible. If this is not the case or ovulation is not desired, the question arises how this control system can be influenced, for example, by hormonal treatments. In order to be able to cover a vast range of external manipulations, the mathematical model must comprise the main components where the processes belonging to the menstrual cycle occur, as well as their interrelations. A system of differential equations serves as the mathematical model, describing the dynamics of hormones, enzymes, receptors, and follicular phases. Since the processes take place in different parts of the body and influence each other with a certain delay, passing over to delay differential equations is deemed a reasonable step. The pulsatile release of the gonadotropin-releasing hormone (GnRH) is controlled by a complex neural network. We choose to model the pulse time points of this GnRH pulse generator by a stochastic process. Focus in this paper is on the model development. This rather elaborate mathematical model is the basis for a detailed analysis and could be helpful for possible drug design.     

Transrectal Doppler sonography of uterine blood flow

Transrectal Doppler ultrasound was used for the noninvasive investigation of uterine blood flow in cows. Both the left and right Aa. uterinae were scanned to obtain blood flow velocity waveforms over 2 consecutive estrous cycles. Blood flow was reflected by the resistance index (RI) and the time-averaged maximum velocity (TAMV). Intra-observer reproducibility of Doppler measurements was evaluated. The intra-class correlation coefficient (Intra-CC) was 0.97 for the RI and 0.95 for TAMV. While RI values did not differ between the left and right A. uterina (P > 0.05), differences in TAMV occurred between both vessels in 2 cows. These differences were not related to the ovary bearing the dominant follicle or to the corpus luteum (P < 0.001). As in all cows, changes of RI and TAMV values between the left and right artery during the estrous cycle were correlated (correlation coefficient r > 0.72; P < 0.0001); the mean values of both sides were used for subsequent analyses. Variance component estimates for the effect of cow on RI and TAMV were 8 and 13% and for the influence of day of estrous cycle they were 70 and 47%, respectively (P <0.0001). Between estrous cycles no significant differences could be measured within cows (P > 0.05). The highest RI and lowest TAMV values occurred on Day 0 (= day of ovulation) and Day 1, while the lowest RI and highest TAMV values were measured between Days -3 and -1 of the estrous cycle, respectively. There was a positive correlation between TAMV and estrogen concentrations and a negative correlation between RI and plasma estrogen levels. Plasma progesterone levels and TAMV were negatively correlated, but no correlation could be measured (P > 0.05) between RI values and plasma progesterone concentrations. While there were no differences in plasma concentrations of estrogens and progesterone between estrous cycles within cows, the levels of these hormones differed between cows. The results show that transrectal Doppler sonography is a useful, noninvasive method for examining uterine blood flows in cows. If there is an influence of uterine perfusion on fertility in cows its role needs further investigation.     

Mitogenic factors of corpora lutea

The mammalian corpus luteum (CL), which plays a central role in the reproductive process because of its production of hormones such as progesterone, appears to be an exceptionally dynamic organ. Its rate of growth and development are extremely rapid and, even when the CL is functionally mature, its rate of cell turnover remains relatively high. Associated with this high rate of cell turnover, the mature CL receives the greatest blood supply per unit tissue of any organ, and also exhibits a relatively high metabolic rate. Although numerous growth factors have been identified in luteal tissue, their role in growth and differentiation of this dynamic organ remains unclear. Recently, while attempting to identify mitogenic factors of ovine and bovine CL, we have found that they produce several mitogens during the estrous cycle as well as pregnancy. The majority of these luteal-derived mitogenic factors are heparin-binding, and although some may represent previously identified factors, several appear to be novel heparin-binding growth factors. Isolation and purification of mitogenic factors produced by the CL will enable us to determine their roles in luteal growth, development and differentiated function, which will contribute to our understanding not only of the regulation of fertility but also of tissue growth and development in general.     

The preovulatory rise of ovarian ornithine decarboxylase is required for progesterone secretion by the corpus luteum

Ovarian progesterone secretion during the diestrus stage of the estrous cycle is produced by luteal cells derived from granulosa and thecal cells after the differentiation process that follows ovulation. Our results show that blockade of the preovulatory rise of ovarian ornithine decarboxylase (ODC), a key enzyme in polyamine biosynthesis, by treatment with the specific inhibitor alpha-difluoromethylornithine (DFMO) leads to a significant decrease in the ovarian progesterone content and a dramatic fall in the plasma levels of this hormone during the following diestrus. The same inhibition was produced in spite of the fact that both luteinizing and follicle stimulating hormones were given concomitantly with DFMO. On the other hand, the acute rise in the plasma progesterone levels observed after administration of human chorionic gonadotropin to mice at different periods of the estrous cycle was not affected by DFMO administration. Our results indicate that although elevated levels of ODC are not required for acute ovarian steroidogenesis, the preovulatory peak of ovarian ODC activity observed in the evening of proestrus may be critical for the establishment of a constitutive steroidogenic pathway and progesterone secretion by the corpus luteum during the diestrus stage of the murine estrous cycle.     

The ability of subordinate follicles of the second follicular wave to become dominant is lost by day 15 of the estrous cycle in cattle

Generally, unilateral ovariectomy before a critical period in the latter part of the estrous cycle induces a transitory increase in plasma FSH, which causes subordinate follicles to develop and maintain ovulation rates characteristic of the species. A limiting period for subordinate follicles to assume dominance and from which ovulation occurs has not been shown for cattle. Growth and/or regression of subordinate follicles were characterized following removal of the dominant follicle at different days of the luteal phase of the estrous cycle in cattle in this study. In the mid-luteal phase (Day 13 or 15), the ovary with the dominant follicle of the second wave was ablated via unilateral ovariectomy; the corpus luteum also was removed. In the late luteal phase (Day 17 or 19), the dominant follicle was ablated with an ultrasonically guided 20 gauge needle. When the dominant follicle was removed on Day 13, the largest subordinate follicle of the second wave of follicular development became dominant and ovulation occurred from this follicle in 4 of 4 animals. However, when the dominant follicle was removed on Day 15, 17 or 19, a new wave of follicular development was induced in 14 of 15 animals. Moreover, the recovered subordinate follicle of the second wave of follicular development had similar growth characteristics to naturally occurring dominant follicles. In conclusion, the subordinate follicle in the second follicular wave in cattle retained the ability to become dominant, but this ability was lost by Day 15 of the estrous cycle. However, cattle then were able to maintain ovulation by developing a new wave of follicular growth.     

Follicular fluid rheology and the duration of the ovulatory process

The fluid dynamics of ovulation were investigated to understand the mechanical role of follicular fluid in oocyte release. A set of equations describing the flow of fluid from an evacuating follicle was derived from basic principles. These equations demonstrate that, subject to assumptions about the available pressure differential and the source of the expulsive force, the size and shape of the ovulatory orifice have the largest influences on the rate of fluid loss, although the viscosity of the fluid is also an important variable. A thorough rheological examination of pig, bovine and human follicular fluids, performed using a cone-plate viscometer, demonstrated that these fluids have complex, non-Newtonian characteristics. The fluids also undergo time-dependent and spontaneous changes in viscosity at constant shear rates; some fluids were subject to coagulation-like events. Viscosity characteristics were unrelated to broad parameters of follicle development. The models used representative viscosity values to demonstrate that variations in the rate and duration of follicle evacuation, as observed by ultrasonography, could be explained largely by variations in fluid viscosity and the characteristics of the ovulatory orifice.     

Plasma concentrations of inhibin a and follicle-stimulating hormone differ between cows with two or three waves of ovarian follicular development in a single estrous cycle

Patterns of ovarian follicle development were monitored daily in Holstein-Friesian cows that had two (n = 4) or three (n = 4) waves of ovarian follicle development during a single estrous cycle. The plasma from daily blood samples was used in assays for inhibin A, FSH, progesterone, and estradiol-17beta. Mean cycle lengths for cows with two and three waves were 21.8 and 25.3 days, respectively (P < 0.02). Although the average number of follicles >3-mm diameter on each pair of ovaries was similar for two- and three-wave cows on Days 2, 3, and 4 (Day 0 = day of ovulation; 8.6 vs. 9.6 follicles), there were more follicles >6-mm diameter on the ovaries of cows with two waves on Days 3 and 4. This difference was associated with a shorter interval from wave emergence to peak concentrations of inhibin A during the first wave in two-wave cows (2.0 vs. 3.8 days; P = 0.03) and with higher peak concentrations (474 vs. 332 pg/ml; P = 0.03). Differences in peak FSH concentrations were not significant (1.7 vs. 1.3 ng/ml; P = 0.10) and were inversely related to inhibin A concentrations. The peak concentrations of inhibin A and FSH in the second nonovulatory wave in the three-wave cows were similar to the low concentrations measured in the first wave (292 vs. 332 pg/ml of inhibin A, 1.3 vs. 1.3 ng/ml of FSH; P > 0.20). Average peak concentrations of inhibin A and FSH were similar during the ovulatory wave for cows with either two or three waves in a cycle (432 vs. 464 pg/ml of inhibin A, 2.3 vs. 2.1 ng/ml of FSH; P > 0.3). The lower concentrations of FSH during the emergence of the first follicular wave in cows with three-wave cycles may have reduced the rate of development of some of the follicles and reduced the concentrations of inhibin A. This pattern of lower concentrations of FSH and inhibin A was repeated in the second nonovulatory wave but not in the ovulatory wave. Subtle differences in the concentrations of these two hormones may underlie the mechanism that influences the number of waves of ovarian follicle development that occur during the bovine estrous cycle.