Masoud Zamani Esteki

Single-cell paired-end genome sequencing reveals structural variation per cell cycle

The nature and pace of genome mutation is largely unknown. Because standard methods sequence DNA from populations of cells, the genetic composition of individual cells is lost, de novo mutations in cells are concealed within the bulk signal and per cell cycle mutation rates and mechanisms remain elusive. Although single-cell genome analyses could resolve these problems, such analyses are error-prone because of whole-genome amplification (WGA) artefacts and are limited in the types of DNA mutation that can be discerned. We developed methods for paired-end sequence analysis of single-cell WGA products that enable (i) detecting multiple classes of DNA mutation, (ii) distinguishing DNA copy number changes from allelic WGA-amplification artefacts by the discovery of matching aberrantly mapping read pairs among the surfeit of paired-end WGA and mapping artefacts and (iii) delineating the break points and architecture of structural variants. By applying the methods, we capture DNA copy number changes acquired over one cell cycle in breast cancer cells and in blastomeres derived from a human zygote after in vitro fertilization. Furthermore, we were able to discover and fine-map a heritable inter-chromosomal rearrangement t(1;16)(p36;p12) by sequencing a single blastomere. The methods will expedite applications in basic genome research and provide a stepping stone to novel approaches for clinical genetic diagnosis.

Nonallelic homologous recombination between retrotransposable elements is a driver of de novo unbalanced translocations

Large-scale analysis of balanced chromosomal translocation breakpoints has shown nonhomologous end joining and microhomology-mediated repair to be the main drivers of interchromosomal structural aberrations. Breakpoint sequences of de novo unbalanced translocations have not yet been investigated systematically. We analyzed 12 de novo unbalanced translocations and mapped the breakpoints in nine. Surprisingly, in contrast to balanced translocations, we identify nonallelic homologous recombination (NAHR) between (retro)transposable elements and especially long interspersed elements (LINEs) as the main mutational mechanism. This finding shows yet another involvement of (retro)transposons in genomic rearrangements and exposes a profoundly different mutational mechanism compared with balanced chromosomal translocations. Furthermore, we show the existence of compound maternal/paternal derivative chromosomes, reinforcing the hypothesis that human cleavage stage embryogenesis is a cradle of chromosomal rearrangements.

Genome-wide copy number profiling of single cells in S-phase reveals DNA-replication domains

Single-cell genomics is revolutionizing basic genome research and clinical genetic diagnosis. However, none of the current research or clinical methods for single-cell analysis distinguishes between the analysis of a cell in G1-, S- or G2/M-phase of the cell cycle. Here, we demonstrate by means of array comparative genomic hybridization that charting the DNA copy number landscape of a cell in S-phase requires conceptually different approaches to that of a cell in G1- or G2/M-phase. Remarkably, despite single-cell whole-genome amplification artifacts, the log2 intensity ratios of single S-phase cells oscillate according to early and late replication domains, which in turn leads to the detection of significantly more DNA imbalances when compared with a cell in G1- or G2/M-phase. Although these DNA imbalances may, on the one hand, be falsely interpreted as genuine structural aberrations in the S-phase cell’s copy number profile and hence lead to misdiagnosis, on the other hand, the ability to detect replication domains genome wide in one cell has important applications in DNA-replication research. Genome-wide cell-type-specific early and late replicating domains have been identified by analyses of DNA from populations of cells, but cell-to-cell differences in DNA replication may be important in genome stability, disease aetiology and various other cellular processes.

Concurrent Whole-Genome Haplotyping and Copy-Number Profiling of Single Cells

Methods for haplotyping and DNA copy-number typing of single cells are paramount for studying genomic heterogeneity and enabling genetic diagnosis. Before analyzing the DNA of a single cell by microarray or next-generation sequencing, a whole-genome amplification (WGA) process is required, but it substantially distorts the frequency and composition of the cells alleles. As a consequence, haplotyping methods suffer from error-prone discrete SNP genotypes (AA, AB, BB) and DNA copy-number profiling remains difficult because true DNA copy-number aberrations have to be discriminated from WGA artifacts. Here, we developed a single-cell genome analysis method that reconstructs genome-wide haplotype architectures as well as the copy-number and segregational origin of those haplotypes by employing phased parental genotypes and deciphering WGA-distorted SNP B-allele fractions via a process we coin haplarithmisis. We demonstrate that the method can be applied as a generic method for preimplantation genetic diagnosis on single cells biopsied from human embryos, enabling diagnosis of disease alleles genome wide as well as numerical and structural chromosomal anomalies. Moreover, meiotic segregation errors can be distinguished from mitotic ones.

Preimplantation genetic diagnosis guided by single-cell genomics

Preimplantation genetic diagnosis (PGD) aims to help couples with heritable genetic disorders to avoid the birth of diseased offspring or the recurrence of loss of conception. Following in vitro fertilization, one or a few cells are biopsied from each human preimplantation embryo for genetic testing, allowing diagnosis and selection of healthy embryos for uterine transfer. Although classical methods, including single-cell PCR and fluorescent in situ hybridization, enable PGD for many genetic disorders, they have limitations. They often require family-specific designs and can be labor intensive, resulting in long waiting lists. Furthermore, certain types of genetic anomalies are not easy to diagnose using these classical approaches, and healthy offspring carrying the parental mutant allele(s) can result. Recently, state-of-the-art methods for single-cell genomics have flourished, which may overcome the limitations associated with classical PGD, and these underpin the development of generic assays for PGD that enable selection of embryos not only for the familial genetic disorder in question, but also for various other genetic aberrations and traits at once. Here, we discuss the latest single-cell genomics methodologies based on DNA microarrays, single-nucleotide polymorphism arrays or next-generation sequence analysis. We focus on their strengths, their validation status, their weaknesses and the challenges for implementing them in PGD.

Tracing the origin of disseminated tumor cells in breast cancer using single-cell sequencing

BackgroundSingle-cell micro-metastases of solid tumors often occur in the bone marrow. These disseminated tumor cells (DTCs) may resist therapy and lay dormant or progress to cause overt bone and visceral metastases. The molecular nature of DTCs remains elusive, as well as when and from where in the tumor they originate. Here, we apply single-cell sequencing to identify and trace the origin of DTCs in breast cancer.ResultsWe sequence the genomes of 63 single cells isolated from six non-metastatic breast cancer patients. By comparing the cells’ DNA copy number aberration (CNA) landscapes with those of the primary tumors and lymph node metastasis, we establish that 53% of the single cells morphologically classified as tumor cells are DTCs disseminating from the observed tumor. The remaining cells represent either non-aberrant “normal” cells or “aberrant cells of unknown origin” that have CNA landscapes discordant from the tumor. Further analyses suggest that the prevalence of aberrant cells of unknown origin is age-dependent and that at least a subset is hematopoietic in origin. Evolutionary reconstruction analysis of bulk tumor and DTC genomes enables ordering of CNA events in molecular pseudo-time and traced the origin of the DTCs to either the main tumor clone, primary tumor subclones, or subclones in an axillary lymph node metastasis.ConclusionsSingle-cell sequencing of bone marrow epithelial-like cells, in parallel with intra-tumor genetic heterogeneity profiling from bulk DNA, is a powerful approach to identify and study DTCs, yielding insight into metastatic processes. A heterogeneous population of CNA-positive cells is present in the bone marrow of non-metastatic breast cancer patients, only part of which are derived from the observed tumor lineages.

Homozygous missense mutation in STYXL1 associated with moderate intellectual disability, epilepsy and behavioural complexities

The introduction of massive parallel sequencing has led to the identification of multiple novel genes for intellectual disability (ID) as well as epilepsy. Whereas dominant de novo mutations have been proven to be a leading cause for these disorders, they do not apply to families suggestive of an autosomal recessive inheritance pattern. In this study, we combined the use of linkage analysis with exome sequencing to elucidate the cause of moderate non-syndromic ID, epilepsy and behavioural problems in a consanguineous Asian family. A founder missense mutation was identified in STYXL1. We propose this as a novel candidate gene involved in ID, accompanied by seizures and behavioural problems. Our findings further confirm the genetic heterogeneity of cognitive disorders and genetic epilepsy.

Association of MOV10L1 gene polymorphisms and male infertility in azoospermic men with complete maturation arrest

PurposeThe present research was undertaken to study probable genetic variations of MOV10L1 in 30 infertile men that had complete maturation arrest in their spermatocyte levels and 70 fertile men as the control group.MethodsWe performed polymerase chain reaction single strand conformation polymorphism (PCR-SSCP) on extracted DNAs and sequencing was used to confirm the results. Identified polymorphisms in the MOV10L1 were further subjected to a haplotype analysis.ResultsWe identified eight single nucleotide polymorphisms (SNPs): one missense (rs2272837) and four nonsense polymorphisms (rs2272836, rs11704548, rs2272838, rs138271) in the exonic sequences and three polymorphisms (rs12170772, rs2272840, rs17248147) in the intronic regions. With the exception of rs2272838, there was a statistically significant association in all polymorphisms between study population (P < 0.05). The result of haplotyping analysis showed ten possible haplotypes, from which five were significantly increased in infertile patients compared with the control group.ConclusionsOur results suggest that MOV10L1 gene polymorphisms in the studied infertile males with complete maturation arrest are linked to infertility.

Copy Number Variation Analysis by Array Analysis of Single Cells Following Whole Genome Amplification

Whole genome amplification is required to ensure the availability of sufficient material for copy number variation analysis of a genome deriving from an individual cell. Here, we describe the protocols we use for copy number variation analysis of non-fixed single cells by array-based approaches following single-cell isolation and whole genome amplification. We are focusing on two alternative protocols, an isothermal and a PCR-based whole genome amplification method, followed by either comparative genome hybridization (aCGH) or SNP array analysis, respectively.

Abstract LB-051: Tumor heterogeneity and dissemination in breast cancer: Deep sequencing of single disseminated cells from bone marrow compared to primary tumor and lymph node metastases

Metastasis is the main cause of death amongst breast cancer patients. Our knowledge of the metastatic cascade and how to inhibit it is limited. Here we dissect the genetic profile of multiple single disseminated tumor cells (DTCs) taken at various time points after diagnosis, and compare them to their matched primary tumors and lymph node metastasis. We have previously published a method for studying CNAs in single DTCs by whole genome sequencing, where we compared two primary breast carcinomas to two corresponding DTCs. Copy number profiles from whole genome sequencing (WGS) from 40 DTCs were analyzed. The single cell whole genome amplified (WGA) DNA was used to generate WGS libraries, and the DTCs were subsequently sequenced on the Illumina HiSeq 2000. The WGS reads were trimmed for WGA adapters and aligned to GRCh37 human reference using Burrows-Wheeler Aligner (BWA). LogR values were calculated for genomic bins and corrected for% GC-bias and segmented using the piecewise constant fitting (PCF) algorithm (the penalty parameter, γ, was set to 25). Copy number was estimated per segment as 2logR × Ψ, where Ψ is the average ploidy. B allele frequency (BAF) was calculated for each known SNP position from dbSNP (dbSNP build 135) and somatic mutations read-outs generated. In this study we compared the mutation spectre and CNAs in six primary tumors, one with corresponding lymph node metastasis and single DTCs. In total, CN profile from 40 DTCs showed evidence of dissemination at both early and late stage of disease progression. At large, the copy number profile of the examined DTCs exhibited either a limited number of alterations, or a pattern similar to the primary tumor and lymph node metastasis suggesting continuous dissemination of single tumor cells throughout the tumor evolution. By demonstrating sub-clonality in the lymph node metastasis we provide novel insight into the metastatic process. Further, the correlation in aberration pattern between the lymph node metastasis and multiple DTCs, implies that cells found in the bone marrow may have originated from the lymph node metastasis. The DTCs exhibited common aberrations typically found in breast carcinomas, and several DTCs had deletion of 16q and17p, and gain of 1q, 8q. Certain DTCs exhibited CNAs not visible in the primary tumor or lymph node including gain of 9q, 14q, 19q and Xq, and loss of 2p, 6p, 8p, 18p and 19p. Two DTCs from time of diagnosis exhibited gain of the whole chromosome 5 that was not observed in the primary tumor or the lymph node. These results reveal the importance of assessing the sub-clonal genetic alterations in the primary tumor, as well as in the lymph node metastasis and DTCs, in order to evaluate patient treatment and prognosis. Citation Format: Elen Moller, Parveen Kumar, Silje Nord, David Wedge, Peter van Loo, April Peterson, Randi R. Mathiesen, Renathe Fjelldal, Masoud Z. Esteki, Jason A. Grundstad, Elin Borgen, Lars O. Baumbusch, Anne-Lise Borresen-Dale, Kevin P. White, Thierry Voet, Bjorn Naume, Vessela N. Kristensen. Tumor heterogeneity and dissemination in breast cancer: Deep sequencing of single disseminated cells from bone marrow compared to primary tumor and lymph node metastases. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr LB-051. doi:10.1158/1538-7445.AM2015-LB-051