Telomeres, aging, and cancer: the big picture.

The role of telomeres in human health and disease is yet to be fully understood. The limitations of mouse models for the study of human telomere biology and difficulties in accurately measuring the length of telomere repeats in chromosomes and cells have diverted attention from many important and relevant observations. The goal of this perspective is to summarize some of these observations and to discuss the antagonistic role of telomere loss in aging and cancer in the context of developmental biology, cell turnover, and evolution. It is proposed that both damage to DNA and replicative loss of telomeric DNA contribute to aging in humans, with the differences in leukocyte telomere length between humans being linked to the risk of developing specific diseases. These ideas are captured in the Telomere Erosion in Disposable Soma theory of aging proposed herein.

Mapping of sister chromatid exchange events and genome alterations in single cells.

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.

A survey of current methods to detect and genotype inversions.

Polymorphic inversions are ubiquitous in humans and they have been linked to both adaptation and disease. Following their discovery in Drosophila more than a century ago, inversions have proved to be more elusive than other structural variants. A wide variety of methods for the detection and genotyping of inversions have recently been developed: multiple techniques based on selective amplification by PCR, short- and long-read sequencing approaches, principal component analysis of small variant haplotypes, template strand sequencing, optical mapping, and various genome assembly methods. Many methods apply complex wet lab protocols or increasingly refined bioinformatic analyses. This review is an attempt to provide a practical summary and comparison of the methods that are in current use, with a focus on metrics such as the maximum size of segmental duplications at inversion breakpoints that each method can tolerate, the size range of inversions that they recover, their throughput, and whether the locations of putative inversions must be known beforehand.

InvertypeR: Bayesian inversion genotyping with Strand-seq data.

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.

breakpointR: an R/Bioconductor package to localize strand state changes in Strand-seq data.

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.

Construction of Whole Genomes from Scaffolds Using Single Cell Strand-Seq Data.

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 directionality 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.

Improved assembly and variant detection of a haploid human genome using single-molecule, high-fidelity long reads.

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 standalone technology for de novo assembly of human genomes.

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.

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.

Guanine quadruplex monoclonal antibody 1H6 cross-reacts with restrained thymidine-rich single stranded DNA.

Previously we reported the production and characterization of monoclonal antibody 1H6 raised against (T4G4)2 intermolecular guanine quadruplex (G4) DNA structures (Henderson A. et al. (2014) Nucleic Acids Res., 42, 860-869; Hoffmann R.F. et al. (2016) Nucleic Acids Res., 44, 152-163). It was shown that 1H6 strongly stains nuclei and has exquisite specificity for heterochromatin by immuno-electron microscopy. Here we extend our studies of 1H6 reactivity using enzyme-linked immunosorbent assay (ELISA) and microscale thermophoresis (MST). As previously reported, 1H6 was found to strongly bind intermolecular G4 structures with a (T4G4)2 sequence motif. However, using both methods we did not detect significant binding to G4 structures without thymidines in their sequence motif or to G4 structures made with (T2G4)2 oligonucleotides. In addition, we observed strong, sequence-specific binding of 1H6 by ELISA to immobilized single stranded poly(T) DNA but not to immobilized poly(C) or poly(A) homo-polymers. Cross-reactivity of 1H6 to poly(T) was not measured in solution using MST. 1H6 was furthermore found to bind to selected areas on DNA fibers but only after DNA denaturation. Based on these observations we propose that 1H6 binds with high affinity to adjacent T's that are restricted in their movement in selected G4 structures and denatured DNA. Cross-reactivity of 1H6 to immobilized single stranded T-rich DNA next to its previously reported specificity for bona fide G4 structures needs to be taken into account in the interpretation of 1H6 binding to (sub-) cellular structures.

Assembling draft genomes using contiBAIT.

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. CONTACT: koneill@bcgsc.ca or mark.hills@stemcell.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

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.

Identification of sister chromatids by DNA template strand sequences.

It is generally assumed that sister chromatids are genetically and functionally identical and that segregation to daughter cells is a random process. However, functional differences between sister chromatids regulate daughter cell fate in yeast and sister chromatid segregation is not random in Escherichia coli. Differentiated sister chromatids, coupled with non-random segregation, have been proposed to regulate cell fate during the development of multicellular organisms. This hypothesis has not been tested because molecular features to reliably distinguish between sister chromatids are not obvious. Here we show that parental 'Watson' and 'Crick' DNA template strands can be identified in sister chromatids of murine metaphase chromosomes using CO-FISH (chromosome orientation fluorescence in situ hybridization) with unidirectional probes specific for centromeric and telomeric repeats. All chromosomes were found to have a uniform orientation with the 5' end of the short arm on the same strand as T-rich major satellite repeats. The invariable orientation of repetitive DNA was used to differentially label sister chromatids and directly study mitotic segregation patterns in different cell types. Whereas sister chromatids appeared to be randomly distributed between daughter cells in cultured lung fibroblasts and embryonic stem cells, significant non-random sister chromatid segregation was observed in a subset of colon crypt epithelial cells, including cells outside positions reported for colon stem cells. Our results establish that DNA template sequences can be used to distinguish sister chromatids and follow their mitotic segregation in vivo.

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.

Strand-seq: a unifying tool for studies of chromosome segregation.

Non random segregation of sister chromatids has been implicated to help specify daughter cell fate (the Silent Sister Hypothesis [1]) or to protect the genome of long-lived stem cells (the Immortal Strand Hypothesis [2]). The idea that sister chromatids are non-randomly segregated into specific daughter cells is only marginally supported by data in sporadic and often contradictory studies. As a result, the field has moved forward rather slowly. The advent of being able to directly label and differentiate sister chromatids in vivo using fluorescence in situ hybridization [3] was a significant advance for such studies. However, this approach is limited by the need for large tracks of unidirectional repeats on chromosomes and the reliance on quantitative imaging of fluorescent probes and rigorous statistical analysis to discern between the two competing hypotheses. A novel method called Strand-seq which uses next-generation sequencing to assay sister chromatid inheritance patterns independently for each chromosome [4] offers a comprehensive approach to test for non-random segregation. In addition Strand-seq enables studies on the deposition of chromatin marks in relation to DNA replication. This method is expected to help unify the field by testing previous claims of non-random segregation in an unbiased way in many model systems in vitro and in vivo.

DNA template strand sequencing of single-cells maps genomic rearrangements at high resolution.

DNA rearrangements such as sister chromatid exchanges (SCEs) are sensitive indicators of genomic stress and instability, but they are typically masked by single-cell sequencing techniques. We developed Strand-seq to independently sequence parental DNA template strands from single cells, making it possible to map SCEs at orders-of-magnitude greater resolution than was previously possible. On average, murine embryonic stem (mES) cells exhibit eight SCEs, which are detected at a resolution of up to 23 bp. Strikingly, Strand-seq of 62 single mES cells predicts that the mm 9 mouse reference genome assembly contains at least 17 incorrectly oriented segments totaling nearly 1% of the genome. These misoriented contigs and fragments have persisted through several iterations of the mouse reference genome and have been difficult to detect using conventional sequencing techniques. The ability to map SCE events at high resolution and fine-tune reference genomes by Strand-seq dramatically expands the scope of single-cell sequencing.

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.