Resurrection of endogenous retroviruses during aging reinforces senescence.

Whether and how certain transposable elements with viral origins, such as endogenous retroviruses (ERVs) dormant in our genomes, can become awakened and contribute to the aging process is largely unknown. In human senescent cells, we found that HERVK (HML-2), the most recently integrated human ERVs, are unlocked to transcribe viral genes and produce retrovirus-like particles (RVLPs). These HERVK RVLPs constitute a transmissible message to elicit senescence phenotypes in young cells, which can be blocked by neutralizing antibodies. The activation of ERVs was also observed in organs of aged primates and mice as well as in human tissues and serum from the elderly. Their repression alleviates cellular senescence and tissue degeneration and, to some extent, organismal aging. These findings indicate that the resurrection of ERVs is a hallmark and driving force of cellular senescence and tissue aging.

Single-Cell Transcriptomic Atlas of Primate Ovarian Aging.

Molecular mechanisms of ovarian aging and female age-related fertility decline remain unclear. We surveyed the single-cell transcriptomic landscape of ovaries from young and aged non-human primates (NHPs) and identified seven ovarian cell types with distinct gene-expression signatures, including oocyte and six types of ovarian somatic cells. In-depth dissection of gene-expression dynamics of oocytes revealed four subtypes at sequential and stepwise developmental stages. Further analysis of cell-type-specific aging-associated transcriptional changes uncovered the disturbance of antioxidant signaling specific to early-stage oocytes and granulosa cells, indicative of oxidative damage as a crucial factor in ovarian functional decline with age. Additionally, inactivated antioxidative pathways, increased reactive oxygen species, and apoptosis were observed in granulosa cells from aged women. This study provides a comprehensive understanding of the cell-type-specific mechanisms underlying primate ovarian aging at single-cell resolution, revealing new diagnostic biomarkers and potential therapeutic targets for age-related human ovarian disorders.

The Transcriptome and DNA Methylome Landscapes of Human Primordial Germ Cells.

Germ cells are vital for transmitting genetic information from one generation to the next and for maintaining the continuation of species. Here, we analyze the transcriptome of human primordial germ cells (PGCs) from the migrating stage to the gonadal stage at single-cell and single-base resolutions. Human PGCs show unique transcription patterns involving the simultaneous expression of both pluripotency genes and germline-specific genes, with a subset of them displaying developmental-stage-specific features. Furthermore, we analyze the DNA methylome of human PGCs and find global demethylation of their genomes. Approximately 10 to 11 weeks after gestation, the PGCs are nearly devoid of any DNA methylation, with only 7.8% and 6.0% of the median methylation levels in male and female PGCs, respectively. Our work paves the way toward deciphering the complex epigenetic reprogramming of the germline with the aim of restoring totipotency in fertilized oocytes.

Reconstructing complex tissues from single-cell analyses.

Techniques for profiling individual cells are rapidly advancing and are providing an unprecedented opportunity for studying the genetic regulation of development and disease. In this issue, Durruthy-Durruthy et al. analyze gene expression at the single-cell level for a simple but highly organized three-dimensional structure, the mouse otocyst.

Programming and inheritance of parental DNA methylomes in mammals.

The reprogramming of parental methylomes is essential for embryonic development. In mammals, paternal 5-methylcytosines (5mCs) have been proposed to be actively converted to oxidized bases. These paternal oxidized bases and maternal 5mCs are believed to be passively diluted by cell divisions. By generating single-base resolution, allele-specific DNA methylomes from mouse gametes, early embryos, and primordial germ cell (PGC), as well as single-base-resolution maps of oxidized cytosine bases for early embryos, we report the existence of 5hmC and 5fC in both maternal and paternal genomes and find that 5mC or its oxidized derivatives, at the majority of demethylated CpGs, are converted to unmodified cytosines independent of passive dilution from gametes to four-cell embryos. Therefore, we conclude that paternal methylome and at least a significant proportion of maternal methylome go through active demethylation during embryonic development. Additionally, all the known imprinting control regions (ICRs) were classified into germ-line or somatic ICRs.

Genome analyses of single human oocytes.

Single-cell genome analyses of human oocytes are important for meiosis research and preimplantation genomic screening. However, the nonuniformity of single-cell whole-genome amplification hindered its use. Here, we demonstrate genome analyses of single human oocytes using multiple annealing and looping-based amplification cycle (MALBAC)-based sequencing technology. By sequencing the triads of the first and second polar bodies (PB1 and PB2) and the oocyte pronuclei from same female egg donors, we phase the genomes of these donors with detected SNPs and determine the crossover maps of their oocytes. Our data exhibit an expected crossover interference and indicate a weak chromatid interference. Further, the genome of the oocyte pronucleus, including information regarding aneuploidy and SNPs in disease-associated alleles, can be accurately deduced from the genomes of PB1 and PB2. The MALBAC-based preimplantation genomic screening in in vitro fertilization (IVF) enables accurate and cost-effective selection of normal fertilized eggs for embryo transfer.

scNanoSeq-CUT&Tag: a single-cell long-read CUT&Tag sequencing method for efficient chromatin modification profiling within individual cells.

Chromatin modifications are fundamental epigenetic marks that determine genome functions, but it remains challenging to profile those of repetitive elements and complex genomic regions. Here, we develop scNanoSeq-CUT&Tag, a streamlined method, by adapting modified cleavage under targets and tagmentation (CUT&Tag) to the nanopore sequencing platform for genome-wide chromatin modification profiling within individual cells. We show that scNanoSeq-CUT&Tag can accurately profile histone marks and transcription factor occupancy patterns at single-cell resolution as well as distinguish different cell types. scNanoSeq-CUT&Tag efficiently maps the allele-specific chromatin modifications and allows analysis of their neighboring region co-occupancy patterns within individual cells. Moreover, scNanoSeq-CUT&Tag can accurately detect chromatin modifications for individual copies of repetitive elements in both human and mouse genomes. Overall, we prove that scNanoSeq-CUT&Tag is a valuable single-cell tool for efficiently profiling histone marks and transcription factor occupancies, especially for previously poorly studied complex genomic regions and blacklist genomic regions.

Human Germline Cell Development: from the Perspective of Single-Cell Sequencing.

Germline cells are the beginning of new individuals in multicellular animals, including humans. Our understanding of these cell types is limited by the difficulty of analyzing the precious and heterogeneous germline tissue samples. The rapid development of single-cell sequencing technologies provides a chance for comprehensive profiling of the omics dynamics of human germline development. In this review, we discuss progress in analyzing the development of human germline cells, including preimplantation and implantation embryos, fetal germ cells (FGCs), and adult spermatogenesis by single-cell transcriptome and epigenome sequencing technologies.

scNanoHi-C: a single-cell long-read concatemer sequencing method to reveal high-order chromatin structures within individual cells.

The high-order three-dimensional (3D) organization of regulatory genomic elements provides a topological basis for gene regulation, but it remains unclear how multiple regulatory elements across the mammalian genome interact within an individual cell. To address this, herein, we developed scNanoHi-C, which applies Nanopore long-read sequencing to explore genome-wide proximal high-order chromatin contacts within individual cells. We show that scNanoHi-C can reliably and effectively profile 3D chromatin structures and distinguish structure subtypes among individual cells. This method could also be used to detect genomic variations, including copy-number variations and structural variations, as well as to scaffold the de novo assembly of single-cell genomes. Notably, our results suggest that extensive high-order chromatin structures exist in active chromatin regions across the genome, and multiway interactions between enhancers and their target promoters were systematically identified within individual cells. Altogether, scNanoHi-C offers new opportunities to investigate high-order 3D genome structures at the single-cell level.

Decoding the development of the human hippocampus.

The hippocampus is an important part of the limbic system in the human brain that has essential roles in spatial navigation and the consolidation of information from short-term memory to long-term memory1,2. Here we use single-cell RNA sequencing and assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis to illustrate the cell types, cell linage, molecular features and transcriptional regulation of the developing human hippocampus. Using the transcriptomes of 30,416 cells from the human hippocampus at gestational weeks 16-27, we identify 47 cell subtypes and their developmental trajectories. We also identify the migrating paths and cell lineages of PAX6+ and HOPX+ hippocampal progenitors, and regional markers of CA1, CA3 and dentate gyrus neurons. Multiomic data have uncovered transcriptional regulatory networks of the dentate gyrus marker PROX1. We also illustrate spatially specific gene expression in the developing human prefrontal cortex and hippocampus. The molecular features of the human hippocampus at gestational weeks 16-20 are similar to those of the mouse at postnatal days 0-5 and reveal gene expression differences between the two species. Transient expression of the primate-specific gene NBPF1 leads to a marked increase in PROX1+ cells in the mouse hippocampus. These data provides a blueprint for understanding human hippocampal development and a tool for investigating related diseases.

Single-cell multiomics analyses of spindle-transferred human embryos suggest a mostly normal embryonic development.

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.

Reconstituting the transcriptome and DNA methylome landscapes of human implantation.

Implantation is a milestone event during mammalian embryogenesis. Implantation failure is a considerable cause of early pregnancy loss in humans1. Owing to the difficulty of obtaining human embryos early after implantation in vivo, it remains unclear how the gene regulatory network and epigenetic mechanisms control the implantation process. Here, by combining an in vitro culture system for the development human embryos after implantation and single-cell multi-omics sequencing technologies, more than 8,000 individual cells from 65 human peri-implantation embryos were systematically analysed. Unsupervised dimensionality reduction and clustering algorithms of the transcriptome data show stepwise implantation routes for the epiblast, primitive endoderm and trophectoderm lineages, suggesting robust preparation for the proper establishment of a mother-to-offspring connection during implantation. Female embryos showed initiation of random X chromosome inactivation based on analysis of parental allele-specific expression of X-chromosome-linked genes during implantation. Notably, using single-cell triple omics sequencing analysis, the re-methylation of the genome in cells from the primitive endoderm lineage was shown to be much slower than in cells of both epiblast and trophectoderm lineages during the implantation process, which indicates that there are distinct re-establishment features in the DNA methylome of the epiblast and primitive endoderm-even though both lineages are derived from the inner cell mass. Collectively, our work provides insights into the complex molecular mechanisms that regulate the implantation of human embryos, and helps to advance future efforts to understanding early embryonic development and reproductive medicine.

Long-read-based single sperm genome sequencing for chromosome-wide haplotype phasing of both SNPs and SVs.

Although localized haploid phasing can be achieved using long read genome sequencing without parental data, reliable chromosome-scale phasing remains a great challenge. Given that sperm is a natural haploid cell, single-sperm genome sequencing can provide a chromosome-wide phase signal. Due to the limitation of read length, current short-read-based single-sperm genome sequencing methods can only achieve SNP haplotyping and come with difficulties in detecting and haplotyping structural variations (SVs) in complex genomic regions. To overcome these limitations, we developed a long-read-based single-sperm genome sequencing method and a corresponding data analysis pipeline that can accurately identify crossover events and chromosomal level aneuploidies in single sperm and efficiently detect SVs within individual sperm cells. Importantly, without parental genome information, our method can accurately conduct de novo phasing of heterozygous SVs as well as SNPs from male individuals at the whole chromosome scale. The accuracy for phasing of SVs was as high as 98.59% using 100 single sperm cells, and the accuracy for phasing of SNPs was as high as 99.95%. Additionally, our method reliably enabled deduction of the repeat expansions of haplotype-resolved STRs/VNTRs in single sperm cells. Our method provides a new opportunity for studying haplotype-related genetics in mammals.

Transcriptome profiling in swine macrophages infected with African swine fever virus at single-cell resolution.

African swine fever virus (ASFV) is the causative agent of African swine fever, a highly contagious and usually fatal disease in pigs. The pathogenesis of ASFV infection has not been clearly elucidated. Here, we used single-cell RNA-sequencing technology to survey the transcriptomic landscape of ASFV-infected primary porcine alveolar macrophages. The temporal dynamic analysis of viral genes revealed increased expression of viral transmembrane genes. Molecular characteristics in the ASFV-exposed cells exhibited the activation of antiviral signaling pathways with increased expression levels of interferon-stimulated genes and inflammatory- and cytokine-related genes. By comparing infected cells with unexposed cells, we showed that the unfolded protein response (UPR) pathway was activated in low viral load cells, while the expression level of UPR-related genes in high viral load cells was less than that in unexposed cells. Cells infected with various viral loads showed signature transcriptomic changes at the median progression of infection. Within the infected cells, differential expression analysis and coregulated virus­host analysis both demonstrated that ASFV promoted metabolic pathways but inhibited interferon and UPR signaling, implying the regulation pathway of viral replication in host cells. Furthermore, our results revealed that the cell apoptosis pathway was activated upon ASFV infection. Mechanistically, the production of tumor necrosis factor alpha (TNF-α) induced by ASFV infection is necessary for cell apoptosis, highlighting the importance of TNF-α in ASFV pathogenesis. Collectively, the data provide insights into the comprehensive host responses and complex virus­host interactions during ASFV infection, which may instruct future research on antiviral strategies.

Integrating single-cell datasets with ambiguous batch information by incorporating molecular network features.

With the rapid development of single-cell sequencing techniques, several large-scale cell atlas projects have been launched across the world. However, it is still challenging to integrate single-cell RNA-seq (scRNA-seq) datasets with diverse tissue sources, developmental stages and/or few overlaps, due to the ambiguity in determining the batch information, which is particularly important for current batch-effect correction methods. Here, we present SCORE, a simple network-based integration methodology, which incorporates curated molecular network features to infer cellular states and generate a unified workflow for integrating scRNA-seq datasets. Validating on real single-cell datasets, we showed that regardless of batch information, SCORE outperforms existing methods in accuracy, robustness, scalability and data integration.

Integrated transcriptomics and epigenomics reveal chamber-specific and species-specific characteristics of human and mouse hearts.

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.

A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex.

The mammalian prefrontal cortex comprises a set of highly specialized brain areas containing billions of cells and serves as the centre of the highest-order cognitive functions, such as memory, cognitive ability, decision-making and social behaviour. Although neural circuits are formed in the late stages of human embryonic development and even after birth, diverse classes of functional cells are generated and migrate to the appropriate locations earlier in development. Dysfunction of the prefrontal cortex contributes to cognitive deficits and the majority of neurodevelopmental disorders; there is therefore a need for detailed knowledge of the development of the prefrontal cortex. However, it is still difficult to identify cell types in the developing human prefrontal cortex and to distinguish their developmental features. Here we analyse more than 2,300 single cells in the developing human prefrontal cortex from gestational weeks 8 to 26 using RNA sequencing. We identify 35 subtypes of cells in six main classes and trace the developmental trajectories of these cells. Detailed analysis of neural progenitor cells highlights new marker genes and unique developmental features of intermediate progenitor cells. We also map the timeline of neurogenesis of excitatory neurons in the prefrontal cortex and detect the presence of interneuron progenitors in early developing prefrontal cortex. Moreover, we reveal the intrinsic development-dependent signals that regulate neuron generation and circuit formation using single-cell transcriptomic data analysis. Our screening and characterization approach provides a blueprint for understanding the development of the human prefrontal cortex in the early and mid-gestational stages in order to systematically dissect the cellular basis and molecular regulation of prefrontal cortex function in humans.

Relationships of tumor differentiation and immune infiltration in gastric cancers revealed by single-cell RNA-seq analyses.

Gastric cancers are highly heterogeneous malignant tumors. To reveal the relationship between differentiation status of cancer cells and tumor immune microenvironments in gastric cancer, single-cell RNA-sequencing was performed on normal mucosa tissue, differentiated gastric cancer (DGC) tissue, poorly differentiated gastric cancer (PDGC) tissue and neuroendocrine carcinoma (NEC) tissue sampled from surgically resected gastric cancer specimens. We identified the signature genes for both DGC and PDGC, and found that signature genes of PDGC strongly enriched in the epithelial-mesenchymal transition (EMT) program. Furthermore, we found that DGC tends to be immune-rich type whereas PDGC tends to be immune-poor type defined according to the density of tumor-infiltrating CD8+ T cells. Additionally, interferon alpha and gamma responding genes were specifically expressed in the immune-rich malignant cells compared with immune-poor malignant cells. Through analyzing the mixed adenoneuroendocrine carcinoma, we identified intermediate state malignant cells during the trans-differentiation process from DGC to NEC, which showed double-negative expressions of both DGC marker genes and NEC marker genes. Interferon-related pathways were gradually downregulated along the DGC to NEC trans-differentiation path, which was accompanied by reduced CD8+ cytotoxic T-cell infiltration. In summary, molecular features of both malignant cells and immune microenvironment cells of DGC, PDGC and NEC were systematically revealed, which may partially explain the strong tumor heterogeneities of gastric cancer. Especially along the DGC to NEC trans-differentiation path, immune-evasion was gradually enhanced with the decreasing activities of interferon pathway responses in malignant cells.

De novo assembly of human genome at single-cell levels.

Genome assembly has been benefited from long-read sequencing technologies with higher accuracy and higher continuity. However, most human genome assembly require large amount of DNAs from homogeneous cell lines without keeping cell heterogeneities, since cell heterogeneity could profoundly affect haplotype assembly results. Herein, using single-cell genome long-read sequencing technology (SMOOTH-seq), we have sequenced K562 and HG002 cells on PacBio HiFi and Oxford Nanopore Technologies (ONT) platforms and conducted de novo genome assembly. For the first time, we have completed the human genome assembly with high continuity (with NG50 of ∼2 Mb using 95 individual K562 cells) at single-cell levels, and explored the impact of different assemblers and sequencing strategies on genome assembly. With sequencing data from 30 diploid individual HG002 cells of relatively high genome coverage (average coverage ∼41.7%) on ONT platform, the NG50 can reach over 1.3 Mb. Furthermore, with the assembled genome from K562 single-cell dataset, more complete and accurate set of insertion events and complex structural variations could be identified. This study opened a new chapter on the practice of single-cell genome de novo assembly.

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.