Detailed analysis of HTT repeat elements in human blood using targeted amplification-free long-read sequencing

Toward understanding the role of genomic repeat elements in neurodegenerative diseases

Neurodegenerative diseases cause great medical and economic burdens for both patients and society; however, the complex molecular mechanisms thereof are not yet well understood. With the development of high-coverage sequencing technology, researchers have started to notice that genomic repeat regions, previously neglected in search of disease culprits, are active contributors to multiple neurodegenerative diseases. In this review, we describe the association between repeat element variants and multiple degenerative diseases through genome-wide association studies and targeted sequencing. We discuss the identification of disease-relevant repeat element variants, further powered by the advancement of long-read sequencing technologies and their related tools, and summarize recent findings in the molecular mechanisms of repeat element variants in brain degeneration, such as those causing transcriptional silencing or RNA-mediated gain of toxic function. Furthermore, we describe how in silico predictions using innovative computational models, such as deep learning language models, could enhance and accelerate our understanding of the functional impact of repeat element variants. Finally, we discuss future directions to advance current findings for a better understanding of neurodegenerative diseases and the clinical applications of genomic repeat elements.

Sequence composition changes in short tandem repeats: heterogeneity, detection, mechanisms and clinical implications

Short tandem repeats (STRs) are a class of repetitive elements, composed of tandem arrays of 1-6 base pair sequence motifs, that comprise a substantial fraction of the human genome. STR expansions can cause a wide range of neurological and neuromuscular conditions, known as repeat expansion disorders, whose age of onset, severity, penetrance and/or clinical phenotype are influenced by the length of the repeats and their sequence composition. The presence of non-canonical motifs, depending on the type, frequency and position within the repeat tract, can alter clinical outcomes by modifying somatic and intergenerational repeat stability, gene expression and mutant transcript-mediated and/or protein-mediated toxicities. Here, we review the diverse structural conformations of repeat expansions, technological advances for the characterization of changes in sequence composition, their clinical correlations and the impact on disease mechanisms.

Long-read sequencing for brain tumors

Brain tumors and genomics have a long-standing history given that glioblastoma was the first cancer studied by the cancer genome atlas. The numerous and continuous advances through the decades in sequencing technologies have aided in the advanced molecular characterization of brain tumors for diagnosis, prognosis, and treatment. Since the implementation of molecular biomarkers by the WHO CNS in 2016, the genomics of brain tumors has been integrated into diagnostic criteria. Long-read sequencing, also known as third generation sequencing, is an emerging technique that allows for the sequencing of longer DNA segments leading to improved detection of structural variants and epigenetics. These capabilities are opening a way for better characterization of brain tumors. Here, we present a comprehensive summary of the state of the art of third-generation sequencing in the application for brain tumor diagnosis, prognosis, and treatment. We discuss the advantages and potential new implementations of long-read sequencing into clinical paradigms for neuro-oncology patients.

The Third-Generation Sequencing Challenge: Novel Insights for the Omic Sciences

The understanding of the human genome has been greatly improved by the advent of next-generation sequencing technologies (NGS). Despite the undeniable advantages responsible for their widespread diffusion, these methods have some constraints, mainly related to short read length and the need for PCR amplification. As a consequence, long-read sequencers, called third-generation sequencing (TGS), have been developed, promising to overcome NGS. Starting from the first prototype, TGS has progressively ameliorated its chemistries by improving both read length and base-calling accuracy, as well as simultaneously reducing the costs/base. Based on these premises, TGS is showing its potential in many fields, including the analysis of difficult-to-sequence genomic regions, structural variations detection, RNA expression profiling, DNA methylation study, and metagenomic analyses. Protocol standardization and the development of easy-to-use pipelines for data analysis will enhance TGS use, also opening the way for their routine applications in diagnostic contexts.

A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Long read sequencing on its way to the routine diagnostics of genetic diseases

The clinical application of technological progress in the identification of DNA alterations has always led to improvements of diagnostic yields in genetic medicine. At chromosome side, from cytogenetic techniques evaluating number and gross structural defects to genomic microarrays detecting cryptic copy number variants, and at molecular level, from Sanger method studying the nucleotide sequence of single genes to the high-throughput next-generation sequencing (NGS) technologies, resolution and sensitivity progressively increased expanding considerably the range of detectable DNA anomalies and alongside of Mendelian disorders with known genetic causes. However, particular genomic regions (i.e., repetitive and GC-rich sequences) are inefficiently analyzed by standard genetic tests, still relying on laborious, time-consuming and low-sensitive approaches (i.e., southern-blot for repeat expansion or long-PCR for genes with highly homologous pseudogenes), accounting for at least part of the patients with undiagnosed genetic disorders. Third generation sequencing, generating long reads with improved mappability, is more suitable for the detection of structural alterations and defects in hardly accessible genomic regions. Although recently implemented and not yet clinically available, long read sequencing (LRS) technologies have already shown their potential in genetic medicine research that might greatly impact on diagnostic yield and reporting times, through their translation to clinical settings. The main investigated LRS application concerns the identification of structural variants and repeat expansions, probably because techniques for their detection have not evolved as rapidly as those dedicated to single nucleotide variants (SNV) identification: gold standard analyses are karyotyping and microarrays for balanced and unbalanced chromosome rearrangements, respectively, and southern blot and repeat-primed PCR for the amplification and sizing of expanded alleles, impaired by limited resolution and sensitivity that have not been significantly improved by the advent of NGS. Nevertheless, more recently, with the increased accuracy provided by the latest product releases, LRS has been tested also for SNV detection, especially in genes with highly homologous pseudogenes and for haplotype reconstruction to assess the parental origin of alleles with de novo pathogenic variants. We provide a review of relevant recent scientific papers exploring LRS potential in the diagnosis of genetic diseases and its potential future applications in routine genetic testing.

Nanopore sequencing method for CTG18.1 expansion in TCF4 in late-onset Fuchs endothelial corneal dystrophy and a comparison of the structural features of cornea with first-degree relatives

BACKGROUND

To evaluate the relationship between the number of trinucleotide repeats (TNR) in late-onset Fuchs corneal endothelial dystrophy (FCED) and to compare the endothelial properties of FCED, first-degree relatives, and controls.

METHODS

Blood samples were collected from FCEDs to determine TNR number. The FCED patients, first-degree relatives, and controls were examined with specular microscopy for central corneal thickness (CCT), endothelial cell density (ECD), pleomorphism and polymegatism, and with corneal topography for specific indicators such as (i) displacement of thinnest point of cornea, (ii) loss of isopachs, (iii) focal posterior surface depression towards anterior chamber.

RESULTS

This study included 92 patients with FCED, 92 first-degree relatives, and 96 controls. CCT was thickest in FCEDs (558.0 μm) (p < 0.05) while there was no difference between relatives (533.0 μm) and controls (530.4 μm) (p = 0.845). ECD was decreased in both FCED (2069.2 mm2) and relatives (2171.4 mm2) than controls (2822.9 mm2) (p < 0.05 in both). The presence of pleomorphism and polymegatism was significant in patients with FCED (93.4% and 93.4%, respectively), relatives (86.9% and 86.04%, respectively), and controls (8.33% and 1.04%, respectively) (p < 0.05). Specific topographic indicators differed among the groups (p < 0.05). The mean repeat number of the FCED patients was 17.48 ± 4.54 (12-25) times. The TNR number of FCED cases correlated with the relative CCT (p < 0.05, R = 0.615) and cell density (p = 0.009, R = -0.499).

CONCLUSIONS

A strong association between the corneal endothelium in relatives and TNR number of FCEDs was defined. Relatives tended to have fewer corneal endothelial cells, even though they did not have clinical findings.

Application of long read sequencing in rare diseases: The longer, the better?

Rare diseases encompass a diverse group of genetic disorders that affect a small proportion of the population. Identifying the underlying genetic causes of these conditions presents significant challenges due to their genetic heterogeneity and complexity. Conventional short-read sequencing (SRS) techniques have been widely used in diagnosing and investigating of rare diseases, with limitations due to the nature of short-read lengths. In recent years, long read sequencing (LRS) technologies have emerged as a valuable tool in overcoming these limitations. This minireview provides a concise overview of the applications of LRS in rare disease research and diagnosis, including the identification of disease-causing tandem repeat expansions, structural variations, and comprehensive analysis of pathogenic variants with LRS.

A novel quantitative targeted analysis of X-chromosome inactivation (XCI) using nanopore sequencing

X-chromosome inactivation (XCI) analyses often assist in diagnostics of X-linked traits, however accurate assessment remains challenging with current methods. We developed a novel strategy using amplification-free Cas9 enrichment and Oxford nanopore technologies sequencing called XCI-ONT, to investigate and rigorously quantify XCI in human androgen receptor gene (AR) and human X-linked retinitis pigmentosa 2 gene (RP2). XCI-ONT measures methylation over 116 CpGs in AR and 58 CpGs in RP2, and separate parental X-chromosomes without PCR bias. We show the usefulness of the XCI-ONT strategy over the PCR-based golden standard XCI technique that only investigates one or two CpGs per gene. The results highlight the limitations of using the golden standard technique when the XCI pattern is partially skewed and the advantages of XCI-ONT to rigorously quantify XCI. This study provides a universal XCI-method on DNA, which is highly valuable in clinical and research framework of X-linked traits.

Applications of long-read sequencing to Mendelian genetics

Advances in clinical genetic testing, including the introduction of exome sequencing, have uncovered the molecular etiology for many rare and previously unsolved genetic disorders, yet more than half of individuals with a suspected genetic disorder remain unsolved after complete clinical evaluation. A precise genetic diagnosis may guide clinical treatment plans, allow families to make informed care decisions, and permit individuals to participate in N-of-1 trials; thus, there is high interest in developing new tools and techniques to increase the solve rate. Long-read sequencing (LRS) is a promising technology for both increasing the solve rate and decreasing the amount of time required to make a precise genetic diagnosis. Here, we summarize current LRS technologies, give examples of how they have been used to evaluate complex genetic variation and identify missing variants, and discuss future clinical applications of LRS. As costs continue to decrease, LRS will find additional utility in the clinical space fundamentally changing how pathological variants are discovered and eventually acting as a single-data source that can be interrogated multiple times for clinical service.

Variation within major internal repeats of KSHV in vivo

Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiologic agent of Kaposi's sarcoma (KS), yet the viral genetic factors that lead to the development of KS in KSHV-infected individuals have not been fully elucidated. Nearly, all previous analyses of KSHV genomic evolution and diversity have excluded the three major internal repeat regions: the two origins of lytic replication, internal repeats 1 and 2 (IR1 and IR2), and the latency-associated nuclear antigen (LANA) repeat domain (LANAr). These regions encode protein domains that are essential to the KSHV infection cycle but have been rarely sequenced due to their extended repetitive nature and high guanine and cytosine (GC) content. The limited data available suggest that their sequences and repeat lengths are more heterogeneous across individuals than in the remainder of the KSHV genome. To assess their diversity, the full-length IR1, IR2, and LANAr sequences, tagged with unique molecular identifiers (UMIs), were obtained by Pacific Biosciences' single-molecule real-time sequencing (SMRT-UMI) from twenty-four tumors and six matching oral swabs from sixteen adults in Uganda with advanced KS. Intra-host single-nucleotide variation involved an average of 0.16 per cent of base positions in the repeat regions compared to a nearly identical average of 0.17 per cent of base positions in the remainder of the genome. Tandem repeat unit (TRU) counts varied by only one from the intra-host consensus in a majority of individuals. Including the TRU indels, the average intra-host pairwise identity was 98.3 per cent for IR1, 99.6 per cent for IR2 and 98.9 per cent for LANAr. More individuals had mismatches and variable TRU counts in IR1 (twelve/sixteen) than in IR2 (two/sixteen). There were no open reading frames in the Kaposin coding sequence inside IR2 in at least fifty-five of ninety-six sequences. In summary, the KSHV major internal repeats, like the rest of the genome in individuals with KS, have low diversity. IR1 was the most variable among the repeats, and no intact Kaposin reading frames were present in IR2 of the majority of genomes sampled.

Targeted long-read sequencing captures CRISPR editing and AAV integration outcomes in brain

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing is an emerging therapeutic modality that shows promise in Huntington's disease and spinocerebellar ataxia (SCA) mouse models. However, advancing CRISPR-based therapies requires methods to fully define in vivo editing outcomes. Here, we use polymerase-free, targeted long-read nanopore sequencing and evaluate single- and dual-gRNA AAV-CRISPR editing of human ATXN2 in transgenic mouse models of SCA type 2 (SCA2). Unbiased high sequencing coverage showed 10%-25% editing. Along with intended edits there was AAV integration, 1%-2% of which contained the entire AAV genome and were largely unmethylated. More than 150 kb deletions at target loci and rearrangements of the transgenic allele (1%) were also found. In contrast, PCR-based nanopore sequencing showed bias for partial AAV fragments and inverted terminal repeats (ITRs) and failed to detect full-length AAV. Cumulatively this work defines the spectrum of outcomes of CRISPR editing in mouse brain after AAV gene transfer using an unbiased long-read sequencing approach.

CRISPR-Cas9 Targeted Enrichment and Next-Generation Sequencing for Mutation Detection

Despite the rapid application of next-generation sequencing (NGS) technologies, target sequencing in regions of the genome is often required to diagnose many genetic diseases. Target enrichment can be an effective factor in reducing the cost of sequencing and the duration of sequencing. Recently, several clustered system regularly interspaced short palindromic repeats (CRISPR)-based methods (amplification-free sequencing) have been developed to target enrichment in combination with one of the NGS platforms. CRISPR-based target enrichment strategies act as an auxiliary tool to improve NGS analytical performance, thereby indirectly facilitating nucleic acid detection. The direct DNA cleavage approach by CRISPR-Cas at genome-specific sites enhances the possibility of separating native large fragments from disease-related genomic regions. The CRISPR-Cas can isolate the target region without any amplification; subsequently, long-read sequencing technologies were also implemented. These methods, as promising tools, have the ability to assess genetic and epigenetic composition for clinical application and treatment responses in cancer precision medicine. By modifying CRISPR-based enrichment protocols, it was possible to identify different types of mutations, including structural variants, short tandem repeats, fusion genes, and mobile elements. The Cas9 can specifically eliminate wild-type sequences, and it also enables the enrichment and detection of small amounts of tumor DNA fragments among the highly heterogeneous fragments of wild-type DNA.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

Gene Therapy, A Potential Therapeutic Tool for Neurological and Neuropsychiatric Disorders: Applications, Challenges and Future Perspective

Neurological and neuropsychiatric disorders are the main risks for the health care system, exhibiting a huge socioeconomic load. The available range of pharmacotherapeutics mostly provides palliative consequences and fails to treat such conditions. The molecular etiology of various neurological and neuropsychiatric disorders is mostly associated with a change in genetic background, which can be inherited/triggered by other environmental factors. To address such conditions, gene therapy is considered a potential approach claiming a permanent cure of the disease primarily by deletion, silencing, or edition of faulty genes and by insertion of healthier genes. In gene therapy, vectors (viral/nonvial) play an important role in delivering the desired gene to a specific region of the brain. Targeted gene therapy has unraveled opportunities for the treatment of many neurological and neuropsychiatric disorders. For improved gene delivery, the current techniques mainly focus on designing a precise viral vector, plasmid transfection, nanotechnology, microRNA, and in vivo clustered regulatory interspaced short palindromic repeats (CRISPR)-based therapy. These latest techniques have great benefits in treating predominant neurological and neurodevelopmental disorders, including Parkinson's disease, Alzheimer's disease, and autism spectrum disorder, as well as rarer diseases. Nevertheless, all these delivery methods have their limitations, including immunogenic reactions, off-target effects, and a deficiency of effective biomarkers to appreciate the effectiveness of therapy. In this review, we present a summary of the current methods in targeted gene delivery, followed by the limitations and future direction of gene therapy for the cure of neurological and neuropsychiatric disorders.

Repeat Detector: versatile sizing of expanded tandem repeats and identification of interrupted alleles from targeted DNA sequencing

Targeted DNA sequencing approaches will improve how the size of short tandem repeats is measured for diagnostic tests and preclinical studies. The expansion of these sequences causes dozens of disorders, with longer tracts generally leading to a more severe disease. Interrupted alleles are sometimes present within repeats and can alter disease manifestation. Determining repeat size mosaicism and identifying interruptions in targeted sequencing datasets remains a major challenge. This is in part because standard alignment tools are ill-suited for repetitive and unstable sequences. To address this, we have developed Repeat Detector (RD), a deterministic profile weighting algorithm for counting repeats in targeted sequencing data. We tested RD using blood-derived DNA samples from Huntington's disease and Fuchs endothelial corneal dystrophy patients sequenced using either Illumina MiSeq or Pacific Biosciences single-molecule, real-time sequencing platforms. RD was highly accurate in determining repeat sizes of 609 blood-derived samples from Huntington's disease individuals and did not require prior knowledge of the flanking sequences. Furthermore, RD can be used to identify alleles with interruptions and provide a measure of repeat instability within an individual. RD is therefore highly versatile and may find applications in the diagnosis of expanded repeat disorders and in the development of novel therapies.

Long-read sequencing for molecular diagnostics in constitutional genetic disorders

Long-read sequencing (LRS) has been around for more than a decade, but widespread adoption of the technology has been slow due to the perceived high error rates and high sequencing cost. This is changing due to the recent advancements to produce highly accurate sequences and the reducing costs. LRS promises significant improvement over short read sequencing in four major areas: (1) better detection of structural variation (2) better resolution of highly repetitive or nonunique regions (3) accurate long-range haplotype phasing and (4) the detection of base modifications natively from the sequencing data. Several successful applications of LRS have demonstrated its ability to resolve molecular diagnoses where short-read sequencing fails to identify a cause. However, the argument for increased diagnostic yield from LRS remains to be validated. Larger cohort studies may be required to establish the realistic boundaries of LRS's clinical utility and analytical validity, as well as the development of standards for clinical applications. We discuss the limitations of the current standard of care, and contrast with the applications and advantages of two major LRS platforms, PacBio and Oxford Nanopore, for molecular diagnostics of constitutional disorders, and present a critical argument about the potential of LRS in diagnostic settings.

Comprehensive Analysis of Fragile X Syndrome: Full Characterization of the FMR1 Locus by Long-Read Sequencing

BACKGROUND

Fragile X syndrome (FXS) is the most frequent cause of inherited X-linked intellectual disability. Conventional FXS genetic testing methods mainly focus on FMR1 CGG expansions and fail to identify AGG interruptions, rare intragenic variants, and large gene deletions.

METHODS

A long-range PCR and long-read sequencing-based assay termed comprehensive analysis of FXS (CAFXS) was developed and evaluated in Coriell and clinical samples by comparing to Southern blot analysis and triplet repeat-primed PCR (TP-PCR).

RESULTS

CAFXS accurately detected the number of CGG repeats in the range of 93 to at least 940 with mass fraction of 0.5% to 1% in the background of normal alleles, which was 2-4-fold analytically more sensitive than TP-PCR. All categories of mutations detected by control methods, including full mutations in 30 samples, were identified by CAFXS for all 62 clinical samples. CAFXS accurately determined AGG interruptions in all 133 alleles identified, even in mosaic alleles. CAFXS successfully identified 2 rare intragenic variants including the c.879A > C variant in exon 9 and a 697-bp microdeletion flanking upstream of CGG repeats, which disrupted primer annealing in TP-PCR assay. In addition, CAFXS directly determined the breakpoints of a 237.1-kb deletion and a 774.0-kb deletion encompassing the entire FMR1 gene in 2 samples.

CONCLUSIONS

Long-read sequencing-based CAFXS represents a comprehensive assay for identifying FMR1 CGG expansions, AGG interruptions, rare intragenic variants, and large gene deletions, which greatly improves the genetic screening and diagnosis for FXS.

Characterization of full-length CNBP expanded alleles in myotonic dystrophy type 2 patients by Cas9-mediated enrichment and nanopore sequencing

Myotonic dystrophy type 2 (DM2) is caused by CCTG repeat expansions in the CNBP gene, comprising 75 to >11,000 units and featuring extensive mosaicism, making it challenging to sequence fully expanded alleles. To overcome these limitations, we used PCR-free Cas9-mediated nanopore sequencing to characterize CNBP repeat expansions at the single-nucleotide level in nine DM2 patients. The length of normal and expanded alleles can be assessed precisely using this strategy, agreeing with traditional methods, and revealing the degree of mosaicism. We also sequenced an entire ~50 kbp expansion, which has not been achieved previously for DM2 or any other repeat-expansion disorders. Our approach precisely counted the repeats and identified the repeat pattern for both short interrupted and uninterrupted alleles. Interestingly, in the expanded alleles, only two DM2 samples featured the expected pure CCTG repeat pattern, while the other seven presented also TCTG blocks at the 3′ end, which have not been reported before in DM2 patients, but confirmed hereby with orthogonal methods. The demonstrated approach simultaneously determines repeat length, structure/motif, and the extent of somatic mosaicism, promising to improve the molecular diagnosis of DM2 and achieve more accurate genotype–phenotype correlations for the better stratification of DM2 patients in clinical trials.

Third-Generation Cytogenetic Analysis: Diagnostic Application of Long-Read Sequencing

Copy number variants (CNVs) play important roles in the pathogenesis of several genetic syndromes. Traditional and molecular karyotyping are considered the first-tier diagnostic tests to detect macroscopic and cryptic deletions/duplications. However, their time-consuming and laborious experimental protocols protract diagnostic times from 3 to 15 days. Nanopore sequencing has the ability to reduce time to results for the detection of CNVs with the same resolution of current state-of-the-art diagnostic tests. Nanopore sequencing was compared to molecular karyotyping for the detection of pathogenic CNVs of seven patients with previously diagnosed causative CNVs of different sizes and cellular fractions. Larger chromosomal anomalies included trisomy 21 and mosaic tetrasomy 12p. Among smaller CNVs, two genomic imbalances of 1.3 Mb, a small deletion of 170 kb, and two mosaic deletions (1.2 Mb and 408 kb) were tested. DNA was sequenced and data generated during runs were analyzed in online mode. All pathogenic CNVs were identified with detection time inversely proportional to size and cellular fraction. Aneuploidies were called after only 30 minutes of sequencing, whereas 30 hours were needed to call small CNVs. These results demonstrate the clinical utility of our approach that allows the molecular diagnosis of genomic disorders within a 30-minute to 30-hour time frame and its easy implementation as a routinary diagnostic tool.

Uncovering Essential Tremor Genetics: The Promise of Long-Read Sequencing

Long-read sequencing (LRS) technologies have been recently introduced to overcome intrinsic limitations of widely-used next-generation sequencing (NGS) technologies, namely the sequencing limited to short-read fragments (150–300 base pairs). Since its introduction, LRS has permitted many successes in unraveling hidden mutational mechanisms. One area in clinical neurology in need of rethinking as it applies to genetic mechanisms is essential tremor (ET). This disorder, among the most common in neurology, is a syndrome often exhibiting an autosomal dominant pattern of inheritance whose large phenotypic spectrum suggest a multitude of genetic etiologies. Exome sequencing has revealed the genetic etiology only in rare ET families (FUS, SORT1, SCN4A, NOS3, KCNS2, HAPLN4/BRAL2, and USP46). We hypothesize that a reason for this shortcoming may be non-classical genetic mechanism(s) underpinning ET, among them trinucleotide, tetranucleotide, or pentanucleotide repeat disorders. In support of this hypothesis, trinucleotide (e.g., GGC repeats in NOTCH2NLC) and pentanucleotide repeat disorders (e.g., ATTTC repeats in STARD7) have been revealed as pathogenic in patients with a past history of what has come to be referred to as “ET plus,” bilateral hand tremor associated with epilepsy and/or leukoencephalopathy. A systematic review of LRS in neurodegenerative disorders showed that 10 of the 22 (45%) genetic etiologies ascertained by LRS include tremor in their phenotypic spectrum, suggesting that future clinical applications of LRS for tremor disorders may uncover genetic subtypes of familial ET that have eluded NGS, particularly those with associated leukoencephalopathy or family history of epilepsy. LRS provides a pathway for potentially uncovering novel genes and genetic mechanisms, helping narrow the large proportion of “idiopathic” ET.

CRISPR-Cas9 induces large structural variants at on-target and off-target sites in vivo that segregate across generations

CRISPR-Cas9 genome editing has potential to cure diseases without current treatments, but therapies must be safe. Here we show that CRISPR-Cas9 editing can introduce unintended mutations in vivo, which are passed on to the next generation. By editing fertilized zebrafish eggs using four guide RNAs selected for off-target activity in vitro, followed by long-read sequencing of DNA from >1100 larvae, juvenile and adult fish across two generations, we find that structural variants (SVs), i.e., insertions and deletions ≥50 bp, represent 6% of editing outcomes in founder larvae. These SVs occur both at on-target and off-target sites. Our results also illustrate that adult founder zebrafish are mosaic in their germ cells, and that 26% of their offspring carries an off-target mutation and 9% an SV. Hence, pre-testing for off-target activity and SVs using patient material is advisable in clinical applications, to reduce the risk of unanticipated effects with potentially large implications.

Detecting MUC1 Variants in Patients Clinicopathologically Diagnosed With Having Autosomal Dominant Tubulointerstitial Kidney Disease

Introduction

Autosomal dominant tubulointerstitial kidney disease (ADTKD)-MUC1 is predominantly caused by frameshift mutations owing to a single-base insertion into the variable number tandem repeat (VNTR) region in MUC1. Because of the complexity of the variant hotspot, identification using short-read sequencers (SRSs) is challenging. Although recent studies have revealed the usefulness of long-read sequencers (LRSs), the prevalence of MUC1 variants in patients with clinically suspected ADTKD remains unknown. We aimed to clarify this prevalence and the genetic characteristics and clinical manifestations of ADTKD-MUC1 in a Japanese population using an SRS and an LRS.

Methods

From January 2015 to December 2019, genetic analysis was performed using an SRS in 48 patients with clinically suspected ADTKD. Additional analyses were conducted using an LRS in patients with negative SRS results.

Results

Short-read sequencing results revealed MUC1 variants in 1 patient harboring a cytosine insertion in the second repeat unit of the VNTR region; however, deeper VNTR regions could not be read by the SRS. Therefore, we conducted long-read sequencing analysis of 39 cases and detected MUC1 VNTR variants in 8 patients (in total, 9 patients from unrelated families). With the inclusion of family-affected patients (n = 31), the median age at the development of end-stage kidney disease (ESKD) was 45 years (95% CI: 40–40 years).

Conclusion

In Japan, the detection rate of MUC1 variants in patients with clinically suspected ADTKD was 18.8%. More than 20% of patients with negative SRS results had MUC1 variants detected by an LRS.

Gene Therapy Approach with an Emphasis on Growth Factors: Theoretical and Clinical Outcomes in Neurodegenerative Diseases

The etiology of many neurological diseases affecting the central nervous system (CNS) is unknown and still needs more effective and specific therapeutic approaches. Gene therapy has a promising future in treating neurodegenerative disorders by correcting the genetic defects or by therapeutic protein delivery and is now an attraction for neurologists to treat brain disorders, like Alzheimer’s disease, Parkinson’s disease, amyotrophic lateral sclerosis, spinal muscular atrophy, spinocerebellar ataxia, epilepsy, Huntington’s disease, stroke, and spinal cord injury. Gene therapy allows the transgene induction, with a unique expression in cells’ substrate. This article mainly focuses on the delivering modes of genetic materials in the CNS, which includes viral and non-viral vectors and their application in gene therapy. Despite the many clinical trials conducted so far, data have shown disappointing outcomes. The efforts done to improve outcomes, efficacy, and safety in the identification of targets in various neurological disorders are also discussed here. Adapting gene therapy as a new therapeutic approach for treating neurological disorders seems to be promising, with early detection and delivery of therapy before the neuron is lost, helping a lot the development of new therapeutic options to translate to the clinic.

Patterns of CAG repeat instability in the central nervous system and periphery in Huntington's disease and in spinocerebellar ataxia type 1

The expanded HTT CAG repeat causing Huntington’s disease (HD) exhibits somatic expansion proposed to drive the rate of disease onset by eliciting a pathological process that ultimately claims vulnerable cells. To gain insight into somatic expansion in humans, we performed comprehensive quantitative analyses of CAG expansion in ~50 central nervous system (CNS) and peripheral postmortem tissues from seven adult-onset and one juvenile-onset HD individual. We also assessed ATXN1 CAG repeat expansion in brain regions of an individual with a neurologically and pathologically distinct repeat expansion disorder, spinocerebellar ataxia type 1 (SCA1). Our findings reveal similar profiles of tissue instability in all HD individuals, which, notably, were also apparent in the SCA1 individual. CAG expansion was observed in all tissues, but to different degrees, with multiple cortical regions and neostriatum tending to have the greatest instability in the CNS, and liver in the periphery. These patterns indicate different propensities for CAG expansion contributed by disease locus-independent trans-factors and demonstrate that expansion per se is not sufficient to cause cell type or disease-specific pathology. Rather, pathology may reflect distinct toxic processes triggered by different repeat lengths across cell types and diseases. We also find that the HTT CAG length-dependent expansion propensity of an individual is reflected in all tissues and in cerebrospinal fluid. Our data indicate that peripheral cells may be a useful source to measure CAG expansion in biomarker assays for therapeutic efforts, prompting efforts to dissect underlying mechanisms of expansion that may differ between the brain and periphery.

Understanding the Potential of Genome Editing in Parkinson's Disease

CRISPR is a simple and cost-efficient gene-editing technique that has become increasingly popular over the last decades. Various CRISPR/Cas-based applications have been developed to introduce changes in the genome and alter gene expression in diverse systems and tissues. These novel gene-editing techniques are particularly promising for investigating and treating neurodegenerative diseases, including Parkinson's disease, for which we currently lack efficient disease-modifying treatment options. Gene therapy could thus provide treatment alternatives, revolutionizing our ability to treat this disease. Here, we review our current knowledge on the genetic basis of Parkinson's disease to highlight the main biological pathways that become disrupted in Parkinson's disease and their potential as gene therapy targets. Next, we perform a comprehensive review of novel delivery vehicles available for gene-editing applications, critical for their successful application in both innovative research and potential therapies. Finally, we review the latest developments in CRISPR-based applications and gene therapies to understand and treat Parkinson's disease. We carefully examine their advantages and shortcomings for diverse gene-editing applications in the brain, highlighting promising avenues for future research.

Straglr: discovering and genotyping tandem repeat expansions using whole genome long-read sequences

Tandem repeat (TR) expansion is the underlying cause of over 40 neurological disorders. Long-read sequencing offers an exciting avenue over conventional technologies for detecting TR expansions. Here, we present Straglr, a robust software tool for both targeted genotyping and novel expansion detection from long-read alignments. We benchmark Straglr using various simulations, targeted genotyping data of cell lines carrying expansions of known diseases, and whole genome sequencing data with chromosome-scale assembly. Our results suggest that Straglr may be useful for investigating disease-associated TR expansions using long-read sequencing.

Deciphering Neurodegenerative Diseases Using Long-Read Sequencing

Neurodegenerative diseases exhibit chronic progressive lesions in the central and peripheral nervous systems with unclear causes. The search for pathogenic mutations in human neurodegenerative diseases has benefited from massively parallel short-read sequencers. However, genomic regions, including repetitive elements, especially with high/low GC content, are far beyond the capability of conventional approaches. Recently, long-read single-molecule DNA sequencing technologies have emerged and enabled researchers to study genomes, transcriptomes, and metagenomes at unprecedented resolutions. The identification of novel mutations in unresolved neurodegenerative disorders, the characterization of causative repeat expansions, and the direct detection of epigenetic modifications on naive DNA by virtue of long-read sequencers will further expand our understanding of neurodegenerative diseases. In this article, we review and compare 2 prevailing long-read sequencing technologies, Pacific Biosciences and Oxford Nanopore Technologies, and discuss their applications in neurodegenerative diseases.

An update on the neurological short tandem repeat expansion disorders and the emergence of long-read sequencing diagnostics

BACKGROUND

Short tandem repeat (STR) expansion disorders are an important cause of human neurological disease. They have an established role in more than 40 different phenotypes including the myotonic dystrophies, Fragile X syndrome, Huntington's disease, the hereditary cerebellar ataxias, amyotrophic lateral sclerosis and frontotemporal dementia.

MAIN BODY

STR expansions are difficult to detect and may explain unsolved diseases, as highlighted by recent findings including: the discovery of a biallelic intronic 'AAGGG' repeat in RFC1 as the cause of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS); and the finding of 'CGG' repeat expansions in NOTCH2NLC as the cause of neuronal intranuclear inclusion disease and a range of clinical phenotypes. However, established laboratory techniques for diagnosis of repeat expansions (repeat-primed PCR and Southern blot) are cumbersome, low-throughput and poorly suited to parallel analysis of multiple gene regions. While next generation sequencing (NGS) has been increasingly used, established short-read NGS platforms (e.g., Illumina) are unable to genotype large and/or complex repeat expansions. Long-read sequencing platforms recently developed by Oxford Nanopore Technology and Pacific Biosciences promise to overcome these limitations to deliver enhanced diagnosis of repeat expansion disorders in a rapid and cost-effective fashion.

CONCLUSION

We anticipate that long-read sequencing will rapidly transform the detection of short tandem repeat expansion disorders for both clinical diagnosis and gene discovery.

Long-read targeted sequencing uncovers clinicopathological associations for C9orf72-linked diseases

To examine the length of a hexanucleotide expansion in C9orf72, which represents the most frequent genetic cause of frontotemporal lobar degeneration and motor neuron disease, we employed a targeted amplification-free long-read sequencing technology: No-Amp sequencing. In our cross-sectional study, we assessed cerebellar tissue from 28 well-characterized C9orf72 expansion carriers. We obtained 3507 on-target circular consensus sequencing reads, of which 814 bridged the C9orf72 repeat expansion (23%). Importantly, we observed a significant correlation between expansion sizes obtained using No-Amp sequencing and Southern blotting (P = 5.0 × 10-4). Interestingly, we also detected a significant survival advantage for individuals with smaller expansions (P = 0.004). Additionally, we uncovered that smaller expansions were significantly associated with higher levels of C9orf72 transcripts containing intron 1b (P = 0.003), poly(GP) proteins (P = 1.3 × 10- 5), and poly(GA) proteins (P = 0.005). Thorough examination of the composition of the expansion revealed that its GC content was extremely high (median: 100%) and that it was mainly composed of GGGGCC repeats (median: 96%), suggesting that expanded C9orf72 repeats are quite pure. Taken together, our findings demonstrate that No-Amp sequencing is a powerful tool that enables the discovery of relevant clinicopathological associations, highlighting the important role played by the cerebellar size of the expanded repeat in C9orf72-linked diseases.

Robust Detection of Somatic Mosaicism and Repeat Interruptions by Long-Read Targeted Sequencing in Myotonic Dystrophy Type 1

Myotonic dystrophy type 1 (DM1) is the most complex and variable trinucleotide repeat disorder caused by an unstable CTG repeat expansion, reaching up to 4000 CTG in the most severe cases. The genetic and clinical variability of DM1 depend on the sex and age of the transmitting parent, but also on the CTG repeat number, presence of repeat interruptions and/or on the degree of somatic instability. Currently, it is difficult to simultaneously and accurately determine these contributing factors in DM1 patients due to the limitations of gold standard methods used in molecular diagnostics and research laboratories. Our study showed the efficiency of the latest PacBio long-read sequencing technology to sequence large CTG trinucleotides, detect multiple and single repeat interruptions and estimate the levels of somatic mosaicism in DM1 patients carrying complex CTG repeat expansions inaccessible to most methods. Using this innovative approach, we revealed the existence of de novo CCG interruptions associated with CTG stabilization/contraction across generations in a new DM1 family. We also demonstrated that our method is suitable to sequence the DM1 locus and measure somatic mosaicism in DM1 families carrying more than 1000 pure CTG repeats. Better characterization of expanded alleles in DM1 patients can significantly improve prognosis and genetic counseling, not only in DM1 but also for other tandem DNA repeat disorders.

What is the Pathogenic CAG Expansion Length in Huntington's Disease?

Huntington's disease (HD) (OMIM 143100) is caused by an expanded CAG repeat tract in the HTT gene. The inherited CAG length is known to expand further in somatic and germline cells in HD subjects. Age at onset of the disease is inversely correlated with the inherited CAG length, but is further modulated by a series of genetic modifiers which are most likely to act on the CAG repeat in HTT that permit it to further expand. Longer repeats are more prone to expansions, and this expansion is age dependent and tissue-specific. Given that the inherited tract expands through life and most subjects develop disease in mid-life, this implies that in cells that degenerate, the CAG length is likely to be longer than the inherited length. These findings suggest two thresholds- the inherited CAG length which permits further expansion, and the intracellular pathogenic threshold, above which cells become dysfunctional and die. This two-step mechanism has been previously proposed and modelled mathematically to give an intracellular pathogenic threshold at a tract length of 115 CAG (95% confidence intervals 70- 165 CAG). Empirically, the intracellular pathogenic threshold is difficult to determine. Clues from studies of people and models of HD, and from other diseases caused by expanded repeat tracts, place this threshold between 60- 100 CAG, most likely towards the upper part of that range. We assess this evidence and discuss how the intracellular pathogenic threshold in manifest disease might be better determined. Knowing the cellular pathogenic threshold would be informative for both understanding the mechanism in HD and deploying treatments.

Approaches to Sequence the HTT CAG Repeat Expansion and Quantify Repeat Length Variation

BACKGROUND

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by the expansion of the HTT CAG repeat. Affected individuals inherit ≥36 repeats and longer alleles cause earlier onset, greater disease severity and faster disease progression. The HTT CAG repeat is genetically unstable in the soma in a process that preferentially generates somatic expansions, the proportion of which is associated with disease onset, severity and progression. Somatic mosaicism of the HTT CAG repeat has traditionally been assessed by semi-quantitative PCR-electrophoresis approaches that have limitations (e.g., no information about sequence variants). Genotyping-by-sequencing could allow for some of these limitations to be overcome.

OBJECTIVE

To investigate the utility of PCR sequencing to genotype large (>50 CAGs) HD alleles and to quantify the associated somatic mosaicism.

METHODS

We have applied MiSeq and PacBio sequencing to PCR products of the HTT CAG repeat in transgenic R6/2 mice carrying ∼55, ∼110, ∼255 and ∼470 CAGs. For each of these alleles, we compared the repeat length distributions generated for different tissues at two ages.

RESULTS

We were able to sequence the CAG repeat full length in all samples. However, the repeat length distributions for samples with ∼470 CAGs were biased towards shorter repeat lengths.

CONCLUSION

PCR sequencing can be used to sequence all the HD alleles considered, but this approach cannot be used to estimate modal allele size or quantify somatic expansions for alleles ⪢250 CAGs. We review the limitations of PCR sequencing and alternative approaches that may allow the quantification of somatic contractions and very large somatic expansions.

Amplification-free long-read sequencing reveals unforeseen CRISPR-Cas9 off-target activity

BACKGROUND

One ongoing concern about CRISPR-Cas9 genome editing is that unspecific guide RNA (gRNA) binding may induce off-target mutations. However, accurate prediction of CRISPR-Cas9 off-target activity is challenging. Here, we present SMRT-OTS and Nano-OTS, two novel, amplification-free, long-read sequencing protocols for detection of gRNA-driven digestion of genomic DNA by Cas9 in vitro.

RESULTS

The methods are assessed using the human cell line HEK293, re-sequenced at 18x coverage using highly accurate HiFi SMRT reads. SMRT-OTS and Nano-OTS are first applied to three different gRNAs targeting HEK293 genomic DNA, resulting in a set of 55 high-confidence gRNA cleavage sites identified by both methods. Twenty-five of these sites are not reported by off-target prediction software, either because they contain four or more single nucleotide mismatches or insertion/deletion mismatches, as compared with the human reference. Additional experiments reveal that 85% of Cas9 cleavage sites are also found by other in vitro-based methods and that on- and off-target sites are detectable in gene bodies where short-reads fail to uniquely align. Even though SMRT-OTS and Nano-OTS identify several sites with previously validated off-target editing activity in cells, our own CRISPR-Cas9 editing experiments in human fibroblasts do not give rise to detectable off-target mutations at the in vitro-predicted sites. However, indel and structural variation events are enriched at the on-target sites.

CONCLUSIONS

Amplification-free long-read sequencing reveals Cas9 cleavage sites in vitro that would have been difficult to predict using computational tools, including in dark genomic regions inaccessible by short-read sequencing.

Xdrop: Targeted sequencing of long DNA molecules from low input samples using droplet sorting

Long-read sequencing can resolve regions of the genome that are inaccessible to short reads, and therefore are ideal for genome-gap closure, solving structural rearrangements and sequencing through repetitive elements. Here we introduce the Xdrop technology: a novel microfluidic-based system that allows for targeted enrichment of long DNA molecules starting from only a few nanograms of DNA. Xdrop is based on the isolation of long DNA fragments in millions of droplets, where the droplets containing a target sequence of interest are fluorescently labeled and sorted using flow cytometry. The final product from the Xdrop procedure is an enriched population of long DNA molecules that can be investigated by sequencing. To demonstrate the capability of Xdrop, we performed enrichment of the human papilloma virus 18 integrated into the genome of human HeLa cells. Analysis of the sequencing reads resolved three HPV18-chr8 integrations at base-pair resolution, and the captured fragments extended up to 30 kb into the human genome at the integration sites. Further, we enriched the complete TP53 locus in a leukemia cell line and could successfully phase coexisting mutations using PacBio sequencing. In summary, our results show that Xdrop is an efficient enrichment technology for studying complex genomic regions.

CRISPR/Cas9 novel therapeutic road for the treatment of neurodegenerative diseases

CRISPR (clustered regularly interspaced short palindromic Repeats)/Cas9 is a new genetic editing technology that can be a beneficial method to advance gene therapy. CRISPR technology is a defense system of some bacteria against invading viruses. Genome editing based on the CRISPR/Cas9 system is an efficient and potential technology that can be a viable alternative to traditional methods. This system is a compound of a short guide RNAs (gRNAs) for identifying the target DNA sequence and Cas9 protein as nuclease for breaking and cutting of DNA. In this review, recent advances in the CRISPR/Cas9-mediated genome editing tools are presented as well as their use in gene therapy strategies for the treatment of neurological disorders including Parkinson's disease, Alzheimer's disease, and Huntington's disease.

TCF4-mediated Fuchs endothelial corneal dystrophy: Insights into a common trinucleotide repeat-associated disease

Fuchs endothelial corneal dystrophy (FECD) is a common cause for heritable visual loss in the elderly. Since the first description of an association between FECD and common polymorphisms situated within the transcription factor 4 (TCF4) gene, genetic and molecular studies have implicated an intronic CTG trinucleotide repeat (CTG18.1) expansion as a causal variant in the majority of FECD patients. To date, several non-mutually exclusive mechanisms have been proposed that drive and/or exacerbate the onset of disease. These mechanisms include (i) TCF4 dysregulation; (ii) toxic gain-of-function from TCF4 repeat-containing RNA; (iii) toxic gain-of-function from repeat-associated non-AUG dependent (RAN) translation; and (iv) somatic instability of CTG18.1. However, the relative contribution of these proposed mechanisms in disease pathogenesis is currently unknown. In this review, we summarise research implicating the repeat expansion in disease pathogenesis, define the phenotype-genotype correlations between FECD and CTG18.1 expansion, and provide an update on research tools that are available to study FECD as a trinucleotide repeat expansion disease. Furthermore, ongoing international research efforts to develop novel CTG18.1 expansion-mediated FECD therapeutics are highlighted and we provide a forward-thinking perspective on key unanswered questions that remain in the field.

•

FECD is a common, age-related corneal dystrophy.

•

The majority of cases are associated with expansion of a CTG repeat (CTG18.1).

•

FECD is the most common trinucleotide repeat expansion disease in humans.

•

Evidence supports multiple molecular mechanisms underlying the pathophysiology.

•

Novel CTG18.1-targeted therapeutics are in development.

Genomic Diagnosis for Pediatric Disorders: Revolution and Evolution

Powerful, recent advances in technologies to analyze the genome have had a profound impact on the practice of medical genetics, both in the laboratory and in the clinic. Increasing utilization of genome-wide testing such as chromosomal microarray analysis and exome sequencing have lead a shift toward a "genotype-first" approach. Numerous techniques are now available to diagnose a particular syndrome or phenotype, and while traditional techniques remain efficient tools in certain situations, higher-throughput technologies have become the de facto laboratory tool for diagnosis of most conditions. However, selecting the right assay or technology is challenging, and the wrong choice may lead to prolonged time to diagnosis, or even a missed diagnosis. In this review, we will discuss current core technologies for the diagnosis of classic genetic disorders to shed light on the benefits and disadvantages of these strategies, including diagnostic efficiency, variant interpretation, and secondary findings. Finally, we review upcoming technologies posed to impart further changes in the field of genetic diagnostics as we move toward "genome-first" practice.

The third generation sequencing: the advanced approach to genetic diseases

Genomic sequencing technologies have revolutionized mutation detection of the genetic diseases in the past few years. In recent years, the third generation sequencing (TGS) has been gaining insight into more genetic diseases owing to the single molecular and real time sequencing technology. This paper reviews the genomic sequencing revolutionary history first and then focuses on the genetic diseases discovered through the TGS and the clinical effects of the TGS, which is followed by the discussion of the improvement in the bioinformatic analysis for the TGS and its limitations. In summary, the TGS has been enhancing the diagnostic accuracy of genetic diseases in molecular level as well as paving a new way for basic researches and therapies.

Structural variant identification and characterization

Structural variant (SV) differences between human genomes can cause germline and mosaic disease as well as inter-individual variation. De-regulation of accurate DNA repair and genomic surveillance mechanisms results in a large number of SVs in cancer. Analysis of the DNA sequences at SV breakpoints can help identify pathways of mutagenesis and regions of the genome that are more susceptible to rearrangement. Large-scale SV analyses have been enabled by high-throughput genome-level sequencing on humans in the past decade. These studies have shed light on the mechanisms and prevalence of complex genomic rearrangements. Recent advancements in both sequencing and other mapping technologies as well as calling algorithms for detection of genomic rearrangements have helped propel SV detection into population-scale studies, and have begun to elucidate previously inaccessible regions of the genome. Here, we discuss the genomic organization of simple and complex SVs, the molecular mechanisms of their formation, and various ways to detect them. We also introduce methods for characterizing SVs and their consequences on human genomes.

Critical assessment of bioinformatics methods for the characterization of pathological repeat expansions with single-molecule sequencing data

A number of studies have reported the successful application of single-molecule sequencing technologies to the determination of the size and sequence of pathological expanded microsatellite repeats over the last 5 years. However, different custom bioinformatics pipelines were employed in each study, preventing meaningful comparisons and somewhat limiting the reproducibility of the results. In this review, we provide a brief summary of state-of-the-art methods for the characterization of expanded repeats alleles, along with a detailed comparison of bioinformatics tools for the determination of repeat length and sequence, using both real and simulated data. Our reanalysis of publicly available human genome sequencing data suggests a modest, but statistically significant, increase of the error rate of single-molecule sequencing technologies at genomic regions containing short tandem repeats. However, we observe that all the methods herein tested, irrespective of the strategy used for the analysis of the data (either based on the alignment or assembly of the reads), show high levels of sensitivity in both the detection of expanded tandem repeats and the estimation of the expansion size, suggesting that approaches based on single-molecule sequencing technologies are highly effective for the detection and quantification of tandem repeat expansions and contractions.

Long-read sequencing for rare human genetic diseases

During the past decade, the search for pathogenic mutations in rare human genetic diseases has involved huge efforts to sequence coding regions, or the entire genome, using massively parallel short-read sequencers. However, the approximate current diagnostic rate is <50% using these approaches, and there remain many rare genetic diseases with unknown cause. There may be many reasons for this, but one plausible explanation is that the responsible mutations are in regions of the genome that are difficult to sequence using conventional technologies (e.g., tandem-repeat expansion or complex chromosomal structural aberrations). Despite the drawbacks of high cost and a shortage of standard analytical methods, several studies have analyzed pathogenic changes in the genome using long-read sequencers. The results of these studies provide hope that further application of long-read sequencers to identify the causative mutations in unsolved genetic diseases may expand our understanding of the human genome and diseases. Such approaches may also be applied to molecular diagnosis and therapeutic strategies for patients with genetic diseases in the future.

CRISPR/Cas9-targeted enrichment and long-read sequencing of the Fuchs endothelial corneal dystrophy-associated TCF4 triplet repeat

PURPOSE

To demonstrate the utility of an amplification-free long-read sequencing method to characterize the Fuchs endothelial corneal dystrophy (FECD)-associated intronic TCF4 triplet repeat (CTG18.1).

METHODS

We applied an amplification-free method, utilizing the CRISPR/Cas9 system, in combination with PacBio single-molecule real-time (SMRT) long-read sequencing, to study CTG18.1. FECD patient samples displaying a diverse range of CTG18.1 allele lengths and zygosity status (n = 11) were analyzed. A robust data analysis pipeline was developed to effectively filter, align, and interrogate CTG18.1-specific reads. All results were compared with conventional polymerase chain reaction (PCR)-based fragment analysis.

RESULTS

CRISPR-guided SMRT sequencing of CTG18.1 provided accurate genotyping information for all samples and phasing was possible for 18/22 alleles sequenced. Repeat length instability was observed for all expanded (≥50 repeats) phased CTG18.1 alleles analyzed. Furthermore, higher levels of repeat instability were associated with increased CTG18.1 allele length (mode length ≥91 repeats) indicating that expanded alleles behave dynamically.

CONCLUSION

CRISPR-guided SMRT sequencing of CTG18.1 has revealed novel insights into CTG18.1 length instability. Furthermore, this study provides a framework to improve the molecular diagnostic accuracy for CTG18.1-mediated FECD, which we anticipate will become increasingly important as gene-directed therapies are developed for this common age-related and sight threatening disease.

Amplification-free long-read sequencing of TCF4 expanded trinucleotide repeats in Fuchs Endothelial Corneal Dystrophy

Amplification of a CAG trinucleotide motif (CTG18.1) within the TCF4 gene has been strongly associated with Fuchs Endothelial Corneal Dystrophy (FECD). Nevertheless, a small minority of clinically unaffected elderly patients who have expanded CTG18.1 sequences have been identified. To test the hypothesis that the CAG expansions in these patients are protected from FECD because they have interruptions within the CAG repeats, we utilized a combination of an amplification-free, long-read sequencing method and a new target-enrichment sequence analysis tool developed by Pacific Biosciences to interrogate the sequence structure of expanded repeats. The sequencing was successful in identifying a previously described interruption within an unexpanded allele and provided sequence data on expanded alleles greater than 2000 bases in length. The data revealed considerable heterogeneity in the size distribution of expanded repeats within each patient. Detailed analysis of the long sequence reads did not reveal any instances of interruptions to the expanded CAG repeats, but did reveal novel variants within the AGG repeats that flank the CAG repeats in two of the five samples from clinically unaffected patients with expansions. This first examination of the sequence structure of CAG repeats in CTG18.1 suggests that factors other than interruptions to the repeat structure account for the absence of disease in some elderly patients with repeat expansions in the TCF4 gene.

Long-Read Sequencing Emerging in Medical Genetics

The wide implementation of next-generation sequencing (NGS) technologies has revolutionized the field of medical genetics. However, the short read lengths of currently used sequencing approaches pose a limitation for the identification of structural variants, sequencing repetitive regions, phasing of alleles and distinguishing highly homologous genomic regions. These limitations may significantly contribute to the diagnostic gap in patients with genetic disorders who have undergone standard NGS, like whole exome or even genome sequencing. Now, the emerging long-read sequencing (LRS) technologies may offer improvements in the characterization of genetic variation and regions that are difficult to assess with the prevailing NGS approaches. LRS has so far mainly been used to investigate genetic disorders with previously known or strongly suspected disease loci. While these targeted approaches already show the potential of LRS, it remains to be seen whether LRS technologies can soon enable true whole genome sequencing routinely. Ultimately, this could allow the de novo assembly of individual whole genomes used as a generic test for genetic disorders. In this article, we summarize the current LRS-based research on human genetic disorders and discuss the potential of these technologies to facilitate the next major advancements in medical genetics.

Tandem-genotypes: robust detection of tandem repeat expansions from long DNA reads

Tandemly repeated DNA is highly mutable and causes at least 31 diseases, but it is hard to detect pathogenic repeat expansions genome-wide. Here, we report robust detection of human repeat expansions from careful alignments of long but error-prone (PacBio and nanopore) reads to a reference genome. Our method is robust to systematic sequencing errors, inexact repeats with fuzzy boundaries, and low sequencing coverage. By comparing to healthy controls, we prioritize pathogenic expansions within the top 10 out of 700,000 tandem repeats in whole genome sequencing data. This may help to elucidate the many genetic diseases whose causes remain unknown.