Mosaic autosomal aneuploidies are detectable from single-cell RNAseq data

Molecular profiling of human blastocysts reveals primitive endoderm defects among embryos of decreased implantation potential

Human embryo implantation is remarkably inefficient, and implantation failure remains among the greatest obstacles in treating infertility. Gene expression data from human embryos have accumulated rapidly in recent years; however, identification of the subset of genes that determine successful implantation remains a challenge. We leverage clinical morphologic grading-known for decades to correlate with implantation potential-and transcriptome analyses of matched embryonic and abembryonic samples to identify factors and pathways enriched and depleted in human blastocysts of good and poor morphology. Unexpectedly, we discovered that the greatest difference was in the state of extraembryonic primitive endoderm (PrE) development, with relative deficiencies in poor morphology blastocysts. Our results suggest that implantation success is most strongly influenced by the embryonic compartment and that deficient PrE development is common among embryos with decreased implantation potential. Our study provides a valuable resource for those investigating the markers and mechanisms of human embryo implantation.

Differential effects of aneuploidy on growth and differentiation in human intestinal stem cells

Aneuploidy, a near ubiquitous genetic feature of tumors, is a context-dependent driver of cancer evolution; however, the mechanistic basis of this role remains unclear. Here, by inducing heterogeneous aneuploidy in non-transformed human colon organoids (colonoids), we investigate how the effects of aneuploidy on cell growth and differentiation may promote malignant transformation. Single-cell RNA sequencing reveals that the gene expression signature across over 100 unique aneuploid karyotypes is enriched with p53 responsive genes. The primary driver of p53 activation is karyotype complexity. Complex aneuploid cells with multiple unbalanced chromosomes activate p53 and undergo G1 cell-cycle arrest, independent of DNA damage and without evidence of senescence. By contrast, simple aneuploid cells with 1-3 chromosomes gained or lost continue to proliferate, demonstrated by single cell tracking in colonoids. Notably, simple aneuploid cells exhibit impaired differentiation when niche factors are withdrawn. These findings suggest that while complex aneuploid cells are eliminated from the normal epithelium due to p53 activation, simple aneuploid cells can escape this checkpoint and may contribute to niche factor-independent growth of cancer-initiating cells.

Identification and Genome-Wide Gene Expression Perturbation of a Trisomy in Chinese Kale (Brassica oleracea var. alboglabra)

Trisomy harbouring an extra copy of the chromosome generally causes a variety of physical and intellectual disabilities in mammals but is an extremely rare and important genetic stock in plants. In this study, a spontaneous trisomy plant in a Chinese kale accession (Brassica oleracea var. alboglabra, CC, 2n = 18) that showed significantly smaller plant architecture when compared to other normal plants was found and subsequently confirmed by cytological analysis in which the chromosome set of 2n = 19 and abnormal chromosome behaviour were observed. Then, based on the gene expression deviation determined by RNA-seq, the extra chromosome copy in this trisomy was identified as chromosome C2 (TC2). Compared to normal plants, TC2 not only showed generally upregulated differentially expressed genes (DEGs) on chromosome C2 (97.21% of 573 DEGs in chromosome C2) but also exhibited a whole-genome expression perturbation, in which 1329 DEGs (69.87% of total DEGs) were observed along two-copy chromosomes (trans-effect). The genes in the high (gene expression value > 100) and medium (100 > gene expression value > 10) groups were more prone to decreased gene expression, but the genes in the low group (10 > gene expression value > 0.1) showed upregulated expression deviation. In addition, GO (Gene ontology) annotation analysis revealed that the upregulated DEGs in the trans-effect group were overrepresented by the genes involved in the response to stress category, while the downregulated DEGs in the trans-effect group were mostly enriched in pathways related to DNA synthesis. In conclusion, we think our results can provide important resources for genetic analysis in B. oleracea and show some novel insights for understanding trisomy plant biology.

Mosaic results after preimplantation genetic testing for aneuploidy may be accompanied by changes in global gene expression

Aneuploidy in preimplantation embryos is a major cause of human reproductive failure. Unlike uniformly aneuploid embryos, embryos diagnosed as diploid-aneuploid mosaics after preimplantation genetic testing for aneuploidy (PGT-A) can develop into healthy infants. However, the reason why these embryos achieve full reproductive competence needs further research. Current RNA sequencing techniques allow for the investigation of the human preimplantation transcriptome, providing new insights into the molecular mechanisms of embryo development. In this prospective study, using euploid embryo gene expression as a control, we compared the transcriptome profiles of inner cell mass and trophectoderm samples from blastocysts with different levels of chromosomal mosaicism. A total of 25 samples were analyzed from 14 blastocysts with previous PGT-A diagnosis, including five low-level mosaic embryos and four high-level mosaic embryos. Global gene expression profiles visualized in cluster heatmaps were correlated with the original PGT-A diagnosis. In addition, gene expression distance based on the number of differentially expressed genes increased with the mosaic level, compared to euploid controls. Pathways involving apoptosis, mitosis, protein degradation, metabolism, and mitochondrial energy production were among the most deregulated within mosaic embryos. Retrospective analysis of the duration of blastomere cell cycles in mosaic embryos revealed several mitotic delays compared to euploid controls, providing additional evidence of the mosaic status. Overall, these findings suggest that embryos with mosaic results are not simply a misdiagnosis by-product, but may also have a genuine molecular identity that is compatible with their reproductive potential.

Metabolic and epigenetic dysfunctions underlie the arrest of in vitro fertilized human embryos in a senescent-like state

Around 60% of in vitro fertilized (IVF) human embryos irreversibly arrest before compaction between the 3- to 8-cell stage, posing a significant clinical problem. The mechanisms behind this arrest are unclear. Here, we show that the arrested embryos enter a senescent-like state, marked by cell cycle arrest, the down-regulation of ribosomes and histones and down-regulation of MYC and p53 activity. The arrested embryos can be divided into 3 types. Type I embryos fail to complete the maternal-zygotic transition, and Type II/III embryos have low levels of glycolysis and either high (Type II) or low (Type III) levels of oxidative phosphorylation. Treatment with the SIRT agonist resveratrol or nicotinamide riboside (NR) can partially rescue the arrested phenotype, which is accompanied by changes in metabolic activity. Overall, our data suggests metabolic and epigenetic dysfunctions underlie the arrest of human embryos.

Molecular contribution to embryonic aneuploidy and karyotypic complexity in initial cleavage divisions of mammalian development

Embryonic aneuploidy is highly complex, often leading to developmental arrest, implantation failure or spontaneous miscarriage in both natural and assisted reproduction. Despite our knowledge of mitotic mis-segregation in somatic cells, the molecular pathways regulating chromosome fidelity during the error-prone cleavage-stage of mammalian embryogenesis remain largely undefined. Using bovine embryos and live-cell fluorescent imaging, we observed frequent micro-/multi-nucleation of mis-segregated chromosomes in initial mitotic divisions that underwent unilateral inheritance, re-fused with the primary nucleus or formed a chromatin bridge with neighboring cells. A correlation between a lack of syngamy, multipolar divisions and asymmetric genome partitioning was also revealed, and single-cell DNA-seq showed propagation of primarily non-reciprocal mitotic errors. Depletion of the mitotic checkpoint protein BUB1B (also known as BUBR1) resulted in similarly abnormal nuclear structures and cell divisions, as well as chaotic aneuploidy and dysregulation of the kinase-substrate network that mediates mitotic progression, all before zygotic genome activation. This demonstrates that embryonic micronuclei sustain multiple fates, provides an explanation for blastomeres with uniparental origins, and substantiates defective checkpoints and likely other maternally derived factors as major contributors to the karyotypic complexity afflicting mammalian preimplantation development.

The impact of maternal age on gene expression during the GV to MII transition in euploid human oocytes

STUDY QUESTION

Are there age-related differences in gene expression during the germinal vesicle (GV) to metaphase II (MII) stage transition in euploid human oocytes?

SUMMARY ANSWER

A decrease in mitochondrial-related transcripts from GV to MII oocytes was observed, with a much greater reduction in MII oocytes with advanced age.

WHAT IS KNOWN ALREADY

Early embryonic development is dependent on maternal transcripts accumulated and stored within the oocyte during oogenesis. Transcriptional activity of the oocyte, which dictates its ultimate developmental potential, may be influenced by age and explain the reduced competence of advanced maternal age (AMA) oocytes compared with the young maternal age (YMA). Gene expression has been studied in human and animal oocytes; however, RNA sequencing could provide further insights into the transcriptome profiling of GV and in vivo matured MII euploid oocytes of YMA and AMA patients.

STUDY DESIGN, SIZE, DURATION

Fifteen women treated for infertility in a single IVF unit agreed to participate in this study. Five GV and 5 MII oocytes from 6, 21-26 years old women (YMA cohort) and 5 GV and 6 MII oocytes from 6, 41-44 years old women (AMA cohort) undergoing IVF treatment were donated. The samples were collected within a time frame of 4 months. RNA was isolated and deep sequenced at the single-cell level. All donors provided either GV or MII oocytes.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Cumulus dissection from donated oocytes was performed 38 h after hCG injection, denuded oocytes were inserted into lysis buffer supplemented with RNase inhibitor. The samples were stored at -80°C until further use. Isolated RNA from GV and MII oocytes underwent library preparation using an oligo deoxy-thymidine (dT) priming approach (SMART-Seq v4 Ultra Low Input RNA assay; Takara Bio, Japan) and Nextera XT DNA library preparation assay (Illumina, USA) followed by deep sequencing. Data processing, quality assessment and bioinformatics analysis were performed using source-software, mainly including FastQC, HISAT2, StringTie and edgeR, along with functional annotation analysis, while scploid R package was employed to determine the ploidy status.

MAIN RESULTS AND THE ROLE OF CHANCE

Following deep sequencing of single GV and MII oocytes in both YMA and AMA cohorts, several hundred transcripts were found to be expressed at significantly different levels. When YMA and AMA MII oocyte transcriptomes were compared, the most significant of these were related to mitochondrial structure and function, including biological processes, mitochondrial respiratory chain complex I assembly and mitochondrial translational termination (false discovery rate (FDR) 6.0E-10 to 1.2E-7). These results indicate a higher energy potential of the YMA MII cohort that is reduced with ageing. Other biological processes that were significantly higher in the YMA MII cohort included transcripts involved in the translation process (FDR 1.9E-2). Lack of these transcripts could lead to inappropriate protein synthesis prior to or upon fertilisation of the AMA MII oocytes.

LARGE SCALE DATA

The RNA sequencing data were deposited in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo), under the accession number: GSE164371.

LIMITATIONS, REASONS FOR CAUTION

The relatively small sample size could be a reason for caution. However, the RNA sequencing results showed homogeneous clustering with low intra-group variation and five to six biological replicates derived from at least three different women per group minimised the potential impact of the sample size.

WIDER IMPLICATIONS OF THE FINDINGS

Understanding the effects of ageing on the oocyte transcriptome could highlight the mechanisms involved in GV to MII transition and identify biomarkers that characterise good MII oocyte quality. This knowledge has the potential to guide IVF regimes for AMA patients.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Medical Research Council (MRC Grant number MR/K020501/1).

SCReadCounts: estimation of cell-level SNVs expression from scRNA-seq data

BACKGROUND

Recent studies have demonstrated the utility of scRNA-seq SNVs to distinguish tumor from normal cells, characterize intra-tumoral heterogeneity, and define mutation-associated expression signatures. In addition to cancer studies, SNVs from single cells have been useful in studies of transcriptional burst kinetics, allelic expression, chromosome X inactivation, ploidy estimations, and haplotype inference.

RESULTS

To aid these types of studies, we have developed a tool, SCReadCounts, for cell-level tabulation of the sequencing read counts bearing SNV reference and variant alleles from barcoded scRNA-seq alignments. Provided genomic loci and expected alleles, SCReadCounts generates cell-SNV matrices with the absolute variant- and reference-harboring read counts, as well as cell-SNV matrices of expressed Variant Allele Fraction (VAFRNA) suitable for a variety of downstream applications. We demonstrate three different SCReadCounts applications on 59,884 cells from seven neuroblastoma samples: (1) estimation of cell-level expression of known somatic mutations and RNA-editing sites, (2) estimation of cell- level allele expression of biallelic SNVs, and (3) a discovery mode assessment of the reference and each of the three alternative nucleotides at genomic positions of interest that does not require prior SNV information. For the later, we applied SCReadCounts on the coding regions of KRAS, where it identified known and novel somatic mutations in a low-to-moderate proportion of cells. The SCReadCounts read counts module is benchmarked against the analogous modules of GATK and Samtools. SCReadCounts is freely available ( https://github.com/HorvathLab/NGS ) as 64-bit self-contained binary distributions for Linux and MacOS, in addition to Python source.

CONCLUSIONS

SCReadCounts supplies a fast and efficient solution for estimation of cell-level SNV expression from scRNA-seq data. SCReadCounts enables distinguishing cells with monoallelic reference expression from those with no gene expression and is applicable to assess SNVs present in only a small proportion of the cells, such as somatic mutations in cancer.

Frequent loss of heterozygosity in CRISPR-Cas9-edited early human embryos

CRISPR-Cas9 genome editing is a promising technique for clinical applications, such as the correction of disease-associated alleles in somatic cells. The use of this approach has also been discussed in the context of heritable editing of the human germ line. However, studies assessing gene correction in early human embryos report low efficiency of mutation repair, high rates of mosaicism, and the possibility of unintended editing outcomes that may have pathologic consequences. We developed computational pipelines to assess single-cell genomics and transcriptomics datasets from OCT4 (POU5F1) CRISPR-Cas9-targeted and control human preimplantation embryos. This allowed us to evaluate on-target mutations that would be missed by more conventional genotyping techniques. We observed loss of heterozygosity in edited cells that spanned regions beyond the POU5F1 on-target locus, as well as segmental loss and gain of chromosome 6, on which the POU5F1 gene is located. Unintended genome editing outcomes were present in ∼16% of the human embryo cells analyzed and spanned 4-20 kb. Our observations are consistent with recent findings indicating complexity at on-target sites following CRISPR-Cas9 genome editing. Our work underscores the importance of further basic research to assess the safety of genome editing techniques in human embryos, which will inform debates about the potential clinical use of this technology.

Depletion of aneuploid cells in human embryos and gastruloids

Chromosomal instability leading to aneuploidy is pervasive in early human embryos^1-3 and is considered as a major cause of infertility and pregnancy wastage^4,5. Here we provide several lines of evidence that blastocysts containing aneuploid cells are worthy of in vitro fertilization transfer. First, we show clinically that aneuploid embryos can lead to healthy births, suggesting the presence of an in vivo mechanism to eliminate aneuploidy. Second, early development and cell specification modelled in micropatterned human 'gastruloids' grown in confined geometry show that aneuploid cells are depleted from embryonic germ layers, but not from extraembryonic tissue, by apoptosis in a bone morphogenetic protein 4 (BMP4)-dependent manner. Third, a small percentage of euploid cells rescues embryonic tissue in mosaic gastruloids when mixed with aneuploid cells. Finally, single-cell RNA-sequencing analysis of early human embryos revealed a decline of aneuploidy beginning on day 3. Our findings challenge two current dogmas: that a single trophectoderm biopsy at blastocyst stage to perform prenatal genetic testing can accurately determine the chromosomal make-up of a human embryo, and that aneuploid embryos should be withheld from embryo transfer in association with in vitro fertilization.

Single-cell sequencing technology in tumor research

Tumor heterogeneity is a key characteristic of malignant tumors and a significant obstacle in cancer treatment and research. Although bulk tissue sequencing has wide coverage and high accuracy, it can only represent the dominant cell signal information of each sample, while masking the unique gene expression of rare cells; therefore it cannot represent genes that are unstable within a subgroup, but unchanged in a majority of cells. With the progress of genomic technology, the emergence of single-cell sequencing (SCS) has effectively solved the above problem. Genetic, transcriptomic and epigenetic sequencing at the single-cell level provides an important basis for us to correctly classify the cell subsets of heterogeneous tumor populations and to reveal the process of complex changes in tumor cells at the molecular level. Single-cell sequencing technology has been applied to the field of cancer, revealing exciting discoveries in the potential mechanisms of tumor driver gene mutation, clonal evolution, invasion and metastasis. It also provides favorable conditions for developing new tumor biomarkers and providing more accurate and individualized targeted tumor therapy. Herein, we review the steps and methods of single-cell sequencing and highlight the application of SCS in tumor diagnosis and clinical treatment.

Single-cell analysis of human embryos reveals diverse patterns of aneuploidy and mosaicism

Less than half of human zygotes survive to birth, primarily due to aneuploidies of meiotic or mitotic origin. Mitotic errors generate chromosomal mosaicism, defined by multiple cell lineages with distinct chromosome complements. The incidence and impacts of mosaicism in human embryos remain controversial, with most previous studies based on bulk DNA assays or comparisons of multiple biopsies of few embryonic cells. Single-cell genomic data provide an opportunity to quantify mosaicism on an embryo-wide scale. To this end, we extended an approach to infer aneuploidies based on dosage-associated changes in gene expression by integrating signatures of allelic imbalance. We applied this method to published single-cell RNA sequencing data from 74 human embryos, spanning the morula to blastocyst stages. Our analysis revealed widespread mosaic aneuploidies, with 59 of 74 (80%) embryos harboring at least one putative aneuploid cell (1% FDR). By clustering copy number calls, we reconstructed histories of chromosome segregation, inferring that 55 (74%) embryos possessed mitotic aneuploidies and 23 (31%) embryos possessed meiotic aneuploidies. We found no significant enrichment of aneuploid cells in the trophectoderm compared to the inner cell mass, although we do detect such enrichment in data from later postimplantation stages. Finally, we observed that aneuploid cells up-regulate immune response genes and down-regulate genes involved in proliferation, metabolism, and protein processing, consistent with stress responses documented in other stages and systems. Together, our work provides a high-resolution view of aneuploidy in preimplantation embryos, and supports the conclusion that low-level mosaicism is a common feature of early human development.

Estimating the Allele-Specific Expression of SNVs From 10× Genomics Single-Cell RNA-Sequencing Data

With the recent advances in single-cell RNA-sequencing (scRNA-seq) technologies, the estimation of allele expression from single cells is becoming increasingly reliable. Allele expression is both quantitative and dynamic and is an essential component of the genomic interactome. Here, we systematically estimate the allele expression from heterozygous single nucleotide variant (SNV) loci using scRNA-seq data generated on the 10×Genomics Chromium platform. We analyzed 26,640 human adipose-derived mesenchymal stem cells (from three healthy donors), sequenced to an average of 150K sequencing reads per cell (more than 4 billion scRNA-seq reads in total). High-quality SNV calls assessed in our study contained approximately 15% exonic and >50% intronic loci. To analyze the allele expression, we estimated the expressed variant allele fraction (VAF_RNA) from SNV-aware alignments and analyzed its variance and distribution (mono- and bi-allelic) at different minimum sequencing read thresholds. Our analysis shows that when assessing positions covered by a minimum of three unique sequencing reads, over 50% of the heterozygous SNVs show bi-allelic expression, while at a threshold of 10 reads, nearly 90% of the SNVs are bi-allelic. In addition, our analysis demonstrates the feasibility of scVAF_RNA estimation from current scRNA-seq datasets and shows that the 3′-based library generation protocol of 10×Genomics scRNA-seq data can be informative in SNV-based studies, including analyses of transcriptional kinetics.