Search

Abstract The successful accomplishment of the first telomere-to-telomere human genome assembly, T2T-CHM13, marked a milestone in achieving completeness of the human reference genome. The upcoming era of genome study will focus on fully phased diploid genome assembly, with an emphasis on genetic differences between individual haplotypes. Most existing sequencing approaches only achieved localized haplotype phasing and relied on additional pedigree information for further whole-chromosome scale phasing. The short-read-based Strand-seq method is able to directly phase single nucleotide polymorphisms (SNPs) at whole-chromosome scale but falls short when it comes to phasing structural variations (SVs). To shed light on this issue, we developed a Nanopore sequencing platform-based Strand-seq approach, which we named NanoStrand-seq. This method allowed for de novo SNP calling with high precision (99.52%) and acheived a superior phasing accuracy (0.02% Hamming error rate) at whole-chromosome scale, a level of performance comparable to Strand-seq for haplotype phasing of the GM12878 genome. Importantly, we demonstrated that NanoStrand-seq can efficiently resolve the MHC locus, a highly polymorphic genomic region. Moreover, NanoStrand-seq enabled independent direct calling and phasing of deletions and insertions at whole-chromosome level; when applied to long genomic regions of SNP homozygosity, it outperformed the strategy that combined Strand-seq with bulk long-read sequencing. Finally, we showed that, like Strand-seq, NanoStrand-seq was also applicable to primary cultured cells. Together, here we provided a novel methodology that enabled interrogation of a full spectrum of haplotype-resolved SNPs and SVs at whole-chromosome scale, with broad applications for species with diploid or even potentially polypoid genomes.

Abstract Time-stamped cross-sectional data, which lack linkage across time points, are commonly generated in single-cell transcriptional profiling. Many previous methods for inferring gene regulatory networks (GRNs) driving cell-state transitions relied on constructing single-cell temporal ordering. Introducing COSLIR (COvariance restricted Sparse LInear Regression), we presented a direct approach to reconstructing GRNs that govern cell-state transitions, utilizing only the first and second moments of samples between two consecutive time points. Simulations validated COSLIR’s perfect accuracy in the oracle case and demonstrated its robust performance in real-world scenarios. When applied to single-cell RT-PCR and RNAseq datasets in developmental biology, COSLIR competed favorably with existing methods. Notably, its running time remained nearly independent of the number of cells. Therefore, COSLIR emerges as a promising addition to GRN reconstruction methods under cell-state transitions, bypassing the single-cell temporal ordering to enhance accuracy and efficiency in single-cell transcriptional profiling.

Abstract DEAD-box helicase 17 (DDX17) is a typical member of the DEAD-box family with transcriptional cofactor activity. Although DDX17 is abundantly expressed in the myocardium, its role in heart is not fully understood. We generated cardiomyocyte-specific Ddx17-knockout mice (Ddx17-cKO), cardiomyocyte-specific Ddx17 transgenic mice (Ddx17-Tg), and various models of cardiomyocyte injury and heart failure (HF). DDX17 is downregulated in the myocardium of mouse models of heart failure and cardiomyocyte injury. Cardiomyocyte-specific knockout of Ddx17 promotes autophagic flux blockage and cardiomyocyte apoptosis, leading to progressive cardiac dysfunction, maladaptive remodeling and progression to heart failure. Restoration of DDX17 expression in cardiomyocytes protects cardiac function under pathological conditions. Further studies showed that DDX17 can bind to the transcriptional repressor B-cell lymphoma 6 (BCL6) and inhibit the expression of dynamin-related protein 1 (DRP1). When DDX17 expression is reduced, transcriptional repression of BCL6 is attenuated, leading to increased DRP1 expression and mitochondrial fission, which in turn leads to impaired mitochondrial homeostasis and heart failure. We also investigated the correlation of DDX17 expression with cardiac function and DRP1 expression in myocardial biopsy samples from patients with heart failure. These findings suggest that DDX17 protects cardiac function by promoting mitochondrial homeostasis through the BCL6-DRP1 pathway in heart failure.

Abstract Single-cell whole-genome sequencing methods have undergone great improvements over the past decade. However, allele dropout, which means the inability to detect both alleles simultaneously in an individual diploid cell, largely restricts the application of these methods particularly for medical applications. Here, we develop a new single-cell whole-genome sequencing method based on third-generation sequencing (TGS) platform named Refresh-seq (restriction fragment ligation-based genome amplification and TGS). It is based on restriction endonuclease cutting and ligation strategy in which two alleles in an individual cell can be cut into equal fragments and tend to be amplified simultaneously. As a new single-cell long-read genome sequencing method, Refresh-seq features much lower allele dropout rate compared with SMOOTH-seq. Furthermore, we apply Refresh-seq to 688 sperm cells and 272 female haploid cells (secondary polar bodies and parthenogenetic oocytes) from F1 hybrid mice. We acquire high-resolution genetic map of mouse meiosis recombination at low sequencing depth and reveal the sexual dimorphism in meiotic crossovers. We also phase the structure variations (deletions and insertions) in sperm cells and female haploid cells with high precision. Refresh-seq shows great performance in screening aneuploid sperm cells and oocytes due to the low allele dropout rate and has great potential for medical applications such as preimplantation genetic diagnosis.

Summary: Our understanding of human fetal cerebellum development during the late second trimester, a critical period for the generation of astrocytes, oligodendrocytes, and unipolar brush cells (UBCs), remains limited. Here, we performed single-cell RNA sequencing (scRNA-seq) in human fetal cerebellum samples from gestational weeks (GWs) 18–25. We find that proliferating UBC progenitors distribute in the subventricular zone of the rhombic lip (RLSVZ) near white matter (WM), forming a layer structure. We also delineate two trajectories from astrogenic radial glia (ARGs) to Bergmann glial progenitors (BGPs) and recognize oligodendrogenic radial glia (ORGs) as one source of primitive oligodendrocyte progenitor cells (PriOPCs). Additionally, our scRNA-seq analysis of the trisomy 21 fetal cerebellum at this stage reveals abnormal upregulated genes in pathways such as the cell adhesion pathway and focal adhesion pathway, which potentially promote neuronal differentiation. Overall, our research provides valuable insights into normal and abnormal development of the human fetal cerebellum.

Summary: Background: The high heterogeneity of tumour and the complexity of tumour microenvironment (TME) greatly impacted the tumour development and the prognosis of cancer in the era of immunotherapy. In this study, we aimed to portray the single cell-characterised landscape of lung adenocarcinoma (LUAD), and develop an integrated signature incorporating both tumour heterogeneity and TME for prognosis stratification. Methods: Single-cell tagged reverse transcription sequencing (STRT-seq) was performed on tumour tissues and matched normal tissues from 14 patients with LUAD for immune landscape depiction and candidate key genes selection for signature construction. Kaplan–Meier survival analyses and in-vitro cell experiments were conducted to confirm the gene functions. The transcriptomic profile of 1949 patients from 11 independent cohorts including nine public datasets and two in-house cohorts were obtained for validation. Findings: We selected 11 key genes closely related to cell-to-cell interaction, tumour development, T cell phenotype transformation, and Ma/Mo cell distribution, including HLA-DPB1, FAM83A, ITGB4, OAS1, FHL2, S100P, FSCN1, SFTPD, SPP1, DBH-AS1, CST3, and established an integrated 11-gene signature, stratifying patients to High-Score or Low-Score group for better or worse prognosis. Moreover, the prognostically-predictive potency of the signature was validated by 11 independent cohorts, and the immunotherapeutic predictive potency was also validated by our in-house cohort treated by immunotherapy. Additionally, the in-vitro cell experiments and drug sensitivity prediction further confirmed the gene function and generalizability of this signature across the entire RNA profile spectrum. Interpretation: This single cell-characterised 11-gene signature might offer insights for prognosis stratification and potential guidance for treatment selection. Funding: Support for the study was provided by National key research and development project (2022YFC2505004, 2022YFC2505000 to Z.W. and J.W.), Beijing Natural Science Foundation (7242114 to J.X.), National Natural Science Foundation of China of China (82102886 to J.X., 81871889 and 82072586 to Z.W.), Beijing Nova Program (20220484119 to J.X.), NSFC general program (82272796 to J.W.), NSFC special program (82241229 to J.W.), CAMS Innovation Fund for Medical Sciences (2021-1-I2M-012, 2022-I2M-1-009 to Z.W. and J.W.), Beijing Natural Science Foundation (7212084 to Z.W.), CAMS Key lab of translational research on lung cancer (2018PT31035 to J.W.), Aiyou Foundation (KY201701 to J.W.). Medical Oncology Key Foundation of Cancer Hospital Chinese Academy of Medical Sciences (CICAMS-MOCP2022003 to J.X.)

Summary: Ovarian endometriosis is characterized by the growth of endometrial tissue within the ovary, causing infertility and chronic pain. However, its pathophysiology remains unclear. Utilizing high-precision single-cell RNA sequencing, we profile the normal, eutopic, and ectopic endometrium from 34 individuals across proliferative and secretory phases. We observe an increased proportion of ciliated cells in both eutopic and ectopic endometrium, characterized by a diminished expression of estrogen sulfotransferase, which likely confers apoptosis resistance. After translocating to ectopic lesions, endometrial epithelium upregulates nicotinamide N-methyltransferase expression that inhibits apoptosis by promoting deacetylation and subsequent nuclear exclusion of transcription factor forkhead box protein O1, thereby leading to the downregulation of the apoptotic gene BIM. Moreover, epithelial cells in ectopic lesions elevate HLA class II complex expression, which stimulates CD4+ T cells and consequently contributes to chronic inflammation. Altogether, our study provides a comprehensive atlas of ovarian endometriosis and highlights potential therapeutic targets for modulating apoptosis and inflammation.

Abstract During embryo development, DNA methylation is established by DNMT3A/3B and subsequently maintained by DNMT1. While much research has been done in this field, the functional significance of DNA methylation in embryogenesis remains unknown. Here, we establish a system of simultaneous inactivation of multiple endogenous genes in zygotes through screening for base editors that can efficiently introduce a stop codon. Embryos with mutations in Dnmts and/or Tets can be generated in one step with IMGZ. Dnmt-null embryos display gastrulation failure at E7.5. Interestingly, although DNA methylation is absent, gastrulation-related pathways are down-regulated in Dnmt-null embryos. Moreover, DNMT1, DNMT3A, and DNMT3B are critical for gastrulation, and their functions are independent of TET proteins. Hypermethylation can be sustained by either DNMT1 or DNMT3A/3B at some promoters, which are related to the suppression of miRNAs. The introduction of a single mutant allele of six miRNAs and paternal IG-DMR partially restores primitive streak elongation in Dnmt-null embryos. Thus, our results unveil an epigenetic correlation between promoter methylation and suppression of miRNA expression for gastrulation and demonstrate that IMGZ can accelerate deciphering the functions of multiple genes in vivo.

Abstract Bone marrow mesenchymal stromal/stem cells (MSCs) are a heterogeneous population that can self-renew and generate stroma, cartilage, fat, and bone. Although a significant progress has been made toward recognizing about the phenotypic characteristics of MSCs, the true identity and properties of MSCs in bone marrow remain unclear. Here, we report the expression landscape of human fetal BM nucleated cells (BMNCs) based on the single-cell transcriptomic analysis. Unexpectedly, while the common cell surface markers such as CD146, CD271, and PDGFRa used for isolating MSCs were not detected, LIFR+PDGFRB+ were identified to be specific markers of MSCs as the early progenitors. In vivo transplantation demonstrated that LIFR+PDGFRB+CD45-CD31-CD235a- MSCs could form bone tissues and reconstitute the hematopoietic microenvironment (HME) effectively in vivo. Interestingly, we also identified a subpopulation of bone unipotent progenitor expressing TM4SF1+CD44+CD73+CD45-CD31-CD235a-, which had osteogenic potentials, but could not reconstitute HME. MSCs expressed a set of different transcription factors at the different stages of human fetal bone marrow, indicating that the stemness properties of MSCs might change during development. Moreover, transcriptional characteristics of cultured MSCs were significantly changed compared with freshly isolated primary MSCs. Our cellular profiling provides a general landscape of heterogeneity, development, hierarchy, microenvironment of the human fetal BM-derived stem cells at single-cell resolution.

High malignancy is a prominent characteristic of epithelial ovarian cancer (EOC), emphasizing the necessity for further elucidation of the potential mechanisms underlying cancer progression. Aneuploidy and copy number variation (CNV) partially contribute to the heightened malignancy observed in EOC; however, the precise features of aneuploidy and their underlying molecular patterns, as well as the relationship between CNV and aneuploidy in EOC, remain unclear. In this study, we employed single-cell sequencing data along with The Cancer Genome Atlas (TCGA) to investigate aneuploidy and CNV in EOC. The technique of fluorescence in situ hybridization (FISH) was employed using specific probes. The copy number variation within the genomic region of chromosome 8 (42754568-47889815) was assessed and utilized as a representative measure for the ploidy status of individual cells in chromosome 8. Differential expression analysis was performed between different subgroups based on chromosome 8 ploidy. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein–protein interaction (PPI), and hub–gene analyses were subsequently utilized to identify crucial genes involved. By classifying enriched tumor cells into distinct subtypes based on chromosome 8 ploidy combined with TCGA data integration, we identified key genes driving chromosome 8 aneuploidy in EOC, revealing that PRKDC gene involvement through the mediated non-homologous end-joining pathway may play a pivotal role in disease progression. Further validation through analysis of the GEO and TCGA database and survival assessment, considering both mRNA expression levels and CNV status of PRKDC, has confirmed its involvement in the progression of EOC. Further functional analysis revealed an upregulation of PRKDC in both ovarian EOC cells and tissues, with its expression showing a significant correlation with the extent of copy number variation (CNV) on chromosome 8. Taken together, CNV amplification and aneuploidy of chromosome 8 are important characteristics of EOC. PRKDC and the mediated NHEJ pathway may play a crucial role in driving aneuploidy on chromosome 8 during the progression of EOC.

Although chromosomal mosaic embryos detected by trophectoderm (TE) biopsy offer healthy embryos available for transfer, high-resolution postnatal karyotyping and chromosome testing of the transferred embryos are insufficient. Here, we applied single-cell multi-omics sequencing for seven infants with blastula chromosomal mosaicism detected by TE biopsy. The chromosome ploidy was examined by single-cell genome analysis, with the cellular identity being identified by single-cell transcriptome analysis. A total of 1616 peripheral leukocytes from seven infants with embryonic chromosomal mosaicism and three control ones with euploid TE biopsy were analyzed. A small number of blood cells showed copy number alterations (CNAs) on seemingly random locations at a frequency of 0%−2.5% per infant. However, none of the cells showed CNAs that were the same as those of the corresponding TE biopsies. The blastula chromosomal mosaicism may be fully self-corrected, probably through the selective loss of the aneuploid cells during development, and the transferred embryos can be born as euploid infants without mosaic CNAs corresponding to the TE biopsies. The results provide a new reference for the evaluations of transferring chromosomal mosaic embryos in certain situations.

Abstract Small-cell lung cancer (SCLC) is the most aggressive and lethal subtype of lung cancer, for which, better understandings of its biology are urgently needed. Single-cell sequencing technologies provide an opportunity to profile individual cells within the tumor microenvironment (TME) and investigate their roles in tumorigenic processes. Here, we performed high-precision single-cell transcriptomic analysis of ~5000 individual cells from primary tumors (PTs) and matched normal adjacent tissues (NATs) from 11 SCLC patients, including one patient with both PT and relapsed tumor (RT). The comparison revealed an immunosuppressive landscape of human SCLC. Malignant cells in SCLC tumors exhibited diverse states mainly related to the cell cycle, immune, and hypoxic properties. Our data also revealed the intratumor heterogeneity (ITH) of key transcription factors (TFs) in SCLC and related gene expression patterns and functions. The non-neuroendocrine (non-NE) tumors were correlated with increased inflammatory gene signatures and immune cell infiltrates in SCLC, which contributed to better responses to immune checkpoint inhibitors. These findings indicate a significant heterogeneity of human SCLC, and intensive crosstalk between cancer cells and the TME at single-cell resolution, and thus, set the stage for a better understanding of the biology of SCLC as well as for developing new therapeutics for SCLC.

Abstract Small bowel adenocarcinomas (SBAs) are rare malignant tumors with a high mortality rate, and their molecular characteristics are still largely unexplored. Here we performed single-cell RNA sequencing for tumor samples from 12 SBA patients and predicted drug candidates for SBA. We identified four prevalent subtypes of malignant cells with distinct signatures including cell cycle program, mitochondria program, metabolism program and epithelial–mesenchymal transition (EMT) program. The progression relationships of these four subtypes of malignant cells were also revealed, which started from the cell cycle program, through the mitochondria program and then progressing into either the metabolism program or the EMT program. Importantly, ligand–receptor interaction pairs were found to be specifically enriched in pairs of EMT-program malignant cells and highly exhausted CD8+ T cells, suggesting that cancer cell subpopulations with EMT features may contribute most to the exhaustion of T cells. We also showed that the duodenal subtype of SBA exhibited molecular features more similar to gastric cancer whereas jejunal subtype of SBA more similar to colorectal cancer. Especially, we predicted specific drugs for SBA based on differential gene expression signatures between malignant cells and normal epithelial cells of SBA, and verified more potent inhibitory effects of volasertib and tozasertib for SBA cancer cells than conventional drugs of SBA at the same concentration, which provides new clues for treatments of SBA. In summary, our study provides a blueprint of the molecular signatures of both tumor cells and tumor microenvironment cells in SBA and reveals potential targets and drug candidates for its clinical treatments.

Abstract Background Colorectal cancer (CRC) ranks as the second-leading cause of cancer-related death worldwide with metastases being the main cause of cancer-related death. Here, we investigated the genomic and transcriptomic alterations in matching adjacent normal tissues, primary tumors, and metastatic tumors of CRC patients. Methods We performed whole genome sequencing (WGS), multi-region whole exome sequencing (WES), simultaneous single-cell RNA-Seq, and single-cell targeted cDNA Sanger sequencing on matching adjacent normal tissues, primary tumors, and metastatic tumors from 12 metastatic colorectal cancer patients (n=84 for genomes, n=81 for exomes, n=9120 for single cells). Patient-derived tumor organoids were used to estimate the anti-tumor effects of a PPAR inhibitor, and self-renewal and differentiation ability of stem cell-like tumor cells. Results We found that the PPAR signaling pathway was prevalently and aberrantly activated in CRC tumors. Blocking of PPAR pathway both suppressed the growth and promoted the apoptosis of CRC organoids in vitro, indicating that aberrant activation of the PPAR signaling pathway plays a critical role in CRC tumorigenesis. Using matched samples from the same patient, distinct origins of the metastasized tumors between lymph node and liver were revealed, which was further verified by both copy number variation and mitochondrial mutation profiles at single-cell resolution. By combining single-cell RNA-Seq and single-cell point mutation identification by targeted cDNA Sanger sequencing, we revealed important phenotypic differences between cancer cells with and without critical point mutations (KRAS and TP53) in the same patient in vivo at single-cell resolution. Conclusions Our data provides deep insights into how driver mutations interfere with the transcriptomic state of cancer cells in vivo at a single-cell resolution. Our findings offer novel knowledge on metastatic mechanisms as well as potential markers and therapeutic targets for CRC diagnosis and therapy. The high-precision single-cell RNA-seq dataset of matched adjacent normal tissues, primary tumors, and metastases from CRCs may serve as a rich resource for further studies.

Abstract Background Lung cancer, one of the most common malignant tumors, exhibits high inter- and intra-tumor heterogeneity which contributes significantly to treatment resistance and failure. Single-cell RNA sequencing (scRNA-seq) has been widely used to dissect the cellular composition and characterize the molecular properties of cancer cells and their tumor microenvironment in lung cancer. However, the transcriptomic heterogeneity among various cancer cells in non-small cell lung cancer (NSCLC) warrants further illustration. Methods To comprehensively analyze the molecular heterogeneity of NSCLC, we performed high-precision single-cell RNA-seq analyses on 7364 individual cells from tumor tissues and matched normal tissues from 19 primary lung cancer patients and 1 pulmonary chondroid hamartoma patient. Results In 6 of 16 patients sequenced, we identified a significant proportion of cancer cells simultaneously expressing classical marker genes for two or even three histologic subtypes of NSCLC—adenocarcinoma (ADC), squamous cell carcinoma (SCC), and neuroendocrine tumor (NET) in the same individual cell, which we defined as mixed-lineage tumor cells; this was verified by both co-immunostaining and RNA in situ hybridization. These data suggest that mixed-lineage tumor cells are highly plastic with mixed features of different types of NSCLC. Both copy number variation (CNV) patterns and mitochondrial mutations clearly showed that the mixed-lineage and single-lineage tumor cells from the same patient had common tumor ancestors rather than different origins. Moreover, we revealed that patients with high mixed-lineage features of different cancer subtypes had worse survival than patients with low mixed-lineage features, indicating that mixed-lineage tumor features were associated with poorer prognosis. In addition, gene signatures specific to mixed-lineage tumor cells were identified, including AKR1B1. Gene knockdown and small molecule inhibition of AKR1B1 can significantly decrease cell proliferation and promote cell apoptosis, suggesting that AKR1B1 plays an important role in tumorigenesis and can serve as a candidate target for tumor therapy of NSCLC patients with mixed-lineage tumor features. Conclusions In summary, our work provides novel insights into the tumor heterogeneity of NSCLC in terms of the identification of prevalent mixed-lineage subpopulations of cancer cells with combined signatures of SCC, ADC, and NET and offers clues for potential treatment strategies in these patients.

Gonadal somatic cells are the main players in gonad development and are important for sex determination and germ cell development. Here, using a time-series single-cell RNA sequencing (scRNA-seq) strategy, we analyzed fetal germ cells (FGCs) and gonadal somatic cells in human embryos and fetuses. Clustering analysis of testes and ovaries revealed several novel cell subsets, including POU5F1+SPARC+ FGCs and KRT19+ somatic cells. Furthermore, our data indicated that the bone morphogenetic protein (BMP) signaling pathway plays cell type-specific and developmental stage-specific roles in testis development and promotes the gonocyte-to-spermatogonium transition (GST) in late-stage testicular mitotic arrest FGCs. Intriguingly, testosterone synthesis function transitioned from fetal Sertoli cells to adult Leydig cells in a stepwise manner. In our study, potential interactions between gonadal somatic cells were systematically explored and we identified cell type-specific developmental defects in both FGCs and gonadal somatic cells in a Turner syndrome embryo (45, XO). Our work provides a blueprint of the complex yet highly ordered development of and the interactions among human FGCs and gonadal somatic cells.

Abstract Background Patient-derived organoid culture is a powerful system for studying the molecular mechanisms of cancers, especially colorectal cancer (CRC), one of the most prevalent cancers worldwide. There are two main types of 3D culture methods for colonic cells, but the similarities and differences between gene expression patterns in different culture media remain largely unexplored. Results Here, we establish patient-derived organoids from colorectal cancer patients and perform single-cell RNA-Seq for pairwise samples from seven patients for both organoids and their corresponding tumor and normal tissues in vivo. We find that organoids derived from tumor tissues faithfully recapitulate the main gene expression signatures of cancer cells in vivo. On the other hand, organoids derived from normal tissues exhibited some tumor-like features at the whole transcriptome level but retained normal genomic features, such as CNVs, point mutations, and normal global DNA methylation levels, for both cultural media. More importantly, we show that conditioned medium outperforms chemical-defined medium in long-term culture of tumor epithelial cells. Finally, we mutually exchange the culture medium for the organoids and find that after interchanging the medium, the organoid cells basically maintain the transcriptome characteristics of the original medium. Conclusions Our work gives a thorough evaluation of both the cultural conditions and the biological features of organoids of CRC patients.

Abstract The epigenomic abnormality of pancreatic ductal adenocarcinoma (PDAC) has rarely been investigated due to its strong heterogeneity. Here, we used single-cell multiomics sequencing to simultaneously analyze the DNA methylome, chromatin accessibility and transcriptome in individual tumor cells of PDAC patients. We identified normal epithelial cells in the tumor lesion, which have euploid genomes, normal patterns of DNA methylation, and chromatin accessibility. Using all these normal epithelial cells as controls, we determined that DNA demethylation in the cancer genome was strongly enriched in heterochromatin regions but depleted in euchromatin regions. There were stronger negative correlations between RNA expression and promoter DNA methylation in cancer cells compared to those in normal epithelial cells. Through in-depth integrated analyses, a set of novel candidate biomarkers were identified, including ZNF667 and ZNF667-AS1, whose expressions were linked to a better prognosis for PDAC patients by affecting the proliferation of cancer cells. Our work systematically revealed the critical epigenomic features of cancer cells in PDAC patients at the single-cell level.

Abstract Background Circulating cell-free DNA (cfDNA) can be released when myocardial damage occurs. Methods Here, we used the methylated CpG tandem amplification and sequencing (MCTA-seq) method for analyzing dynamic changes in heart-derived DNA in plasma samples from myocardial infarction (MI) patients. Results We identified six CGCGCGG loci showing heart-specific hypermethylation patterns. MCTA-seq deconvolution analysis combining these loci detected heart-released cfDNA in MI patients at hospital admission, and showed that the prominently elevated total cfDNA level after percutaneous coronary intervention (PCI) was derived from both the heart and white blood cells. Furthermore, for the top marker CORO6, we developed a digital droplet PCR (ddPCR) assay that clearly detected heart damage signals in cfDNA of MI patients at hospital admission. Conclusions Our study provides insights into MI pathologies and developed a new ddPCR assay for detecting myocardial damage in clinical applications.

Non-coding RNAs exert diverse functions in many cell types. In addition to transcription factors from coding genes, non-coding RNAs may also play essential roles in shaping and directing the fate of germ cells. The presence of many long non-coding RNAs (lncRNAs) which are specifically expressed in the germ cells during human gonadal development were reported and one divergent lncRNA, LNC1845, was functionally characterized. Comprehensive bioinformatic analysis of these lncRNAs indicates that divergent lncRNAs occupied the majority of female and male germ cells. Integrating lncRNA expression into the bioinformatic analysis also enhances the cell-type classification of female germ cells. Functional dissection using in vitro differentiation of human pluripotent stem cells to germ cells revealed the regulatory role of LNC1845 on a transcription factor essential for ovarian follicle development, LHX8, by modulating the levels of histone modifications, H3K4me3 and H3K27Ac. Hence, bioinformatical analysis and experimental verification provide a comprehensive analysis of lncRNAs in developing germ cells and elucidate how an lncRNA function as a cis regulator during human germ cell development.

The full-term developmental profile of male germ cells remains undefined. Here, the authors use single-cell sequencing to investigate the transcriptome landscapes of mouse male germ cells throughout development and find several critical regulators for prenatal cell-fate determination.

Abstract There is no effective way to detect structure variations (SVs) and extra-chromosomal circular DNAs (ecDNAs) at single-cell whole-genome level. Here, we develop a novel third-generation sequencing platform-based single-cell whole-genome sequencing (scWGS) method named SMOOTH-seq (single-molecule real-time sequencing of long fragments amplified through transposon insertion). We evaluate the method for detecting CNVs, SVs, and SNVs in human cancer cell lines and a colorectal cancer sample and show that SMOOTH-seq reliably and effectively detects SVs and ecDNAs in individual cells, but shows relatively limited accuracy in detection of CNVs and SNVs. SMOOTH-seq opens a new chapter in scWGS as it generates high fidelity reads of kilobases long.

Mitochondrial DNA (mtDNA) mutations are often associated with incurable diseases and lead to detectable pathogenic variants in 1 out of 200 babies. Uncoupling of the inheritance of mtDNA and the nuclear genome by spindle transfer (ST) can potentially prevent the transmission of mtDNA mutations from mother to offspring. However, no well-established studies have critically assessed the safety of this technique. Here, using single-cell triple omics sequencing method, we systematically analyzed the genome (copy number variation), DNA methylome, and transcriptome of ST and control blastocysts. The results showed that, compared to that in control embryos, the percentage of aneuploid cells in ST embryos did not significantly change. The epiblast, primitive endoderm, and trophectoderm (TE) of ST blastocysts presented RNA expression profiles that were comparable to those of control blastocysts. However, the DNA demethylation process in TE cells of ST blastocysts was slightly slower than that in the control blastocysts. Collectively, our results suggest that ST seems generally safe for embryonic development, with a relatively minor delay in the DNA demethylation process at the blastocyst stage.

Editor’s summary_NCOMMS-19-41732B The anterior pituitary gland controls body growth and reproduction but how early development is dynamically regulated is unclear. Here, the authors use scRNA-seq of human fetal pituitaries to identify different developmental routes and state transitions of five hormone-producing cell lineages, and a hybrid epithelial/mesenchymal state of pituitary stem cells.

Arterial degeneration, closely associated with cardiovascular diseases, is driven by aging-related vascular cell-specific transcriptomics changes. This study provides a single-cell transcriptomic atlas for senile aortic and coronary arteries and underscores FOXO3A-based the transcriptional network in vasoprotection during aging.

Summary: Inhibitors of Mek1/2 and Gsk3β, known as 2i, and, together with leukemia inhibitory factor, enhance the derivation of embryonic stem cells (ESCs) and promote ground-state pluripotency (2i/L-ESCs). However, recent reports show that prolonged Mek1/2 suppression impairs developmental potential of ESCs, and is rescued by serum (S/L-ESCs). Here, we show that culturing ESCs in Activin A and BMP4, and in the absence of MEK1/2 inhibitor (ABC/L medium), establishes advanced stem cells derived from ESCs (esASCs). We demonstrate that esASCs contributed to germline lineages, full-term chimeras and generated esASC-derived mice by tetraploid complementation. We show that, in contrast to 2i/L-ESCs, esASCs display distinct molecular signatures and a stable hypermethylated epigenome, which is reversible and similar to serum-cultured ESCs. Importantly, we also derived novel ASCs (blASCs) from blastocysts in ABC/L medium. Our results provide insights into the derivation of novel ESCs with DNA hypermethylation from blastocysts in chemically defined medium. : In this article, Bao and colleagues show that ASCs cultured in Activin A and BMP4, and in the absence of MEK1/2 inhibitor, possess reversible epigenetic changes that extend their developmental potency, distinct transcriptome pattern, and a stable hypermethylated epigenome. Importantly, they also derived novel hypermethylated ASCs from blastocysts in ABC/L medium. Keywords: mouse, embryonic stem cells, epigenetics, hypermethylation, development, chemically defined, genomic imprinting, blastocyst, pluripotency, implantation

DNA methylation, chromatin accessibility, and gene expression represent different levels information in biological process, but a comprehensive multiomics analysis of the mammalian heart is lacking. Here, we applied nucleosome occupancy and methylome sequencing, which detected DNA methylation and chromatin accessibility simultaneously, as well as RNA-seq, for multiomics analysis of the 4 chambers of adult and fetal human hearts, and adult mouse hearts. Our results showed conserved region-specific patterns in the mammalian heart at transcriptome and DNA methylation level. Adult and fetal human hearts showed distinct features in DNA methylome, chromatin accessibility, and transcriptome. Novel long noncoding RNAs were identified in the human heart, and the gene expression profiles of major cardiovascular diseases associated genes were displayed. Furthermore, cross-species comparisons revealed human-specific and mouse-specific differentially expressed genes between the atria and ventricles. We also reported the relationship among multiomics and found there was a bell-shaped relationship between gene-body methylation and expression in the human heart. In general, our study provided comprehensive spatiotemporal and evolutionary insights into the regulation of gene expression in the heart.

The development of next generation sequencing (NGS) platform-based single-cell RNA sequencing (scRNA-seq) techniques has tremendously changed biological researches, while there are still many questions that cannot be addressed by them due to their short read lengths. We developed a novel scRNA-seq technology based on third-generation sequencing (TGS) platform (single-cell amplification and sequencing of full-length RNAs by Nanopore platform, SCAN-seq). SCAN-seq exhibited high sensitivity and accuracy comparable to NGS platform-based scRNA-seq methods. Moreover, we captured thousands of unannotated transcripts of diverse types, with high verification rate by reverse transcription PCR (RT-PCR)-coupled Sanger sequencing in mouse embryonic stem cells (mESCs). Then, we used SCAN-seq to analyze the mouse preimplantation embryos. We could clearly distinguish cells at different developmental stages, and a total of 27,250 unannotated transcripts from 9,338 genes were identified, with many of which showed developmental stage-specific expression patterns. Finally, we showed that SCAN-seq exhibited high accuracy on determining allele-specific gene expression patterns within an individual cell. SCAN-seq makes a major breakthrough for single-cell transcriptome analysis field.

Epigenetic dynamics, such as DNA methylation and chromatin accessibility, have been extensively explored in human preimplantation embryos. However, the active demethylation process during this crucial period remains largely unexplored. In this study, we use single-cell chemical-labeling-enabled C-to-T conversion sequencing (CLEVER-seq) to quantify the DNA 5-formylcytosine (5fC) levels of human preimplantation embryos. We find that 5-formylcytosine phosphate guanine (5fCpG) exhibits genomic element-specific distribution features and is enriched in L1 and endogenous retrovirus-K (ERVK), the subfamilies of repeat elements long interspersed nuclear elements (LINEs) and long terminal repeats (LTRs), respectively. Unlike in mice, paired pronuclei in the same zygote present variable difference of 5fCpG levels, although the male pronuclei experience stronger global demethylation. The nucleosome-occupied regions show a higher 5fCpG level compared with nucleosome-depleted ones, suggesting the role of 5fC in organizing nucleosome position. Collectively, our work offers a valuable resource for ten-eleven translocation protein family (TET)-dependent active demethylation-related study during human early embryonic development.

Summary: We show that a human pluripotent stem cell (hPSC) population cultured on a low-adhesion substrate developed two hPSC subtypes with different colony morphologies: flat and domed. Notably, the dome-like cells showed higher active proliferation capacity and increased several pluripotent genes’ expression compared with the flat monolayer cells. We further demonstrated that cell-matrix adhesion mediates the interaction between cell morphology and expression of KLF4 and KLF5 through a serum response factor (SRF)-based regulatory double loop. Our results provide a mechanistic view on the coupling among adhesion, stem cell morphology, and pluripotency, shedding light on the critical role of cell-matrix adhesion in the induction and maintenance of hPSC. : When culturing hPSCs on low-adhesion substrate (gelatin nanofiber), Dr. Liu Li and her colleagues found two subtypes with different colony morphologies. The dome-like cells showed higher proliferation capacity and KLF4/5 and NANOG expression than the monolayer cells. A serum response factor-based regulatory double loop was proposed to explain how cell-matrix adhesion mediates the interaction between cell morphology and pluripotency genes. Keywords: human pluripotent stem cells, cell-matrix adhesion, stem cell morphology, single-cell culture, heterogeneity, serum response factor, nanofiber, mathematical model

Summary: Transforming growth factor β1 (TGF-β1) plays a role in the maintenance of quiescent hematopoietic stem cells (HSCs) in vivo. We asked whether TGF-β1 controls the cell cycle status of HSCs in vitro to enhance the reconstitution activity. To examine the effect of TGF-β1 on the HSC function, we used an in vitro culture system in which single HSCs divide with the retention of their short- and long-term reconstitution ability. Extensive single-cell analyses showed that, regardless of its concentration, TGF-β1 slowed down the cell cycle progression of HSCs but consequently suppressed their self-renewal potential. Cycling HSCs were not able to go back to quiescence with TGF-β1. This study revealed a negative role of TGF-β1 in the regulation of the HSC number and reconstitution activity. : Dr. Ema and colleagues report in vitro effect of TGF-β1 on single hematopoietic stem cells (HSCs). TGF-β1 slowed down the cell cycle progression of HSCs but consequently suppressed their self-renewal potential. Cycling HSCs were unable to return to quiescence with TGF-β1. This study revealed a negative role of TGF-β1 in the regulation of the HSC number and reconstitution activity. Keywords: hematopoietic stem cells, transforming growth factor β1, cell cycle, quiescence, self-renewal, apoptosis, differentiation, G0 phase

Abstract Background Organogenesis is crucial for proper organ formation during mammalian embryonic development. However, the similarities and shared features between different organs and the cellular heterogeneity during this process at single-cell resolution remain elusive. Results We perform single-cell RNA sequencing analysis of 1916 individual cells from eight organs and tissues of E9.5 to E11.5 mouse embryos, namely, the forebrain, hindbrain, skin, heart, somite, lung, liver, and intestine. Based on the regulatory activities rather than the expression patterns, all cells analyzed can be well classified into four major groups with epithelial, mesodermal, hematopoietic, and neuronal identities. For different organs within the same group, the similarities and differences of their features and developmental paths are revealed and reconstructed. Conclusions We identify mutual interactions between epithelial and mesenchymal cells and detect epithelial cells with prevalent mesenchymal features during organogenesis, which are similar to the features of intermediate epithelial/mesenchymal cells during tumorigenesis. The comprehensive transcriptome at single-cell resolution profiled in our study paves the way for future mechanistic studies of the gene-regulatory networks governing mammalian organogenesis.

Summary: Trimethylation of histone H3 at lysine-4 (H3K4me3) is associated with eukaryotic gene promoters and poises their transcriptional activation during development. To examine the in vivo function of H3K4me3 in the absence of DNA replication, we deleted CXXC finger protein 1 (CFP1), the DNA-binding subunit of the SETD1 histone H3K4 methyltransferase, in developing oocytes. We find that CFP1 is required for H3K4me3 accumulation and the deposition of histone variants onto chromatin during oocyte maturation. Decreased H3K4me3 in oocytes caused global downregulation of transcription activity. Oocytes lacking CFP1 failed to complete maturation and were unable to gain developmental competence after fertilization, due to defects in cytoplasmic lattice formation, meiotic division, and maternal-zygotic transition. Our study highlights the importance of H3K4me3 in continuous histone replacement for transcriptional regulation, chromatin remodeling, and normal developmental progression in a non-replicative system. : Yu et al. define CFP1 as a key regulator of oocyte chromatin changes and developmental potential. Deletion of CFP1 in oocytes causes decreased H3K4me3 levels and transcription, in turn, leading to cytoskeletal defects, meiotic division, maternal-zygotic transition, and, ultimately, infertility. Keywords: CxxC1 finger protein 1, histone H3K4 trimethylation, oocyte, maternal-zygotic transition, histone exchange, female fertility, zygotic genome activation

In vitro production of human oocytes for the treatment of female infertility is a goal in reproductive medicine. Here, the authors establish in vitroconditions to generate human ovarian follicle-like cells from human embryonic stem cells.

The adult zebrafish heart is capable of regeneration but the molecular mechanisms are poorly understood. Here the authors show that chromatin remodeling factor Brg1 represses cyclin-dependent kinase inhibitors to promote myocardial regeneration.

The developmental pathway of the neural retina (NR) and retinal pigment epithelium (RPE) has been revealed by extensive research in mice. However, the molecular mechanisms underlying the development of the human NR and RPE, as well as the interactions between these two tissues, have not been well defined. Here, we analyzed 2,421 individual cells from human fetal NR and RPE using single-cell RNA sequencing (RNA-seq) technique and revealed the tightly regulated spatiotemporal gene expression network of human retinal cells. We identified major cell classes of human fetal retina and potential crucial transcription factors for each cell class. We dissected the dynamic expression patterns of visual cycle- and ligand-receptor interaction-related genes in the RPE and NR. Moreover, we provided a map of disease-related genes for human fetal retinal cells and highlighted the importance of retinal progenitor cells as potential targets of inherited retinal diseases. Our findings captured the key in vivo features of the development of the human NR and RPE and offered insightful clues for further functional studies.

Cellular senescence is a driver of various aging-associated disorders, including osteoarthritis. Here, we identified a critical role for Yes-associated protein (YAP), a major effector of Hippo signaling, in maintaining a younger state of human mesenchymal stem cells (hMSCs) and ameliorating osteoarthritis in mice. Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR associated protein 9 nuclease (Cas9)-mediated knockout (KO) of YAP in hMSCs resulted in premature cellular senescence. Mechanistically, YAP cooperated with TEA domain transcriptional factor (TEAD) to activate the expression of forkhead box D1 (FOXD1), a geroprotective protein. YAP deficiency led to the down-regulation of FOXD1. In turn, overexpression of YAP or FOXD1 rejuvenated aged hMSCs. Moreover, intra-articular administration of lentiviral vector encoding YAP or FOXD1 attenuated the development of osteoarthritis in mice. Collectively, our findings reveal YAP-FOXD1, a novel aging-associated regulatory axis, as a potential target for gene therapy to alleviate osteoarthritis.

Summary: CBX4, a component of polycomb repressive complex 1 (PRC1), plays important roles in the maintenance of cell identity and organ development through gene silencing. However, whether CBX4 regulates human stem cell homeostasis remains unclear. Here, we demonstrate that CBX4 counteracts human mesenchymal stem cell (hMSC) aging via the maintenance of nucleolar homeostasis. CBX4 protein is downregulated in aged hMSCs, whereas CBX4 knockout in hMSCs results in destabilized nucleolar heterochromatin, enhanced ribosome biogenesis, increased protein translation, and accelerated cellular senescence. CBX4 maintains nucleolar homeostasis by recruiting nucleolar protein fibrillarin (FBL) and heterochromatin protein KRAB-associated protein 1 (KAP1) at nucleolar rDNA, limiting the excessive expression of rRNAs. Overexpression of CBX4 alleviates physiological hMSC aging and attenuates the development of osteoarthritis in mice. Altogether, our findings reveal a critical role of CBX4 in counteracting cellular senescence by maintaining nucleolar homeostasis, providing a potential therapeutic target for aging-associated disorders. : Ren et al. identify a geroprotective role for CBX4 in human mesenchymal stem cells (hMSCs) by maintaining nucleolar homeostasis. Overexpression of CBX4 alleviates hMSC aging and attenuates the development of osteoarthritis, highlighting a potential avenue for the use of CBX4 gene therapy vector in treating aging and aging-related disorders. Keywords: CBX4, stem cell, aging, nucleolus, rDNA, epigenetics, heterochromatin, osteoarthritis, CRISPR/Cas9, gene editing

Summary: The heart is the central organ of the circulatory system, and its proper development is vital for maintaining human life. Here, we used single-cell RNA sequencing to profile the gene expression landscapes of ∼4,000 cardiac cells from human embryos and identified four major types of cells: cardiomyocytes (CMs), cardiac fibroblasts, endothelial cells (ECs), and valvar interstitial cells (VICs). Atrial and ventricular CMs acquired distinct features early in heart development. Furthermore, both CMs and fibroblasts show stepwise changes in gene expression. As development proceeds, VICs may be involved in the remodeling phase, and ECs display location-specific characteristics. Finally, we compared gene expression profiles between humans and mice and identified a series of unique features of human heart development. Our study lays the groundwork for elucidating the mechanisms of in vivo human cardiac development and provides potential clues to understand cardiac regeneration. : Cui et al. reveal the transcriptional landscapes of human fetal heart development at single-cell resolution and identify critical biological features of different cell types, providing insights into the molecular mechanisms of human heart development. Keywords: human fetal heart development, single-cell RNA-seq, cardiomyocytes, cardiac fibroblasts, endothelial cells, valvar interstitial cells, extracellular matrix, cross-species comparison, single-cell transcriptome analysis, human embryonic development

Summary: Healthy renal function depends on normal nephrogenesis during embryonic development. However, a comprehensive gene expression profile of human fetal kidney development remains largely unexplored. Here, using a single-cell RNA-sequencing technique, we analyzed >3,000 human fetal renal cells spanning 4 months of development in utero. Unsupervised analysis identified two progenitor subtypes during cap mesenchyme development, suggesting a mechanism for sustaining their progenitor states. Furthermore, we identified critical transcriptional regulators and signaling pathways involved in the segmentation of nephron tubules. We explored the development of the highly heterogeneous collecting duct epithelia and dissected the metabolic gene repertoire and the extracellular matrix composition of the glomerular mesangium. The results provide insights on the molecular basis and regulatory events in human renal development. Moreover, the cell-type-specific expression features of causal genes in congenital renal diseases may be helpful in the treatment of these diseases. : To investigate the gene expression profiles of human fetal kidney development, Wang et al. performed single-cell RNA-seq for >3,000 renal cells. They highlight the heterogeneity of cap mesenchyme and regulatory factors in nephron formation. The results reveal cell-type-specific risk genes in congenital renal diseases. Keywords: human fetal kidney development, single-cell RNA-seq, nephron progenitors, nephron tubule segmentation, congenital renal diseases, single cell transcriptome analysis, STRT

Methylation of cytosine at the 5 position of the pyrimidine ring is the most prevalent and significant epigenetic modifications in mammalian DNA. The CpG methylation level shows a bimodal distribution but the bimodality can be overestimated due to the heterogeneity of per-base depth. Here, we developed an algorithm to eliminate the effect of per-base depth inhomogeneity on the bimodality and obtained a random CpG methylation distribution. By quantifying the deviation of the observed methylation distribution and the random one using the information formula, we find that in tetranucleotides 5'-N5CGN3-3' (N5, N3 = A, C, G or T), GCGN3 and CCGN3 show less apparent deviation than ACGN3 and TCGN3, indicating that GCGN3 and CCGN3 are less variant in their level of methylation. The methylation variation of N5CGN3 are conserved among different cells, tissues and species, implying common features in the mechanisms of methylation and demethylation, presumably mediated by DNMTs and TETs in mammalians, respectively. Sequence dependence of DNA methylation variation also relates to gene regulatory and promotes the reexamination of the role of DNA sequence in fundamental biological processes.

Recently, Multiple Annealing and Looping-Based Amplification Cycles (MALBAC) has been developed for whole genome amplification of an individual cell, relying on quasilinear instead of exponential amplification to achieve high coverage. Here we adapt MALBAC for single-cell transcriptome amplification, which gives consistently high detection efficiency, accuracy and reproducibility. With this newly developed technique, we successfully amplified and sequenced single cells from 3 germ layers from mouse embryos in the early gastrulation stage, and examined the epithelial-mesenchymal transition (EMT) program among cells in the mesoderm layer on a single-cell level.

Stochastic and deterministic allele specific gene expression (ASE) might influence single cell phenotype, but the extent and nature of the phenomenon at the onset of early mouse development is unknown. Here we performed single cell RNA-Seq analysis of single blastomeres of mouse embryos, which revealed significant changes in the transcriptome. Importantly, over half of the transcripts with detectable genetic polymorphisms exhibit ASE, most notably, individual blastomeres from the same two-cell embryo show similar pattern of ASE. However, about 6% of them exhibit stochastic expression, indicated by altered expression ratio between the two alleles. Thus, we demonstrate that ASE is both deterministic and stochastic in early blastomeres. Furthermore, we also found that 1,718 genes express two isoforms with different lengths of 3'UTRs, with the shorter one on average 5-6 times more abundant in early blastomeres compared to the transcripts in epiblast cells, suggesting that microRNA mediated regulation of gene expression acquires increasing importance as development progresses.

MicroRNAs (miRNAs) are critical regulators of transcriptional and post-transcriptional gene silencing, which are involved in multiple developmental processes in many organisms. Apart from miRNAs, mouse germ cells express another type of small RNA, piwi-interacting RNAs (piRNAs). Although it has been clear that piRNAs play a role in repression of retrotransposons during spermatogenesis, the function of miRNA in mouse germ cells has been unclear.In this study, we first revealed the expression pattern of miRNAs by using a real-time PCR-based 220-plex miRNA expression profiling method. During development of germ cells, miR-17-92 cluster, which is thought to promote cell cycling, and the ES cell-specific cluster encoding miR-290 to -295 (miR-290-295 cluster) were highly expressed in primordial germ cells (PGCs) and spermatogonia. A set of miRNAs was developmentally regulated. We next analysed function of miRNA biogenesis in germ cell development by using conditional Dicer-knockout mice in which Dicer gene was deleted specifically in the germ cells. Dicer-deleted PGCs and spermatogonia exhibited poor proliferation. Retrotransposon activity was unexpectedly suppressed in Dicer-deleted PGCs, but not affected in the spermatogonia. In Dicer-deleted testis, spermatogenesis was retarded at an early stage when proliferation and/or early differentiation. Additionally, we analysed spermatogenesis in conditional Argonaute2-deficient mice. In contrast to Dicer-deficient testis, spermatogenesis in Argonaute2-deficient testis was indistinguishable from that in wild type.These results illustrate that miRNAs are important for the proliferation of PGCs and spermatogonia, but dispensable for the repression of retrotransposons in developing germ cells. Consistently, miRNAs promoting cell cycling are highly expressed in PGCs and spermatogonia. Furthermore, based on normal spermatogenesis in Argonaute2-deficient testis, the critical function of Dicer in spermatogenesis is independent of Argonaute2.

BackgroundSpecification of primordial germ cells (PGCs) results in the conversion of pluripotent epiblast cells into monopotent germ cell lineage. Blimp1/Prmt5 complex plays a critical role in the specification and maintenance of the early germ cell lineage. However, PGCs can be induced to dedifferentiate back to a pluripotent state as embryonic germ (EG) cells when exposed to exogenous signaling molecules, FGF-2, LIF and SCF.Methodology and principal findingsHere we show that Trichostatin A (TSA), an inhibitor of histone deacetylases, is a highly potent agent that can replace FGF-2 to induce dedifferentiation of PGCs into EG cells. A key early event during dedifferentiation of PGCs in response to FGF-2 or TSA is the down-regulation of Blimp1, which reverses and apparently relieves the cell fate restriction imposed by it. Notably, the targets of Blimp1, which include c-Myc and Klf-4, which represent two of the key factors known to promote reprogramming of somatic cells to pluripotent state, are up-regulated. We also found early activation of the LIF/Stat-3 signaling pathway with the translocation of Stat-3 into the nucleus. By contrast, while Prmt5 is retained in EG cells, it translocates from the nucleus to the cytoplasm where it probably has an independent role in regulating pluripotency.Conclusions/significanceWe propose that dedifferentiation of PGCs into EG cells may provide significant mechanistic insights on early events associated with reprogramming of committed cells to a pluripotent state.

Hypertrophic cardiomyopathy (HCM) is a common inherited cardiovascular disease, which can cause heart failure and lead to death. In this study, we performed highresolution single-cell RNA-sequencing of 2115 individual cardiomyocytes obtained from HCM patients and normal controls. Signature up- and down-regulated genes in HCM were identified by integrative analysis across 37 patients and 41 controls from our data and published human single-cell and single-nucleus RNA-seq datasets, which were further classified into gene modules by single-cell co-expression analysis. Using our high-resolution dataset, we also investigated the heterogeneity among individual cardiomyocytes and revealed five distinct clusters within HCM cardiomyocytes. Interestingly, we showed that some extracellular matrix (ECM) genes were up-regulated in the HCM cardiomyocytes, suggesting that they play a role in cardiac remodelling. Taken together, our study comprehensively profiled the transcriptomic programs of HCM cardiomyocytes and provided insights into molecular mechanisms underlying the pathogenesis of HCM. [ABSTRACT FROM AUTHOR]