Single-cell template strand sequencing by Strand-seq enables the characterization of individual homologs

The ability to distinguish between genome sequences of homologous chromosomes in single cells is important for studies of copy-neutral genomic rearrangements (such as inversions and translocations), building chromosome-length haplotypes, refining genome assemblies, mapping sister chromatid exchange events and exploring cellular heterogeneity. Strand-seq is a single-cell sequencing technology that resolves the individual homologs within a cell by restricting sequence analysis to the DNA template strands used during DNA replication. This protocol, which takes up to 4 d to complete, relies on the directionality of DNA, in which each single strand of a DNA molecule is distinguished based on its 5'-3' orientation. Culturing cells in a thymidine analog for one round of cell division labels nascent DNA strands, allowing for their selective removal during genomic library construction. To preserve directionality of template strands, genomic preamplification is bypassed and labeled nascent strands are nicked and not amplified during library preparation. Each single-cell library is multiplexed for pooling and sequencing, and the resulting sequence data are aligned, mapping to either the minus or plus strand of the reference genome, to assign template strand states for each chromosome in the cell. The major adaptations to conventional single-cell sequencing protocols include harvesting of daughter cells after a single round of BrdU incorporation, bypassing of whole-genome amplification, and removal of the BrdU⁺ strand during Strand-seq library preparation. By sequencing just template strands, the structure and identity of each homolog are preserved.

Deletion of the MAD2L1 spindle assembly checkpoint gene is tolerated in mouse models of acute T-cell lymphoma and hepatocellular carcinoma

Chromosome instability (CIN) is deleterious to normal cells because of the burden of aneuploidy. However, most human solid tumors have an abnormal karyotype implying that gain and loss of chromosomes by cancer cells confers a selective advantage. CIN can be induced in the mouse by inactivating the spindle assembly checkpoint. This is lethal in the germline but we show here that adult T cells and hepatocytes can survive conditional inactivation of the Mad2l1 SAC gene and resulting CIN. This causes rapid onset of acute lymphoblastic leukemia (T-ALL) and progressive development of hepatocellular carcinoma (HCC), both lethal diseases. The resulting DNA copy number variation and patterns of chromosome loss and gain are tumor-type specific, suggesting differential selective pressures on the two tumor cell types.

Centrosome Amplification Is Sufficient to Promote Spontaneous Tumorigenesis in Mammals

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.

Defects in lymphocyte telomere homeostasis contribute to cellular immune phenotype in patients with cartilage-hair hypoplasia

BACKGROUND

Mutations in the long noncoding RNA RNase component of the mitochondrial RNA processing endoribonuclease (RMRP) give rise to the autosomal recessive condition cartilage-hair hypoplasia (CHH). The CHH disease phenotype has some overlap with dyskeratosis congenita, a well-known "telomere disorder." RMRP binds the telomerase reverse transcriptase (catalytic subunit) in some cell lines, raising the possibility that RMRP might play a role in telomere biology.

OBJECTIVE

We sought to determine whether a telomere phenotype is present in immune cells from patients with CHH and explore mechanisms underlying these observations.

METHODS

We assessed proliferative capacity and telomere length using flow-fluorescence in situ hybridization (in situ hybridization and flow cytometry) of primary lymphocytes from patients with CHH, carrier relatives, and control subjects. The role of telomerase holoenzyme components in gene expression and activity were assessed by using quantitative PCR and the telomere repeat amplification protocol from PBMCs and enriched lymphocyte cultures.

RESULTS

Lymphocyte cultures from patients with CHH display growth defects in vitro, which is consistent with an immune deficiency cellular phenotype. Here we show that telomere length and telomerase activity are impaired in primary lymphocyte subsets from patients with CHH. Notably, telomerase activity is affected in a gene dose-dependent manner when comparing heterozygote RMRP carriers with patients with CHH. Telomerase deficiency in patients with CHH is not mediated by abnormal telomerase gene transcript levels relative to those of endogenous genes.

CONCLUSION

These findings suggest that telomere deficiency is implicated in the CHH disease phenotype through an as yet unidentified mechanism.

Direct chromosome-length haplotyping by single-cell sequencing

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.

Characterizing polymorphic inversions in human genomes by single-cell sequencing

Identifying genomic features that differ between individuals and cells can help uncover the functional variants that drive phenotypes and disease susceptibilities. For this, single-cell studies are paramount, as it becomes increasingly clear that the contribution of rare but functional cellular subpopulations is important for disease prognosis, management, and progression. Until now, studying these associations has been challenged by our inability to map structural rearrangements accurately and comprehensively. To overcome this, we coupled single-cell sequencing of DNA template strands (Strand-seq) with custom analysis software to rapidly discover, map, and genotype genomic rearrangements at high resolution. This allowed us to explore the distribution and frequency of inversions in a heterogeneous cell population, identify several polymorphic domains in complex regions of the genome, and locate rare alleles in the reference assembly. We then mapped the entire genomic complement of inversions within two unrelated individuals to characterize their distinct inversion profiles and built a nonredundant global reference of structural rearrangements in the human genome. The work described here provides a powerful new framework to study structural variation and genomic heterogeneity in single-cell samples, whether from individuals for population studies or tissue types for biomarker discovery.

Single-cell sequencing reveals karyotype heterogeneity in murine and human malignancies

BACKGROUND

Chromosome instability leads to aneuploidy, a state in which cells have abnormal numbers of chromosomes, and is found in two out of three cancers. In a chromosomal instable p53 deficient mouse model with accelerated lymphomagenesis, we previously observed whole chromosome copy number changes affecting all lymphoma cells. This suggests that chromosome instability is somehow suppressed in the aneuploid lymphomas or that selection for frequently lost/gained chromosomes out-competes the CIN-imposed mis-segregation.

RESULTS

To distinguish between these explanations and to examine karyotype dynamics in chromosome instable lymphoma, we use a newly developed single-cell whole genome sequencing (scWGS) platform that provides a complete and unbiased overview of copy number variations (CNV) in individual cells. To analyse these scWGS data, we develop AneuFinder, which allows annotation of copy number changes in a fully automated fashion and quantification of CNV heterogeneity between cells. Single-cell sequencing and AneuFinder analysis reveals high levels of copy number heterogeneity in chromosome instability-driven murine T-cell lymphoma samples, indicating ongoing chromosome instability. Application of this technology to human B cell leukaemias reveals different levels of karyotype heterogeneity in these cancers.

CONCLUSION

Our data show that even though aneuploid tumours select for particular and recurring chromosome combinations, single-cell analysis using AneuFinder reveals copy number heterogeneity. This suggests ongoing chromosome instability that other platforms fail to detect. As chromosome instability might drive tumour evolution, karyotype analysis using single-cell sequencing technology could become an essential tool for cancer treatment stratification.

Bromodeoxyuridine does not contribute to sister chromatid exchange events in normal or Bloom syndrome 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.

Guanine quadruplex structures localize to heterochromatin

Increasing amounts of data support a role for guanine quadruplex (G4) DNA and RNA structures in various cellular processes. We stained different organisms with monoclonal antibody 1H6 specific for G4 DNA. Strikingly, immuno-electron microscopy showed exquisite specificity for heterochromatin. Polytene chromosomes from Drosophila salivary glands showed bands that co-localized with heterochromatin proteins HP1 and the SNF2 domain-containing protein SUUR. Staining was retained in SUUR knock-out mutants but lost upon overexpression of SUUR. Somatic cells in Macrostomum lignano were strongly labeled, but pluripotent stem cells labeled weakly. Similarly, germline stem cells in Drosophila ovaries were weakly labeled compared to most other cells. The unexpected presence of G4 structures in heterochromatin and the difference in G4 staining between somatic cells and stem cells with germline DNA in ciliates, flatworms, flies and mammals point to a conserved role for G4 structures in nuclear organization and cellular differentiation.

The developmental potential of iPSCs is greatly influenced by reprogramming factor selection

Induced pluripotent stem cells (iPSCs) are commonly generated by transduction of Oct4, Sox2, Klf4, and Myc (OSKM) into cells. Although iPSCs are pluripotent, they frequently exhibit high variation in terms of quality, as measured in mice by chimera contribution and tetraploid complementation. Reliably high-quality iPSCs will be needed for future therapeutic applications. Here, we show that one major determinant of iPSC quality is the combination of reprogramming factors used. Based on tetraploid complementation, we found that ectopic expression of Sall4, Nanog, Esrrb, and Lin28 (SNEL) in mouse embryonic fibroblasts (MEFs) generated high-quality iPSCs more efficiently than other combinations of factors including OSKM. Although differentially methylated regions, transcript number of master regulators, establishment of specific superenhancers, and global aneuploidy were comparable between high- and low-quality lines, aberrant gene expression, trisomy of chromosome 8, and abnormal H2A.X deposition were distinguishing features that could potentially also be applicable to human.

How to count chromosomes in a cell: An overview of current and novel technologies

Aneuploidy, an aberrant number of chromosomes in a cell, is a feature of several syndromes associated with cognitive and developmental defects. In addition, aneuploidy is considered a hallmark of cancer cells and has been suggested to play a role in neurodegenerative disease. To better understand the relationship between aneuploidy and disease, various methods to measure the chromosome numbers in cells have been developed, each with their own advantages and limitations. While some methods rely on dividing cells and thus bias aneuploidy rates to that population, other, more unbiased methods can only detect the average aneuploidy rates in a cell population, cloaking cell-to-cell heterogeneity. Furthermore, some techniques are more prone to technical artefacts, which can result in over- or underestimation of aneuploidy rates. In this review, we provide an overview of several "traditional" karyotyping methods as well as the latest high throughput next generation sequencing karyotyping protocols with their respective advantages and disadvantages.

Resolution of telomere associations by TRF1 cleavage in mouse embryonic stem cells

Telomere associations have been observed during key cellular processes such as mitosis, meiosis, and carcinogenesis and must be resolved before cell division to prevent genome instability. Here we establish that telomeric repeat-binding factor 1 (TRF1), a core component of the telomere protein complex, is a mediator of telomere associations in mammalian cells. Using live-cell imaging, we show that expression of TRF1 or yellow fluorescent protein (YFP)-TRF1 fusion protein above endogenous levels prevents proper telomere resolution during mitosis. TRF1 overexpression results in telomere anaphase bridges and aggregates containing TRF1 protein and telomeric DNA. Site-specific protein cleavage of YFP-TRF1 by tobacco etch virus protease resolves telomere aggregates, indicating that telomere associations are mediated by TRF1. This study provides novel insight into the formation and resolution of telomere associations.

Detection of G-quadruplex DNA in mammalian cells

It has been proposed that guanine-rich DNA forms four-stranded structures in vivo called G-quadruplexes or G4 DNA. G4 DNA has been implicated in several biological processes, but tools to study G4 DNA structures in cells are limited. Here we report the development of novel murine monoclonal antibodies specific for different G4 DNA structures. We show that one of these antibodies designated 1H6 exhibits strong nuclear staining in most human and murine cells. Staining intensity increased on treatment of cells with agents that stabilize G4 DNA and, strikingly, cells deficient in FANCJ, a G4 DNA-specific helicase, showed stronger nuclear staining than controls. Our data strongly support the existence of G4 DNA structures in mammalian cells and indicate that the abundance of such structures is increased in the absence of FANCJ. We conclude that monoclonal antibody 1H6 is a valuable tool for further studies on the role of G4 DNA in cell and molecular biology.

Cord blood T cells retain early differentiation phenotype suitable for immunotherapy after TCR gene transfer to confer EBV specificity

Adoptive T cell therapy can be effective for Epstein-Barr virus (EBV)-associated posttransplant lymphoproliferative disease and melanoma. Transducing high-affinity TCR genes into T lymphocytes is an emerging method to improve potency and specificity of tumor-specific T cells. However, both methods necessitate in vitro lymphocyte proliferation, generating highly differentiated effector cells that display reduced survival and antitumor efficacy postinfusion. TCR-transduction of naive lymphocytes isolated from peripheral blood is reported to provide superior in vivo survival and function. We utilized cord blood (CB) lymphocytes, which comprise mainly naive cells, for transducing EBV-specific TCR. Comparable TCR expression was achieved in adult and CB cells, but the latter expressed an earlier differentiation profile. Further antigen-driven stimulation skewed adult lymphocytes to a late differentiation phenotype associated with immune exhaustion. In contrast, CB T cells retained a less differentiated phenotype after antigen stimulation, remaining CD57-negative but were still capable of antigen-specific polyfunctional cytokine expression and cytotoxicity in response to EBV antigen. CB T cells also retained longer telomeres and in general possessed higher telomerase activity indicative of greater proliferative potential. CB lymphocytes therefore have qualities indicating prolonged survival and effector function favorable to immunotherapy, especially in settings where donor lymphocytes are unavailable such as in solid organ and CB transplantation.

The mammalian proteins MMS19, MIP18, and ANT2 are involved in cytoplasmic iron-sulfur cluster protein assembly

Iron-sulfur (Fe-S) clusters are essential cofactors of proteins with a wide range of biological functions. A dedicated cytosolic Fe-S cluster assembly (CIA) system is required to assemble Fe-S clusters into cytosolic and nuclear proteins. Here, we show that the mammalian nucleotide excision repair protein homolog MMS19 can simultaneously bind probable cytosolic iron-sulfur protein assembly protein CIAO1 and Fe-S proteins, confirming that MMS19 is a central protein of the CIA machinery that brings Fe-S cluster donor proteins and the receiving apoproteins into proximity. In addition, we show that mitotic spindle-associated MMXD complex subunit MIP18 also interacts with both CIAO1 and Fe-S proteins. Specifically, it binds the Fe-S cluster coordinating regions in Fe-S proteins. Furthermore, we show that ADP/ATP translocase 2 (ANT2) interacts with Fe-S apoproteins and MMS19 in the CIA complex but not with the individual proteins. Together, these results elucidate the composition and interactions within the late CIA complex.

Epigenetic differences between sister chromatids?

Semi-conservative replication ensures that the DNA sequence of sister chromatids is identical except for replication errors and variation in the length of telomere repeats resulting from replicative losses and variable end processing. What happens with the various epigenetic marks during DNA replication is less clear. Many chromatin marks are likely to be copied onto both sister chromatids in conjunction with DNA replication, whereas others could be distributed randomly between sister chromatids. Epigenetic differences between sister chromatids could also emerge in a more predictable manner, for example, following processes that are associated with lagging strand DNA replication. The resulting epigenetic differences between sister chromatids could result in different gene expression patterns in daughter cells. This possibility has been difficult to test because techniques to distinguish between parental sister chromatids require analysis of single cells and are not obvious. Here, we briefly review the topic of sister chromatid epigenetics and discuss how the identification of sister chromatids in cells could change the way we think about asymmetric cell divisions and stochastic variation in gene expression between cells in general and paired daughter cells in particular.

Collapse of telomere homeostasis in hematopoietic cells caused by heterozygous mutations in telomerase genes

Telomerase activity is readily detectable in extracts from human hematopoietic stem and progenitor cells, but appears unable to maintain telomere length with proliferation in vitro and with age in vivo. We performed a detailed study of the telomere length by flow FISH analysis in leukocytes from 835 healthy individuals and 60 individuals with reduced telomerase activity. Healthy individuals showed a broad range in average telomere length in granulocytes and lymphocytes at any given age. The average telomere length declined with age at a rate that differed between age-specific breakpoints and between cell types. Gender differences between leukocyte telomere lengths were observed for all cell subsets studied; interestingly, this trend could already be detected at birth. Heterozygous carriers for mutations in either the telomerase reverse transcriptase (hTERT) or the telomerase RNA template (hTERC) gene displayed striking and comparable telomere length deficits. Further, non-carrier relatives of such heterozygous individuals had somewhat shorter leukocyte telomere lengths than expected; this difference was most profound for granulocytes. Failure to maintain telomere homeostasis as a result of partial telomerase deficiency is thought to trigger cell senescence or cell death, eventually causing tissue failure syndromes. Our data are consistent with these statements and suggest that the likelihood of similar processes occurring in normal individuals increases with age. Our work highlights the essential role of telomerase in the hematopoietic system and supports the notion that telomerase levels in hematopoietic cells, while limiting and unable to prevent overall telomere shortening, are nevertheless crucial to maintain telomere homeostasis with age.

Human blood cells all originate from a common precursor, the hematopoietic stem cell. Telomerase, the enzyme responsible for adding telomere repeats to chromosome ends, is active in human hematopoietic stem cells but appears unable to maintain a constant telomere length with age. We first document the telomere length of different blood cell subsets from 835 healthy individuals between birth and 100 years, to delineate the normal rate of telomere attrition with age. Telomere lengths of blood cells were found to be slightly longer in women than in men, from birth and throughout life. We then compared this reference data to the telomere length in similar blood cell subsets from individuals with reduced telomerase activity as a result of a mutation in one of the genes encoding telomerase and from their direct relatives. Strikingly short telomeres were found in telomerase-deficient individuals, consistent with their cellular pathology and disease susceptibility, and somewhat shorter telomeres than expected were found in cells of relatives with normal telomerase maintenance. Our data can be used as a reference for blood cell telomere length in studies of normal and accelerated aging.

RTEL1 contributes to DNA replication and repair and telomere maintenance

Telomere maintenance and DNA repair are important processes that protect the genome. The essential helicase mRtel1 functions in homologous recombination repair and replication. In addition, telomeres in mRtel-deficient ES cells appear relatively stable in length, suggesting that mRtel1 is required to allow extension by telomerase.

Telomere maintenance and DNA repair are important processes that protect the genome against instability. mRtel1, an essential helicase, is a dominant factor setting telomere length in mice. In addition, mRtel1 is involved in DNA double-strand break repair. The role of mRtel1 in telomere maintenance and genome stability is poorly understood. Therefore we used mRtel1-deficient mouse embryonic stem cells to examine the function of mRtel1 in replication, DNA repair, recombination, and telomere maintenance. mRtel1-deficient mouse embryonic stem cells showed sensitivity to a range of DNA-damaging agents, highlighting its role in replication and genome maintenance. Deletion of mRtel1 increased the frequency of sister chromatid exchange events and suppressed gene replacement, demonstrating the involvement of the protein in homologous recombination. mRtel1 localized transiently at telomeres and is needed for efficient telomere replication. Of interest, in the absence of mRtel1, telomeres in embryonic stem cells appeared relatively stable in length, suggesting that mRtel1 is required to allow extension by telomerase. We propose that mRtel1 is a key protein for DNA replication, recombination, and repair and efficient elongation of telomeres by telomerase.

Telomere length is associated with disease severity and declines with age in dyskeratosis congenita

BACKGROUND

Dyskeratosis congenita is a cancer-prone bone marrow failure syndrome caused by aberrations in telomere biology.

DESIGN AND METHODS

We studied 65 patients with dyskeratosis congenita and 127 unaffected relatives. Telomere length was measured by automated multicolor flow fluorescence in situ hybridization in peripheral blood leukocyte subsets. We age-adjusted telomere length using Z-scores (standard deviations from the mean for age).

RESULTS

We confirmed that telomere lengths below the first percentile for age are very sensitive and specific for the diagnosis of dyskeratosis congenita. We provide evidence that lymphocytes alone and not granulocytes may suffice for clinical screening, while lymphocyte subsets may be required for challenging cases, including identification of silent carriers. We show for the first time using flow fluorescence in situ hybridization that the shortest telomeres are associated with severe variants (Hoyeraal-Hreidarsson and Revesz syndromes), mutations in DKC1, TINF2, or unknown genes, and moderate or severe aplastic anemia. In the first longitudinal follow up of dyskeratosis congenita patients, we demonstrate that telomere lengths decline with age, in contrast to the apparent stable telomere length observed in cross-sectional data.

CONCLUSIONS

Telomere length by flow fluorescence in situ hybridization is an important diagnostic test for dyskeratosis congenita; age-adjusted values provide a quantitative measure of disease severity (clinical subset, mutated gene, and degree of bone marrow failure). Patients with dyskeratosis congenita have accelerated telomere shortening. This study is registered at www.clinicaltrials.gov (identifier: NCT00027274).

Analysis of repetitive DNA in chromosomes by flow cytometry

We developed a flow cytometry method, chromosome flow fluorescence in situ hybridization (FISH), called CFF, to analyze repetitive DNA in chromosomes using FISH with directly labeled peptide nucleic acid (PNA) probes. We used CFF to measure the abundance of interstitial telomeric sequences in Chinese hamster chromosomes and major satellite sequences in mouse chromosomes. Using CFF we also identified parental homologs of human chromosome 18 with different amounts of repetitive DNA.

Ancestral mutation in telomerase causes defects in repeat addition processivity and manifests as familial pulmonary fibrosis

The telomerase reverse transcriptase synthesizes new telomeres onto chromosome ends by copying from a short template within its integral RNA component. During telomere synthesis, telomerase adds multiple short DNA repeats successively, a property known as repeat addition processivity. However, the consequences of defects in processivity on telomere length maintenance are not fully known. Germline mutations in telomerase cause haploinsufficiency in syndromes of telomere shortening, which most commonly manifest in the age-related disease idiopathic pulmonary fibrosis. We identified two pulmonary fibrosis families that share two non-synonymous substitutions in the catalytic domain of the telomerase reverse transcriptase gene hTERT: V791I and V867M. The two variants fell on the same hTERT allele and were associated with telomere shortening. Genealogy suggested that the pedigrees shared a single ancestor from the nineteenth century, and genetic studies confirmed the two families had a common founder. Functional studies indicated that, although the double mutant did not dramatically affect first repeat addition, hTERT V791I-V867M showed severe defects in telomere repeat addition processivity in vitro. Our data identify an ancestral mutation in telomerase with a novel loss-of-function mechanism. They indicate that telomere repeat addition processivity is a critical determinant of telomere length and telomere-mediated disease.