Decoding molecular features of bovine oocyte fate during antral follicle growth via single-cell multi-omics analysis†

Antral follicle size is a useful predictive marker of the competency of enclosed oocytes for yielding an embryo following in vitro maturation and fertilization. However, the molecular mechanisms underpinning oocyte developmental potential during bovine antral follicle growth are still unclear. Here, we used a modified single-cell multi-omics approach to analyze the transcriptome, DNA methylome, and chromatin accessibility in parallel for oocytes and cumulus cells collected from bovine antral follicles of different sizes. Transcriptome profiling identified three types of oocytes (small, medium, and large) that underwent different developmental trajectories, with large oocytes exhibiting the largest average follicle size and characteristics resembling metaphase-II oocytes. Differential expression analysis and real-time polymerase chain reaction assay showed that most replication-dependent histone genes were highly expressed in large oocytes. The joint analysis of multi-omics data revealed that the transcription of 20 differentially expressed genes in large oocytes was associated with both DNA methylation and chromatin accessibility. In addition, oocyte-cumulus interaction analysis showed that inflammation, DNA damage, and p53 signaling pathways were active in small oocytes, which had the smallest average follicle sizes. We further confirmed that p53 pathway inhibition in the in vitro maturation experiments using oocytes obtained from small antral follicles could improve the quality of oocytes and increased the blastocyte rate after in vitro fertilization and culture. Our work provides new insights into the intricate orchestration of bovine oocyte fate determination during antral folliculogenesis, which is instrumental for optimizing in vitro maturation techniques to optimize oocyte quality.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Influences of Supplementing Selective Members of the Interleukin-6 Cytokine Family on Bovine Oocyte Competency

This work explored whether supplementing selective members of the interleukin-6 (IL6) cytokine family during in vitro bovine oocyte maturation affects maturation success, cumulus-oocyte complex (COC) gene expression, fertilization success, and embryo development potential. Human recombinant proteins for IL6, IL11, and leukemia inhibitory factor (LIF) were supplemented to COCs during the maturation period, then fertilization and embryo culture commenced without further cytokine supplementation. The first study determined that none of these cytokines influenced the rate that oocytes achieved arrest at meiosis II. The second study identified that LIF and IL11 supplementation increases AREG transcript abundance. Supplementation with IL6 supplementation did not affect AREG abundance but reduced HAS2 transcript abundance. Several other transcriptional markers of oocyte competency were not affected by any of the cytokines. The third study determined that supplementing these cytokines during maturation did not influence fertilization success, but either LIF or IL11 supplementation increased blastocyst development. No effect of IL6 supplementation on subsequent blastocyst development was detected. The fourth experiment explored whether each cytokine treatment affects the post-thaw survivability of cryopreserved IVP blastocysts. None of the cytokines supplemented during oocyte maturation produced any positive effects on post-thaw blastocyst re-expansion and hatching. In conclusion, these outcomes implicate IL11 and LIF as potentially useful supplements for improving bovine oocyte competency.

A multi-omics analysis identifies molecular features associated with fertility in heifers (Bos taurus)

Infertility or subfertility is a critical barrier to sustainable cattle production, including in heifers. The development of heifers that do not produce a calf within an optimum window of time is a critical factor for the profitability and sustainability of the cattle industry. In parallel, heifers are an excellent biomedical model for understanding the underlying etiology of infertility because well-nourished heifers can still be infertile, mostly because of inherent physiological and genetic causes. Using a high-density single nucleotide polymorphism (SNP) chip, we collected genotypic data, which were analyzed using an association analysis in PLINK with Fisher's exact test. We also produced quantitative transcriptome data and proteome data. Transcriptome data were analyzed using the quasi-likelihood test followed by the Wald's test, and the likelihood test and proteome data were analyzed using a generalized mixed model and Student's t-test. We identified two SNPs significantly associated with heifer fertility (rs110918927, chr12: 85648422, P = 6.7 × 10-7; and rs109366560, chr11:37666527, P = 2.6 × 10-5). We identified two genes with differential transcript abundance (eFDR ≤ 0.002) between the two groups (Fertile and Sub-Fertile): Adipocyte Plasma Membrane Associated Protein (APMAP, 1.16 greater abundance in the Fertile group) and Dynein Axonemal Intermediate Chain 7 (DNAI7, 1.23 greater abundance in the Sub-Fertile group). Our analysis revealed that the protein Alpha-ketoglutarate-dependent dioxygenase FTO was more abundant in the plasma collected from Fertile heifers relative to their Sub-Fertile counterparts (FDR < 0.05). Lastly, an integrative analysis of the three datasets identified a series of molecular features (SNPs, gene transcripts, and proteins) that discriminated 21 out of 22 heifers correctly based on their fertility category. Our multi-omics analyses confirm the complex nature of female fertility. Very importantly, our results also highlight differences in the molecular profile of heifers associated with fertility that transcend the constraints of breed-specific genetic background.

In vivo and in vitro matured bovine oocytes present a distinct pattern of single-cell gene expression

Oocyte gene expression is a well controlled event that promotes gamete competence to undergo maturation, fertilization, and to support early embryo development, directly affecting reproductive outcomes. Considering that in vivo controlled ovarian stimulation or in vitro maturation (IVM) for the acquisition of mature oocytes has distinct implications for gene expression, we sought to evaluate the effects of these procedures on the expression of competence-related genes in single-cell oocytes. Healthy Nelore cows of reproductive age were synchronized to harvest in vivo matured oocytes; ovaries from slaughtered animals were used to obtain cumulus-oocyte complexes that were in vitro matured. Single-cell gene expression was performed using TaqMan Low-Density Arrays and 42 genes were evaluated. In silico analysis of protein interactions and Gene Ontology (GO) analysis was performed. Reduced gene expression was observed for 24 targets in IVM oocytes when compared with those of in vivo matured oocytes (P < 0.05). Differences ranged from 1.5-fold to 4.8-fold higher in in vivo oocytes and the BMP15 (5.28), GDF9 (6.23), NOBOX (7.25), HSPA8 (7.85) and MSX1 (11.00) showed the greatest fold increases. The strongest score of functional interactions was observed between the CDC20 and CKS2, with the differentially expressed gene CDC20 being the main marker behind GO enrichment. IVM negatively affected the expression of important genes related to oocyte competency, and showed higher expression levels in in vivo matured oocytes. In vivo controlled ovarian stimulation may be a better strategy to achieve proper oocyte competence and increase the success of assisted reproductive technologies.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Cleavage kinetics is a better indicator of embryonic developmental competency than brilliant cresyl blue staining of oocytes

In vitro production of embryos (IVP) is a valuable technology to produce embryos of high genetic value. Despite advances in IVP, the efficiency of culture systems remains low. One method to increase IVP success is the early selection of oocytes or embryos that may have greater developmental potential. Here, we investigated two methods of selection, namely BCB staining and cleavage kinetics, both individually and in conjunction, for improved developmental outcomes in vitro. We hypothesized that a synergistic use of both BCB staining and cleavage kinetics would result in identification of embryos of greater developmental potential. The selection of oocytes by BCB staining does select for those oocytes with higher developmental potential, as noted by a greater blastocyst development between BCB positive (32.6%) and BCB negative (22.0%) on day 8 post-fertilization. However, the utilization of BCB staining and cleavage kinetics in tandem resulted in a complete masking of the effect observed when using BCB alone. We obtained the highest proportion of blastocyst development per selection group using cleavage kinetics alone, in which 53.1% of embryos grouped as Fast produced a blastocyst, which was significantly different from the three other groups (Fast+, Slow, not cleaved). We observed, however, that the separation of embryos by cleavage kinetics did not predict their survival to cryopreservation. In conclusion, in standard culture systems, cleavage kinetics is an effective method for the selection of embryos with increased developmental potential to develop blastocysts, however, it may not be effective to select healthy embryos for transfer following cryopreservation.

Tight gene co-expression in BCB positive cattle oocytes and their surrounding cumulus cells

BACKGROUND

Cytoplasmic and nuclear maturation of oocytes, as well as interaction with the surrounding cumulus cells, are important features relevant to the acquisition of developmental competence.

METHODS

Here, we utilized Brilliant cresyl blue (BCB) to distinguish cattle oocytes with low activity of the enzyme Glucose-6-Phosphate Dehydrogenase, and thus separated fully grown (BCB positive) oocytes from those in the growing phase (BCB negative). We then analyzed the developmental potential of these oocytes, mitochondrial DNA (mtDNA) copy number in single oocytes, and investigated the transcriptome of single oocytes and their surrounding cumulus cells of BCB positive versus BCB negative oocytes.

RESULTS

The BCB positive oocytes were twice as likely to produce a blastocyst in vitro compared to BCB- oocytes (P < 0.01). We determined that BCB negative oocytes have 1.3-fold more mtDNA copies than BCB positive oocytes (P = 0.004). There was no differential transcript abundance of genes expressed in oocytes, however, 172 genes were identified in cumulus cells with differential transcript abundance (FDR < 0.05) based on the BCB staining of their oocyte. Co-expression analysis between oocytes and their surrounding cumulus cells revealed a subset of genes whose co-expression in BCB positive oocytes (n = 75) and their surrounding cumulus cells (n = 108) compose a unique profile of the cumulus-oocyte complex.

CONCLUSIONS

If oocytes transition from BCB negative to BCB positive, there is a greater likelihood of producing a blastocyst, and a reduction of mtDNA copies, but there is no systematic variation of transcript abundance. Cumulus cells present changes in transcript abundance, which reflects in a dynamic co-expression between the oocyte and cumulus cells.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Oocyte quality following in vitro follicle development†

In vitro follicle development (IVFD) is an adequate model to obtain basic knowledge of folliculogenesis and provides a tool for ovarian toxicity screening. IVFD yielding competent oocytes may also offer an option for fertility and species preservation. To promote follicle growth and oocyte maturation in vitro, various culture systems are utilized for IVFD in rodents, domestic animals, wild animals, nonhuman primates, and humans. Follicle culture conditions have been improved by optimizing gonadotropin levels, regulatory factors, nutrient supplements, oxygen concentration, and culture matrices. This review summarizes quality assessment of oocytes generated from in vitro-developed antral follicles from the preantral stage, including oocyte epigenetic and genetic profile, cytoplasmic and nuclear maturation, preimplantation embryonic development following in vitro fertilization, as well as pregnancy and live offspring after embryo transfer. The limitations of oocyte quality evaluation following IVFD and the gaps in our knowledge of IVFD to support proper oocyte development are also discussed. The information may advance our understanding of the requirements for IVFD, with a goal of producing competent oocytes with genetic integrity to sustain embryonic development resulting in healthy offspring.

Isolation of high-quality total RNA and RNA sequencing of single bovine oocytes

Studying individual mammalian oocytes has been extremely valuable for the understanding of the molecular composition of oocytes including RNA storage. Here, a detailed protocol for isolation of oocytes, extraction of total RNA from single oocytes followed by full-length cDNA amplification, and library preparation is presented. The procedure permits the production of cost-effective and high-quality sequencing libraries. This protocol can be adapted for transcriptome analysis of oocytes from other species and be used to generate high-quality data from single embryos.

For complete details on the use and execution of this protocol, please refer to Biase and Kimble (2018).

•

Isolation of high-quality total RNA from single bovine oocytes

•

Detailed procedures for amplification of complementary DNA for library preparation

•

A bioinformatic pipeline for the quantitation of transcript abundance

•

The protocol also enables high-quality data production from single embryos

Folic acid supplementation during oocytes maturation influences in vitro production and gene expression of bovine embryos

Embryos that are produced in vitro frequently present epigenetic modifications. However, maternal supplementation with folic acid (FA) may improve oocyte maturation and embryo development, preventing epigenetic errors in the offspring. We sought to evaluate the influence of FA supplementation during in vitro maturation of grade I (GI) and grade III (GIII) bovine oocytes on embryo production rate and the expression of IGF2 and KCNQ1OT1 genes. The oocytes were matured in vitro with different concentrations of FA (0, 10, 30 and 100 μM), followed by in vitro fertilization and embryo culture. On the seventh day (D7) of culture, embryo production was evaluated and gene expression was measured using real-time qPCR. Supplementation with 10 μM of FA did not affect embryo production for GI and GIII oocytes. Moderate supplementation (30 μM) seemed to be a positive influence, increasing embryo production for GIII (P = 0.012), while the highest dose (100 μM) reduced embryo production (P = 0.010) for GI, and IGF2 expression was not detected. In GIII, only embryos whose oocyte maturation was not supplemented with FA demonstrated detected IGF2 expression. The lowest concentration of FA (10 μM) reduced KCNQ1OT1 expression (P = 0.05) on embryos from GIII oocytes. Different FA concentrations induced different effects on bovine embryo production and gene expression that was related to oocyte quality. Despite the epigenetic effects of FA, supplementation seems to be a promising factor to improve bovine embryo production if used carefully, as concentration is an important factor, especially in oocytes with impaired quality.

Meiotic arrest with roscovitine and sexual maturity improve competence of mouse oocytes by regulating expression of competence-related genes

We have studied the mechanisms by which meiotic arrest maintenance (MAM) with roscovitine, female sexual maturity, and the surrounded nucleoli (SN) chromatin configuration improve the competence of mouse oocytes by observing the expression of oocyte competence-related genes in non-surrounded nucleoli (NSN) and SN oocytes from prepubertal and adult mice following maturation with or without MAM. The results demonstrated that MAM with roscovitine significantly improved the developmental potential of adult SN and prepubertal NSN oocytes, but had no effect on that of prepubertal SN oocytes. Without MAM, while 40% of the 2-cell embryos derived from prepubertal SN oocytes developed into 4-cell embryos, none of the 2-cell embryos derived from prepubertal NSN oocytes did, and while 42% of the 4-cell embryos derived from adult SN oocytes developed into blastocysts, only 1% of the 4-cell embryos derived from prepubertal SN oocytes developed into blastocysts. Furthermore, MAM with roscovitine, SN configuration, and female sexual maturity significantly increased the mRNA levels of competence-beneficial genes and decreased those of competence-detrimental genes. In conclusion, our results suggest that MAM with roscovitine, SN chromatin configuration, and female sexual maturity improve oocyte competence by regulating the expression of competence-related genes, suggesting that Oct4, Stella, Mater, Zar1, Mapk8, and Bcl2 are oocyte competence-beneficial genes, whereas Foxj2, Ship1, and Bax are competence-detrimental genes.