Preimplantation genetic testing: A remarkable history of pioneering, technical challenges, innovations, and ethical considerations.

Preimplantation genetic testing (PGT) has emerged as a powerful companion to assisted reproduction technologies. The origins and history of PGT are reviewed here, along with descriptions of advances in molecular assays and sampling methods, their capabilities, and their applications in preventing genetic diseases and enhancing pregnancy outcomes. Additionally, the potential for increasing accuracy and genome coverage is considered, as well as some of the emerging ethical and legislative considerations related to the expanding capabilities of PGT.

Predictive modeling of oocyte maternal mRNA features for five mammalian species reveals potential shared and species-restricted regulators during maturation.

Oocyte maturation is accompanied by changes in abundances of thousands of mRNAs, many degraded and many preferentially stabilized. mRNA stability can be regulated by diverse features including GC content, codon bias, and motifs within the 3'-untranslated region (UTR) interacting with RNA binding proteins (RBPs) and miRNAs. Many studies have identified factors participating in mRNA splicing, bulk mRNA storage, and translational recruitment in mammalian oocytes, but the roles of potentially hundreds of expressed factors, how they regulate cohorts of thousands of mRNAs, and to what extent their functions are conserved across species has not been determined. We performed an extensive in silico cross-species analysis of features associated with mRNAs of different stability classes during oocyte maturation (stable, moderately degraded, and highly degraded) for five mammalian species. Using publicly available RNA sequencing data for germinal vesicle (GV) and MII oocyte transcriptomes, we determined that 3'-UTR length and synonymous codon usage are positively associated with stability, while greater GC content is negatively associated with stability. By applying machine learning and feature selection strategies, we identified RBPs and miRNAs that are predictive of mRNA stability, including some across multiple species and others more species-restricted. The results provide new insight into the mechanisms regulating maternal mRNA stabilization or degradation.NEW & NOTEWORTHY Conservation across species of mRNA features regulating maternal mRNA stability during mammalian oocyte maturation was analyzed. 3'-Untranslated region length and synonymous codon usage are positively associated with stability, while GC content is negatively associated. Just three RNA binding protein motifs were predicted to regulate mRNA stability across all five species examined, but associated pathways and functions are shared, indicating oocytes of different species arrive at comparable physiological destinations via different routes.

Ovarian stimulation with excessive FSH doses causes cumulus cell and oocyte dysfunction in small ovarian reserve heifers.

Excessive FSH doses during ovarian stimulation in the small ovarian reserve heifer (SORH) cause premature cumulus expansion and follicular hyperstimulation dysgenesis (FHD) in nearly all ovulatory-size follicles with predicted disruptions in cell-signaling pathways in cumulus cells and oocytes (before ovulatory hCG stimulation). These observations support the hypothesis that excessive FSH dysregulates cumulus cell function and oocyte maturation. To test this hypothesis, we determined whether excessive FSH-induced differentially expressed genes (DEGs) in cumulus cells identified in our previously published transcriptome analysis were altered independent of extreme phenotypic differences observed amongst ovulatory-size follicles, and assessed predicted roles of these DEGs in cumulus and oocyte biology. We also determined if excessive FSH alters cumulus cell morphology, and oocyte nuclear maturation before (premature) or after an ovulatory hCG stimulus or during IVM. Excessive FSH doses increased expression of 17 cumulus DEGs with known roles in cumulus cell and oocyte functions (responsiveness to gonadotrophins, survival, expansion, and oocyte maturation). Excessive FSH also induced premature cumulus expansion and oocyte maturation but inhibited cumulus expansion and oocyte maturation post-hCG and diminished the ability of oocytes with prematurely expanded cumulus cells to undergo IVF or nuclear maturation during IVM. Ovarian stimulation with excessive FSH is concluded to disrupt cumulus cell and oocyte functions by inducing premature cumulus expansion and dysregulating oocyte maturation without an ovulatory hCG stimulus yielding poor-quality cumulus-oocyte complexes that may be incorrectly judged morphologically as suitable for IVF during ART.

Preimplantation embryo gene expression: 56 years of discovery, and counting.

The biology of preimplantation embryo gene expression began 56 years ago with studies of the effects of protein synthesis inhibition and discovery of changes in embryo metabolism and related enzyme activities. The field accelerated rapidly with the emergence of embryo culture systems and progressively evolving methodologies that have allowed early questions to be re-addressed in new ways and in greater detail, leading to deeper understanding and progressively more targeted studies to discover ever more fine details. The advent of technologies for assisted reproduction, preimplantation genetic testing, stem cell manipulations, artificial gametes, and genetic manipulation, particularly in experimental animal models and livestock species, has further elevated the desire to understand preimplantation development in greater detail. The questions that drove enquiry from the earliest years of the field remain drivers of enquiry today. Our understanding of the crucial roles of oocyte-expressed RNA and proteins in early embryos, temporal patterns of embryonic gene expression, and mechanisms controlling embryonic gene expression has increased exponentially over the past five and a half decades as new analytical methods emerged. This review combines early and recent discoveries on gene regulation and expression in mature oocytes and preimplantation stage embryos to provide a comprehensive understanding of preimplantation embryo biology and to anticipate exciting future advances that will build upon and extend what has been discovered so far.

Effects of early lactation body condition loss in dairy cows on serum lipid profiles and on oocyte and cumulus cell transcriptomes.

The objective of this study was to determine the effect of early lactation body condition (BC) loss in multiparous dairy cows on serum lipids and the effect of these changes on oocyte and cumulus cell transcriptomes. Body condition loss in dairy cattle after parturition is associated with reduced fertility and increased pregnancy loss. The complex interplay between BC, nutrition, dry matter intake, milk production, and time of calving has presented a barrier to understanding mechanisms leading to reduced fertility. We identified cows that lost BC (L group; n = 10) or maintained or gained BC (M/G group; n = 8) during the first 27 to 33 d in milk and investigated changes in serum fatty acids and oocyte and cumulus cell transcriptomes at 75 to 81 d in milk. The L group had increased serum levels of nonesterified fatty acids and mead acid, and reduced serum levels of petroselaidic acid and behenic acid. Transcriptome analyses revealed 38 differentially expressed genes (DEG) in oocytes and 71 DEG in cumulus cells of L (n = 3) compared with M/G group (n = 3). Network analysis connected serum fatty acid changes to downstream effects including reduced inflammatory response and mitochondrial membrane depolarization, increased production of reactive oxygen species, and functions related to fatty acid metabolism and cytoplasmic organization in oocytes. These effects were associated with predicted effects on signaling in oocytes through calcium, insulin, O-GlcNAcase (OGA), fibroblast growth factor receptor 4 (FGF4R), peroxisome proliferator activated receptor gamma coactivator 1 α (PPARGC1A), and phospholipase D2 (PLD2) pathways, with a connection to the cumulus cell via calcium signaling. These results connect BC loss following parturition to changes in serum lipid levels, and changes potentially affecting oocyte quality; thus, these results provide new insight into mechanism of reduced fertility.

Follicular Hyperstimulation Dysgenesis: New Explanation for Adverse Effects of Excessive FSH in Ovarian Stimulation.

High follicle-stimulating hormone (FSH) doses during ovarian stimulation protocols for assisted reproductive technologies (ART) are detrimental to ovulatory follicle function and oocyte quality. However, the mechanisms are unclear. In a small ovarian reserve heifer model, excessive FSH doses lead to phenotypic heterogeneity of ovulatory size follicles, with most follicles displaying signs of premature luteinization and a range in severity of abnormalities. By performing whole transcriptome analyses of granulosa cells, cumulus cells, and oocytes from individual follicles of animals given standard or excessive FSH doses, we identified progressive changes in the transcriptomes of the 3 cell types, with increasing severity of follicular abnormality with the excessive doses. The granulosa and cumulus cells each diverged progressively from their normal phenotypes and became highly similar to each other in the more severely affected follicles. Pathway analysis indicates a possible dysregulation of the final stages of folliculogenesis, with processes characteristic of ovulation and luteinization occurring concurrently rather than sequentially in the most severely affected follicles. These changes were associated with disruptions in key pathways in granulosa and cumulus cells, which may account for previously reported reduced estradiol production, enhanced progesterone and oxytocin production and diminished ovulation rates. Predicted deficiencies in oocyte survival, stress response, and fertilization suggest likely reductions in oocyte health, which could further compromise oocyte quality and ART outcomes.

Limitations in use of ovarian reserve biomarkers to predict the superovulation response in small ovarian reserve heifers.

High FSH doses during superovulation of heifers with a small ovarian reserve increase the number of dysfunctional ovulatory-size follicles that do not ovulate in response to human chorionic gonadotropin (hCG). Thus, anti-Müllerian hormone (AMH) and antral follicle count (AFC), two well-established biomarkers of responsiveness of individuals to superovulation, are hypothesized to be positively linked to number of dysfunctional ovulatory-size follicles developing in response to superovulation with high FSH doses. To test this hypothesis, heifers with a small ovarian reserve were stimulated beginning on Day 1 of the estrous cycle with twice daily treatments for 4 days with each of four Folltropin-V (FSH) doses (35 IU, 70 IU (industry standard), 140 IU, or 210 IU) followed by prostaglandin F2α to regress corpora lutea (CL) from the previous estrous cycle and hCG to induce ovulation. Ovulatory-size follicles were classified as functional or dysfunctional based on whether they ovulated and formed CL in response to hCG. FSH dose did not impact the relationship between AMH, AFC and the number of functional or dysfunctional ovulatory-size follicles developing in response to superovulation. Thus, data from the four superovulations were averaged for each heifer. AMH and AFC were positively associated with the subsequent number of functional and dysfunctional ovulatory-size follicles and the proportion of ovulatory-size follicles that are dysfunctional after superovulation. Because measurements of AMH concentration and AFC predict the number but not functionality of ovulatory-size follicles, which may also impact oocyte quality, these ovarian reserve biomarkers are concluded to be unlikely useful to improve IVF or embryo transfer outcomes in heifers with a small ovarian reserve.

Excessive follicle-stimulating hormone during ovarian stimulation of cattle may induce premature luteinization of most ovulatory-size follicles†.

High follicle-stimulating hormone (FSH) doses during ovarian stimulation are detrimental to ovulatory follicle function and decrease live birth rate in cattle and women. However, the mechanism whereby excessive FSH causes ovarian dysfunction is unknown. This study tested the hypothesis that excessive FSH during ovarian stimulation induces premature luteinization of ovulatory-size follicles. Small ovarian reserve heifers were injected twice daily for 4 days with 70 IU (N = 7 heifers) or 210 IU (N = 6 heifers) Folltropin-V [commercial FSH-enriched preparation of porcine pituitary glands with minor (<1%) luteinizing hormone (LH) contamination, cpFSH]. Ovulatory-size (≥10 mm) follicles were excised from ovaries after the last cpFSH injection and hormone concentrations in follicular fluid (FF) were determined using ELISA. Luteinization was monitored by assessing cumulus cell-oocyte complex (COC) morphology and measuring concentrations of estradiol (E), progesterone (P), and oxytocin (O) in FF. COCs were classified as having compact (cCOC) or expanded (eCOC) cumulus cell layers, and as estrogen-active (E:P in FF ≥1), estrogen-inactive (EI, E:P in FF ≤1 > 0.1), or extreme-estrogen-inactive (EEI, E:P in FF ≤0.1). A high proportion (72%) of ovulatory-size follicles in 210 IU, but not 70 IU, dose heifers displayed eCOCs. The high doses also produced higher proportions of EI or EEI follicles which had lower E:P ratio and/or E but higher P and/or O concentrations compared with the 70 IU dose heifers. In conclusion, excessive cpFSH doses during ovarian stimulation may induce premature luteinization of most ovulatory-size follicles in heifers with small ovarian reserves.

Cross-species meta-analysis of transcriptome changes during the morula-to-blastocyst transition: metabolic and physiological changes take center stage.

The morula-to-blastocyst transition (MBT) culminates with formation of inner cell mass (ICM) and trophectoderm (TE) lineages. Recent studies identified signaling pathways driving lineage specification, but some features of these pathways display significant species divergence. To better understand evolutionary conservation of the MBT, we completed a meta-analysis of RNA sequencing data from five model species and ICMTE differences from four species. Although many genes change in expression during the MBT within any given species, the number of shared differentially expressed genes (DEGs) is comparatively small, and the number of shared ICMTE DEGs is even smaller. DEGs related to known lineage determining pathways (e.g., POU5F1) are seen, but the most prominent pathways and functions associated with shared DEGs or shared across individual species DEG lists impact basic physiological and metabolic activities, such as TCA cycle, unfolded protein response, oxidative phosphorylation, sirtuin signaling, mitotic roles of polo-like kinases, NRF2-mediated oxidative stress, estrogen receptor signaling, apoptosis, necrosis, lipid and fatty acid metabolism, cholesterol biosynthesis, endocytosis, AMPK signaling, homeostasis, transcription, and cell death. We also observed prominent differences in transcriptome regulation between ungulates and nonungulates, particularly for ICM- and TE-enhanced mRNAs. These results extend our understanding of shared mechanisms of the MBT and formation of the ICM and TE and should better inform the selection of model species for particular applications.

Essential shared and species-specific features of mammalian oocyte maturation-associated transcriptome changes impacting oocyte physiology.

Oogenesis is a complex process resulting in the production of a truly remarkable cell-the oocyte. Oocytes execute many unique processes and functions such as meiotic segregation of maternal genetic material, and essential life-generating functions after fertilization including posttranscriptional support of essential homeostatic and metabolic processes, and activation and reprogramming of the embryonic genome. An essential goal for understanding female fertility and infertility in mammals is to discover critical features driving the production of quality oocytes, particularly the complex regulation of oocyte maternal mRNAs. We report here the first in-depth meta-analysis of oocyte maturation-associated transcriptome changes, using eight datasets encompassing 94 RNAseq libraries for human, rhesus monkey, mouse, and cow. A majority of maternal mRNAs are regulated in a species-restricted manner, highlighting considerable divergence in oocyte transcriptome handling during maturation. We identified 121 mRNAs changing in relative abundance similarly across all four species (92 of high homology), and 993 (670 high homology) mRNAs regulated similarly in at least three of the four species, corresponding to just 0.84% and 6.9% of mRNAs analyzed. Ingenuity Pathway Analysis (IPA) revealed an association of these shared mRNAs with many shared pathways and functions, most prominently oxidative phosphorylation and mitochondrial function. These shared functions were reinforced further by primate-specific and species-specific differentially expressed genes (DEGs). Thus, correct downregulation of mRNAs related to oxidative phosphorylation and mitochondrial function is a major shared feature of mammalian oocyte maturation.

Shared aspects of mRNA expression associated with oocyte maturation failure in humans and rhesus monkeys indicating compromised oocyte quality.

Oocyte maturation failure observed in assisted reproduction technology (ART) cycles can limit the number of quality oocytes obtained and present a pronounced barrier for some patients. The potential exists to use unmatured oocytes for ART through in vitro maturation. Understanding the molecular basis of oocyte maturation failure is pertinent to minimizing this loss of oocytes and considerations of whether such oocytes can be used safely for ART. We identified shared transcriptome abnormalities for rhesus monkey and human failed-to-mature (FTM) oocytes relative to healthy matured MII stage oocytes. We discovered that, although the number of shared affected genes was comparatively small, FTM oocytes in both species shared effects for several pathways and functions, including predicted activation of oxidative phosphorylation (OxPhos) with additional effects on mitochondrial function, lipid metabolism, transcription, nucleotide excision repair, endoplasmic reticulum stress, unfolded protein response, and cell viability. RICTOR emerged as a prominent upstream regulator with predicted inhibition across all analyses. Alterations in KDM5A, MTOR, MTORC1, INSR, CAB39L, and STK11 activities were implicated along with RICTOR in modulating mitochondrial activity and OxPhos. Defects in cell cycle progression were not a prominent feature of FTM oocytes. These results identify a common set of transcriptome abnormalities associated with oocyte maturation failure. While our results do not demonstrate causality, they indicate that fundamental aspects of cellular function are abnormal in FTM oocytes and raise significant concerns about the potential risks of using FTM oocytes for ART.

Negative impact of high doses of follicle-stimulating hormone during superovulation on the ovulatory follicle function in small ovarian reserve dairy heifers†.

When women with small ovarian reserves are subjected to assisted reproductive technologies, high doses of gonadotropins are linked to high oocyte and embryo wastage and low live birth rates. We hypothesized that excessive follicle-stimulating hormone (FSH) doses during superovulation are detrimental to ovulatory follicle function in individuals with a small ovarian reserve. To test this hypothesis, heifers with small ovarian reserves were injected twice daily for 4 days, beginning on Day 1 of the estrous cycle with 35, 70, 140, or 210 IU doses of Folltropin-V (FSH). Each heifer (n = 8) was superovulated using a Williams Latin Square Design. During each superovulation regimen, three prostaglandin F2α injections were given at 12-h interval, starting at the seventh FSH injection to regress the newly formed corpus luteum (CL). Human chorionic gonadotropin was injected 12 h after the last (8th) FSH injection to induce ovulation. Daily ultrasonography and blood sampling were used to determine the number and size of follicles and corpora lutea, uterine thickness, and circulating concentrations of estradiol, progesterone, and anti-Müllerian hormone (AMH). The highest doses of FSH did not increase AMH, progesterone, number of ovulatory-size follicles, uterine thickness, or number of CL. However, estradiol production and ovulation rate were lower for heifers given high FSH doses compared to lower doses, indicating detrimental effects on ovulatory follicle function.

Probing lasting cryoinjuries to oocyte-embryo transcriptome.

Clinical applications of oocytes cryopreservation include preservation of future fertility of young cancer patients, substitution of embryo freezing to avoid associated legal and ethical issues, and delaying childbearing years. While the outcome of oocyte cryopreservation has recently been improved, currently used vitrification method still suffer from increased biosafety risk and handling issues while slow freezing techniques yield overall low success. Understanding better the mechanism of cryopreservation-induced injuries may lead to development of more reliable and safe methods for oocyte cryopreservation. Using the mouse model, a microarray study was conducted on oocyte cryopreservation to identify cryoinjuries to transcriptionally active genome. To this end, metaphase II (MII) oocytes were subjected to standard slow freezing, and then analyzed at the four-cell stage after embryonic genome activation. Non-frozen four-cell embryos served as controls. Differentially expressed genes were identified and validated using RT-PCR. Embryos produced from the cryopreserved oocytes displayed 200 upregulated and 105 downregulated genes, associated with the regulation of mitochondrial function, protein ubiquitination and maintenance, cellular response to stress and oxidative states, fatty acid and lipid regulation/metabolism, and cell cycle maintenance. These findings reveal previously unrecognized effects of standard slow oocyte freezing on embryonic gene expression, which can be used to guide improvement of oocyte cryopreservation methods.

Listening to mother: Long-term maternal effects in mammalian development.

The oocyte is a complex cell that executes many crucial and unique functions at the start of each life. These functions are fulfilled by a unique collection of macromolecules and other factors, all of which collectively support meiosis, oocyte activation, and embryo development. This review focuses on the effects of oocyte components on developmental processes that occur after the initial stages of embryogenesis. These include long-term effects on genome function, metabolism, lineage allocation, postnatal progeny health, and even subsequent generations. Factors that regulate chromatin structure, genome programming, and mitochondrial function are elements that contribute to these oocyte functions.

Genotypic divergence in mouse oocyte transcriptomes: possible pathways to hybrid vigor impacting fertility and embryogenesis.

What causes hybrid vigor phenotypes in mammalian oocytes and preimplantation embryos? Answering this question should provide new insight into determinants of oocyte and embryo quality and infertility. Hybrid vigor could arise through a variety of mechanisms, many of which must operate through posttranscriptional mechanisms affecting oocyte mRNA accumulation, stability, translation, and degradation. The differential regulation of such mRNAs may impact essential pathways and functions within the oocyte. We conducted in-depth transcriptome comparisons of immature and mature oocytes of C57BL/6J and DBA/2J inbred strains and C57BL/6J × DBA/2J F1 (BDF1) hybrid oocytes with RNA sequencing, combined with novel computational methods of analysis. We observed extensive differences in mRNA expression and regulation between parental inbred strains and between inbred and hybrid genotypes, including mRNAs encoding proposed markers of oocyte quality. Unique BDF1 oocyte characteristics arise through a combination of additive dominance and incomplete dominance features in the transcriptome, with a lesser degree of transgressive mRNA expression. Special features of the BDF1 transcriptome most prominently relate to histone expression, mitochondrial function, and oxidative phosphorylation. The study reveals the major underlying mechanisms that contribute to superior properties of hybrid oocytes in a mouse model.

A New Role for SMCHD1 in Life's Master Switch and Beyond.

Structural maintenance of chromosomes flexible hinge-domain containing protein 1 (SMCHD1) has emerged as a key regulator of embryonic genome function. Its functions have now extended well beyond the initial findings of effects on X chromosome inactivation associated with lethality in female embryos homozygous for a null allele. Autosomal dominant effects impact stem cell properties as well as postnatal health. Recent studies have revealed that SMCHD1 plays an important role as a maternal effect gene that regulates the master switch of life, namely embryonic genome activation, as well as subsequent preimplantation development and term viability. These discoveries mark SMCHD1 as a major regulator linking developmental processes to adult disorders including a form of muscular dystrophy.

SMCHD1 terminates the first embryonic genome activation event in mouse two-cell embryos and contributes to a transcriptionally repressive state.

Embryonic genome activation (EGA) in mammals begins with transient expression of a large group of genes (EGA1). Importantly, entry into and exit from the 2C/EGA state is essential for viability. Dux family member genes play an integral role in EGA1 by activating other EGA marker genes such as Zscan4 family members. We previously reported that structural maintenance of chromosomes flexible hinge domain-containing protein 1 (Smchd1) is expressed at the mRNA and protein levels in mouse oocytes and early embryos and that elimination of Smchd1 expression inhibits inner cell mass formation, blastocyst formation and hatching, and term development. We extend these observations here by showing that siRNA knockdown of Smchd1 in zygotes results in overexpression of Dux and Zscan4 in two-cell embryos, with continued overexpression of Dux at least through the eight-cell stage as well as prolonged expression of Zscan4. These results are consistent with a role for SMCHD1 in promoting exit from the EGA1 state and establishing SMCHD1 as a maternal effect gene and the first chromatin regulatory factor identified with this role. Additionally, bioinformatics analysis reveals that SMCHD1 also contributes to the creation of a transcriptionally repressive state to allow correct gene regulation.

Inosine RNA modifications are enriched at the codon wobble position in mouse oocytes and eggs†.

Mammalian oocytes and eggs are transcriptionally quiescent and depend on post-transcriptional mechanisms for proper maturation. Post-transcriptional mRNA modifications comprise an important regulatory mechanism that can alter protein and miRNA recognition sites, splicing, stability, secondary structure, and protein coding. We discovered that fully grown mouse germinal vesicle oocytes and metaphase II eggs display abundant inosine mRNA modifications compared to growing oocytes from postnatal day 12 oocytes. These inosines were enriched in mRNA protein coding regions (CDS) and specifically located at the third codon base, or wobble position. Inosines, observed at lower frequencies in CDS of somatic tissues, were similarly enriched at the codon wobble position. In oocytes and eggs, inosine modifications lead primarily to synonymous changes in mRNA transcripts. Inosines may ultimately affect maternal mRNA stability by changing codon usage, thereby altering translational efficiency and translationally coupled mRNA degradation. These important observations advance our understanding of post-transcriptional mechanisms contributing to mammalian oocyte maturation.

Undernutrition and hyperandrogenism during pregnancy: Role in programming of cardiovascular disease and infertility.

Maternal nutritional status programs the development of several systems in female offspring, with effects that depend on the severity, duration, and window of development when the nutritional perturbation is imposed. On the basis of the developmental origins of health and disease concept, we hypothesize that gestational low caloric intake may induce maternal subclinical hyperandrogenism during early pregnancy and compromise cardiovascular health and fertility in the female offspring. To examine this possibility, a literature search for human and animal studies was conducted using two electronic databases, PubMed and Cochrane until April 2019 to address the following questions: (a) Do androgens have a developmental role in cardiovascular and ovarian development? (b) Is excess maternal testosterone linked to cardiovascular disease and infertility? and (c) Could early pregnancy undernutrition enhance maternal androgen production and compromise health and fertility in female offspring? The observations reviewed, establish a potential causative link between maternal undernutrition and subclinical hyperandrogenism with hypertension and reduced ovarian reserve in the progeny. Further studies in appropriate models are needed to better understand whether low energy intake and subclinical maternal hyperandrogenism during early pregnancy can negatively affect the health of the female offspring.

Temporal patterns of gene regulation and upstream regulators contributing to major developmental transitions during Rhesus macaque preimplantation development.

The preimplantation period of life in mammals encompasses a tremendous amount of restructuring and remodeling of the embryonic genome and reprogramming of gene expression. These vast changes support metabolic activation and cellular processes that drive early cleavage divisions and enable the creation of the earliest primitive cell lineages. A major question in mammalian embryology is how such vast, sweeping changes in gene expression are orchestrated, so that changes in gene expression are exactly appropriate to meet the developmental needs of the embryo over time. Using the rhesus macaque as an experimentally tractable model species closely related to the human, we combined high quality RNA-seq libraries, in-depth sequencing and advanced systems analysis to discover the underlying mechanisms that drive major changes in gene regulation during preimplantation development. We identified the major changes in mRNA population and the biological pathways and processes impacted by those changes. Most importantly, we identified 24 key upstream regulators that are themselves modulated during development and that are associated with the regulation of over 1000 downstream genes. Through their roles in extensive gene networks, these 24 upstream regulators are situated to either drive major changes in target gene expression or modify the cellular environment in which other genes function, thereby directing major developmental transitions in the preimplantation embryo. The data presented here highlight some of the specific molecular features that likely drive preimplantation development in a nonhuman primate species and provides an extensive database for novel hypothesis-driven studies.