NPGREAT: assembly of 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 mapped ultralong nanopore reads spanning subtelomeric segmental duplications to potentially overcome these difficulties. Linked-Read sets of DNA sequences 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). Mapped 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. Our method was tested for a subset of representative subtelomeres with ultralong nanopore read coverage in the haploid human cell line CHM13. A 10X Linked-Read dataset from CHM13 was combined with ultralong nanopore reads from the same genome to provide improved subtelomere assemblies. Comparison of Nanopore-only assemblies using SHASTA with our NPGREAT assemblies in the distal-most subtelomere regions showed that NPGREAT produced higher-quality and more complete assemblies than SHASTA alone when these regions had low ultralong nanopore coverage (such as cases where large segmental duplications were immediately adjacent to (TTAGGG) tracts). CONCLUSION: In genomic regions with large segmental duplications adjacent to telomeres, NPGREAT offers an alternative economical approach to improving assembly accuracy and coverage using linked-read datasets when more expensive HiFi datasets of 10-20 kb reads are unavailable.

Unraveling novel survivin mRNA transcripts in cancer cells using an in-house developed targeted high-throughput sequencing approach.

The human baculoviral IAP repeat containing 5 (BIRC5), also known as survivin, is a conserved member of the inhibitor of apoptosis protein (IAPs) family, which is normally expressed during embryonic and fetal development. Although the expression levels of survivin are low in terminally differentiated cells and/or tissues, they can be found notably increased in certain pathological conditions including malignant tumors. Conventional cloning and sequencing techniques have already confirmed that alternative splicing events of the survivin pre-mRNA result in five distinct alternative transcript variants. In the present study, however, we implemented an innovative, in-house developed, targeted DNA-seq assay to identify novel survivin alternative transcript variants with increased depth and coverage that high-throughput sequencing approaches offer. Bioinformatics analysis of the derived NGS datasets unveiled several novel splice junctions between annotated exons of survivin gene as well as the existence of a novel exon of 117 nt, spanning between the annotated exons 3 and 3B. Validation of the NGS findings with PCR-based assays, using variant-specific primers, led to the identification of fourteen novel survivin alternative splice variants (BIRC5 v.4 - v.17), which demonstrate wide expression profiles in a broad established panel of human cell lines. Although the presented novel findings provide a crystal-clear overview of the survivin mRNAs that are actually generated from the pre-mRNA, future studies should focus on the impending necessity of characterizing the biological function of all novel alternative transcript variants as well as the putative protein isoforms. Such studies will further contribute to our understanding of how the balance between survivin isoforms regulate malignant cell proliferation and apoptosis, providing novel diagnostic, prognostic and predictive biomarkers as well as therapeutic targets.

Nanopore Guided Assembly of Segmental Duplications near Telomeres.

Human subtelomere regions are highly enriched in large segmental duplications and structural variants, leading to many gaps and misassemblies in these regions. We develop a novel method, NPGREAT (NanoPore Guided REgional Assembly Tool), which combines Nanopore ultralong read datasets and short-read assemblies derived from 10x linked-reads to efficiently assemble these subtelomere regions into a single continuous sequence. We show that with the use of ultralong Nanopore reads as a guide, the highly accurate shorter linked-read sequence contigs are correctly oriented, ordered, spaced and extended. In the rare cases where a linked-read sequence contig contains inaccurately assembled segments, the use of Nanopore reads allows for detection and correction of this error. We tested NPGREAT on four representative subtelomeres of the NA12878 human genome (10p, 16p, 19q and 20p). The results demonstrate that the final computed assembly of each subtelomere is accurate and complete.