Single-cell RNA-seq reveals mRNAs and lncRNAs important for oocytes in vitro matured in pigs

Meiotic and developmental competence of growing pig oocytes derived from small antral follicles is enhanced in culture medium containing FGF2, LIF, and IGF1 (FLI medium)

BACKGROUND

Oocytes of large animal species isolated from small ovarian follicles (< 2 mm) are less competent to support early embryonic development after in vitro maturation and fertilization than their counterparts isolated from medium-sized and preovulatory follicles. This study aimed to assess the effect of a new maturation medium containing FGF2, LIF, and IGF1 (FLI medium) on the meiotic and developmental competence of pig cumulus-oocytes complexes (COCs) derived from the small and medium-sized follicles.

METHODS

The growing oocytes were isolated from 1 to 2 (small follicle; SF) and the fully-grown ones from 3 to 6 (large follicle; LF) mm follicles and matured in a control M199 medium with gonadotropins and EGF and the FLI medium enriched by the triplet of growth factors. The matured oocytes were parthenogenetically activated and cultured to the blastocyst stage. Chromatin configuration before and during the culture and MAP kinase activity were assessed in the oocytes. Finally, the expression of cumulus cell genes previously identified as markers of oocyte quality was assessed.

RESULTS

The maturation and blastocyst rates of oocytes gained from LF were significantly higher than that from SF in the control medium. In contrast, similar proportions of oocytes from LF and SF completed meiosis and developed to blastocysts when cultured in FLI. Most of the oocytes freshly isolated from SF possessed germinal vesicles with fine filaments of chromatin (GV0) or chromatin surrounding the nucleolus (GVI; 30%); the oocytes from LF were mainly in GVI (or GVII) exhibiting a few small lumps of chromatin beneath the nuclear membrane. When cultured in the FLI medium for 16 h, an acceleration of the course of maturation in oocytes both from SF and LF compared to the control medium was observed and a remarkable synchrony in the course of chromatin remodeling was noticed in oocytes from SF and LF.

CONCLUSIONS

This work demonstrates that the enrichment of culture medium by FGF2, LIF, and IGF1 can enhance the meiotic and developmental competence of not only fully-grown, but also growing pig oocytes and significantly thus expanding the number of oocytes available for various assisted reproductive technology applications.

Transcriptome analyses reveal transcriptional profiles of horse oocytes before and after in vitro maturation

Oocyte in vitro maturation is necessary for the study and application of animal-assisted reproduction technology in animal reproduction and breeding. The comprehensive transcriptional profile of equine oocyte maturated in vitro has not been fully mined yet, which makes many key transcriptional events still unidentified. Here, Smart-seq2 was performed to analyse the gene expression pattern and the underlying regulatory mechanism of horse germinal vesicle (GV) and in vitro metaphase II (MII) oocytes. The results showed that 6402 genes (2640 up-regulated and 3762 down-regulated in MII samples compared to GV) and 4021 lncRNA transcripts (1210 up-regulated and 2811 down-regulated in MII samples compared to GV) were differentially expressed in GV and MII oocytes. Further, GO and KEGG analysis found that differentially expressed mRNAs and lncRNAs were mainly enriched in the pathways related to energy and lipid metabolism. In addition, LGALS3 was found a key gene in mediating the regulation of oocyte meiosis recovery and fertilization ability. This study provides novel knowledge about gene expression and energy metabolism during equine oocyte maturation and a reference for the further study and application of assisted reproductive technology in horse reproduction and breeding.

Porcine oocytes matured in a chemically defined medium are transcriptionally active

The statement that fully-grown porcine oocytes (oocytes from follicles with diameter from 3 to 6 mm) are transcriptionally quiescent is not as strongly supported as it was before. Currently, we know that there is a difference between the transcription profile of germinal vesicle (GV) and metaphase II (MII) oocytes. The goal of our study was to compare the transcription profile of GV, germinal vesicle breakdown (GVBD), metaphase I (MI), and MII oocytes matured in the chemically defined medium FLI. Oocytes were sequenced, and the results were subsequently validated using quantitative reverse transcription polymerase chain reaction (RT-qPCR). We detected multiple differentially transcribed mRNAs, of which many were upregulated. Among them we found mRNAs necessary for protein production, mitochondrial functions and cytoplasmic maturation. Collectively, these data support the hypothesis that transcription activity in fully-grown porcine oocytes is necessary for key processes during their successful maturation in vitro in a chemically defined maturation medium.

Advances in Oocyte Maturation In Vivo and In Vitro in Mammals

The quality and maturation of an oocyte not only play decisive roles in fertilization and embryo success, but also have long-term impacts on the later growth and development of the fetus. Female fertility declines with age, reflecting a decline in oocyte quantity. However, the meiosis of oocytes involves a complex and orderly regulatory process whose mechanisms have not yet been fully elucidated. This review therefore mainly focuses on the regulation mechanism of oocyte maturation, including folliculogenesis, oogenesis, and the interactions between granulosa cells and oocytes, plus in vitro technology and nuclear/cytoplasm maturation in oocytes. Additionally, we have reviewed advances made in the single-cell mRNA sequencing technology related to oocyte maturation in order to improve our understanding of the mechanism of oocyte maturation and to provide a theoretical basis for subsequent research into oocyte maturation.

Transcriptome analysis reveals mRNAs and long non-coding RNAs associated with fecundity in the hypothalamus of high-and low-fecundity goat

As an important organ that coordinates the neuroendocrine system, the hypothalamus synthesizes and secretes reproductive hormones that act on the goat organism, thereby precisely regulating follicular development and reproductive processes in goats. However, it is still elusive to explore the mechanism of hypothalamic effects on goat fertility alone. Therefore, RNA-seq was used to analyze the gene expression in hypothalamic tissues of goats in high fertility group (HFG: litter size per litter ≥2) and low fertility group (LFG: litter size per litter = 1), and identified the differential lncRNAs and mRNAs and their associated pathways related to their fertility. The results showed that a total of 23 lncRNAs and 57 mRNAs were differentially expressed in the hypothalamic tissue of high and low fertility goats. GO terms and KEGG functional annotation suggest that DE lncRNAs and DE mRNAs were significantly enriched in hormone-related pathways regulating ovarian development, hormone synthesis and secretion, regulation of reproductive processes, Estrogen signaling pathway, Oxytocin signaling pathway and GnRH signaling pathway. And we constructed a co-expression network of lncRNAs and target genes, and identified reproduction-related genes such as NMUR2, FEZF1, and WT1. The sequencing results of the hypothalamic transcriptome have broadened our understanding of lncRNA and mRNA in goat hypothalamic tissue and provided some new insights into the molecular mechanisms of follicle development and regulation of its fertility in goats.

Supplementation of SDF1 during Pig Oocyte In Vitro Maturation Improves Subsequent Embryo Development

The quality of in vitro matured oocytes is inferior to that of in vivo matured oocytes, which translates to low developmental capacity of embryos derived from in vitro matured oocytes. The developmental potential of in vitro matured oocytes is usually impaired due to oxidative stress. Stromal cell-derived factor-l (SDF1) can reduce oxidative stress and inhibit apoptosis. The aim of this study was to investigate the effects of SDF1 supplementation during pig oocyte in vitro maturation (IVM) on subsequent embryo development, and to explore the acting mechanisms of SDF1 in pig oocytes. We found that the IVM medium containing 20 ng/mL SDF1 improved the maturation rate of pig oocytes, as well as the cleavage rate and blastocyst rate of embryos generated by somatic cell nuclear transfer, in vitro fertilization, and parthenogenesis. Supplementation of 20 ng/mL SDF1 during IVM decreased the ROS level, increased the mitochondrial membrane potential, and altered the expression of apoptosis-related genes in the pig oocytes. The porcine oocyte transcriptomic data showed that SDF1 addition during IVM altered the expression of genes enriched in the purine metabolism and TNF signaling pathways. SDF1 supplementation during pig oocyte IVM also upregulated the mRNA and protein levels of YY1 and TET1, two critical factors for oocyte development. In conclusion, supplementation of SDF1 during pig oocyte IVM reduces oxidative stress, changes expression of genes involved in regulating apoptosis and oocyte growth, and enhances the ability of in vitro matured pig oocytes to support subsequent embryo development. Our findings provide a theoretical basis and a new method for improving the developmental potential of pig in vitro matured oocytes.

Network analysis reveals different hub genes and molecular pathways for pig in vitro fertilized early embryos and parthenogenotes

Maternal-to-zygotic transition (MZT) occurs when maternal transcripts decay and zygotic genome is activated gradually at early stage of embryo development. Previously, single cell RNA-seq (scRNA-seq) has helped us to uncover the MZT-associated mRNA dynamics of in vitro produced pig early embryos. Here, to further investigate functional modules and hub genes associated with MZT process, the weighted gene-coexpression network analysis (WGCNA) was performed on our previously generated 45 scRNA-seq datasets. For the in vitro fertilized embryo (IVF) group, 5 significant modules were identified (midnightblue/black/red and blue/brown modules, positively correlated with 1-cell (IVF1) and 8-cell (IVF8), respectively), containing genes mainly enriched in signaling pathways such as Wnt, regulation of RNA transcription, fatty acid metabolic process, poly(A) RNA binding and lysosome. For the parthenogenetically activated embryo (PA) group, 9 significant modules were identified (black/purple/red, brown/turquoise/yellow, and magenta/blue/green modules, positively correlated with MII oocytes, 1-cell (PA1), and 8-cell (PA8), respectively), mainly enriched in extracellular exosome, poly(A) RNA binding, mitochondrion, transcription factor activity. Moreover, some of identified hub genes within 3 IVF and 9 PA significant modules, including ADCY2, DHX34, KDM4A, GDF10, ABCC10, PAFAH2, HEXIM2, COQ9, DCAF11, SGK1, ESRRB etc., have been reported to play vital roles in different biological processes. Our findings provide information and resources for subsequent in-depth study on the regulation and function of MZT in pig embryos.

Transcriptome Analysis of mRNAs and Long Non-Coding RNAs During Subsequent Embryo Development of Porcine Cloned Zygotes After Vitrification

Cryopreservation of porcine cloned zygotes has important implications for biotechnology and biomedicine research; however, lower embryo developmental potential remains an urgent problem to be resolved. For exploring the sublethal cryodamages during embryo development, this study was designed to acquire the mRNA and long non-coding RNA (lncRNA) profiles of 2-cells, 4-cells and blastocysts derived from vitrified porcine cloned zygotes using transcriptome sequencing. We identified 167 differentially expressed (DE) mRNAs and 516 DE lncRNAs in 2-cell stage, 469 DE mRNAs and 565 lncRNAs in 4-cell stage, and 389 DE mRNAs and 816 DE lncRNAs in blastocyst stage. Functional enrichment analysis revealed that the DE mRNAs during embryo development were involved in many regulatory mechanisms related to cell cycle, cell proliferation, apoptosis, metabolism and others. Moreover, the target genes of DE lncRNAs in the three embryonic stages were also enriched in many key GO terms or pathways such as “defense response”, “linoleic acid metabolic process”, “embryonic axis specification”, “negative regulation of protein neddylation”, etc., In conclusion, the present study provided comprehensive transcriptomic data about mRNAs and lncRNAs for the vitrified porcine cloned zygotes during different developmental stages, which contributed to further understand the potential cryodamage mechanisms responsible for impaired embryo development.

Identification of lncRNAs involved in maternal-to-zygotic transition of in vitro-produced porcine embryos by single-cell RNA-seq

Long non-coding RNAs (lncRNAs) function through multiple tiers of molecular circuits and are vital to gamete maturation and early embryo development. However, in pig early embryos, identification and expression dynamics of lncRNAs remain less studied. Here, we systematically analysed the expression dynamics of lncRNAs based on our previously published single-cell RNA-seq data from pig mature oocytes (GSE160334), and single blastomeres biopsied from pig in vitro fertilized (IVF) and early parthenogenetically activated (PA) embryos (1- to 8-cell stages; GSE164812). With the progression of embryo development, the total number of expressed lncRNAs gradually decreased and showed great variation at each developmental stage for both IVF and PA groups. Consecutive stage pairwise comparison of MII oocytes, 1-cell zygotes, 2-cell, 4-cell and 8-cell IVF embryos identified 151, 245, 1119 and 188 differentially expressed (DE) lncRNAs, including 119, 80, 867, 77 up-regulated and 32, 165, 252, 111 down-regulated, while 289, 437, 895 and 495 DE lncRNAs (141, 89, 768, 97 up-regulated and 148, 348, 127, 398 down-regulated) were identified in PA embryos at the same stages. The DE lncRNAs identified within IVF embryos were much different from that identified within PA embryos, showing embryo type-specific manner. Further cross-comparison between PA and IVF embryos identified 184, 656, 2502 and 266 DE lncRNAs for the 1- to 8-cell embryo stages, respectively. Further GO and KEGG enrichment analysis of DE mRNAs targeted by DELs indicated that different signalling pathways were involved in maternal-only and bi-parental embryo development. Collectively, comparative profiling of lncRNA expression dynamics between pig IVF and PA embryos provides a valuable resource, to investigate further regulatory mechanisms of lncRNAs associated with ZGA and maternal RNA decay during early embryo development.

Global 3'-untranslated region landscape mediated by alternative polyadenylation during meiotic maturation of pig oocytes

Alternative polyadenylation affects the length and composition of 3'-untranslated region (3'-UTR) and regulates mRNA stability or translational activity to affect important biological processes. However, global 3'-UTR landscape and its relationship with gamete maturation remain less studied. Here, we analysed our previously reported single-cell RNA-seq data of germinal vesicle and metaphase II stage oocytes in pigs to systematically catalogue the 3'-UTR dynamics during oocyte maturation. Two softwares (DaPars and APAtrap) were employed and identified 110 and 228 mRNAs with significantly different 3'-UTRs (adjusted p ≤ .05), respectively. Gene enrichment analyses found signalling pathways related with biological processes of female gametophyte production, methyltransferase activity and mRNA surveillance (DaPars) and cell cycle process, regulation of ERK1 and ERK2 cascade, regulation of translation, spindle organization, kinetochore, condensed chromosome and progesterone-mediated oocyte maturation (APAtrap), respectively. Moreover, 18 of 110 mRNAs (|△PDUI| ≥ 0.25 and |log₂ PDUI ratio| ≥ 0.59) and 15 of 228 mRNAs (Perc. diff. ≥ 0.5) were with greater difference of 3'-UTR length or abundance, and integrative genomics viewer analysis further identified 4 (Alg10, Hadhb, Hsd17b4 and Sbds) of 18 mRNAs to be with 3'-UTR length differed ≥150 bp and 6 (Gcc1, Hnrnpa2b1, Lsm6, Prpf18, Sfr1 and Ust) of 15 mRNAs to be with 3'-UTR abundance extremely differed. Furthermore, the location, sequences and number of cis-elements were predicted, which were shown to derange cytoplasmic polyadenylation element, poly(A) site and microRNA binding sites within 3'-UTRs of Alg10, Hadhb, Hsd17b4 and Sbds mRNAs. Taken together, global 3'-UTR landscape changes dynamically with oocyte meiotic maturation, potentially involved in regulating oocyte meiotic process in pigs.

The Application of Single-Cell RNA Sequencing in Mammalian Meiosis Studies

Meiosis is a cellular division process that produces gametes for sexual reproduction. Disruption of complex events throughout meiosis, such as synapsis and homologous recombination, can lead to infertility and aneuploidy. To reveal the molecular mechanisms of these events, transcriptome studies of specific substages must be conducted. However, conventional methods, such as bulk RNA-seq and RT-qPCR, are not able to detect the transcriptional variations effectively and precisely, especially for identifying cell types and stages with subtle differences. In recent years, mammalian meiotic transcriptomes have been intensively studied at the single-cell level by using single-cell RNA-seq (scRNA-seq) approaches, especially through two widely used platforms, Smart-seq2 and Drop-seq. The scRNA-seq protocols along with their downstream analysis enable researchers to accurately identify cell heterogeneities and investigate meiotic transcriptomes at a higher resolution. In this review, we compared bulk RNA-seq and scRNA-seq to show the advantages of the scRNA-seq in meiosis studies; meanwhile, we also pointed out the challenges and limitations of the scRNA-seq. We listed recent findings from mammalian meiosis (male and female) studies where scRNA-seq applied. Next, we summarized the scRNA-seq analysis methods and the meiotic marker genes from spermatocytes and oocytes. Specifically, we emphasized the different features of the two scRNA-seq protocols (Smart-seq2 and Drop-seq) in the context of meiosis studies and discussed their strengths and weaknesses in terms of different research purposes. Finally, we discussed the future applications of scRNA-seq in the meiosis field.

Ovarian Aging: Role of Pituitary-Ovarian Axis Hormones and ncRNAs in Regulating Ovarian Mitochondrial Activity

The number of mitochondria in the oocyte along with their functions (e.g., energy production, scavenger activity) decline with age progression. Such multifaceted functions support several processes during oocyte maturation, ranging from energy supply to synthesis of the steroid hormones. Hence, it is hardly surprising that their impairment has been reported in both physiological and premature ovarian aging, wherein they are crucial players in the apoptotic processes that arise in aged ovaries. In any form, ovarian aging implies the progressive damage of the mitochondrial structure and activities as regards to ovarian germ and somatic cells. The imbalance in the circulating hormones and peptides (e.g., gonadotropins, estrogens, AMH, activins, and inhibins), active along the pituitary-ovarian axis, represents the biochemical sign of ovarian aging. Despite the progress accomplished in determining the key role of the mitochondria in preserving ovarian follicular number and health, their modulation by the hormonal signalling pathways involved in ovarian aging has been poorly and randomly explored. Yet characterizing this mechanism is pivotal to molecularly define the implication of mitochondrial dysfunction in physiological and premature ovarian aging, respectively. However, it is fairly difficult considering that the pathways associated with ovarian aging might affect mitochondria directly or by altering the activity, stability and localization of proteins controlling mitochondrial dynamics and functions, either unbalancing other cellular mediators, released by the mitochondria, such as non-coding RNAs (ncRNAs). We will focus on the mitochondrial ncRNAs (i.e., mitomiRs and mtlncRNAs), that retranslocate from the mitochondria to the nucleus, as active players in aging and describe their role in the nuclear-mitochondrial crosstalk and its modulation by the pituitary-ovarian hormone dependent pathways. In this review, we will illustrate mitochondria as targets of the signaling pathways dependent on hormones and peptides active along the pituitary/ovarian axis and as transducers, with a particular focus on the molecules retrieved in the mitochondria, mainly ncRNAs. Given their regulatory function in cellular activities we propose them as potential diagnostic markers and/or therapeutic targets.