A map of canine sequence variation relative to a Greenland wolf outgroup.

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

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.

Genome-wide methylation patterns from canine nanopore assemblies.

Recent advances in long-read sequencing have enabled the creation of reference-quality genome assemblies for multiple individuals within a species. In particular, 8 long-read genome assemblies have recently been published for the canine model (dogs and wolves). These assemblies were created using a range of sequencing and computational approaches, with only limited comparisons described among subsets of the assemblies. Here we present 3 high-quality de novo reference assemblies based upon Oxford Nanopore long-read sequencing: 2 Bernese Mountain Dogs (BD & OD) and a Cairn terrier (CA611). These breeds are of particular interest due to the enrichment of unresolved genetic disorders. Leveraging advancement in software technologies, we utilized published data of Labrador Retriever (Yella) to generate a new assembly, resulting in a ∼280-fold increase in continuity (N50 size of 91 kbp vs 25.75 Mbp). In conjunction with these 4 new assemblies, we uniformly assessed 8 existing assemblies for generalized quality metrics, sequence divergence, and a detailed BUSCO assessment. We identified a set of ∼400 conserved genes during the BUSCO analysis missing in all assemblies. Genome-wide methylation profiles were generated from the nanopore sequencing, resulting in broad concordance with existing whole-genome and reduced-representation bisulfite sequencing, while highlighting superior overage of mobile elements. These analyses demonstrate the ability of Nanopore sequencing to resolve the sequence and epigenetic profile of canine genomes.

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.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Identification of a Hypomorphic FANCG Variant in Bernese Mountain Dogs.

Bernese mountain dogs (BMDs), have an overall cancer incidence of 50%, half of which is comprised of an otherwise rare tumor, histiocytic sarcoma (HS). While recent studies have identified driver mutations in the MAPK pathway, identification of key predisposing genes has been elusive. Studies have identified several loci to be associated with predisposition to HS in BMDs, including near the MTAP/CDKN2A region, but no causative coding variant has been identified. Here we report the presence of a coding polymorphism in the gene encoding FANCG, near the MTAP/CDKN2A locus. This variant is in a conserved region of the protein and appears to be specific to BMDs. Canine fibroblasts derived from dogs homozygous for this variant are hypersensitive to cisplatin. We show this canine FANCG variant and a previously defined hypomorphic FANCG allele in humans impart similar defects in DNA repair. However, our data also indicate that this variant is neither necessary nor sufficient for the development of HS. Furthermore, BMDs homozygous for this FANCG allele display none of the characteristic phenotypes associated with Fanconi anemia (FA) such as anemia, short stature, infertility, or an earlier age of onset for HS. This is similar to findings in FA deficient mice, which do not develop overt FA without secondary genetic mutations that exacerbate the FA deficit. In sum, our data suggest that dogs with deficits in the FA pathway are, like mice, innately resistant to the development of FA.

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.

Small Molecule 20S Proteasome Enhancer Regulates MYC Protein Stability and Exhibits Antitumor Activity in Multiple Myeloma.

Despite the addition of several new agents to the armamentarium for the treatment of multiple myeloma (MM) in the last decade and improvements in outcomes, the refractory and relapsing disease continues to take a great toll, limiting overall survival. Therefore, additional novel approaches are needed to improve outcomes for MM patients. The oncogenic transcription factor MYC drives cell growth, differentiation and tumor development in many cancers. MYC protein levels are tightly regulated by the proteasome and an increase in MYC protein expression is found in more than 70% of all human cancers, including MM. In addition to the ubiquitin-dependent degradation of MYC by the 26S proteasome, MYC levels are also regulated in a ubiquitin-independent manner through the REGγ activation of the 20S proteasome. Here, we demonstrate that a small molecule activator of the 20S proteasome, TCH-165, decreases MYC protein levels, in a manner that parallels REGγ protein-mediated MYC degradation. TCH-165 enhances MYC degradation and reduces cancer cell growth in vitro and in vivo models of multiple myeloma by enhancing apoptotic signaling, as assessed by targeted gene expression analysis of cancer pathways. Furthermore, 20S proteasome enhancement is well tolerated in mice and dogs. These data support the therapeutic potential of small molecule-driven 20S proteasome activation for the treatments of MYC-driven cancers, especially MM.

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.

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.

Activating Mutations in PTPN11 and KRAS in Canine Histiocytic Sarcomas.

While the genetic contributions to the predisposition of Bernese mountain dogs (BMDs) to histiocytic sarcoma (HS) remains unclear, some insights into key genetic drivers have been gained. Our group recently reported a mutation in the PTPN11 gene (E76K). We have now identified a second missense mutation in PTPN11 (G503V), and a mutation in KRAS (Q61H) present in HS cell lines. These mutations are associated with malignancies in humans, and known to be gain-of-function mutations that result in activation of the mitogen-activated protein kinase (MAPK) pathway. The goal of the present study was to evaluate the prevalence of these mutations in a large sample of HS cases from BMDs and golden retrievers, and in lymphoma cases, from a cohort of BMDs. Mutations in PTPN11 were present in HS in 41/96 (43%) BMDs, and in 3/13 (23%) golden retrievers. PTPN11 mutations E76K and G503V did not coexist in the same neoplasm. The KRAS mutation was much less frequent, with a prevalence of 3.1% (3/96). We did not identify either PTPN11 nor KRAS mutations in any of the lymphoma samples. These results point out the potential relevance of PTPN11 and KRAS mutations as activators of the oncogenic MAPK pathway for canine HS, particularly in BMDs.

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.

Targeting MEK in a Translational Model of Histiocytic Sarcoma.

Histiocytic sarcoma in humans is an aggressive orphan disease with a poor prognosis as treatment options are limited. Dogs are the only species that spontaneously develops histiocytic sarcoma with an appreciable frequency, and may have value as a translational model system. In the current study, high-throughput drug screening utilizing histiocytic sarcoma cells isolated from canine neoplasms identified these cells as particularly sensitive to a MEK inhibitor, trametinib. One of the canine cell lines carries a mutation in PTPN11 (E76K), and another one in KRAS (Q61H), which are associated with the activation of oncogenic MAPK signaling. Both mutations were previously reported in human histiocytic sarcoma. Trametinib inhibited sensitive cell lines by promoting cell apoptosis, indicated by a significant increase in caspase 3/7. Furthermore, in vitro findings were successfully recapitulated in an intrasplenic orthotopic xenograft mouse model, which represents a disseminated aggressive form of histiocytic sarcoma. Mice with histiocytic sarcoma xenograft neoplasms that were treated with trametinib had significantly longer survival times. Target engagement was validated as activity of ERK, downstream of MEK, was significantly downregulated in neoplasms of treated mice. Additionally, trametinib was found in plasma and neoplastic tissues within projected therapeutic levels. These findings demonstrate that in dogs, histiocytic sarcoma may be associated with a dysfunctional MAPK pathway, at least in some cases, and may be effectively targeted through MEK inhibition. Clinical trials to test safety and efficacy of trametinib in dogs with histiocytic sarcoma are warranted, and may provide valuable translational information to similar diseases in humans. Mol Cancer Ther; 17(11); 2439-50. ©2018 AACR.

Transcriptome analysis of rhesus monkey failed-to-mature oocytes: deficiencies in transcriptional regulation and cytoplasmic maturation of the oocyte mRNA population.

STUDY QUESTION: Which different pathways and functions are altered in rhesus monkey oocytes that fail to mature after an ovulatory stimulus? SUMMARY ANSWER: Failed to mature (FTM) oocytes complete a large portion of the transition in transcriptome composition associated with normal maturation, but also manifest numerous differences that indicate incomplete transcriptional repression and cytoplasmic maturation affecting multiple processes. WHAT IS KNOWN ALREADY: Oocyte maturation defects contribute to unexplained female infertility. Failure of some oocytes to undergo germinal vesicle breakdown or progress to second meiotic metaphase in response to an ovulatory stimulus can limit the number of high quality oocytes available for ART. STUDY DESIGN, SIZE, DURATION: The transcriptome of rhesus monkey oocytes that failed to mature (FTM; n = 11, 5 donors) in response to an ovulatory stimulus in vivo was compared to those of normal germinal vesicle stage (GV, n = 7, 2 donors) and metaphase II stage (MII, n = 7, 5 donors) oocytes by RNA-sequencing (RNAseq). PARTICIPANTS/MATERIALS, SETTING, METHODS: Female rhesus monkeys of normal breeding age (6-12 years old) and with regular menstrual cycles were used. Animals underwent a controlled ovarian stimulation protocol for the collection of oocytes by ultrasound-guided needle aspiration of follicles. MAIN RESULTS AND THE ROLE OF CHANCE: We obtained a high quality RNAseq dataset consisting of n = 7, n = 7, and n = 11 libraries for normal GV, normal MII and FTM oocytes, respectively. Total reads acquired were an average of 34 million for each GV sample, 41 million for each FTM sample and 59 million for each MII oocyte sample. Approximately 44% of the total reads were exonic reads that successfully aligned to the rhesus monkey genome as unique non-rRNA gene transcript sequences, providing high depth of coverage. Approximately 44% of the mRNAs that undergo changes in abundance during normal maturation display partial modulations to intermediate abundances, and 9.2% fail to diverge significantly from GV stage oocytes. Additionally, a small group of mRNAs are grossly mis-regulated in the FTM oocyte. Differential expression was seen for mRNAs associated with mitochondrial functions, fatty acid beta oxidation, lipid accumulation, meiosis, zona pellucida formation, Hippo pathway signaling, and maternal mRNA regulation. A deficiency DNA methyltransferase one mRNA expression indicates a potential defect in transcriptional silencing. LARGE SCALE DATA: All RNAseq data are published in the Gene Expression Omnibus Database (GSE112536). LIMITATIONS, REASONS FOR CAUTION: These results do not establish cause of maturation failure but reveal novel correlates of incompetence to mature. Transcriptome studies likely do not capture all post-transcriptional or post-translational events that inhibit maturation, but do reveal mRNA expression changes that lie downstream of such events or that are related to effects on upstream regulators. The use of an animal model allows the study of oocyte maturation failure independent of covariates and confounders, such as pre-existing conditions of the female, which is a significant concern in human studies. Depending on the legislation, it may not be possible to collect and study oocytes from healthy women; and using surplus oocytes from patients undergoing ART may introduce confounders that vary from case to case. FTM oocytes were at various stages of meiotic progression, so correlates of specific times of arrest are not revealed. All the FTM oocytes failed to respond appropriately to an ovulatory stimulus in vivo. Therefore, this analysis informs us about common transcriptome features associated with meiotic incompetence. WIDER IMPLICATIONS OF THE FINDINGS: These results reveal that some diagnostic markers of oocyte quality may not reflect developmental competence because even meiotically incompetent oocytes display many normal gene expression features. The results also reveal potential mechanisms by which maternal and environmental factors may impact transcriptional repression and cytoplasmic maturation, and prevent oocyte maturation. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by grants from the National Institutes of Health Office of Research Infrastructure Programs Division of Comparative Medicine Grants R24 [OD012221 to K.E.L., OD011107/RR00169 (California National Primate Research Center), and OD010967/RR025880 to C.A.V.]; the Eunice Kennedy Shriver National Institute of Child Health and Human Development of the National Institutes of Health under the award number T32HD087166; MSU AgBioResearch, Michigan State University. Authors have nothing to disclose.