Dynamic changes in gene expression of growing nonhuman primate antral follicles.

The growth of the ovarian antral follicle is a complex process that is difficult to study, especially in human and nonhuman primates. Understanding the antral stage of development is key to new approaches to regulating reproduction. This study analyzed cohorts of three sizes of developing antral follicles obtained from adult rhesus macaque females using RNA sequencing of oocytes and cumulus and granulosa cells. The overall objective of this study was to identify key developmental changes in gene expression in oocytes, granulosa, and cumulus cells, as nonhuman primate antral stage follicles transition through progressively larger sizes in the absence of exogenous hormonal stimulation. Only a relatively small number of genes displayed altered mRNA expression levels in any of the three cell types during this period. Most of the identified differentially expressed genes (DEGs) decreased in the granulosa cells or increased in the cumulus cells. Although the number of DEGs observed was small, these DEGs indicate predicted effects on distinct upstream regulators in the cumulus and granulosa cells. This study is particularly important because it shows for the first time the gene expression changes during antral follicle growth in a medically relevant model.

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.

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-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.

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.