Search

Mammalian genomes encode over a hundred different helicases, many of which are implicated in the repair of DNA lesions by acting on DNA structures arising during DNA replication, recombination or transcription. Defining the in vivo substrates of such DNA helicases is a major challenge given the large number of helicases in the genome, the breadth of potential substrates in the genome and the degree of genetic pleiotropy among DNA helicases in resolving diverse substrates. Helicases such as WRN, BLM and RECQL5 are implicated in the resolution of error-free recombination events known as sister chromatid exchange events (SCEs). Single cell Strand-seq can be used to map the genomic location of individual SCEs at a resolution that exceeds that of classical cytogenetic techniques by several orders of magnitude. By mapping the genomic locations of SCEs in the absence of different helicases, it should in principle be possible to infer the substrate specificity of specific helicases. Here we describe how the genome can be interrogated for such DNA repair events using single-cell template strand sequencing (Strand-seq) and bioinformatic tools. SCEs and copy-number alterations were mapped to genomic locations at kilobase resolution in haploid KBM7 cells. Strategies, possibilities, and limitations of Strand-seq to study helicase function are illustrated using these cells before and after CRISPR/Cas9 knock out of WRN, BLM and/or RECQL5.

BACKGROUND: Circulating tumour cells (CTCs) can be used to monitor cancer longitudinally, but their use in non-small cell lung cancer (NSCLC) is limited due to low numbers in the peripheral blood. Through diagnostic leukapheresis (DLA) CTCs can be obtained from larger blood volumes.METHODS: Patients with all stages of NSCLC were selected. One total body blood volume was screened by DLA before and after treatment. Peripheral blood was drawn pre- and post DLA for CTC enumeration by CellSearch. CTCs were detected in the DLA product (volume equalling 2 × 108 leucocytes) and after leucocyte depletion (RosetteSep, 9 mL DLA product). Single-cell, whole-genome sequencing was performed on isolated CTCs.RESULTS: Fifty-six patients were included. Before treatment, CTCs were more often detected in DLA (32/55, 58%) than in the peripheral blood (pre-DLA: 18/55, 33%; post DLA: 13/55, 23%, both at p CONCLUSIONS: DLA detected nine times more CTCs than in the peripheral blood. The sustained presence of CTCs in DLA after treatment was associated with therapy failure and shortened PFS.TRIAL REGISTRATION: The study was approved by the Medical Ethical Committee (NL55754.042.15) and was registered in the Dutch trial register (NL5423).

Following DNA replication, equal amounts of chromatin proteins are distributed over sister chromatids by re-deposition of parental chromatin proteins and deposition of newly synthesized chromatin proteins. Molecular mechanisms balancing the allocation of new and old chromatin proteins remain largely unknown. Here, we studied the genome-wide distribution of new chromatin proteins relative to parental DNA template strands and replication initiation zones using the double-click-seq. Under control conditions, new chromatin proteins were preferentially found on DNA replicated by the lagging strand machinery. Strikingly, replication stress induced by hydroxyurea or curaxin treatment and inhibition of ataxia telangiectasia and Rad3-related protein (ATR) or p53 inactivation inverted the observed chromatin protein deposition bias to the strand replicated by the leading strand polymerase in line with previously reported effects on replication protein A occupancy. We propose that asymmetric deposition of replication protein occupancy. propose asymmetric deposition newly synthesized chromatin proteins onto sister chromatids reflects differences in the processivity of leading and lagging strand synthesis.

Dense local haplotypes can now readily be extracted from long-read or droplet-based sequence data. However, these methods struggle to combine subchromosomal haplotype blocks into global chromosome-length haplotypes. Strand-seq is a single cell sequencing technique that uses read orientation to capture sparse global phase information by sequencing only one of two DNA strands for each parental homolog. In combination with dense local haplotypes from other technologies, Strand-seq data can be used to obtain complete chromosome-length phase information. In this chapter, we run the R package StrandPhaseR to phase SNVs using publicly available sequence data for sample HG005 of the Genome in a Bottle project.

Aneuploidy and chromosomal instability are both commonly found in cancer. Chromosomal instability leads to karyotype heterogeneity in tumors and is associated with therapy resistance, metastasis and poor prognosis. It has been hypothesized that aneuploidy per se is sufficient to drive CIN, however due to limited models and heterogenous results, it has remained controversial which aspects of aneuploidy can drive CIN. In this study we systematically tested the impact of different types of aneuploidies on the induction of CIN. We generated a plethora of isogenic aneuploid clones harboring whole chromosome or segmental aneuploidies in human p53-deficient RPE-1 cells. We observed increased segregation errors in cells harboring trisomies that strongly correlated to the number of gained genes. Strikingly, we found that clones harboring only monosomies do not induce a CIN phenotype. Finally, we found that an initial chromosome breakage event and subsequent fusion can instigate breakage-fusion-bridge cycles. By investigating the impact of monosomies, trisomies and segmental aneuploidies on chromosomal instability we further deciphered the complex relationship between aneuploidy and CIN.

Single-cell Strand-seq generates directional genomic information to study DNA repair, assemble genomes, and map structural variation onto chromosome-length haplotypes. We report a nanoliter-volume, one-pot (OP) Strand-seq library preparation protocol in which reagents are added cumulatively, DNA purification steps are avoided, and enzymes are inactivated with a thermolabile protease. OP-Strand-seq libraries capture 10%-25% of the genome from a single-cell with reduced costs and increased throughput.

SummarySister chromatid exchanges (SCEs) are products of joint DNA molecule resolution, and are considered to form through homologous recombination (HR). Indeed, upon generation of irradiation-induced DNA breaks, SCE induction was compromised in cells deficient for canonical HR factors BRCA1, BRCA2 and RAD51. Contrarily, replication-blocking agents, including PARP inhibitors, induced SCEs independently of BRCA1, BRCA2 and RAD51. PARP inhibitor-induced SCEs were enriched at common fragile sites (CFSs), and were accompanied by post-replicative single-stranded DNA (ssDNA) gaps. Moreover, PARP inhibitor-induced replication lesions were transmitted into mitosis, suggesting that SCEs originate from mitotic processing of under-replicated DNA. We found that DNA polymerase theta (POLQ) was recruited to mitotic DNA lesions, and loss of POLQ resulted in reduced SCE numbers and severe chromosome fragmentation upon PARP inhibition in HR-deficient cells. Combined, our data show that PARP inhibition generates under-replicated DNA, which is transferred into mitosis and processed into SCEs, independently of canonical HR factors.

Aneuploidy and chromosomal instability are both commonly found in cancer. Chromosomal instability leads to karyotype heterogeneity in tumors and is associated with therapy resistance, metastasis and poor prognosis. It has been hypothesized that aneuploidy per se is sufficient to drive CIN, however due to limited models and heterogenous results, it has remained controversial which aspects of aneuploidy can drive CIN. In this study we systematically tested the impact of different types of aneuploidies on the induction of CIN. We generated a plethora of isogenic aneuploid clones harboring whole chromosome or segmental aneuploidies in human p53-deficient RPE-1 cells. We observed increased segregation errors in cells harboring trisomies that strongly correlated to the number of gained genes. Strikingly, we found that clones harboring only monosomies do not induce a CIN phenotype. Finally, we found that an initial chromosome breakage event and subsequent fusion can instigate breakage-fusion-bridge cycles in segmental aneuploidies. This suggests that monosomies, trisomies and segmental aneuploidies have fundamentally different effects on chromosomal instability and these results help us to decipher the complex relationship between aneuploidy and CIN.

Background Single cell Strand-seq is a unique tool for the discovery and phasing of genomic inversions. Conventional methods to discover inversions with Strand-seq data are blind to known inversion locations, limiting their statistical power for the detection of inversions smaller than 10 Kb. Moreover, the methods rely on manual inspection to separate false and true positives. Results Here we describe “InvertypeR”, a method based on a Bayesian binomial model that genotypes inversions using fixed genomic coordinates. We validated InvertypeR by re-genotyping inversions reported for three trios by the Human Genome Structural Variation Consortium. Although 6.3% of the family inversion genotypes in the original study showed Mendelian discordance, this was reduced to 0.5% using InvertypeR. By applying InvertypeR to published inversion coordinates and predicted inversion hotspots (n = 3701), as well as coordinates from conventional inversion discovery, we furthermore genotyped 66 inversions not previously reported for the three trios. Conclusions InvertypeR discovers, genotypes, and phases inversions without relying on manual inspection. For greater accessibility, results are presented as phased chromosome ideograms with inversions linked to Strand-seq data in the genome browser. InvertypeR increases the power of Strand-seq for studies on the role of inversions in phenotypic variation, genome instability, and human disease.

Chromosome segregation errors cause aneuploidy and genomic heterogeneity, which are hallmarks of cancer in humans. A persistent high frequency of these errors (chromosomal instability (CIN)) is predicted to profoundly impact tumor evolution and therapy response. It is unknown, however, how prevalent CIN is in human tumors. Using three-dimensional live-cell imaging of patient-derived tumor organoids (tumor PDOs), we show that CIN is widespread in colorectal carcinomas regardless of background genetic alterations, including microsatellite instability. Cell-fate tracking showed that, although mitotic errors are frequently followed by cell death, some tumor PDOs are largely insensitive to mitotic errors. Single-cell karyotype sequencing confirmed heterogeneity of copy number alterations in tumor PDOs and showed that monoclonal lines evolved novel karyo-types over time in vitro. We conclude that ongoing CIN is common in colorectal cancer organoids, and propose that CIN levels and the tolerance for mitotic errors shape aneuploidy landscapes and karyotype heterogeneity.

Following DNA replication, equal amounts of histones are distributed over sister chromatids by re-deposition of parental histones and deposition of newly synthesized histones. Molecular mechanisms balancing the allocation of new and old histones remain largely unknown. Here, we studied the genome-wide distribution of new histones relative to parental DNA template strands and replication initiation zones using double-click-seq. In control conditions, new histones were preferentially found on DNA replicated by the lagging strand machinery. Strikingly, replication stress induced by hydroxyurea or curaxin treatment, and inhibition of ATR or p53 inactivation, inverted the observed histone deposition bias to the strand replicated by the leading strand polymerase in line with previously reported effects on RPA occupancy. We propose that asymmetric deposition of newly synthesized histones onto sister chromatids reflects differences in the processivity of leading and lagging strand synthesis.

Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality1, 2. Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes, by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.

Human genomes are typically assembled as consensus sequences that lack information on parental haplotypes. Here we describe a reference-free workflow for diploid de novo genome assembly that combines the chromosome-wide phasing and scaffolding capabilities of single-cell strand sequencing1,2 with continuous long-read or high-fidelity3 sequencing data. Employing this strategy, we produced a completely phased de novo genome assembly for each haplotype of an individual of Puerto Rican descent (HG00733) in the absence of parental data. The assemblies are accurate (quality value > 40) and highly contiguous (contig N50 > 23 Mbp) with low switch error rates (0.17%), providing fully phased single-nucleotide variants, indels and structural variants. A comparison of Oxford Nanopore Technologies and Pacific Biosciences phased assemblies identified 154 regions that are preferential sites of contig breaks, irrespective of sequencing technology or phasing algorithms., Assembly of haplotype-resolved human genomes is achieved by combining short and long reads.

Circulating tumor cells (CTCs) detected by CellSearch are prognostic in non-small-cell lung cancer (NSCLC), but rarely found. CTCs can be extracted from the blood together with mononuclear cell populations by diagnostic leukapheresis (DLA), therefore concentrating them. However, CellSearch can only process limited DLA volumes (&asymp, 2 mL). Therefore, we established a protocol to enumerate CTCs in DLA products with Isolation by SizE of Tumor cells (ISET), and compared CTC counts between CellSearch&reg, and ISET. DLA was performed in NSCLC patients who started a new therapy. With an adapted protocol, ISET could process 10 mL of DLA. CellSearch detected CTCs in a volume equaling 2 &times, 108 leukocytes (mean 2 mL). CTC counts per mL were compared. Furthermore, the live cell protocol of ISET was tested in eight patients. ISET successfully processed all DLA products&mdash, 16 with the fixed cell protocol and 8 with the live cell protocol. In total, 10&ndash, 20 mL of DLA was processed. ISET detected CTCs in 88% (14/16), compared to 69% (11/16, p &lt, 0.05) with CellSearch. ISET also detected higher number of CTCs (ISET median CTC/mL = 4, interquartile range [IQR] = 2&ndash, 6, CellSearch median CTC/mL = 0.9, IQR = 0&ndash, 1.8, p &lt, 0.01). Cells positive for the epithelial cell adhesion molecule (EpCAM+) per mL were detected in similar counts by both methods. Eight patients were processed with the live cell protocol. All had EpCAM+, CD45&minus, CD235- cells isolated by fluorescence-activated cell sorting (FACS). Overall, ISET processed larger volumes and detected higher CTC counts compared to CellSearch. EpCAM+ CTCs were detected in comparable rates.

The diploid nature of the human genome is neglected in many analyses done today, where a genome is perceived as a set of unphased variants with respect to a reference genome. This lack of haplotype-level analyses can be explained by a lack of methods that can produce dense and accurate chromosome-length haplotypes at reasonable costs. Here we introduce an integrative phasing strategy that combines global, but sparse haplotypes obtained from strand-specific single-cell sequencing (Strand-seq) with dense, yet local, haplotype information available through long-read or linked-read sequencing. We provide comprehensive guidance on the required sequencing depths and reliably assign more than 95% of alleles (NA12878) to their parental haplotypes using as few as 10 Strand-seq libraries in combination with 10-fold coverage PacBio data or, alternatively, 10X Genomics linked-read sequencing data. We conclude that the combination of Strand-seq with different technologies represents an attractive solution to chart the genetic variation of diploid genomes., Haplotype information is important in investigating many biological phenomena. Here, Porubsky et al. combine Strand-seq with long-read or linked-read sequencing to obtain complete and genome-wide haplotypes of a single individual genome at manageable costs.

Summary: The presence of an abnormal karyotype has been shown to be profoundly detrimental at the cellular and organismal levels but is an overt hallmark of cancer. Aneuploidy can lead to p53 activation and thereby prevents proliferation, but the exact trigger for p53 activation has remained controversial. Here, we have used a system to induce aneuploidy in untransformed human cells to explore how cells deal with different segregation errors. We show that p53 is activated only in a subset of the cells with altered chromosome content. Importantly, we find that at least a subset of whole-chromosome aneuploidies can be propagated in p53-proficient cells, indicating that aneuploidy does not always lead to activation of p53. Finally, we demonstrate that propagation of structural aneuploidies (gain or loss of part of a chromosome) induced by segregation errors is limited to p53-deficient cells. : Chromosome segregation errors can result in whole-chromosome aneuploidies or structural aneuploidies (involving only chromosome fragments). Here, Soto et al. show that whole-chromosome aneuploidies do not always lead to p53 activation and can therefore be propagated in a p53-proficient setting, whereas structural imbalances that result from segregation errors cannot. Keywords: aneuploidy, p53, chromosome missegregation, single-cell sequencing

Centrosome amplification is a common feature of human tumors, but whether this is a cause or a consequence of cancer remains unclear. Here, we test the consequence of centrosome amplification by creating mice in which centrosome number can be chronically increased in the absence of additional genetic defects. We show that increasing centrosome number elevated tumor initiation in a mouse model of intestinal neoplasia. Most importantly, we demonstrate that supernumerary centrosomes are sufficient to drive aneuploidy and the development of spontaneous tumors in multiple tissues. Tumors arising from centrosome amplification exhibit frequent mitotic errors and possess complex karyotypes, recapitulating a common feature of human cancer. Together, our data support a direct causal relationship among centrosome amplification, genomic instability, and tumor development.

Motivation Strand-seq is a specialized single-cell DNA sequencing technique centered around the directionality of single-stranded DNA. Computational tools for Strand-seq analyses must capture the strand-specific information embedded in these data. Results Here we introduce breakpointR, an R/Bioconductor package specifically tailored to process and interpret single-cell strand-specific sequencing data obtained from Strand-seq. We developed breakpointR to detect local changes in strand directionality of aligned Strand-seq data, to enable fine-mapping of sister chromatid exchanges, germline inversion and to support global haplotype assembly. Given the broad spectrum of Strand-seq applications we expect breakpointR to be an important addition to currently available tools and extend the accessibility of this novel sequencing technique. Availability and implementation R/Bioconductor package https://bioconductor.org/packages/breakpointR. Supplementary information Supplementary data are available at Bioinformatics online.

The prevailing genome assembly paradigm is to produce consensus sequences that “collapse” parental haplotypes into a consensus sequence. Here, we leverage the chromosome-wide phasing and scaffolding capabilities of single-cell strand sequencing (Strand-seq)1,2 and combine them with high-fidelity (HiFi) long sequencing reads3, in a novel reference-free workflow for diploid de novo genome assembly. Employing this strategy, we produce completely phased de novo genome assemblies separately for each haplotype of a single individual of Puerto Rican origin (HG00733) in the absence of parental data. The assemblies are accurate (QV > 40), highly contiguous (contig N50 > 25 Mbp) with low switch error rates (0.4%) providing fully phased single-nucleotide variants (SNVs), indels, and structural variants (SVs). A comparison of Oxford Nanopore and PacBio phased assemblies identifies 150 regions that are preferential sites of contig breaks irrespective of sequencing technology or phasing algorithms.

Guanine quadruplex (G4) structures are among the most stable secondary DNA structures that can form in vitro, and evidence for their existence in vivo has been steadily accumulating. Originally described mainly for their deleterious effects on genome stability, more recent research has focused on (potential) functions of G4 structures in telomere maintenance, gene expression, and other cellular processes. The combined research on G4 structures has revealed that properly regulating G4 DNA structures in cells is important to prevent genome instability and disruption of normal cell function. In this short review we provide some background and historical context of our work resulting in the identification of FANCJ, RTEL1 and BLM as helicases that act on G4 structures in vivo. Taken together these studies highlight important roles of different G4 DNA structures and specific G4 helicases at selected genomic locations and telomeres in regulating gene expression and maintaining genome stability.

Purpose: Tumor cells from patients with lung cancer are expelled from the primary tumor into the blood, but difficult to detect in the peripheral circulation. We studied the release of circulating tumor cells (CTCs) during surgery to test the hypothesis that CTC counts are influenced by hemodynamic changes (caused by surgical approach) and manipulation. Experimental Design: Patients undergoing video-assisted thoracic surgery (VATS) or open surgery for (suspected) primary lung cancer were included. Blood samples were taken before surgery (T0) from the radial artery (RA), from both the RA and pulmonary vein (PV) when the PV was located (T1) and when either the pulmonary artery (T2 open) or the PV (T2 VATS) was dissected. The CTCs were enumerated using the CellSearch system. Single-cell whole-genome sequencing was performed on isolated CTCs for aneuploidy. Results: CTCs were detected in 58 of 138 samples (42%) of 31 patients. CTCs were more often detected in the PV (70%) compared with the RA (22%, P < 0.01) and in higher counts (P < 0.01). After surgery, the RA but not the PV showed less often CTCs (P = 0.02). Type of surgery did not influence CTC release. Only six of 496 isolated CTCs showed aneuploidy, despite matched primary tumor tissue being aneuploid. Euploid so-called CTCs had a different morphology than aneuploid. Conclusions: CTCs defined by CellSearch were identified more often and in higher numbers in the PV compared with the RA, suggesting central clearance. The majority of cells in the PV were normal epithelial cells and outnumbered CTCs. Release of CTCs was not influenced by surgical approach.

The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6 kbp, CLR 191.5 kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective stand-alone technology for de novo assembly of human genomes.

The incomplete identification of structural variants (SVs) from whole-genome sequencing data limits studies of human genetic diversity and disease association. Here, we apply a suite of long-read, short-read, strand-specific sequencing technologies, optical mapping, and variant discovery algorithms to comprehensively analyze three trios to define the full spectrum of human genetic variation in a haplotype-resolved manner. We identify 818,054 indel variants (, Structural variants (SVs) in human genomes contribute diversity and diseases. Here, the authors use a multi-platform strategy to generate haplotype-resolved SVs for three human parent–child trios.

Introduction: Circulating tumor cells (CTC) can be used to monitor malignant disease longitudinally, but their use in non-small cell lung cancer (NSCLC) is limited due to low numbers in peripheral blood. Through Diagnostic leukapheresis (DLA) CTC can be obtained from larger blood volumes. We studied CTC in DLA product of NSCLC patients before and after treatment. Methods: One total blood volume was screened by DLA before and 1-3 months after treatment. Peripheral blood was drawn pre and post DLA for CTC enumeration by CellSearch. CTC were detected in DLA product directly (volume equaling 2×10^8 leukocytes) and after leukocyte depletion (RosetteSep, 9mL DLA product). Single cell whole genome sequencing was performed on isolated CTC. Results: Before treatment, CTC were more often detected in DLA (32/55, 58%) compared to blood (pre: 18/55, 33%, p

High-throughput next generation sequencing karyotyping has emerged as a powerful tool for the detection of genomic heterogeneity in normal tissues and cancers. Here we describe a single-cell whole genome sequencing (scWGS) platform to assess whole-chromosome aneuploidy, structural aneuploidies involving only chromosome fragments and more local small copy number alterations in individual cells. We provide a detailed protocol for the isolation, library preparation, low coverage sequencing and data analysis of single cells. Since our approach does not involve a whole-genome preamplification step, our method allows for acquisition of reliable high-resolution single-cell copy number profiles. Moreover, the protocol allows multiplexing of 384 single-cell libraries in one sequencing run, thereby significantly reducing sequencing costs and can be completed in 3-4 days starting from single cell isolation to analysis of sequencing data.

Haplotypes are fundamental to fully characterize the diploid genome of an individual, yet methods to directly chart the unique genetic makeup of each parental chromosome are lacking. Here we introduce single-cell DNA template strand sequencing (Strand-seq) as a novel approach to phasing diploid genomes along the entire length of all chromosomes. We demonstrate this by building a complete haplotype for a HapMap individual (NA12878) at high accuracy (concordance 99.3%), without using generational information or statistical inference. By use of this approach, we mapped all meiotic recombination events in a family trio with high resolution (median range ∼14 kb) and phased larger structural variants like deletions, indels, and balanced rearrangements like inversions. Lastly, the single-cell resolution of Strand-seq allowed us to observe loss of heterozygosity regions in a small number of cells, a significant advantage for studies of heterogeneous cell populations, such as cancer cells. We conclude that Strand-seq is a unique and powerful approach to completely phase individual genomes and map inheritance patterns in families, while preserving haplotype differences between single cells.

Sister chromatid exchanges (SCEs) are considered sensitive indicators of genome instability. Detection of SCEs typically requires cells to incorporate bromodeoxyuridine (BrdU) during two rounds of DNA synthesis. Previous studies have suggested that SCEs are induced by DNA replication over BrdU-substituted DNA and that BrdU incorporation alone could be responsible for the high number of SCE events observed in cells from patients with Bloom syndrome (BS), a rare genetic disorder characterized by marked genome instability and high SCE frequency. Here we show using Strand-seq, a single cell DNA template strand sequencing technique, that the presence of variable BrdU concentrations in the cell culture medium and in DNA template strands has no effect on SCE frequency in either normal or BS cells. We conclude that BrdU does not induce SCEs and that SCEs detected in either normal or BS cells reflect DNA repair events that occur spontaneously.

Summary Massively parallel sequencing is now widely used, but data interpretation is only as good as the reference assembly to which it is aligned. While the number of reference assemblies has rapidly expanded, most of these remain at intermediate stages of completion, either as scaffold builds, or as chromosome builds (consisting of correctly ordered, but not necessarily correctly oriented scaffolds separated by gaps). Completion of de novo assemblies remains difficult, as regions that are repetitive or hard to sequence prevent the accumulation of larger scaffolds, and create errors such as misorientations and mislocalizations. Thus, complementary methods for determining the orientation and positioning of fragments are important for finishing assemblies. Strand-seq is a method for determining template strand inheritance in single cells, information that can be used to determine relative genomic distance and orientation between scaffolds, and find errors within them. We present contiBAIT, an R/Bioconductor package which uses Strand-seq data to repair and improve existing assemblies. Availability and Implementation contiBAIT is available on Bioconductor. Source files available from GitHub. Supplementary information Supplementary data are available at Bioinformatics online.

High-throughput next generation sequencing karyotyping has emerged as a powerful tool for the detection of genomic heterogeneity in normal tissues and cancers. Here we describe a single-cell whole genome sequencing (scWGS) platform to assess whole-chromosome aneuploidy, structural aneuploidies involving only chromosome fragments and more local small copy number alterations in individual cells. We provide a detailed protocol for the isolation, library preparation, low coverage sequencing and data analysis of single cells. Since our approach does not involve a whole-genome preamplification step, our method allows for acquisition of reliable high-resolution single-cell copy number profiles. Moreover, the protocol allows multiplexing of 384 single-cell libraries in one sequencing run, thereby significantly reducing sequencing costs and can be completed in 3-4 days starting from single cell isolation to analysis of sequencing data.