Nanopore Guided Regional Assembly

The telomeres are the “caps” of the chromosomes and their vital role is to protect them. Possible telomere dysfunction caused by telomere rearrangements can be fatal for the cell and result in age-related diseases, including cancer. The telomeres and subtelomeres are regions that are hard to investigate. The current technology cannot provide their complete sequence, instead the DNA is given in multiple pieces. Current methods of assembling the pieces of these regions are not accurate enough due to the region’s high variability and complex repeated patterns. We propose a hybrid assembly method, the NPGREAT, which utilizes two of the latest available data: Linked-Reads and ultralong Nanopore reads. It consists of five main steps: (i) The input selection of the data, (ii) the Orientation, Order and Enhanced Correction of the short contigs by using the long reads as scaffolds, upon which the short contigs are mapped to. Particularly, the Enhanced Correction step allows for the correction of potential misassemblies within the short contigs due to deletions in tandem repeat regions. The nanopore sequence is used to fill the missing portion, representing the tandem repeat region accurately, a region which is highly variable from one human to another. Next, in the (iii) Region Extraction step, the segments of the multiple long reads that can be used to connect the short contigs, are extracted. Then, in the (iv) Gap Filling step, all possible segments are taken into account and one is selected to fill each gap. Finally, in the (v) Combination step, the corrected short pieces are combined with the connector segments. The output is the subtelomere region of the chromosome. NPGREAT is evaluated with the use of the QUAST tool and the resulting assemblies are of high quality.https://digitalcommons.odu.edu/gradposters2021_sciences/1003/thumbnail.jp

NPGreat: Assembly of the Human Subtelomere Regions with the Use of Ultralong Nanopore Reads and Linked Reads

Background: Human subtelomeric DNA regulates the length and stability of adjacent telomeres that are critical for cellular function, and contains many gene/pseudogene families. Large evolutionarily recent segmental duplications and associated structural variation in human subtelomeres has made complete sequencing and assembly of these regions difficult to impossible for many loci, complicating or precluding a wide range of genetic analyses to investigate their function. Results: We present a hybrid assembly method, NanoPore Guided REgional Assembly Tool (NPGREAT), which combines Linked-Read data with ultralong nanopore reads spanning subtelomeric segmental duplications to potentially overcome these difficulties. Linked-Read sets identified by matches with 1-copy subtelomere sequence adjacent to segmental duplications are assembled and extended into the segmental duplication regions using Regional Extension of Assemblies using Linked-Reads (REXTAL). Telomere-containing ultralong nanopore reads are then used to provide contiguity and correct orientation for matching REXTAL sequence contigs as well as identification/correction of any misassemblies (associated primarily with tandem repeats). While we focus on subtelomeres, the method is generally applicable to assembly of segmental duplications and other complex genome regions. Our method was tested for a subset of representative subtelomeres with ultralong nanopore read coverage in GM12878. 10X Linked-Read datasets with high depth of coverage and a TELL-seq Linked-Read dataset with lower depth of coverage were each combined with the ultralong nanopore reads from the same genome to provide improved assemblies. Tandem repeat regions of the short-read assemblies, which are especially prone to misassembly due to collapse of matching tandemly repeated reads, were readily identified and properly sized by comparison with the nanopore reads. Conclusion: The NPGREAT method resulted in extension of high-quality assemblies into otherwise inaccessible segmental duplication regions near telomeres, enhancing our ability to accurately assemble human subtelomere DNA. This information will enable improved analyses of the structure, function, and evolution of these key regions

Metaphase and Interphase Cytogenetics with Alu-PCR-amplified Yeast Artificial Chromosome Clones Containing the BCR Gene and the Protooncogenes c-raf-1, c-fms, and c-erbB-21

A human yeast artificial chromosome (YAC) library was screened by polymerase chain reaction with oligonucleotide primers defined for DNA sequences of the BCR gene and the protooncogenes c-raf-1, c-fms, and c-erB-2. Alu-PCR-generated human DNA sequences were obtained from the respective YAC clones and used for fluorescence in situ hybridization experiments under suppression conditions. After chromosomal in situ suppression hybridization to GTG-banded human prometaphase chromosomes, seven of nine initially isolated YAC clones yielded strong signals exclusively in the chromosome bands containing the respective genes. Two clones yielded additional signals on other chromosomes and were excluded from further tests. The band-specific YACs were successfully applied to visualize specific structural chromosome aberrations in peripheral blood cells from patients with myelodysplasia exhibiting del(5)(q13q34), chronic myeloid leukemia and acute lymphocytic leukemia with t(9;22)(q34;q11), acute promyelocytic leukemia (M3) with t(15;17)(q22;q21), and in a cell line established from a proband with the constitutional translocation t(3;8)(p14.2;q24). In addition to the analysis of metaphase spreads, we demonstrate the particular usefulness of these YAC clones in combination with whole chromosome painting to analyze specific chromosome aberrations directly in the interphase nucleus

Chromosomal bar codes produced by multicolor fluorescence in situ hybridization with multiple YAC clones and whole chromosome painting probes

Colored chromosome staining patterns, termed chromosomal ‘bar codes’ (CBCs), were obtained on human chromosomes by fluorescence in situ hybridization (FISH) with pools of Alu-PCR products from YAC dones containing human DNA inserts ranging from 100 kbp to 1 Mbp. In contrast to conventional G- or R-bands, the chromosomal position, extent, Individual color and relative signal intensity of each ‘bar’ could be modified depending on probe selection and labeling procedures. Alu-PCR amplification products were generated from 31 YAC clones which mapped to 37 different chromosome bands. For multiple color FISH, Alu-PCR amplification products from various clones were either biotinylated or labeled with digoxigenin. Probes from up to twenty YAC clones were used simultaneously to produce CBCs on selected human chromosomes. Evaluation using a cooled CCD camera and digital image analysis confirmed the high reproducibility of the bars from one metaphase spread to another. Combinatorial FISH with mixtures of whole chromosome paint probes was applied to paint seven chromosomes simultaneously in different colors along with a set of YAC clones which map to these chromosomes. We discuss the potential to construct analytical chromosomal bar codes adapted to particular needs of cytogenetic investigations and automated image analysis

Comprehensive Analysis of Human Subtelomeres by Whole Genome Mapping

Detailed comprehensive knowledge of the structures of individual long-range telomere-terminal haplotypes are needed to understand their impact on telomere function, and to delineate the population structure and evolution of subtelomere regions. However, the abundance of large evolutionarily recent segmental duplications and high levels of large structural variations have complicated both the mapping and sequence characterization of human subtelomere regions. Here, we use high throughput optical mapping of large single DNA molecules in nanochannel arrays for 154 human genomes from 26 populations to present a comprehensive look at human subtelomere structure and variation. The results catalog many novel long-range subtelomere haplotypes and determine the frequencies and contexts of specific subtelomeric duplicons on each chromosome arm, helping to clarify the currently ambiguous nature of many specific subtelomere structures as represented in the current reference sequence (HG38). The organization and content of some duplicons in subtelomeres appear to show both chromosome arm and population-specific trends. Based upon these trends we estimate a timeline for the spread of these duplication blocks

Analysis of Subtelomeric REXTAL Assemblies Using QUAST

Genomic regions of high segmental duplication content and/or structural variation have led to gaps and misassemblies in the human reference sequence, and are refractory to assembly from whole-genome short-read datasets. Human subtelomere regions are highly enriched in both segmental duplication content and structural variations, and as a consequence are both impossible to assemble accurately and highly variable from individual to individual. Recently, we developed a pipeline for improved region-specific assembly called Regional Extension of Assemblies Using Linked-Reads (REXTAL). In this study, we evaluate REXTAL and genome-wide assembly (Supernova) approaches on 10X Genomics linked-reads data sets partitioned and barcoded using the Gel Bead in Emulsion (GEM) microfluidic method. Our results describe the accuracy and relative performance of these two approaches using the reference-based assessment module of QUAST. We show that REXTAL dramatically outperforms the Supernova whole genome assembler in subtelomeric segmental duplication regions, and results in highly accurate assemblies. Nearly all of the REXTAL misassemblies identified using default QUAST parameters simply pinpoint locations of tandem repeat arrays in the reference sequence where the repeat array length differs from that in the cognate REXTAL assembly by \u3e 1000 bp

Single-Molecule Analysis of Subtelomeres and Telomeres in Alternative Lengthening of Telomeres (ALT) Cells

BACKGROUND: Telomeric DNA is typically comprised of G-rich tandem repeat motifs and maintained by telomerase (Greider CW, Blackburn EH; Cell 51:887-898; 1987). In eukaryotes lacking telomerase, a variety of DNA repair and DNA recombination based pathways for telomere maintenance have evolved in organisms normally dependent upon telomerase for telomere elongation (Webb CJ, Wu Y, Zakian VA; Cold Spring Harb Perspect Biol 5:a012666; 2013); collectively called Alternative Lengthening of Telomeres (ALT) pathways. By measuring (TTAGGG) n tract lengths from the same large DNA molecules that were optically mapped, we simultaneously analyzed telomere length dynamics and subtelomere-linked structural changes at a large number of specific subtelomeric loci in the ALT-positive cell lines U2OS, SK-MEL-2 and Saos-2. RESULTS: Our results revealed loci-specific ALT telomere features. For example, while each subtelomere included examples of single molecules with terminal (TTAGGG) n tracts as well as examples of recombinant telomeric single molecules, the ratio of these molecules was subtelomere-specific, ranging from 33:1 (19p) to 1:25 (19q) in U2OS. The Saos-2 cell line shows a similar percentage of recombinant telomeres. The frequency of recombinant subtelomeres of SK-MEL-2 (11%) is about half that of U2OS and Saos-2 (24 and 19% respectively). Terminal (TTAGGG) n tract lengths and heterogeneity levels, the frequencies of telomere signal-free ends, and the frequency and size of retained internal telomere-like sequences (ITSs) at recombinant telomere fusion junctions all varied according to the specific subtelomere involved in a particular cell line. Very large linear extrachromosomal telomere repeat (ECTR) DNA molecules were found in all three cell lines; these are in principle capable of templating synthesis of new long telomere tracts via break-induced repair (BIR) long-tract DNA synthesis mechanisms and contributing to the very long telomere tract length and heterogeneity characteristic of ALT cells. Many of longest telomere tracts (both end-telomeres and linear ECTRs) displayed punctate CRISPR/Cas9-dependent (TTAGGG) n labeling patterns indicative of interspersion of stretches of non-canonical telomere repeats. CONCLUSION: Identifying individual subtelomeres and characterizing linked telomere (TTAGGG) n tract lengths and structural changes using our new single-molecule methodologies reveals the structural consequences of telomere damage, repair and recombination mechanisms in human ALT cells in unprecedented molecular detail and significant differences in different ALT-positive cell lines