scholarly journals Time and space dimensions of gene dosage imbalance of aneuploidies revealed by single cell transcriptomes

2018 ◽  
Author(s):  
Georgios Stamoulis ◽  
Marco Garieri ◽  
Periklis Makrythanasis ◽  
Audrey Letourneau ◽  
Michel Guipponi ◽  
...  

AbstractThe mechanisms underlying cellular and organismal phenotypes due to copy number alterations (CNA) are not fully understood. Aneuploidy is a major source of gene dosage imbalance due to CNA and viable human trisomies are model disorders of altered gene expression. To understand the cellular impact of gene dosage imbalance, we studied gene and allele specific expression (ASE) of 9668 single-cell fibroblasts in trisomies T21, T18, T13 and T8. To limit the bias of interindividual noise, all comparisons between euploid and trisomic single-cells were performed on an isogenic setting for all trisomies studied. Initially we examined 928 single cells with deep RNA-Seq. For T21 we used fibroblasts from one pair of monozygotic twins discordant for T21 and from mosaic T21. For T18, T13 and T8 we analyzed single cells from mosaic individuals. Single-cell analyses revealed inconsistencies concerning the overexpression of some genes observed in differential trisomic vs euploid bulk RNAseq while this imbalance was not detectable in trisomic vs. euploid single cells. Moreover, ASE profiling of all single cells uncovered a substantial monoallelic pattern of expression in the trisomic fraction of the genome. By classifying genes according to the level of mono and bi-allelic transcription, we have observed that, for genes with monoallelic and low-to-average expression, the altered gene dosage is mainly due to the higher fraction of cells simultaneously expressing these genes in the trisomic samples. These results were confirmed in a further experiment of 8740 single fibroblasts from the monozygotic twins discordant for T21 samples. We conclude that gene dosage imbalance is of bidimensional nature: over time (simultaneous expression of all alleles resulting in increased accumulation of RNA of copy altered genes in each single cell) as previously stated, and over space (increased fraction of cells simultaneously expressing copy altered genes). These results strongly suggest that each class of genes contributes to the phenotypic variability of trisomies according to its temporal and spatial behavior and propose an improved model to understand the effects of copy number alterations.

2019 ◽  
Vol 10 (1) ◽  
Author(s):  
Georgios Stamoulis ◽  
Marco Garieri ◽  
Periklis Makrythanasis ◽  
Audrey Letourneau ◽  
Michel Guipponi ◽  
...  

Abstract Aneuploidy is a major source of gene dosage imbalance due to copy number alterations (CNA), and viable human trisomies are model disorders of altered gene expression. We study gene and allele-specific expression (ASE) of 9668 single-cell fibroblasts from trisomy 21 (T21) discordant twins and from mosaic T21, T18, T13 and T8. We examine 928 single cells with deep scRNAseq. Expected and observed overexpression of trisomic genes in trisomic vs. diploid bulk RNAseq is not detectable in trisomic vs. diploid single cells. Instead, for trisomic genes with low-to-average expression, their altered gene dosage is mainly due to the higher fraction of trisomic cells simultaneously expressing these genes, in agreement with a stochastic 2-state burst-like model of transcription. These results, confirmed in a further analysis of 8740 single fibroblasts with shallow scRNAseq, suggest that the specific transcriptional profile of each gene contributes to the phenotypic variability of trisomies. We propose an improved model to understand the effects of CNA and, generally, of gene regulation on gene dosage imbalance.


2021 ◽  
Author(s):  
Lu Liu ◽  
He Chen ◽  
Cheng Sun ◽  
Jianyun Zhang ◽  
Juncheng Wang ◽  
...  

Genomic-scale somatic copy number alterations in healthy humans are difficult to investigate because of low occurrence rates and the structural variations' stochastic natures. Using a Tn5-transposase assisted single-cell whole genome sequencing method, we sequenced over 20,000 single lymphocytes from 16 individuals. Then, with the scale increased to a few thousand single cells per individual, we found that about 7.5% of the cells had large-size copy number alterations. Trisomy 21 was the most prevalent aneuploid event among all autosomal copy number alterations, while monosomy X occurred most frequently in over-30-year-old females. In the monosomy X single cells from individuals with phased genomes and identified X- inactivation ratios in bulk, the inactive X Chromosomes were lost more often than were the active ones.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 122-122
Author(s):  
Nicola E Potter ◽  
Luca Ermini ◽  
Elli Papaemmanuil ◽  
Gowri Vijayaraghavan ◽  
Ian Titley ◽  
...  

Abstract Abstract 122 Cancer clone development is widely regarded as an evolutionary or Darwinian process of genetic diversification and natural (or therapeutic) selection within tissue ecosystems. Emerging studies are providing strong evidence that dynamic and complex branching sub-clonal genetic architectures are a common feature of cancer (Greaves M and Maley CC Nature 2012). This complexity may underpin the intransigence of advanced cancer to therapeutic control, particularly as the critical 'driver' cells – cancer or leukaemic stem cells, also appear to be genetically diverse within individual patients (Anderson K et al Nature 2011, Notta F et al Nature 2011). Sub-clonal architecture can only be fully determined through the study of large numbers of single cells uniformly sampled from the individual cancer of interest and assessed for composite genotype. Various technologies and approaches from fluorescent in situ hybridisation (FISH) to whole-genome sequencing of single cells have been applied to cancer and leukaemic cells but each approach has limitations. We have developed a novel multiplex microfluidic Q-PCR approach that allows unbiased single cell sampling, high throughput analysis of hundreds of individual cells and simultaneous detection of multiple genetic alterations in a single cell, including fusion genes, DNA copy number alterations (CNAs) and sequence-based mutations. As a proof of principle study we have applied this technique to REH, an acute lymphoblastic leukaemia (ALL) cell line that harbors the ETV6-RUNX1 fusion and a SNP in the EPO receptor gene, which we used as a surrogate mutation. We further determined a detailed sub-clonal genetic architecture for two ETV6-RUNX1 positive ALL patient samples with multiple point mutations and copy number alterations (determined by whole-genome sequencing) by interrogating approximately 400 flow cytometry sorted single cells with validation by FISH and standard sequencing. Briefly, single cells were lysed prior to multiplex specific (DNA) target amplification (STA) and Q-PCR using the 96.96 dynamic microfluidic array and the BioMarkï HD (Fluidigm, UK). Phylogenetic trees were constructed using maximum parsimony with PAUP analysis software. Interrogation of REH revealed that all single cells registered the ETV6-RUNX1 fusion and EPO receptor SNP, but 42% of cells gained either 1 or 2 additional copies of chromosome 21. Patient sample data revealed branching sub-clonal architectures in Case A in which all leukaemic cells harbored the fusion with additional point mutations but only sub-clones showed CNAs. In contrast, the sub-clonal architecture of Case B showed that whilst the ETV6-RUNX1 fusion was the earliest (or universal) genomic event, CNAs were relatively early events preceding the acquisition of point mutations (Figure 1). In both cases, the numerically predominant sub-clone harbored both point mutations and CNAs in addition to the presumptive initiating lesion, ETV6-RUNX1. These detailed and complex sub-clonal architectures would be masked by other genetic techniques. Single cell genetics coupled with deep genome sequencing is now technically feasible and provides an accurate portrait of the dynamic clonal complexity in leukaemia (and other cancers). Variegated genetics and clonal complexity in individual leukaemias has important implications for our understanding of molecular pathogenesis and for therapeutic targeting. Figure 1. This sub-clonal genetic architecture depicts the branching structure found for Case B, illustrating that in this case the ETV6-RUNX1 fusion was the earliest genomic event, followed by CNAs and the acquisition of point mutations. Those populations highlighted grey are within the experimental error rate but potentially true populations. Figure 1. This sub-clonal genetic architecture depicts the branching structure found for Case B, illustrating that in this case the ETV6-RUNX1 fusion was the earliest genomic event, followed by CNAs and the acquisition of point mutations. Those populations highlighted grey are within the experimental error rate but potentially true populations. Disclosures: No relevant conflicts of interest to declare.


2021 ◽  
Author(s):  
Lu Liu ◽  
He Chen ◽  
Cheng Sun ◽  
Jianyun Zhang ◽  
Juncheng Wang ◽  
...  

Genomic-scale somatic copy number alterations in healthy humans are difficult to investigate because of low occurrence rates and the structural variations’ stochastic natures. Using a Tn5-transposase-assisted single-cell whole-genome sequencing method, we sequenced over 20,000 single lymphocytes from 16 individuals. Then, with the scale increased to a few thousand single cells per individual, we found that about 7.5% of the cells had large-size copy number alterations. Trisomy 21 was the most prevalent aneuploid event among all autosomal copy number alterations, whereas monosomy X occurred most frequently in over-30-yr-old females. In the monosomy X single cells from individuals with phased genomes and identified X-inactivation ratios in bulk, the inactive X Chromosomes were lost more often than the active ones.


2020 ◽  
Vol 22 (Supplement_3) ◽  
pp. iii408-iii408
Author(s):  
Marina Danilenko ◽  
Masood Zaka ◽  
Claire Keeling ◽  
Stephen Crosier ◽  
Rafiqul Hussain ◽  
...  

Abstract Medulloblastomas harbor clinically-significant intra-tumoral heterogeneity for key biomarkers (e.g. MYC/MYCN, β-catenin). Recent studies have characterized transcriptional heterogeneity at the single-cell level, however the underlying genomic copy number and mutational architecture remains to be resolved. We therefore sought to establish the intra-tumoural genomic heterogeneity of medulloblastoma at single-cell resolution. Copy number patterns were dissected by whole-genome sequencing in 1024 single cells isolated from multiple distinct tumour regions within 16 snap-frozen medulloblastomas, representing the major molecular subgroups (WNT, SHH, Group3, Group4) and genotypes (i.e. MYC amplification, TP53 mutation). Common copy number driver and subclonal events were identified, providing clear evidence of copy number evolution in medulloblastoma development. Moreover, subclonal whole-arm and focal copy number alterations covering important genomic loci (e.g. on chr10 of SHH patients) were detected in single tumour cells, yet undetectable at the bulk-tumor level. Spatial copy number heterogeneity was also common, with differences between clonal and subclonal events detected in distinct regions of individual tumours. Mutational analysis of the cells allowed dissection of spatial and clonal heterogeneity patterns for key medulloblastoma mutations (e.g. CTNNB1, TP53, SMARCA4, PTCH1) within our cohort. Integrated copy number and mutational analysis is underway to establish their inter-relationships and relative contributions to clonal evolution during tumourigenesis. In summary, single-cell analysis has enabled the resolution of common mutational and copy number drivers, alongside sub-clonal events and distinct patterns of clonal and spatial evolution, in medulloblastoma development. We anticipate these findings will provide a critical foundation for future improved biomarker selection, and the development of targeted therapies.


2019 ◽  
Author(s):  
Haoyun Lei ◽  
Bochuan Lyu ◽  
E. Michael Gertz ◽  
Alejandro A. Schaeffer ◽  
Xulian Shi ◽  
...  

2019 ◽  
Vol 47 (19) ◽  
pp. e122-e122
Author(s):  
Ramya Viswanathan ◽  
Elsie Cheruba ◽  
Lih Feng Cheow

Abstract Genome-wide profiling of copy number alterations and DNA methylation in single cells could enable detailed investigation into the genomic and epigenomic heterogeneity of complex cell populations. However, current methods to do this require complex sample processing and cleanup steps, lack consistency, or are biased in their genomic representation. Here, we describe a novel single-tube enzymatic method, DNA Analysis by Restriction Enzyme (DARE), to perform deterministic whole genome amplification while preserving DNA methylation information. This method was evaluated on low amounts of DNA and single cells, and provides accurate copy number aberration calling and representative DNA methylation measurement across the whole genome. Single-cell DARE is an attractive and scalable approach for concurrent genomic and epigenomic characterization of cells in a heterogeneous population.


Blood ◽  
2021 ◽  
Vol 138 (Supplement 1) ◽  
pp. 1590-1590
Author(s):  
Mehmet K. Samur ◽  
Anil Aktas-Samur ◽  
Romain Lannes ◽  
Jill Corre ◽  
Anjan Thakurta ◽  
...  

Abstract New generation immunotherapies in Multiple Myeloma (MM) targeting BCMA, have shown remarkable clinical benefits. However relapse still occurs due to tumor intrinsic and extrisic resistance mechanisms including antigen loss related to mutation, deletion and splicing pattern changes. Two recent case reports including ours highlighted biallelic loss of BCMA as a cause for resistance to anti-BCMA targeting therapy. In both studies BCMA locus at 16p was deleted bringing in focus importance of del16p. Here, we have evaluated 2883 MM patients at diagnosis and relapse to understand frequency characteristics of somatic events targeting BCMA. We first evaluated the frequency of deletion involving the BCMA locus (16p13.13) in MM patients from multiple studies using WGS sequencing data as well as using Affymetrix Cytoscan HD and SNP 6.0 arrays. We observed del16p in 8.58 % (7.6% to 14.6% in individual studies) of newly-diagnosed patients (n=2458). Similar frequency was observed in relapsed MM patients not previously exposed to BCMA targeting therapy. Next, we evaluated genome wide copy number alterations (CNAs) in all patients with loss of BCMA locus and observed similar frequency of loss in both hyperdiploid MM (HMM) and non-HMM suggesting its independence from cytogentic subtypes of MM. Overall copy number loss was significantly higher in patients with BCMA loss compared to rest of the MM patients. Patients with loss of BCMA locus have increased mutational load (8202 with 95% HDI 6921 and 9535) compared to those without BCMA locus loss (6975 with 95% HDI 6626 - 7343); probability of difference greater than 0 was 96.8% and difference of the means were 1222 [95% CI -112 - 2589] We next evaluated co-occurrence of BCMA loss with other high risk events and observed del1p and del17p as being significantly associated with loss of BCMA locus [Odds ratio 19.37 (13.13-25.80), FDR = 1.57e-65; and 8.8 (6.39-12.15), FDR = 5.57E-39, respectively)]. Furthermore, we observed that when both BCMA and TP53 loss are present, they have same log ratio (sequencing) or smoothed copy numbers (SNP array). Similarly, we used CDKN2C as a proxy to chromosome 1p loss and observed that when both BCMA and CKDN2C loss are present in the same patient they tend to show similar copy number values. These data suggested a possibility of co-occurrence of these events in the same cell. To further investigate this observation, we used single cell DNA sequencing data from patients with sub clonal and clonal BCMA locus loss. scDNA sequencing showed that almost all cells with BCMA deletion also had TP53 deletion (95%). Interestingly, almost all cells with BCMA loss also had p53 loss, while not all cells with p53 loss had BCMA loss suggesting that the chronology of this copy number alternation may suggest first p53 loss followed by BCMA loss. We further investigated whether a bi-allelic BCMA loss was observed after anti-BCMA targeted CAR-T cell therapy by imputing the copy number alterations using single cell RNA sequencing data. Our data from this case also indicated that BCMA loss tend to co-occur with TP53 deletions (OR=5.67 [95% CI 4.12-7.84], p value < 0.0001). Moreover, TP53 mutations were also more frequent in patients with del16p and del17p, compared to patients who only had del16p or del17p. In summary, our data from large scale copy number profiles at the diagnosis and relapse showed that monoallelic BCMA deletions are frequent events, patients with these events show increased aneuploidy, mostly deletions, potentially making these cells vulnerable for biallelic loss of genes, especially under the pressure of targeted therapy. Our results also highlight that BCMA expressions in bulk sample may not detect the presence or absence of cells with target loss and therefore combining strategies at bulk and single cell level are necessary to understand the disease status. These results suggest the need to study del16p in patients being targeted for BCMA-directed therapy and its association with other risk factors in MM. Disclosures Thakurta: Bristol Myers Squibb: Current Employment, Current equity holder in publicly-traded company. Anderson: Celgene: Membership on an entity's Board of Directors or advisory committees; Janssen: Membership on an entity's Board of Directors or advisory committees; Gilead: Membership on an entity's Board of Directors or advisory committees; Sanofi-Aventis: Membership on an entity's Board of Directors or advisory committees; Millenium-Takeda: Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Membership on an entity's Board of Directors or advisory committees; Pfizer: Membership on an entity's Board of Directors or advisory committees; Scientific Founder of Oncopep and C4 Therapeutics: Current equity holder in publicly-traded company, Current holder of individual stocks in a privately-held company; AstraZeneca: Membership on an entity's Board of Directors or advisory committees; Mana Therapeutics: Membership on an entity's Board of Directors or advisory committees. Munshi: Takeda: Consultancy; Adaptive Biotechnology: Consultancy; Amgen: Consultancy; Karyopharm: Consultancy; Celgene: Consultancy; Abbvie: Consultancy; Oncopep: Consultancy, Current equity holder in publicly-traded company, Other: scientific founder, Patents & Royalties; Novartis: Consultancy; Legend: Consultancy; Pfizer: Consultancy; Janssen: Consultancy; Bristol-Myers Squibb: Consultancy.


Sign in / Sign up

Export Citation Format

Share Document