Transcriptome analyses of inner cell mass and trophectoderm cells isolated by magnetic-activated cell sorting from bovine blastocysts using single cell RNA-seq

Single-cell RNA sequencing in osteoarthritis

Osteoarthritis is a progressive and heterogeneous joint disease with complex pathogenesis. The various phenotypes associated with each patient suggest that better subgrouping of tissues associated with genotypes in different phases of osteoarthritis may provide new insights into the onset and progression of the disease. Recently, single-cell RNA sequencing was used to describe osteoarthritis pathogenesis on a high-resolution view surpassing traditional technologies. Herein, this review summarizes the microstructural changes in articular cartilage, meniscus, synovium and subchondral bone that are mainly due to crosstalk amongst chondrocytes, osteoblasts, fibroblasts and endothelial cells during osteoarthritis progression. Next, we focus on the promising targets discovered by single-cell RNA sequencing and its potential applications in target drugs and tissue engineering. Additionally, the limited amount of research on the evaluation of bone-related biomaterials is reviewed. Based on the pre-clinical findings, we elaborate on the potential clinical values of single-cell RNA sequencing for the therapeutic strategies of osteoarthritis. Finally, a perspective on the future development of patient-centred medicine for osteoarthritis therapy combining other single-cell multi-omics technologies is discussed. This review will provide new insights into osteoarthritis pathogenesis on a cellular level and the field of applications of single-cell RNA sequencing in personalized therapeutics for osteoarthritis in the future.

Pre-Implantation Bovine Embryo Evaluation-From Optics to Omics and Beyond

Approximately 80% of the ~1.5 million bovine embryos transferred in 2021 were in vitro produced. However, only ~27% of the transferred IVP embryos will result in live births. The ~73% pregnancy failures are partly due to transferring poor-quality embryos, a result of erroneous stereomicroscopy-based morphological evaluation, the current method of choice for pre-transfer embryo evaluation. Numerous microscopic (e.g., differential interference contrast, electron, fluorescent, time-lapse, and artificial-intelligence-based microscopy) and non-microscopic (e.g., genomics, transcriptomics, epigenomics, proteomics, metabolomics, and nuclear magnetic resonance) methodologies have been tested to find an embryo evaluation technique that is superior to morphologic evaluation. Many of these research tools can accurately determine embryo quality/viability; however, most are invasive, expensive, laborious, technically sophisticated, and/or time-consuming, making them futile in the context of in-field embryo evaluation. However accurate they may be, using complex methods, such as RNA sequencing, SNP chips, mass spectrometry, and multiphoton microscopy, at thousands of embryo production/collection facilities is impractical. Therefore, future research is warranted to innovate field-friendly, simple benchtop tests using findings already available, particularly from omics-based research methodologies. Time-lapse monitoring and artificial-intelligence-based automated image analysis also have the potential for accurate embryo evaluation; however, further research is warranted to innovate economically feasible options for in-field applications.

The Evolution of Single-Cell RNA Sequencing Technology and Application: Progress and Perspectives

As an emerging sequencing technology, single-cell RNA sequencing (scRNA-Seq) has become a powerful tool for describing cell subpopulation classification and cell heterogeneity by achieving high-throughput and multidimensional analysis of individual cells and circumventing the shortcomings of traditional sequencing for detecting the average transcript level of cell populations. It has been applied to life science and medicine research fields such as tracking dynamic cell differentiation, revealing sensitive effector cells, and key molecular events of diseases. This review focuses on the recent technological innovations in scRNA-Seq, highlighting the latest research results with scRNA-Seq as the core technology in frontier research areas such as embryology, histology, oncology, and immunology. In addition, this review outlines the prospects for its innovative application in traditional Chinese medicine (TCM) research and discusses the key issues currently being addressed by scRNA-Seq and its great potential for exploring disease diagnostic targets and uncovering drug therapeutic targets in combination with multiomics technologies.

Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

Pluripotent Core in Bovine Embryos: A Review

Early development in mammals is characterized by the ability of each cell to produce a complete organism plus the extraembryonic, or placental, cells, defined as pluripotency. During subsequent development, pluripotency is lost, and cells begin to differentiate to a particular cell fate. This review summarizes the current knowledge of pluripotency features of bovine embryos cultured in vitro, focusing on the core of pluripotency genes (OCT4, NANOG, SOX2, and CDX2), and main chemical strategies for controlling pluripotent networks during early development. Finally, we discuss the applicability of manipulating pluripotency during the morula to blastocyst transition in cattle species.

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.

The global gene expression outline of the bovine blastocyst: reflector of environmental conditions and predictor of developmental capacity

Background

Morphological evaluation of embryos has been used to screen embryos for transfer. However, the repeatability and accuracy of this method remains low. Thus, evaluation of an embryo’s gene expression signature with respect to its developmental capacity could provide new opportunities for embryo selection. Since the gene expression outline of an embryo is considered as an aggregate of its intrinsic characteristics and culture conditions, we have compared transcriptome profiles of in vivo and in vitro derived blastocysts in relation to pregnancy outcome to unravel the discrete effects of developmental competence and environmental conditions on bovine embryo gene expression outlines. To understand whether the gene expression patterns could be associated with blastocyst developmental competency, the global transcriptome profile of in vivo (CVO) and in vitro (CVT) derived competent blastocysts that resulted in pregnancy was investigated relative to that of in vivo (NVO) and in vitro (NVT) derived blastocysts which did not establish initial pregnancy, respectively while to unravel the effects of culture condition on the transcriptome profile of embryos, the transcriptional activity of the CVO group was compared to the CVT group and the NVO group was compared to the NVT ones.

Results

A total of 700 differentially expressed genes (DEGs) were identified between CVO and NVO blastocysts. These gene transcripts represent constitutive regions, indel variants, 3′-UTR sequence variants and novel transcript regions. The majority (82%) of these DEGs, including gene clusters like ATP synthases, eukaryotic translation initiation factors, ribosomal proteins, mitochondrial ribosomal proteins, NADH dehydrogenase and cytochrome c oxidase subunits were enriched in the CVO group. These DEGs were involved in pathways associated with glycolysis/glycogenesis, citrate acid cycle, pyruvate metabolism and oxidative phosphorylation. Similarly, a total of 218 genes were differentially expressed between CVT and NVT groups. Of these, 89%, including TPT1, PDIA6, HSP90AA1 and CALM, were downregulated in the CVT group and those DEGs were overrepresented in pathways related to protein processing, endoplasmic reticulum, spliceasome, ubiquitone mediated proteolysis and steroid biosynthesis. On the other hand, although both the CVT and CVO blastocyst groups resulted in pregnancy, a total of 937 genes were differential expressed between the two groups. Compared to CVO embryos, the CVT ones exhibited downregulation of gene clusters including ribosomal proteins, mitochondrial ribosomal protein, eukaryotic translation initiation factors, ATP synthases, NADH dehydrogenase and cytochrome c oxidases. Nonetheless, downregulation of these genes could be associated with pre and postnatal abnormalities observed after transfer of in vitro embryos.

Conclusion

The present study provides a detailed inventory of differentially expressed gene signatures and pathways specifically reflective of the developmental environment and future developmental capacities of bovine embryos suggesting that transcriptome activity observed in blastocysts could be indicative of further pregnancy success but also adaptation to culture environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07693-0.

Pluripotency and Growth Factors in Early Embryonic Development of Mammals: A Comparative Approach

The regulation of early events in mammalian embryonic development is a complex process. In the early stages, pluripotency, cellular differentiation, and growth should occur at specific times and these events are regulated by different genes that are expressed at specific times and locations. The genes related to pluripotency and cellular differentiation, and growth factors that determine successful embryonic development are different (or differentially expressed) among mammalian species. Some genes are fundamental for controlling pluripotency in some species but less fundamental in others, for example, Oct4 is particularly relevant in bovine early embryonic development, whereas Oct4 inhibition does not affect ovine early embryonic development. In addition, some mechanisms that regulate cellular differentiation do not seem to be clear or evolutionarily conserved. After cellular differentiation, growth factors are relevant in early development, and their effects also differ among species, for example, insulin-like growth factor improves the blastocyst development rate in some species but does not have the same effect in mice. Some growth factors influence genes related to pluripotency, and therefore, their role in early embryo development is not limited to cell growth but could also involve the earliest stages of development. In this review, we summarize the differences among mammalian species regarding the regulation of pluripotency, cellular differentiation, and growth factors in the early stages of embryonic development.

Establishment of bovine expanded potential stem cells

Embryonic stem cells (ESCs) and induced pluripotent stem cells have the potential to differentiate to all cell types of an adult individual and are useful for studying development and for translational research. However, extrapolation of mouse and human ESC knowledge to deriving stable ESC lines of domestic ungulates and large livestock species has been challenging. In contrast to ESCs that are usually established from the blastocyst, mouse expanded potential stem cells (EPSCs) are derived from four-cell and eight-cell embryos. We have recently used the EPSC approach and established stem cells from porcine and human preimplantation embryos. EPSCs are molecularly similar across species and have broader developmental potential to generate embryonic and extraembryonic cell lineages. We further explore the EPSC technology for mammalian species refractory to the standard ESC approaches and report here the successful establishment of bovine EPSCs (bEPSCs) from preimplantation embryos of both wild-type and somatic cell nuclear transfer. bEPSCs express high levels of pluripotency genes, propagate robustly in feeder-free culture, and are genetically stable in long-term culture. bEPSCs have enriched transcriptomic features of early preimplantation embryos and differentiate in vitro to cells of the three somatic germ layers and, in chimeras, contribute to both the embryonic (fetal) and extraembryonic cell lineages. Importantly, precise gene editing is efficiently achieved in bEPSCs, and genetically modified bEPSCs can be used as donors in somatic cell nuclear transfer. bEPSCs therefore hold the potential to substantially advance biotechnology and agriculture.

Non-invasive Embryo Assessment: Altered Individual Protein Profile in Spent Culture Media from Embryos Transferred at Day 5

In order to improve ART outcome, non-invasive embryo assessment is gaining more and more attention. Therefore, the aim of this study is to determine the consecutive implantation potential via the secretome between blastocysts with or without implantation and to analyse possible interactions between these differentially expressed proteins. In this prospective study, 69 spent culture media from blastocysts transferred at day 5 were collected from patients undergoing IVF/ICSI treatment in a single IVF centre between April 2015 and November 2018 after informed consent and analysed individually. Exclusion criteria were the absence of informed consent, PCOS, endometriosis and maternal age > 42 years. Dependent on the treatment outcome, media were subsequently divided into two groups: from embryos who implanted successfully (n = 37) and from embryos without implantation (n = 32). Ninety-two proteins were measured simultaneously using the proximity extension assay (PEA) technology with the Olink® CVD III panel employing oligonucleotide-labelled antibodies. Statistical analysis was performed using the Kolmogorov-Smirnov test, Student's t test, the Mann-Whitney U test and Fisher's exact test. Media from implanted blastocysts showed significantly higher expression of EPHB4, ALCAM, CSTB, BMH, TIMP4, CCL24, SELE, FAS, JAM-A, PON3, PDGF-A, vWF and PECAM-1 compared with media from blastocysts without subsequent implantation. The highest relative expression change could be demonstrated for PECAM-1 and TIMP4. PECAM-1, SELE and vWF were co-expressed. Especially EPHB 4, SELE, ALCAM, MCP-1, CCL24, FAS, JAM-A and PDGF-A have already been described in early embryonic development and metabolism. Therefore, these proteins together with PECAM-1 indicate possible biomarkers for non-invasive embryo assessment in the future. However, due to the innovative methodology, defining a threshold for the use as biomarkers remains to be assessed.

Isolation and Identification of Bovine Preadipocytes and Screening of MicroRNAs Associated with Adipogenesis

Simple Summary

Promoting fat deposition in beef cattle has been a focus of modern animal breeding research. However, previous researchers have not examined the mechanism of adipogenesis in much detail. MicroRNAs (miRNAs) are small noncoding RNAs that play a pivotal role in adipogenesis. In this study, to explore the molecular regulatory mechanism of adipocyte differentiation and formation, bovine preadipocytes were isolated and induced into adipocytes, and then the expression patterns of miRNAs between preadipocytes and adipocytes were detected through RNA sequencing. Deep sequence analysis has identified 78, 71, and 48 novel miRNAs and 497, 491, and 524 known miRNAs in the preadipocytes, and 44, 54, and 47 novel miRNAs and 519, 522, and 504 known miRNAs in the adipocytes. Among the annotated miRNAs, 131 bovine miRNAs were upregulated in adipocytes, and 119 bovine miRNAs were downregulated in adipocytes, such as bta-miR-3604, bta-miR-23b-3p, bta-miR-26a, and bta-miR-129-3p. Bovine target gene prediction results of these miRNAs show that numerous genes are associated with lipid metabolism. These results can provide both technical support and a research basis for promoting bovine adipocyte fat deposition.

Abstract

The elucidation of the mechanisms of preadipocyte differentiation and fat accumulation in adipocytes is a major work in beef cattle breeding. As important post-transcriptional regulators, microRNAs (miRNAs) take part in cell proliferation, differentiation, apoptosis, and fat metabolism through binding seed sites of targeting mRNAs. The aim of this study was to isolate and identify bovine preadipocytes and screen miRNAs associated with adipogenesis. Bovine preadipocytes were isolated from subcutaneous fatty tissue and induced to differentiate into adipocytes. Verification of preadipocytes and adipocytes was performed by qRT-PCR (real-time quantitative reverse transcription PCR), Oil Red O staining, and immunofluorescence staining. Total RNA was extracted for small RNA sequencing. The sequencing data showed that 131 miRNAs were highly expressed in adipocytes, and 119 miRNAs were highly expressed in preadipocytes. Stem–loop qPCR (stem–loop quantitative real-time PCR) results showed that the expression patterns of 11 miRNAs were consistent with the sequencing results (miR-149-5p, miR-24-3p, miR-199a-5p, miR-33a, etc.). According to KEGG pathway and Gene Ontology (GO) analyses, multiple predicted target genes were associated with lipid metabolism. In summary, this study provides a protocol of isolating bovine preadipocytes and screening various differently expressed miRNAs during preadipocyte differentiation.

Genome-wide histone modification profiling of inner cell mass and trophectoderm of bovine blastocysts by RAT-ChIP

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) has revolutionized our understanding of chromatin-related biological processes. The method, however, requires thousands of cells and has therefore limited applications in situations where cell numbers are limited. Here we describe a novel method called Restriction Assisted Tagmentation Chromatin Immunoprecipitation (RAT-ChIP) that enables global histone modification profiling from as few as 100 cells. The method is simple, cost-effective and takes a single day to complete. We demonstrate the sensitivity of the method by deriving the first genome-wide maps of histone H3K4me3 and H3K27me3 modifications of inner cell mass and trophectoderm of bovine blastocyst stage embryos.

The expression level of SOX2 at the blastocyst stage regulates the developmental capacity of bovine embryos up to day-13 of in vitro culture

Quality of in vitro-produced embryos is influenced by changes in gene expression in response to adverse conditions. Gene markers for predicting 'good embryos' do not exist at present. We propose that the expression of pluripotency markers OCT4-SOX2-NANOG in D9 (day 9) bovine demi-embryos correlated with development at D13 (day 13). Day 8 in vitro-produced blastocysts were split in two cloned halves, one half (D9) was subjected to analysis of pluripotency markers and the other was kept in culture until D13 of development. Embryo development was scored and correlated with its own status at D9 and assigned to one of two categories: G1, arrested/dead; or G2, development up to D13. SOX2 and NANOG expression levels were significantly higher in embryos from G1 and there was also negative correlation between SOX2 and embryo survival to D13 (G3; r = -0.37; P = 0.03). We observed a significant reduction in the expression of the three studied genes from D9 to D13. Furthermore, there was a correlation between the expression of pluripotency markers at D9 and embryo diameter and the expression of trophoblastic markers at D13 (TP1-EOMES-FGF4-CDX2-TKDP1). Finally, the quotient between the relative expression of SOX2 and OCT4 in the D9 blastocysts from G1 and G2 showed that embryos that were considered as competent (G2) had a quotient close to one, while the other group had a quotient of 2.3 due to a higher expression of SOX2. These results might indicate that overexpression of SOX2 at the blastocyst stage had a negative effect on the control of embryonic developmental potential.

Specification of trophoblast from embryonic stem cells exposed to BMP4

Trophoblast (TB) comprises the outer cell layers of the mammalian placenta that make direct contact with the maternal uterus and, in species with a highly invasive placenta, maternal blood. It has its origin as trophectoderm, a single epithelial layer of extra-embryonic ectoderm that surrounds the embryo proper at the blastocyst stage of development. Here, we briefly compare the features of TB specification and determination in the mouse and the human. We then review research on a model system that has been increasingly employed to study TB emergence, namely the BMP4 (bone morphogenetic protein-4)-directed differentiation of human embryonic stem cells (ESCd), and discuss why outcomes using it have proved so uneven. We also examine the controversial aspects of this model, particularly the issue of whether or not the ESCd represents TB at all. Our focus here has been to explore similarities and potential differences between the phenotypes of ESCd, trophectoderm, placental villous TB, and human TB stem cells. We then explore the role of BMP4 in the differentiation of human pluripotent cells to TB and suggest that it converts the ESC into a totipotent state that is primed for TB differentiation when self-renewal is blocked. Finally we speculate that the TB formed from ESC is homologous to the trophectoderm-derived, invasive TB that envelopes the implanting conceptus during the second week of pregnancy.

Role of chemokine (C-C motif) ligand 24 in spatial arrangement of the inner cell mass of the bovine embryo

The process of spatial rearrangement of cells of the inner cell mass (ICM) that are destined to become hypoblast is not well understood. The observation that the chemokine (C-C motif) ligand 24 (CCL24) and several other genes involved in chemokine signaling are expressed more in the ICM than in the trophectoderm of the bovine embryo resulted in the hypothesis that CCL24 participates in spatial organization of the ICM. Temporally, expression of CCL24 in the bovine embryo occurs coincidently with blastocyst formation: transcript abundance was low until the late morula stage, peaked in the blastocyst at Day 7 of development and declined by Day 9. Treatment of embryos with two separate antagonists of C-C motif chemokine receptor 3 (the prototypical receptor for CCL24) decreased the percent of GATA6+ cells (hypoblast precursors) that were located in the outside of the ICM. Similarly, injection of zygotes with a CCL24-specific morpholino decreased the percent of GATA6+ cells in the outside of the ICM. In conclusion, CCL24 assists in spatial arrangement of the ICM in the bovine embryo. This experiment points to new functions of chemokine signaling in the bovine embryo and is consistent with the idea that cell migration is involved in the spatial organization of hypoblast cells in the blastocyst.

Bovine lineage specification revealed by single-cell gene expression analysis from zygote to blastocyst

Preimplantation embryos undergo zygotic genome activation and lineage specification resulting in three distinct cell types in the late blastocyst. The molecular mechanisms underlying this progress are largely unknown in bovines. Here, we sought to analyze an extensive set of regulators at the single-cell level to define the events involved in the development of the bovine blastocyst. Using a quantitative microfluidics approach in single cells, we analyzed mRNA levels of 96 genes known to function in early embryonic development and maintenance of stem cell pluripotency in parallel in 384 individual cells from bovine preimplantation embryos. The developmental transitions can be distinguished by distinctive gene expression profiles and we identified NOTCH1, expressed in early developmental stages, while T-box 3 (TBX3) and fibroblast growth factor receptor 4 (FGFR4), expressed in late developmental stages. Three lineages can be segregated in bovine expanded blastocysts based on the expression patterns of lineage-specific genes such as disabled homolog 2 (DAB2), caudal type homeobox 2 (CDX2), ATPase H+/K+ transporting non-gastric alpha2 subunit (ATP12A), keratin 8 (KRT8), and transcription factor AP-2 alpha (TFAP2A) for trophectoderm; GATA binding protein 6 (GATA6) and goosecoid homeobox (GSC) for primitive endoderm; and Nanog homeobox (NANOG), teratocarcinoma-derived growth factor 1 (TDGF1), and PR/SET domain 14 (PRDM14) for epiblast. Moreover, some lineage-specific genes were coexpressed in blastomeres from the morula. The commitment to trophectoderm and inner cell mass lineages in bovines occurs later than in the mouse, and KRT8 might be an earlier marker for bovine trophectoderm cells. We determined that TDGF1 and PRDM14 might play pivotal roles in the primitive endoderm and epiblast specification of bovine blastocysts. Our results shed light on early cell fate determination in bovine preimplantation embryos and offer theoretical support for deriving bovine embryonic stem cells.

Single-cell gene expression of the bovine blastocyst

The first two differentiation events in the embryo result in three cell types - epiblast, trophectoderm (TE) and hypoblast. The purpose here was to identify molecular markers for each cell type in the bovine and evaluate the differences in gene expression among individual cells of each lineage. The cDNA from 67 individual cells of dissociated blastocysts was used to determine transcript abundance for 93 genes implicated as cell lineage markers in other species or potentially involved in developmental processes. Clustering analysis indicated that the cells belonged to two major populations (clades A and B) with two subpopulations of clade A and four of clade B. Use of lineage-specific markers from other species indicated that the two subpopulations of clade A represented epiblast and hypoblast respectively while the four subpopulations of clade B were TE. Among the genes upregulated in epiblast were AJAP1, DNMT3A, FGF4, H2AFZ, KDM2B, NANOG, POU5F1, SAV1 and SLIT2 Genes overexpressed in hypoblast included ALPL, FGFR2, FN1, GATA6, GJA1, HDAC1, MBNL3, PDGFRA and SOX17, while genes overexpressed in all four TE populations were ACTA2, CDX2, CYP11A1, GATA2, GATA3, IFNT, KRT8, RAC1 and SFN The subpopulations of TE varied among each other for multiple genes including the prototypical TE marker IFNT. New markers for each cell type in the bovine blastocyst were identified. Results also indicate heterogeneity in gene expression among TE cells. Further studies are needed to confirm whether subpopulations of TE cells represent different stages in the development of a committed TE phenotype.