This dissertation investigates synchronization properties of slow-fast oscillators inspired from neuroendocrinology and neuronal dynamics, focusing on the effects of canard phenomena and dynamic bifurcations on the collective behavior. We start from a 4-dimensional system which accounts for the qualitative and quantitative dynamical features of the secretion pattern of the neurohormone GnRH (gonadotropin releasing hormone) along a whole ovarian cycle. This model involves 2 FitzHugh-Nagumo oscillators with different timescales. Unidirectional coupling from the slow oscillator (representing the mean-field activity of a population of regulating neurons) to the fast oscillator (representing the mean-field activity of a population of the secreting neurons) gives a three timescale structure. The behavior of the fast oscillator is characterized by an alternation between a relaxation cycle and a quasi-stationary state which introduces canard-mediated transitions in the model; these transitions have a strong impact on the secretion pattern of the 4-dimensional system. We make a first step forward in multiscale modeling (in space) of the GnRH system, namely, we extend the original system to 6 dimensions by considering two distinct subpopulations of secreting neurons receiving the same signal from the regulating neurons. This step allows us to enrich further the GnRH secretion pattern while keeping a compact dynamic framework and preserving the sequence of neurosecretory events captured by the 4-dimensional model, both qualitatively and quantitatively. An initial analysis of the extended 6-dimensional GnRH model is presented in Chapter 2, where we prove using a 5D minimal model the existence of canard trajectories in coupled systems with folded singularities. Coupling causes separation of trajectories corresponding to each secretor by driving them to different sides of the maximal canard (associated with either a folded-node or a folded-saddle singularity). We explore the impact of the relationship between canard structures and coupling on the collective secretion pattern of the 6-dimensional model. We identify two different sources of canard-mediated (de)synchronization in the secretory events, which depend on the type of underlying folded singularity. In Chapter 3, we attempt to model complex behaviors of the GnRH secretion not captured by the 4-dimensional model, namely, a surge with 2 bumps and partial desynchronization before the surge, by using the 6-dimensional model previously constructed. Regulatory-dependent asymmetric coupling functions and heterogeneity in the secretor subpopulations are essential for obtaining such a 2-bump surge. During the pulsatile regime, we find that the slowly varying regulatory signal causes a dynamic bifurcation, which is responsible for loss of synchrony in asymmetrically coupled nonidentical secretors. We introduce analytic and numerical tools to shape and quantify the additional features embedded within the whole secretion pattern.