Transcriptomic profiles reveal the characteristics of oocytes and cumulus cells at GV, MI, and MII in follicles before ovulation.

BACKGROUND: The oocyte and its surrounding cumulus cells (CCs) exist as an inseparable entity. The maturation of the oocyte relies on communication between the oocyte and the surrounding CCs. However, oocyte evaluation is primarily based on morphological parameters currently, which offer limited insight into the quality and competence of the oocyte. Here, we conducted transcriptomic profiling of oocytes and their CCs from 47 patients undergoing preimplantation genetic testing for aneuploidy (PGT-A). We aimed to investigate the molecular events occurring between oocytes and CCs at different stages of oocyte maturation (germinal vesicle [GV], metaphase I [MI], and metaphase II [MII]). Our goal is to provide new insights into in vitro oocyte maturation (IVM). RESULTS: Our findings indicate that oocyte maturation is a complex and dynamic process and that MI oocytes can be further classified into two distinct subtypes: GV-like-MI oocytes and MII-like-MI oocytes. Human oocytes and cumulus cells at three different stages of maturation were analyzed using RNA-seq, which revealed unique transcriptional machinery, stage-specific genes and pathways, and transcription factor networks that displayed developmental stage-specific expression patterns. We have also identified that both lipid and cholesterol metabolism in cumulus cells is active during the late stage of oocyte maturation. Lipids may serve as a more efficient energy source for oocytes and even embryogenesis. CONCLUSIONS: Overall, our study provides a relatively comprehensive overview of the transcriptional characteristics and potential interactions between human oocytes and cumulus cells at various stages of maturation before ovulation. This study may offer novel perspectives on IVM and provide a reliable reference data set for understanding the transcriptional regulation of follicular maturation.

Telomeres cooperate in zygotic genome activation by affecting DUX4/Dux transcription.

Zygotic genome activation (ZGA) is initiated once the genome chromatin state is organized in the newly formed zygote. Telomeres are specialized chromatin structures at the ends of chromosomes and are reset during early embryogenesis, while the details and significance of telomere changes in preimplantation embryos remain unclear. We demonstrated that the telomere length was shortened in the minor ZGA stage and significantly elongated in the major ZGA stage of human and mouse embryos. Expression of the ZGA pioneer factor DUX4/Dux was negatively correlated with the telomere length. ATAC sequencing data revealed that the chromatin accessibility peaks on the DUX4 promoter region (i.e., the subtelomere of chromosome 4q) were transiently augmented in human minor ZGA. Reduction of telomeric heterochromatin H3K9me3 in the telomeric region also synergistically activated DUX4 expression with p53 in human embryonic stem cells. We propose herein that telomeres regulate the expression of DUX4/Dux through chromatin remodeling and are thereby involved in ZGA.

Translation regulatory factor BZW1 regulates preimplantation embryo development and compaction by restricting global non-AUG Initiation.

Protein synthesis is an essential step in gene expression during the development of mammalian preimplantation embryos. This is a complex and highly regulated process. The accuracy of the translation initiation codon is important in various gene expression programs. However, the mechanisms that regulate AUG and non-AUG codon initiation in early embryos remain poorly understood. BZW1 is a key factor in determining the mRNA translation start codon. Here, we show that BZW1 is essential for early embryonic development in mice. Bzw1-knockdown embryos fail to undergo compaction, and show decreased blastocyst formation rates. We also observe defects in the differentiation capacity and implantation potential after Bzw1 interference. Further investigation revealed that Bzw1 knockdown causes the levels of translation initiation with CUG as the start codon to increase. The decline in BZW1 levels result in a decrease in protein synthesis in preimplantation embryos, whereas the total mRNA levels are not altered. Therefore, we concluded that BZW1 contributes to protein synthesis during early embryonic development by restricting non-AUG translational initiation.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development.

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

The Dynamic Changes of Transcription Factors During the Development Processes of Human Biparental and Uniparental Embryos.

Previous studies have revealed that transcription factors (TFs) play important roles in biparental (BI) early human embryogenesis. However, the contribution of TFs during early uniparental embryo development is still largely unknown. Here we systematically studied the expression profiles of transcription factors in early embryonic development and revealed the dynamic changes of TFs in human biparental and uniparental embryogenesis by single-cell RNA sequencing (scRNA-seq). In general, the TF expression model of uniparental embryos showed a high degree of conformity with biparental embryos. The detailed network analysis of three different types of embryos identified that 10 out of 17 hub TFs were shared or specifically owned, such as ZNF480, ZNF581, PHB, and POU5F1, were four shared TFs, ZFN534, GTF3A, ZNF771, TEAD4, and LIN28A, were androgenic (AG) specific TFs, and ZFP42 was the only one parthenogenetic (PG) specific TF. All the four shared TFs were validated using human embryonic stem cell (hESC) differentiation experiments; most of their target genes are responsible for stem cell maintenance and differentiation. We also found that Zf-C2H2, HMG, and MYB were three dominant transcription factor families that appeared in early embryogenesis. Altogether, our work provides a comprehensive regulatory framework and better understanding of TF function in human biparental and uniparental embryogenesis.

ß-estradiol promotes the growth of primary human fetal spermatogonial stem cells via the induction of stem cell factor in Sertoli cells.

BACKGROUND: Mammalian spermatogenesis is responsible for male fertility and is supported by the self-renewal and differentiation of spermatogonial stem cells (SSCs). Sertoli cells provide a supportive microenvironment for SSCs, in part by the production of stem cell factor (SCF), which is a potent regulator of spermatogonia proliferation and survival. METHODS: We investigated the novel role of ß-estradiol in modulating the proliferation and apoptosis of fetal SSCs via the regulation of SCF secretion in Sertoli cells isolated from human fetal testes. The proliferation of SSCs in the co-culture system was determined by colony formation and BrdU incorporation assays. TUNEL assay was used to measure SSC apoptosis in co-culture in response to treatment with control, ß-estradiol, or the combination of ß-estradiol and the estrogen receptor inhibitor ICI 182780. RESULTS: In the system with purified human fetal Sertoli cells (MIS+/c-Kit-/AP-), ß-estradiol upregulated the production of SCF in a dose- and time-dependent manner. In the co-culture system of primary human fetal SSCs (c-Kit+/SSEA-4+/Oct-4+/AP+) and Sertoli cells (MIS+), ß-estradiol markedly increased the proliferation of SSCs. Moreover, SSC apoptosis was significantly inhibited by ß-estradiol and was completely reversed by the combination of ß-estradiol and ICI 182780. CONCLUSION: Here we report, for the first time, that ß-estradiol can induce the increase of SCF expression in human fetal Sertoli cells and regulates the growth and survival of human fetal SSCs. These novel findings provide new perspectives on the current understanding of the role of estrogen in human spermatogenesis.

CDCA8 regulates meiotic spindle assembly and chromosome segregation during human oocyte meiosis.

As a member of the chromosomal passenger complex, CDCA8 (cell division cycle associated 8) plays an important role in human mitosis, but its roles in human meiosis are unknown. Here, we show that CDCA8 expression is increased and its encoded protein has dynamic localization in human oocytes from germinal vesicle breakdown (GVBD) to metaphase Ⅱ (MⅡ), and that there are multipolar spindles, disordered chromosomes, and that microtubule assembly is affected after CDCA8 RNA interference (RNAi) in GV-stage oocytes. The GVBD and polar body extrusion (PBE) rates were not affected following CDCA8 depletion, but the PBE time was extended. There was no statistical difference between CDCA8 expression of oocytes from older and younger women, but the first polar body from older women was prone to chromosome abnormalities, and oocytes with such abnormalities had lower CDCA8 expression than oocytes with normal polar bodies. These results indicate that CDCA8 is associated with bipolar spindle formation, chromosome segregation, PBE during human oocyte meiosis, and that it may affect the incidence of aneuploidy embryos in older women.

Single-Cell Transcriptome Analysis of Uniparental Embryos Reveals Parent-of-Origin Effects on Human Preimplantation Development.

To investigate the contribution of parental genomes to early embryogenesis, we profiled the single-cell transcriptomes of human biparental and uniparental embryos systematically from the 1-cell to the morula stage. We observed that uniparental embryos exhibited variable and decreased embryonic genome activation (EGA). Comparative transcriptome analysis identified 807 maternally biased expressed genes (MBGs) and 581 paternally biased expressed genes (PBGs) in the preimplantation stages. MBGs became apparent at the 4-cell stage and contributed to the initiation of EGA, whereas PBGs preferentially appeared at the 8-cell stage and might affect embryo compaction and trophectoderm specification. Regulatory network analysis revealed that DUX4, EGR2, and DUXA are key transcription factors in MBGs' expression; ZNF263 and KLF3 are important for PBGs' expression. We demonstrated that parent-specific DNA methylation might account for the expression of most PBGs. Our results provide a valuable resource to understand parental genome activation and might help to elucidate parent-of-origin effects in early human development.

Single-cell RNA-seq reveals distinct dynamic behavior of sex chromosomes during early human embryogenesis.

Several animal and human studies have demonstrated that sex affects kinetics and metabolism during early embryo development. However, the mechanism governing these differences at the molecular level before the expression of the sex-determining gene SRY is unknown. We performed a systematic profiling of gene expression comparing male and female embryos using available single-cell RNA-sequencing data of 1607 individual cells from 99 human preimplantation embryos, covering development stages from 4-cell to late blastocyst. We observed consistent chromosome-wide transcription of autosomes, whereas expression from sex chromosomes exhibits significant differences after embryonic genome activation (EGA). Activation of the Y chromosome is initiated by expression of two genes, RPS4Y1 and DDX3Y, whereas the X chromosome is widely activated, with both copies in females being activated after EGA. In contrast to the stable activation of the Y chromosome, expression of X-linked genes in females declines at the late blastocyst stage, especially in trophectoderm cells, revealing a rapid process of dosage compensation. This dynamic behavior results in a dosage imbalance between male and female embryos, which influences genes involved in cell cycle, protein translation and metabolism. Our results reveal the dynamics of sex chromosomes expression and silencing during early embryogenesis. Studying sex differences during human embryogenesis, as well as understanding the process of X chromosome inactivation and their effects on the sex bias development of in vitro fertilized embryos, will expand the capabilities of assisted reproductive technology and possibly improve the treatment of infertility and enhance reproductive health.

An integrated chromatin accessibility and transcriptome landscape of human pre-implantation embryos.

Human pre-implantation embryonic development involves extensive changes in chromatin structure and transcriptional activity. Here, we report on LiCAT-seq, a technique that enables simultaneous profiling of chromatin accessibility and gene expression with ultra-low input of cells, and map the chromatin accessibility and transcriptome landscapes for human pre-implantation embryos. We observed global difference in chromatin accessibility between sperm and all stages of embryos, finding that the accessible regions in sperm tend to occur in gene-poor genomic regions. Integrative analyses between the two datasets reveals strong association between the establishment of accessible chromatin and embryonic genome activation (EGA), and uncovers transcription factors and endogenous retrovirus (ERVs) specific to EGA. In particular, a large proportion of the early activated genes and ERVs are bound by DUX4 and become accessible as early as the 2- to 4-cell stages. Our results thus offer mechanistic insights into the molecular events inherent to human pre-implantation development.

Deconvolution of single-cell multi-omics layers reveals regulatory heterogeneity.

Integrative analysis of multi-omics layers at single cell level is critical for accurate dissection of cell-to-cell variation within certain cell populations. Here we report scCAT-seq, a technique for simultaneously assaying chromatin accessibility and the transcriptome within the same single cell. We show that the combined single cell signatures enable accurate construction of regulatory relationships between cis-regulatory elements and the target genes at single-cell resolution, providing a new dimension of features that helps direct discovery of regulatory patterns specific to distinct cell identities. Moreover, we generate the first single cell integrated map of chromatin accessibility and transcriptome in early embryos and demonstrate the robustness of scCAT-seq in the precise dissection of master transcription factors in cells of distinct states. The ability to obtain these two layers of omics data will help provide more accurate definitions of "single cell state" and enable the deconvolution of regulatory heterogeneity from complex cell populations.

Identification of biparental and diploid blastocysts from monopronuclear zygotes with the use of a single-nucleotide polymorphism array.

OBJECTIVE: To select normal fertilized diploid blastocysts in patients who had only monopronucleated (1PN) embryos for transfer. DESIGN: Experimental study. SETTING: University-affiliated center. PATIENT(S): Couples who were undergoing intracytoplasmic sperm injection treatment and had 1PN blastocysts. INTERVENTION(S): In a preliminary test, limited cells of parthenogenetic human embryonic stem cells (phESCs) and normal fertilized blastocysts were analyzed with the use of a low-density single-nucleotide polymorphism (SNP) array to identify the distribution pattern and rate of heterozygosity. In the clinical application, 1PN blastocysts were analyzed with the use of the SNP array. Only diagnosed normal blastocysts were transferred. The diagnosed uniparental blastocysts were validated by imprinted gene expression. MAIN OUTCOME MEASURE(S): Distribution pattern and rate of heterozygosity between parthenogenesis and normal fertilization. RESULT(S): In the pretest, phESCs exhibited distinct distribution pattern and lower rate of heterozygosity, compared with normal fertilized blastocysts after SNP analysis. In particular, homozygous hESCs showed a panhomozygosity distribution pattern, hybrid phESCs showed a partial homozygosity distribution pattern, and normal fertilized blastocysts exhibited a panheterozygosity distribution pattern with an average of 20.21% heterozygosity rate; 13.6% was found to be the minimum cutoff to predict normal fertilized samples. In the clinical application, 24 1PN blastocysts were analyzed; 10/24 showed chromosomal abnormalities, 3/24 showed panhomozygosity with 0.45%-0.8% heterozygosity, and 1/24 showed partial homozygosity with 6.54% heterozygosity. The remaining 10 blastocysts, with a panheterozygosity distribution pattern and higher genomic heterozygosity rate, were diagnosed as normal-fertilization diploid embryos; three were transferred and resulted in two healthy newborns. CONCLUSION(S): The low-density SNP array might serve as a cost-effective method to identify biparental origin and diploid 1PN blastocysts for transfer.

Self-diploidization of human haploid parthenogenetic embryos through the Rho pathway regulates endomitosis and failed cytokinesis.

A diploid genome is necessary for normal mammalian development, thus haploid parthenogenetic embryos undergo frequent self-diploidization during preimplantation development; however, the underlying mechanism is unclear. In this study, time-lapse recording revealed that human haploid parthenotes (HPs) undergo self-diploidization via failed cytokinesis (FC) and endomitosis (EM). The frequencies of FC/EM were significantly higher in HPs than in normal fertilized embryos (26.3% vs. 1.6%, P < 0.01; 19.7% vs. 0, P < 0.01), and above 90% of FC/EM occurred at the first cell cycle in HPs. Fluorescent in situ hybridization of chromosome 16,18 and X in HPs identified diploid recovery after the appearance of FC/EM, and FC/EM HPs showed improved blastocyst formation compared with non-FC/EM HPs (18.8% and 40.0% vs. 15.4%, P > 0.05). In 66.7% of the 1-cell stage HPs, furrow ingression was not observed during the time for normal cleavage, and both immunostaining and gene expression analysis of 1-cell stage HPs revealed the absence or down-regulation of several key genes of the Rho pathway, which regulates cytomitosis. Our results suggested that the major mechanism for self-diploidization is Rho pathway inhibition leading to FC/EM in the first cell cycle, and fine-tuning of this signalling pathway may help to generate stable haploid embryos for stem cell biology studies.

The dynamic changes of X chromosome inactivation during early culture of human embryonic stem cells.

X chromosome inactivation (XCI) is required for dosage compensation of X-linked genes in human female cells. Several previous reports have described the promiscuous XCI status in long-term cultured female human embryonic stem cells (hESCs), and the majority of them exhibit non-random XCI. However, when and how such female hESCs acquire the aberrant XCI states during culture is unknown. Herein, through comparing the XCI states in 18 paired hES cell lines throughout early culture, we revealed a uniform dynamic change during this culture period under a widely used culture condition. The female initial hESCs (ihESCs, P4-P9) expressed XIST RNA, H3K27me3 punctate enrichment and displayed random XCI pattern. By further culturing, the female early hESCs (ehESCs, P20-P30) lost the expression of XIST RNA, H3K27me3 punctate enrichment and exhibited a completely skewed XCI pattern. Importantly, a subset of X-linked genes was up-regulated in ehESCs, including some cancer-related genes. At last, we found 5% physiological oxygen was beneficial for the expression of XIST and H3K27me3 punctate enrichment, but not for the XCI pattern. We conclude that the XCI dynamic change is a frequent epigenetic instability event during early culture, which is accompanied by the up-regulation of some X-linked genes. Furthermore, we emphasize that physiological oxygen is beneficial for XCI fidelity.

Differentiation of primordial germ cells from induced pluripotent stem cells of primary ovarian insufficiency.

STUDY QUESTION: Can the induced pluripotent stem cells (iPSCs) derived from women with primary ovarian insufficiency (POI) differentiate into germ cells for potential disease modeling in vitro? SUMMARY ANSWER: The iPSC lines derived from POI patients with 46, X, del(X)(q26) or 46, X, del(X)(q26)9qh+ could differentiate into germ cells and expressed lower levels of genes in the deletion region of the X chromosome. WHAT IS KNOWN ALREADY: iPSC technology has been envisioned as an approach for generating patient-specific stem cells for disease modeling and for developing novel therapies. It has also been confirmed that iPSCs differentiate into germ cells. STUDY DESIGN, SIZE, DURATION: We compared the differentiation ability of germ cells and the gene expression level of germ cell-related genes in the X chromosome deletion region of iPSC lines derived from POI patients (n = 2) with an iPSC line derived from normal fibroblasts (n = 1). PARTICIPANTS/MATERIALS, SETTING, METHODS: We established three iPSC lines from two patients with partial Xq deletion-induced POI and normal fibroblasts by overexpressing four factors: octamer-binding transcription factor 4 (OCT4), sex-determining region Y-box 2 (SOX2), Nanog homeobox (NANOG), and lin-28 homolog (LIN28), using lentiviral vectors. We then generated stable-transfected fluorescent reporter cell lines under the control of the Asp-Glu-Ala-Asp box polypeptide 4 (DDX4, also called VASA) promoter, and selected clonal derived sublines. We induced subline differentiation into germ cells by adding Wnt3a (30 ng/ml) and bone morphogenetic protein 4 (100 ng/ml). After 12 days of differentiation, green fluorescent protein (GFP)-positive and GFP-negative cells were isolated via fluorescence-activated cell sorting and analyzed for endogenous VASA protein (immunostaining) and for germ cell markers and genes expressed in the deleted region of the X chromosome (quantitative RT-PCR). MAIN RESULTS AND THE ROLE OF CHANCE: The POI- and normal fibroblast-derived iPSCs had typical self-renewal and pluripotency characteristics. After stable transfection with the VASA-GFP construct, the sublines POI1-iPS-V.1, POI2-iPS-V.1 and hEF-iPS-V.1 produced green fluorescent cells in the differentiated cultures, and the percentage of GFP-positive cells increased over the 12 days of differentiation to a maximum of 6.9 ± 0.33%, 5.3 ± 0.57% and 8.5 ± 0.29%, respectively, of the total cell population. Immunohistochemical analysis confirmed that endogenous VASA was enriched in the GFP-positive cells. Quantitative reverse transcription-PCR revealed significantly higher expression of germ cell markers [PR domain containing 1, with ZNF domain (PRDM1, BLIMP1), developmental pluripotency-associated 3 (DPPA3, STELLA), deleted in azoospermia-like (DAZL), and VASA (DDX4)] in GFP-positive cells than in GFP-negative cells. Moreover, the GFP-positive cells from POI-iPSCs had reduced expression of the family with sequence similarity 122C (FAM122C), inhibitor of kappa light polypeptide gene enhancer in B-cells, kinase gamma (IKBKG), and RNA binding motif protein, X-linked (RBMX), genes located in the deleted region of the X chromosome and that are highly expressed in differentiated germ cells, compared with cells from normal iPSCs. LIMITATIONS, REASONS FOR CAUTION: Gene expression profiling indicated that the germ cells differentiated from POI-iPSCs were pre-meiotic. Therefore, how the differentiated primordial germ cells could progress further to meiosis and form follicles remains to be determined in the study of POI. WIDER IMPLICATIONS OF THE FINDINGS: Our results might provide an in vitro model for studying germ cell development in patients with POI. STUDY FUNDING/COMPETING INTERESTS: This work was supported by grants from the Major State Basic Research Development Program of China (No. 2012CB944901), the National Science Foundation of China (No. 81222007 and 81471432), the Program for New Century Excellent Talents in University and the Fundamental Research Funds for Central Universities (No. 721500003). The authors have no competing interests to declare. TRIAL REGISTRATION NUMBER: Not applicable.